since vision loss from glaucoma can not be recovered , improved screening and detection methods for the angle closure glaucoma are essential to preserve vision and maintain a good quality of life .
although glaucoma commonly progresses to blindness over many years , acute angle closure can result in permanent blindness in a matter of hours .
studies suggest that glaucomatous optic neuropathology can be prevented when effective prophylactic treatment such as laser peripheral iridotomy is performed at an early and appropriate time for the eyes with anatomically narrow angles .
angle closure glaucoma is characterized by obstruction of aqueous fluid drainage through the trabecular meshwork from the eye 's anterior chamber .
a wide unobstructed iridocorneal angle allows sufficient drainage of aqueous humor , whereas a narrow angle may impede the drainage system and leave the patient susceptible to angle closure glaucoma .
the imaging of angle between the iris and the cornea is the key for open angle and closed angle glaucoma diagnosis
gonioscopy , is considered the gold standard technique to examine the iridocorneal angle , but it is subjective and dependent on the operator .
oct is a high resolution , rapid , and noninvasive screening tool for angle closure glaucoma and its technology has evolved rapidly from time - domain to spectral - domain .
swept source oct ( casia ss-1000 )   is a newly developed novel imaging technology which provides a detailed examination of the structures of anterior chamber ( ac ) .
it is a fourier - domain system and is designed specifically for imaging the anterior segment . with a substantial improvement in scan speed ( 30,000 a - scans per second ) , the anterior chamber angles can be imaged in 128 cross - sections ( each with 512 a - scans ) 360 degree around the anterior segment in 1.2 seconds .
high - definition ss - oct ( ss - oct ) imaging has the potential to become an important tool for the assessment of the anterior chamber angle and detection of angle closure [ 4 , 5 ] .
anterior segment ss - oct imaging has significantly altered the diagnosis and evaluation of acg .
the information gained with new imaging modalities provides clinicians with both qualitative and quantitative information about anatomical structure of the anterior chamber .
hu et al .   compared gonioscopy with visante and cirrus optical coherence tomography ( oct ) for identifying angle structures and the presence of angle closure in patients with glaucoma .
several techniques to assess the anterior chamber are based on manual or automatic detection of some landmarks such as scleral spur ( ss )   and schwalbe 's line ( sl ) [ 8 , 9 ] .
however , the existing methods use a single distance or area , for example , aod500 , tisa500 , and aodsl , for measurement of iridocorneal angle , without consideration of the whole angle profile .
cheung et al .   showed that an irregular iris surface made the aodsl very inaccurate .
our previous study   measured two new parameters , maod and atsl , based on the continuous aod to assess the anterior chamber angle .
these two new parameters tried to overcome the limitations in the single measurement of aod when the iris surface is irregular and the angle is more occludable .
the iris curvature of anterior segment oct ( as - oct ) is usually calculated by drawing a line from the most peripheral to the most central points of the iris pigment epithelium , and then a perpendicular line is extended from this line to the iris pigment epithelium at the point of the greatest convexity .
however , the iris pigment epithelium can not be exactly detected , only the stoma .
the segmentation of the iris lower boundary is challenging due to the low contrast in this region .
the nature of the angle dynamics which has not been considered yet in the previous studies , still involved in different forms of angle closure or open angle glaucoma .
hence , we propose to measure the iris curvature denoted by hiris from the extracted iris image . besides
, we assume that the iridocorneal angle is the irregular shape and the complexity or changes of the iridocorneal angle can be measured quantitatively and qualitatively with shape descriptors .
the shape analysis might give more freedom in computation and is less sensitive to the accurate detection of the landmarks  the scleral spur ( ss ) and schwalbe 's line ( sl ) . using this approach ,
the shape analysis may be successfully applied both quantitatively and objectively to characterize angle shape of anterior chamber of ss - oct images .
the complex and erratic shape description in terms of self - similarity was introduced by mandelbrot .
the concept of fractal dimension ( fd ) is useful in the measurement , analysis , and classification of shape and texture .
the fractal dimension , therefore , might serve as a sensitive descriptor of the iridocorneal angle shape .
fractals provide a new field for characterization of irregularity and complexity , yet self - similar structures in nature . however , although fd has been used extensively in characterizing self - affinity in various kinds of biomedical research [ 13 , 14 ] , little attention has been paid in automated glaucoma subtype classification using the feasibility of fractal and multifractal theory on retinal nerve fiber layer ( rnfl ) and optic disc .
the measurement of moment invariants like the triangularity ( ti ) and ellipticity ( ei )   was successfully used in several applications like the classification of mammographic masses and lung field boundaries .
both shape analyses are a translation , rotation , and scale of the object and general affine transformations .
therefore , the paper proposes a new strategy to analyze the iridocorneal angle by morphological shape analysis as well as biometric angle parameters . in this paper
, we present an automatic angle closure glaucoma detection system based on machine learning and image analysis for the estimation of quantitative parameters to classify ss - oct images into two classes : open angle and angle closure glaucoma .
we propose new strategies of feature extraction methods : ( 1 ) based on measurement of biometric parameters and ( 2 ) based on the shape analysis of iridocorneal angle to capture much more information .
section 2 describes the overview of the algorithm for preprocessing , feature extractions , and classification method of angle closure glaucoma .
it mainly consists of four steps : preprocessing , segmentation of anterior chamber , anterior chamber angle analysis for feature extractions , and classification of angle closure and open angle glaucoma .
this section starts by defining the notation that will be used throughout the paper and figure 2 shows the region of interest for iridocorneal angle analysis .
firstly , the original image denoted by i(x , y ) is preprocessed to remove the vertical saturation artifacts as in   and the processed image is denoted by i^(x , y ) , where ( x , y)[1 , m][1 , n ] are the pixels coordinates and m  n is the size of the image ; then i^(x , y ) is segmented to detect the cornea , the iris , and the anterior chamber .
the lower boundaries of the cornea ee and inner boundary of the iris ei are detected using the same approach as the method in , and the iris image denoted by ii(x , y ) is extracted .
the six features are extracted to quantify the anterior chamber ; the anterior chamber assessment parameters are maod , atsl , the mean curvature of the iris hiris , and shape descriptors fd , ti , and ei .
the definitions of the variables used in this paper are listed here : original gray scale image : i(x , y ) , where ( x , y)[1 , m][1 , n];the preprocessed image : i^(x , y ) , where ( x , y)[1 , m][1 , n];the segmented corneal endothelium , the anterior surface of the iris , and anterior chamber mask in i(x , y ) : ee , ei , and mac;the extracted iris image : ii(x , y);automatically detected schwalbe 's line : sl;angle opening distance : aod;mean of the continuous measurement of the angle opening distance : maod;area of the trapezoidal profile of the iridocorneal angle centered at schwalbe 's line : atsl;mean of the iris curvature from the extracted iris image ; hiris;the extracted iridocorneal angle image : iira(x , y ) , where ( x , y)[1,m][1,n];complex shape descriptor to quantify the complexity or changes of iridocorneal angle : fractal dimension ( fd);ellipticity moment shape descriptor : ei;triangularity moment shape descriptor : ti .
original gray scale image : i(x , y ) , where ( x , y)[1 , m][1 , n ] ; the preprocessed image : i^(x , y ) , where ( x , y)[1 , m][1 , n ] ; the segmented corneal endothelium , the anterior surface of the iris , and anterior chamber mask in i(x , y ) : ee , ei , and mac ; the extracted iris image : ii(x , y ) ; automatically detected schwalbe 's line : sl ; angle opening distance : aod ; mean of the continuous measurement of the angle opening distance : maod ; area of the trapezoidal profile of the iridocorneal angle centered at schwalbe 's line : atsl ; mean of the iris curvature from the extracted iris image ; hiris ; the extracted iridocorneal angle image : iira(x , y ) , where ( x , y)[1,m][1,n ] ; complex shape descriptor to quantify the complexity or changes of iridocorneal angle : fractal dimension ( fd ) ; ellipticity moment shape descriptor : ei ; triangularity moment shape descriptor : ti .
it is observed that some ss - oct images contain vertical saturation artifacts , which hinder the accurate interpretation of the image if thresholding and component labeling are used for the segmentation .
hence , reducing and removing the artifacts are performed prior to the segmentation step . as illustrated in figure 3(a ) , the vertical saturation artifacts are marked by sudden increases in intensity as compared to surrounding area and are caused by the saturation of intensity due to the strong reflection in certain locations .
therefore , we detected such artifacts by searching for abrupt changes in the average intensity of each column as shown in figure 3(b ) .
after removing the vertical saturation artifacts , the segmentation was performed in the ss - oct images .
the method consists of segmenting the anterior chamber , the cornea , and the iris and extracting their edges .
the anterior chamber region , the lower cornea boundary , and the upper iris boundary are extracted by image segmentation method .
this section describes the analysis of aca for the extraction of six features by the biometric parameter measurement and shape analysis .
the biometric parameters , previously proposed in , from the iridocorneal angle based on the continuous measurement of aod are quantified as follows.maod : the average of continuous serial aod measured every 25  m away from the sl till 500  m in both directions ( anterior and posterior to / from schwalbe 's line ) as shown in figure 4(a).atsl : the trapezoidal area of iridocorneal angle bounded by the angle recess , aod_psl line , corneal endothelium , and anterior surface of the iris as shown in figure 4(a).hiris : the mean iris curvature measured from the whole iris image as shown in figure 4(b ) .
maod : the average of continuous serial aod measured every 25  m away from the sl till 500  m in both directions ( anterior and posterior to / from schwalbe 's line ) as shown in figure 4(a ) .
atsl : the trapezoidal area of iridocorneal angle bounded by the angle recess , aod_psl line , corneal endothelium , and anterior surface of the iris as shown in figure 4(a ) .
hiris : the mean iris curvature measured from the whole iris image as shown in figure 4(b ) . after calculating the biometric parameters , we performed anterior chamber shape analysis on the selected region of iridocorneal angle image as shown in figure 5 .
we utilized the fractal complex shape descriptor on the iridocorneal angle image denoted by iira(x , y ) , where ( x , y)[1,m][1,n ] are the pixels coordinates and mn is the size of the extracted angle image .
moment shape descriptors such as the triangularity and ellipticity are used to compare the performance of fractal dimension analysis . the angle structure configuration is then quantified by the following three features : complex shape descriptor to quantify the complexity or changes of iridocorneal angle , fractal dimension ( fd);ellipticity moment shape descriptor ( ei);triangularity moment shape descriptor ( ti ) .
complex shape descriptor to quantify the complexity or changes of iridocorneal angle , fractal dimension ( fd ) ; ellipticity moment shape descriptor ( ei ) ; triangularity moment shape descriptor ( ti ) .
the ability of each feature extraction method to separate open angle and angle closure glaucoma is quantified by the fuzzy k - nearest neighbour classifier .
the basis of the fuzzy k - nn algorithm is to assign membership as a function of the vector 's distance from its k - nearest neighbors and those neighbors ' memberships in the possible classes .
the singaporean chinese population study by the singapore national eye centre recruited 148 subjects ( 90 females and 45 males ) with average age of 59.48  8.97 .
all subjects underwent a standard examination of dark room gonioscopy and the anterior segment imaging by ss - oct on the same day .
the casia ss-1000 oct ( tomey , nagoya , japan ) was used as the imaging modality to visualize the anterior segment of the eye .
we used the 2d angle high - definition ( hd ) mode of ss - oct imaging with the ( 8  mm , 8  mm ) scan dimension . in the hd scan mode , the identification of both the sl and ss was possible in over 98% of ss - oct images .
the aca was graded using the modified shaffer grading system in each ss - oct image .
the doctor examined four different quadrants of the eye , namely , inferior , superior , nasal , and temporal ( i , s , n , and t ) scanning .
a subset of 29 subjects ( 23.4% ) was imaged bilaterally to assess the differences between eyes .
we selected the images with two criteria : ( 1 ) sl could be identified automatically and ( 2 ) only one image from the four scans of each eye , for less bias for classifier .
so , in total 264 ss - oct hd images in which sl could be seen were selected from the dataset for further analysis .
they were 132 nasal scan images , 70 temporal scan images , 29 superior scan images , and 33 inferior scan images .
using gonioscopic shaffer grading as the gold standard , 135 images with closed angle ( 1.131  0.72 grade ) and 129 images with open angle ( 3.321  0.469 grade ) were evaluated for the analysis .
the vertical saturation artifacts were presented in 77 images ( ~46% ) out of 264 images in the dataset on various regions of image .
so , it is necessary to perform artifact removal for affected images prior to image segmentation .  
figure 6 shows the segmentation results of with and without artifact removal in ss - oct image . from the results
as shown in figure 6 , we found that the preprocessing steps relatively improved the segmentation of the anterior chamber region . 

figure 7(a ) shows the segmentation results of the anterior chamber , the upper cornea , the lower cornea , and the inner boundary of the iris .
the extracted iris image is shown in figure 7(b ) which yields the iris curvature of 0.294 .  
figure 7(c ) shows an example of the resulting angle profile consisting of 40 continuous aod centered on sl which reveals 0.262  mm for maod and 0.276  mm for atsl , respectively .
then , we find the region of interest for the fractal dimension and the moment shape analysis as shown in figure 7(d ) .
the results of the triangularity and the ellipticity moment shape analysis are 0.899 and 0.165 , respectively , as shown in figure 7(e ) .  
figure 7(f ) shows the fractal image of the region of interest of figure 7(d ) and its estimated average fd is 1.862 .
the comparison of the iridocorneal angle features between open angle and angle closure glaucoma ss - oct images using gonioscopic grading as a reference is illustrated in table 1 .
the fractal dimensions of open and closed angle are 1.87 and 1.84 , respectively , as shown in table 1 .
the results of fractal dimension depend slightly on the specific images of open angle and angle closure glaucoma
. all the features of closed angle ss - oct images are smaller than those of the open angle ss - oct images .
the coefficients of variance ( covs ) of fd are smaller in both open and closed angle images .
the covs of hiris are the largest in both open and closed angle ss - oct images .
so , hiris results are considered to be high variance and do not necessarily correspond to closed and open angle images in this study .  
figure 8 shows the scatter plot of fd , maod , and atsl features for both open angle and closed angle ss - oct images .
after the feature vectors were computed , training , cross - validation , and testing sets were formed by 1476 vectors ( 264  images  6  features ) : 3 biometric features ( maod , atsl , and hiris ) and 3 shape descriptors ( fd , ti , and ei ) .
then , we perform the feature selections from all features for improving classification accuracy or decreasing the size of the structure without significantly decreasing classification accuracy of the classifier which is built using only the selected features . reducing the number of irrelevant / redundant features
drastically reduces the running time of a learning algorithm and yields a more general concept .
this helps in getting a better insight into the underlying concept of a real - world classification problem .
feature selection methods try to pick a subset of features that are relevant to the target concept and these features are used as the input of classifiers .
the association between the gonioscopic grading and the measured iridocorneal angle features was evaluated for feature selection using spearman correlation coefficient (  ) .
there is a high correlation between fd , maod , and atsl and gonioscopic grading as shown in table 2 .
so , we grouped those features for classification and compared the consistency with other separate features . in the system , these features were firstly normalized to mean zero and variance of one to improve the classification process .
then , the classification was performed based on the normalized features by fuzzy k - nearest neighbor classifier ( fknn ) .
the performance of the fknn was evaluated by comparison with widely used machine learning algorithms , namely , linear discriminant analysis ( lda ) , k - nearest neighbour ( knn ) , and support vector machines ( svm ) , to verify the effectiveness of the proposed model .
the classification was performed on the 5-fold cross - validation ( 5-fold cv ) of various combinations of features .
to verify the effectiveness of the proposed fknn classifier , we firstly find the relationship between the classification performance and the fuzzy strength parameter m which varies from 1 to 2 , with the step size of 0.1 .
it can be observed that the accuracy was achieved between 93% and 99% and auc value fluctuates between 87% and 99% based on the feature groups as shown in figures 9(a ) and 9(b ) .
it reveals that fuzzy strength parameter m has a big impact on the performance of fknn classifier .
the average 5-fold cv accuracy and the corresponding standard deviation of fknn are shown in table 3 .
the classification using fd , maod , and atsl features group yields the best performance accuracy 99.11%  0.76% and auc 0.98  0.012 .
these findings suggest that the proposed method is highly promising in providing accurate diagnosis tools for acg and will make a greater clinical impact if the study can be done on a larger image database .
we also evaluated the performance of fuzzy knn by using only three features ( fd , maod , and atsl ) with other classification methods ( lda , knn , and svm ) .
svm slightly outperforms lda with an accuracy of 83.26%  4.74% ( 87.78%  7.92% sensitivity and 68.4%  9.29% specificity ) and an auc of 0.78  0.03 was reached .
knn was found to outperform lda as well as svm with an accuracy of 85.80%  5.16% ( 90.10%  5.77% sensitivity and 71.82%  5.41% specificity ) and 0.80  0.04 auc value .
the best classification performance was provided by the fuzzy knn classifier with fuzzy membership of the samples ( ~18% ) better than other classifers ( lda , knn , and svm ) as in table 4 .
additionally , fknn has been shown to be an efficient and robust classification method to diagnose angle closure glaucoma in ss - oct images .
the proposed fractal shape descriptor also showed better accuracy in comparison to the moment shape descriptor method .
the moment shape descriptor method is shown in the literature as an efficient technique to obtain shape descriptors . however , the results presented in this work suggest that the fractal analysis is a worthy option for providing shape descriptors for classification tasks , as it is invariant to rotation , translation , and scale . besides , our fractal dimension analysis is performed for each pixel on gray level images that can capture more detailed information of the angle structure .
in addition , that fractal shape analysis gives more freedom in computation and is less sensitive to the accurate detection of the landmarks  the scleral spur ( ss ) and schwalbe 's line ( sl ) . moreover , our classification results are comparable to those of xu et al . .
they achieved the classification accuracies of 0.921  0.036 auc and 84.0%  5.7% balanced accuracy at an 85% specificity using histogram equalized pixel ( hep ) values and svm classification in oct images .
we also observed a significant advantage in terms of classification performance of using the fuzzy knn algorithms in glaucoma diagnosis . in summary
, the proposed framework has been shown to be a useful tool in screening for acg . while our system can provide useful classification and support to the medical experts through identification of features , human intervention to exploit the extracted knowledge
we believe that our findings in this study can serve as a basic grading system for angle closure glaucoma diagnosis in the future .
we evaluated several techniques for extracting useful information from ss - oct images such as traditional biometric parameter measurement and complex shape descriptors using fractal dimension analysis .
the experimental results demonstrated that the proposed technique which is the combination of biometric parameter , fractal dimension analysis , and classification by fuzzy knn method achieved great accuracy in the classification of the open and closed angle glaucoma images with an accuracy rate of 99.11%  0.76% and 0.98  0.012 auc value .
the performance of the fully automatic system presented here is comparable to medical experts in detecting glaucomatous eyes and could argue clinicians ' diagnosis of angle closure glaucoma .
 optical coherence tomography is a high resolution , rapid , and noninvasive diagnostic tool for angle closure glaucoma . in this paper 
 , we present a new strategy for the classification of the angle closure glaucoma using morphological shape analysis of the iridocorneal angle . 
 the angle structure configuration is quantified by the following six features : ( 1 ) mean of the continuous measurement of the angle opening distance ; ( 2 ) area of the trapezoidal profile of the iridocorneal angle centered at schwalbe 's line ; ( 3 ) mean of the iris curvature from the extracted iris image ; ( 4 ) complex shape descriptor , fractal dimension , to quantify the complexity , or changes of iridocorneal angle ; ( 5 ) ellipticity moment shape descriptor ; and ( 6 ) triangularity moment shape descriptor . 
 then , the fuzzy k nearest neighbor ( fknn ) classifier is utilized for classification of angle closure glaucoma . 
 two hundred and sixty - four swept source optical coherence tomography ( ss - oct ) images from 148 patients were analyzed in this study . from the experimental results , 
 the fknn reveals the best classification accuracy ( 99.11  0.76% ) and auc ( 0.98  0.012 ) with the combination of fractal dimension and biometric parameters . 
 it showed that the proposed approach has promising potential to become a computer aided diagnostic tool for angle closure glaucoma ( acg ) disease . 
 
1. Introduction
2. Methods
3. Experiments and Results
4. Conclusions