# Anterior Chamber Angle Shape Analysis and Classification of Glaucoma in SS-OCT Images

Authors: Soe Ni Ni; J. Tian; Pina Marziliano; Hong-Tym Wong
Journal: Journal of Ophthalmology
 (2014)

---

## Abstract

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 fuzzyk 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.


---

## Body

## 1. Introduction

The detection of the angle closure glaucoma is important for preventing irreversible blindness. Since vision loss from glaucoma cannot 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 [1]. Angle closure glaucoma is characterized by obstruction of aqueous fluid drainage through the trabecular meshwork from the eye’s anterior chamber. The width of the angle is one factor affecting the drainage of aqueous humor. 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. Early detection of ACG using imaging technology has recently gained much clinical interest. Gonioscopy [2], 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) [3] 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. [6] 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) [7] 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. [10] showed that an irregular iris surface made the AODsl very inaccurate. Our previous study [11] 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 [10].In addition to these two parameters, the iris curvature must also be considered. 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 [12], 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 cannot 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. An important property of shapes is their complexity. The complex and erratic shape description in terms of self-similarity was introduced by Mandelbrot [13]. 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 [15]. Another approach in shape analysis is founded on convexity [16]. The measurement of moment invariants like the triangularity (TI) and ellipticity (EI) [17] was successfully used in several applications like the classification of mammographic masses and lung field boundaries [18]. Both shape analyses are a translation, rotation, and scale of the object and general affine transformations [19]. 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.The rest of the paper is organized as follows. Section2 describes the overview of the algorithm for preprocessing, feature extractions, and classification method of angle closure glaucoma. The experiment and results are presented in Section 3. Finally Section 4 concludes the paper with the future work.


## 2. Methods

The architecture of the overall proposed system is shown in Figure1. 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.Figure 1
Overview of the automatic classification system for open angle and angle closure glaucoma using SS-OCT images.Figure 2
Anterior chamber imaging showing Schwalbe’s line and yellow bounding box for region of interest (ROI) for angle analysis.Firstly, the original image denoted byI(x,y) is preprocessed to remove the vertical saturation artifacts as in [11] 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 [9], 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. Lastly, the classification step is performed. The short descriptions of each step are presented in the following subsequent sections. The extended description of each step can be found in the appendices.


### 2.1. Notation

The definitions of the variables used in this paper are listed here:(1)
original gray scale image:I(x,y), where (x,y)∈[1,M]×[1,N];(2)
the preprocessed image:I^(x,y), where (x,y)∈[1,M]×[1,N];(3)
the segmented corneal endothelium, the anterior surface of the iris, and anterior chamber mask inI(x,y): Ee, Ei, and MAC⁡;(4)
the extracted iris image:II(x,y);(5)
automatically detected Schwalbe’s line: SL;(6)
angle opening distance: AOD;(7)
mean of the continuous measurement of the angle opening distance: mAOD;(8)
area of the trapezoidal profile of the iridocorneal angle centered at Schwalbe’s line: ATsl;(9)
mean of the iris curvature from the extracted iris image;Hiris;(10)
the extracted iridocorneal angle image:IIRA(x,y), where (x,y)∈[1,M¯]×[1,N¯];(11)
complex shape descriptor to quantify the complexity or changes of iridocorneal angle: fractal dimension (FD);(12)
ellipticity moment shape descriptor:EI;(13)
triangularity moment shape descriptor:TI.


### 2.2. Preprocessing of SS-OCT Images

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 Figure3(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).(a) The original closed angle SS-OCT image with contaminated vertical saturation artifact; (b) detection of the artifact by searching the abrupt changes in intensities of the SS-OCT images where the vertical dotted lines denote the region of detected vertical artifact and the horizontal dotted lines denote the mean intensity of the SS-OCT images.
(a)
(b)


### 2.3. Segmentation of Anterior Chamber and Iris

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. The details of segmentation can be found in appendices.


### 2.4. Feature Extractions

This section describes the analysis of ACA for the extraction of six features by the biometric parameter measurement and shape analysis.


#### 2.4.1. Measurement of Biometric Parameters

The biometric parameters, previously proposed in [11], from the iridocorneal angle based on the continuous measurement of AOD are quantified as follows.(i)
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).(ii)
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).(iii)
Hiris: the mean iris curvature measured from the whole iris image as shown in Figure 4(b).(a) Biometric parameter measurement in SS-OCT image. (b) Extracted region of iris to compute iris curvature. SL: Schwalbe’s line, ATsl: area of trapezoidal profile of iridocorneal angle centered at SL, AOD_psl: angle opening distance 500 μm posterior from SL, Hiris: the mean curvature of iris.
(a)
(b)


#### 2.4.2. Anterior Chamber Shape Analysis

After calculating the biometric parameters, we performed anterior chamber shape analysis on the selected region of iridocorneal angle image as shown in Figure5. 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 M¯×N¯ 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:(1)
complex shape descriptor to quantify the complexity or changes of iridocorneal angle, fractal dimension (FD);(2)
ellipticity moment shape descriptor (EI);(3)
triangularity moment shape descriptor (TI).Figure 5
Shape analysis on the region of interest in SS-OCT image where SL is Schwalbe’s line and AOD_psl is the angle opening distance 500 μm posterior from SL.


### 2.5. Classification

The ability of each feature extraction method to separate open angle and angle closure glaucoma is quantified by the fuzzyk-nearest neighbour classifier [20]. The basis of the fuzzy k-NN algorithm is to assign membership as a function of the vector’s distance from itsk-nearest neighbors and those neighbors’ memberships in the possible classes.


## 2.1. Notation

The definitions of the variables used in this paper are listed here:(1)
original gray scale image:I(x,y), where (x,y)∈[1,M]×[1,N];(2)
the preprocessed image:I^(x,y), where (x,y)∈[1,M]×[1,N];(3)
the segmented corneal endothelium, the anterior surface of the iris, and anterior chamber mask inI(x,y): Ee, Ei, and MAC⁡;(4)
the extracted iris image:II(x,y);(5)
automatically detected Schwalbe’s line: SL;(6)
angle opening distance: AOD;(7)
mean of the continuous measurement of the angle opening distance: mAOD;(8)
area of the trapezoidal profile of the iridocorneal angle centered at Schwalbe’s line: ATsl;(9)
mean of the iris curvature from the extracted iris image;Hiris;(10)
the extracted iridocorneal angle image:IIRA(x,y), where (x,y)∈[1,M¯]×[1,N¯];(11)
complex shape descriptor to quantify the complexity or changes of iridocorneal angle: fractal dimension (FD);(12)
ellipticity moment shape descriptor:EI;(13)
triangularity moment shape descriptor:TI.


## 2.2. Preprocessing of SS-OCT Images

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 Figure3(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).(a) The original closed angle SS-OCT image with contaminated vertical saturation artifact; (b) detection of the artifact by searching the abrupt changes in intensities of the SS-OCT images where the vertical dotted lines denote the region of detected vertical artifact and the horizontal dotted lines denote the mean intensity of the SS-OCT images.
(a)
(b)


## 2.3. Segmentation of Anterior Chamber and Iris

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. The details of segmentation can be found in appendices.


## 2.4. Feature Extractions

This section describes the analysis of ACA for the extraction of six features by the biometric parameter measurement and shape analysis.


### 2.4.1. Measurement of Biometric Parameters

The biometric parameters, previously proposed in [11], from the iridocorneal angle based on the continuous measurement of AOD are quantified as follows.(i)
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).(ii)
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).(iii)
Hiris: the mean iris curvature measured from the whole iris image as shown in Figure 4(b).(a) Biometric parameter measurement in SS-OCT image. (b) Extracted region of iris to compute iris curvature. SL: Schwalbe’s line, ATsl: area of trapezoidal profile of iridocorneal angle centered at SL, AOD_psl: angle opening distance 500 μm posterior from SL, Hiris: the mean curvature of iris.
(a)
(b)


### 2.4.2. Anterior Chamber Shape Analysis

After calculating the biometric parameters, we performed anterior chamber shape analysis on the selected region of iridocorneal angle image as shown in Figure5. 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 M¯×N¯ 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:(1)
complex shape descriptor to quantify the complexity or changes of iridocorneal angle, fractal dimension (FD);(2)
ellipticity moment shape descriptor (EI);(3)
triangularity moment shape descriptor (TI).Figure 5
Shape analysis on the region of interest in SS-OCT image where SL is Schwalbe’s line and AOD_psl is the angle opening distance 500 μm posterior from SL.


## 2.4.1. Measurement of Biometric Parameters

The biometric parameters, previously proposed in [11], from the iridocorneal angle based on the continuous measurement of AOD are quantified as follows.(i)
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).(ii)
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).(iii)
Hiris: the mean iris curvature measured from the whole iris image as shown in Figure 4(b).(a) Biometric parameter measurement in SS-OCT image. (b) Extracted region of iris to compute iris curvature. SL: Schwalbe’s line, ATsl: area of trapezoidal profile of iridocorneal angle centered at SL, AOD_psl: angle opening distance 500 μm posterior from SL, Hiris: the mean curvature of iris.
(a)
(b)


## 2.4.2. Anterior Chamber Shape Analysis

After calculating the biometric parameters, we performed anterior chamber shape analysis on the selected region of iridocorneal angle image as shown in Figure5. 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 M¯×N¯ 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:(1)
complex shape descriptor to quantify the complexity or changes of iridocorneal angle, fractal dimension (FD);(2)
ellipticity moment shape descriptor (EI);(3)
triangularity moment shape descriptor (TI).Figure 5
Shape analysis on the region of interest in SS-OCT image where SL is Schwalbe’s line and AOD_psl is the angle opening distance 500 μm posterior from SL.


## 2.5. Classification

The ability of each feature extraction method to separate open angle and angle closure glaucoma is quantified by the fuzzyk-nearest neighbour classifier [20]. The basis of the fuzzy k-NN algorithm is to assign membership as a function of the vector’s distance from itsk-nearest neighbors and those neighbors’ memberships in the possible classes.


## 3. Experiments and Results

The Singaporean Chinese population study by the Singapore National Eye Centre recruited 148 subjects (90 females and 45 males) with average age of59.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 [21]. 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.(a) Segmentation of anterior chamber without artifact removal; (b) segmentation of anterior chamber, the cornea, and the anterior of iris surface with artifact removal in SS-OCT image. Green region is extracted anterior chamber and red lines are edges of cornea and anterior of iris.
(a)
(b)Figure7(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 mm2 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 results of segmentation and feature extractions on SS-OCT image; (a) segmentation of anterior chamber; (b) extracted iris image for curvature analysis; (c) the resulting angle profile for biometric parameter measurement; (d) region of interest for fractal dimension analysis; (e) extracted shape for moment shape analysis; and (f) fractal image of extracted region of SS-OCT image.
(a)
(b)
(c)
(d)
(e)
(f)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 Table1. 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).Table 1
Comparison of iridocorneal angle features between open angle and angle closure SS-OCT images.
Features
Open angle images (n=129)
Closed angle images    (n=135)
Mean (SD)
COV
Mean (SD)
COV
FD
1.87 (0.04)
0.021
1.84 (0.06)
0.03
mAOD
0.35 (0.13)
0.371
0.21 (0.11)
0.523
ATsl
0.38 (0.15)
0.395
0.18 (0.04)
0.222
T
I
0.89 (0.09)
0.101
0.85 (0.16)
0.188
E
I
0.15 (0.03)
0.2
0.17 (0.056)
0.329
H
iris
−0.03 (0.04)
1.333
−0.018 (0.039)
2.16
SD: standard deviation, FD: fractal dimension, mAOD: mean of continuous AOD, ATsl: area of the trapezoidal profile of iridocorneal angle, TI: triangularity,  EI: ellipticity, Hiris: mean curvature of iris, and COV: coefficient of variance.(a) Scatter plot of FD features versus mAOD for 264 images and (b) scatter plot of FD features versus ATsl for 264 images.
(a)
(b)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 [22]. 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 [23]. 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. We also explored the effect of combining different kinds of features on classifier performance. 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 fuzzyk-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.Table 2
The association between the gonioscopic grading and the measured iridocorneal angle features.
Gonio
FD
mAOD
ATsl
T
I
E
I
H
iris
Gonio
1
FD
0.517
1
mAOD
0.582
0.843
1
ATsl
0.53
0.84
0.979
T
I
0.152
0.185
0.193
0.277
1
E
I
−0.16
−0.21
−0.212
−0.273
−0.254
1
H
iris
−0.09
−0.021
−0.013
−0.088
−0.065
−0.064
1
Gonio: gonioscopic Shaffer grading.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 parameterm 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 best classification performance was achieved with the m parameter of 1.2.(a) The relationship between the fuzzy strength parameterm and the classification accuracy and (b) the relationship between the fuzzy strength parameter m and the AUC value.
(a)
(b)The average 5-fold CV accuracy and the corresponding standard deviation of fkNN are shown in Table3. The AUC ranges from 0.88±0.02 to 0.98±0.012. The biometric features and the fractal feature are performed with a similar accuracy. 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.Table 3
The average classification accuracy and the standard deviation of fuzzy kNN classifier with various feature groups (m=1.2,  k=5).
Features
Accuracy (%)
Error rate (%)
Sensitivity (%)
Specificity (%)
AUC
Mean
SD
Mean
SD
Mean
SD
Mean
SD
Mean
SD
T
I
,
E
I
,
H
iris
94.69
0.71
5.24
0.71
99.9
0.24
77.71
3.39
0.88
0.02
mAOD, ATsl
98.21
0.82
1.79
0.82
1.00
0.00
92.14
0.38
0.96
0.01
FD
98.93
0.62
1.07
0.62
1.00
0.00
95.43
2.69
0.98
0.013
FD, mAOD, ATsl
99.11
0.76
0.89
0.76
99.9
0.24
96.11
0.24
0.98
0.012
T
I: triangularity, EI: ellipticity, Hiris: mean curvature of iris, mAOD: mean of continuous AOD, FD: fractal dimension, ATsl: area of trapezoidal profile of iridocorneal angle, AUC: area under the curve, and SD: standard deviation.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). All classification algorithms provided a correct classification rate higher than 81%. 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. Minimum differences can be appreciated among LDA, kNN, and SVM classifiers. 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. This is an indication of how important the fuzziness of the membership function is. Additionally, fkNN has been shown to be an efficient and robust classification method to diagnose angle closure glaucoma in SS-OCT images.Table 4
Accuracy of fuzzy kNN versus the accuracy of common classification models.
Classifier
Accuracy (%)
Error rate (%)
Sensitivity (%)
Specificity (%)
AUC
Mean
SD
Mean
SD
Mean
SD
Mean
SD
Mean
SD
LDA
81.09
4.55
18.90
4.55
80.94
4.84
81.59
8.16
0.812
0.05
kNN (k=5)
85.80
5.16
14.19
5.16
90.10
5.77
71.82
5.41
0.80
0.04
SVM
83.26
4.74
16.73
4.74
87.78
7.92
68.40
9.29
0.78
0.03
fkNN (k=5)
99.11
0.76
0.89
0.76
99.99
0.24
96.11
0.24
0.98
0.012The 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. [24]. 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 is strongly recommended. We believe that our findings in this study can serve as a basic grading system for angle closure glaucoma diagnosis in the future.


## 4. Conclusions

We proposed a novel automated angle classification system using SS-OCT images. 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 of99.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.



---

*Source: 942367-2014-08-05.xml*