Joint Optic Disc and Cup Segmentation Based on Multi-label Deep Network and Polar Transformation

01/03/2018 ∙ by Huazhu Fu, et al. ∙ CVTE Agency for Science, Technology and Research 0

Glaucoma is a chronic eye disease that leads to irreversible vision loss. The cup to disc ratio (CDR) plays an important role in the screening and diagnosis of glaucoma. Thus, the accurate and automatic segmentation of optic disc (OD) and optic cup (OC) from fundus images is a fundamental task. Most existing methods segment them separately, and rely on hand-crafted visual feature from fundus images. In this paper, we propose a deep learning architecture, named M-Net, which solves the OD and OC segmentation jointly in a one-stage multi-label system. The proposed M-Net mainly consists of multi-scale input layer, U-shape convolutional network, side-output layer, and multi-label loss function. The multi-scale input layer constructs an image pyramid to achieve multiple level receptive field sizes. The U-shape convolutional network is employed as the main body network structure to learn the rich hierarchical representation, while the side-output layer acts as an early classifier that produces a companion local prediction map for different scale layers. Finally, a multi-label loss function is proposed to generate the final segmentation map. For improving the segmentation performance further, we also introduce the polar transformation, which provides the representation of the original image in the polar coordinate system. The experiments show that our M-Net system achieves state-of-the-art OD and OC segmentation result on ORIGA dataset. Simultaneously, the proposed method also obtains the satisfactory glaucoma screening performances with calculated CDR value on both ORIGA and SCES datasets.

READ FULL TEXT VIEW PDF
POST COMMENT

Comments

There are no comments yet.

Authors

page 1

page 3

page 4

page 7

This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.

I Introduction

Glaucoma is the second leading cause of blindness worldwide (only second to cataracts), as well as the foremost cause of irreversible blindness [1]. Since vision loss from glaucoma cannot be reversed, early screening and detection methods are essential to preserve vision and life quality. One major glaucoma screening technique is optic nerve head (ONH) assessment, which employs a binary classification to identify the glaucomatous and healthy subjects [2]. However, the manual assessment by trained clinicians is time consuming and costly, and not suitable for population screening.

Fig. 1: Structure of the optic nerve head. The region enclosed by the green dotted circle is the optic disc (OD); the central bright zone enclosed by the blue dotted circle is the optic cup (OC); and the region between them is the neuroretinal rim. The vertical cup to disc ratio (CDR) is calculated by the ratio of vertical cup diameter (VCD) to vertical disc diameter (VDD). PPA: Peripapillary Atrophy.

For large-scale screening, automatic ONH assessment methods are needed. Some clinical measurements are proposed, such as the vertical cup to disc ratio (CDR) [3], rim to disc area ratio (RDAR), and disc diameter [4]. In them, CDR is well accepted and commonly used by clinicians. In color fundus image, the optic disc (OD) appears as a bright yellowish elliptical region, and can be divided into two distinct zones: a central bright zone as optic cup (OC) and a peripheral region as the neuroretinal rim, as shown in Fig. 1. The CDR is calculated by the ratio of vertical cup diameter (VCD) to vertical disc diameter (VDD). In general, a larger CDR suggests a higher risk of glaucoma and vice versa. Accurate segmentations of OD and OC are essential for CDR measurement. Some methods automatically measure the disc and cup from 3-D optical coherence tomography (OCT) [5, 6, 7, 8]. However, OCT is not easily available due to its high cost, the fundus image is still referred to by most clinicians. A number of works have been proposed to segment the OD and/or OC from the fundus image [9, 10, 11, 12]. The main segmentation techniques include color and contrast thresholding, boundary detection, and region segmentation methods [12]. In these methods, the pixels or patches of fundus images are determined as background, disc and cup regions, through a learned classifier with various visual features. However, most existing methods are based on hand-crafted features (e.g., RGB color, texture, Gabor filter, and gradient), which lack sufficiently discriminative representations and are easily affected by pathological regions and low contrast quality. In addition, most methods segment the OD and OC separately, i.e., segmenting OD first, followed by the OC without considering the mutual relation of them. In this paper, we consider OD and OC together, and provide a one-stage framework based on deep learning technique.

Deep learning techniques have been recently demonstrated to yield highly discriminative representations that have aided in many computer vision tasks. For example, Convolutional Neural Networks (CNNs) have brought heightened performance in image classification 

[13] and segmentation [14]. For retinal image, Gulshan et al. have demonstrated that the deep learning system could obtain a high sensitivity and specificity for detecting referable diabetic retinopathy [15]. In fundus vessel segmentation, the deep learning systems [16, 17, 18] also achieve the state-of-the-art performances. These successes have motivated our investigation of deep learning for disc and cup segmentation from fundus images.

In our paper, we address OD and OC segmentation as a multi-label task and solve it using a novel end-to-end deep network. The main contributions of our work include:

  1. We propose a fully automatic method for joint OD and OC segmentation using a multi-label deep network, named M-Net. Our M-Net is an end-to-end deep learning system, which contains a multi-scale U-shape convolutional network with the side-output layer to learn discriminative representations and produces segmentation probability map.

  2. For joint OD and OC segmentation, a multi-label loss function based on Dice coefficient is proposed, which deals well with the multi-label and imbalance data of pixel-wise segmentation for fundus image.

  3. Moreover, a polar transformation is utilized in our method to transfer the fundus image into a polar coordinate system, which introduces the advantages of spatial constraint, equivalent augmentation, and balancing cup proportion, and improves the segmentation performance.

  4. At last, we evaluate the effectiveness and generalization capability of the proposed M-Net on ORIGA dataset. Our M-Net achieves state-of-the-art segmentation performance, with the average overlapping error of and for OD and OC segmentation, respectively.

  5. Furthermore, the CDR is calculated based on segmented OD and OC for glaucoma screening. Our proposed method obtains highest performances with areas under curve (AUC) of and on ORIGA and SCES datasets.

The remainders of this paper are organized as follows. We begin by reviewing techniques related to OD/OC segmentation in Section II. The details of our system and its components are presented in Section III. To verify the efficacy of our method, extensive experiments are conducted in Section IV, and then we conclude with final remarks in Section V.

Ii Related Works

Optic Disc Segmentation: The OD is the location where ganglion cell axons exit the eye to form the optic nerve through which visual information of the photo-receptors is transmitted to the brain. Earlier, the template based methods are proposed firstly to obtain the OD boundary. For example, Lowell et al. employed the active contour model [19] to detect the contour based on image gradient. In [20, 21], the Circular-based Transformation techniques are employed to obtain the OD boundary. In [9], the local texture features around each point of interest in multidimensional feature space are utilized to provide robustness against variations in OD region. Recently, the pixel classification based method is proposed to transfer the boundary detection problem into a pixel classification task, which obtains a satisfactory performance. Cheng et al. [10] utilizes superpixel classifier to segment the OD and OC, which exploits using various hand-crafted visual features at superpixel level to enhance the detection accuracy. In [22], the disparity values extracted from stereo image pairs are introduced to distinguish the OD and background. However, reliance on hand-crafted features make these methods susceptible to low quality images and pathological regions.

Optic Cup Segmentation: The OC is restricted to the region inside OD. Segmenting OC from fundus images is a more challenging task due to the low contrast boundary. In [23], an automatic OC segmentation algorithm based on a variational level set is proposed. Later, the blood vessel kinks are found to be useful for OC segmentation [24], and a similar concept but named as vessel bend is utilized in [9]. The main challenge in detecting kinks or vessel bending is that it is often affected by natural vessel bending that does not lie on the OC boundary. Moreover, the pixel classification based methods similar to OD segmentation [10] is also introduced to OC segmentation. The various hand-crafted visual features (e.g., center surround statistics, color histogram, and low-rank superpixel representation) are employed in [25, 10, 26] to represent the pixel/superpixel for OC segmentation. A common limitation of these algorithms is that they highly relied on hand-crafted visual feature, which are mainly based on contrast between the neuroretinal rim and the cup.

Fig. 2: Illustration of our segmentation framework, which mainly includes fundus polar transformation and M-Net segmentation. Firstly, the optic disc is localized, and a polar transformation generates the representation of original fundus image in the polar coordinate system based on the detected disc center. Than our M-Net produces the multi-label prediction maps for disc and cup regions. Our M-Net architecture consists of multi-scale input layer, U-shape convolutional network, side-output layer, and multi-label loss function. The (De)Convolutional layer parameters are denoted as “(De)Conv

kernel size, stride

”. Finally, the inverse polar transformation recovers the segmentation result back to the Cartesian coordinate.

Joint OD and OC Segmentation: Most existing methods only focus on the single region segmentation (i.e., OC or OD). Especially, for the cup segmentation, the OD boundary could provide some useful prior informations, e.g., shape constraint and structure constraint [27]. The works in [9, 10] deal with the OD and OC by two separate stages with different features. Zheng et al. integrated the OD and OC segmentation within a graph-cut framework [11]. However, they consider the OD and OC as two mutual labels, which means for any pixel in fundus, it can only belong to one label (i.e., background, OD, and OC). Moreover, the method [11]

only employs color feature within a Gaussian Mixture Model to decide a posterior probability of the pixel, which makes it unsuitable for fundus image with low contrast. In 

[28], a modified U-Net deep network is introduced to segment the OD and OC. However, it still separates OD and OC segmentation in a sequential way. In [29], an ensemble learning method is proposed to extract OC and OD based on the CNN architecture. An entropy sampling technique is used to select informative points, and then graph cut algorithm is employed to obtain the final segmentation result. However this multiple step deep system limits its effectiveness in the training phase.

Iii Proposed Method

Fig. 2 illustrates the overall flowchart of our OD and OC segmentation method, which contains M-Net deep network and the fundus image polar transformation. In our method, we firstly localize the disc center by using the existing automatic disc detection method [30], and then transfers the original fundus image into polar coordinate system based on the detected disc center. Then the transferred image is fed into our M-Net, and generates the multi-label probability maps for OD and OC regions. Finally, the inverse polar transformation recovers the segmentation map back to the Cartesian coordinate.

Iii-a M-Net Architecture

Our M-Net is an end-to-end multi-label deep network, which consists of four main parts. The first is a multi-scale layer used to construct an image pyramid input and achieve multi-level receptive field fusion. The second is a U-shape convolutional network, which is employed as the main body structure to learn a rich hierarchical representation. The third part is side-output layer that works on the early convolutional layers to support deep layer supervision. Finally, a multi-label loss function is proposed to guarantee OD and OC segmentation jointly.

Iii-A1 U-shape Convolutional Network

In our paper, we modify the U-shape convolutional network (U-Net) in [31]

as the main body in our deep architecture. U-Net is an efficient fully convolutional neural network for the biomedical image segmentation. Similar to the original U-Net architecture, our method consists of the encoder path (left side) and decoder path (right side). Each encoder path performs convolution layer with a filter bank to produce a set of encoder feature maps, and the element-wise rectified-linear non-linearity (ReLU) activation function is utilized. The decoder path also utilizes the convolution layer to output decoder feature map. The skip connections transfer the corresponding feature map from encoder path and concatenate them to up-sampled decoder feature maps.

Finally, the high dimensional feature representation at the output of the final decoder layer is fed to a trainable multi-label classifier. In our method, the final classifier utilizes convolutional layer with Sigmoid activation as the pixel-wise classification to produce the probability map. For multi-label segmentation, the output is a channel probability map, where is the class number ( for OD and OC in our work). The predicted probability map corresponds to the class with maximum probability at each pixel.

Iii-A2 Multi-scale Input Layer

The multi-scale input or image pyramid has been demonstrated to improve the quality of segmentation effectively. Different from the other works, which fed the multi-scale images to multi-scream networks separately and combine the final output map in the last layer [32, 33], our M-Net employs the average pooling layer to downsample the image naturally and construct a multi-scale input in the encoder path. Our multi-scale input layer has following advantages: 1) integrating multi-scale inputs into the decoder layers to avoid the large growth of parameters; 2) increasing the network width of decoder path.

Iii-A3 Side-output Layer

In our M-Net, we also introduce the side-output layer, which acts as a classifier that produces a companion local output map for early layers [34]. Let denote the parameters of all the standard convolutional layers, and there are side-output layers in the network, where the corresponding weights are denoted as . The objective function of the side-output layer is given as:

(1)

where is the loss function fusion-weight for each side-output layer ( in our paper), is the side-output number, and denotes the multi-label loss of the -th side-output layer. To directly utilize side-output prediction map, we employ an average layer to combine all side-output maps as the final prediction map. The main advantages of the side-output layer are: first, the side-output layer back-propagates the side-output loss to the early layer in the decoder path with the final layer loss, which could relieve gradient vanishing problem and help the early layer training. It can be treated as a special bridge link between the loss and early layer; second, the multi-scale fusion has been demonstrated to achieve a high performance, and the side-output layer supervises the output map of each scale to output the better result.

Iii-A4 Multi-label Loss Function

In our work, we formulate the OD and OC segmentation as a multi-label problem. The existing segmentation methods usually belong to the multi-class setting, which assign each instance to one unique label of multiple classes. By contrast, multi-label method learns an independent binary classifier for each class, and assigns each instance to multiple binary labels. Especially for OD and OC segmentation, the disc region overlays the cup pixels, which means the pixel marked as cup also has the label as disc. Moreover, for the glaucoma cases, the disc pixels excluded cup region shapes as a thin ring, which makes the disc label extremely imbalance to background label under the multi-class setting. Thus, multi-label method, considering OD and OC as two independent binary classifiers, is more suitable for addressing these issues. In our method, we propose a multi-label loss function based on Dice coefficient. The Dice coefficient is a measure of overlap widely used to assess segmentation performance when the ground truth is available [35]. Our multi-label loss function is defined as:

(2)

where is the pixel number, and denote predicted probability and binary ground truth label for class , respectively. is the class number, and are the class weights. Our multi-label loss function in Eq. (2) is equivalent to the traditional Dice coefficient by setting . In our method, we set for OD and OC segmentation. Note that the Dice loss function indicates the foreground mask overlapping ratio, and can deals with the imbalance issue in the pixels of foreground (i.e., OD or OC) region and background. Under our multi-label setting, the pixel can be labeled as OD or/and OC independently. Thus, the imbalance issue does not exist between OD and OC. in Eq. (2) is the trade-off weight to control the contribution of OD and OC. For glaucoma screening, both the OD and OC are important, thus we set . Our multi-label loss function can be differentiated yielding the gradient as:

(3)

This loss is efficiently integrated into back-propagation via standard stochastic gradient descent.

Iii-B Polar Transformation for Fundus Image

Fig. 3: Illustration of the mapping from Cartesian coordinate system (A) to the polar coordinate system (C) by using the polar transformation. The point in Cartesian coordinate corresponds to the point in polar coordinate. (B) and (D) are the corresponding ground truth, where yellow, red, and black regions denote the optic cup, optic disc, and background, respectively.

In our method, we introduce a polar transformation for improving the OD and OC segmentation performance. The pixel-wise polar transformation transfers the original fundus image to the polar coordinate system. Let denotes the point on fundus image plane, where the origin is set as the disc center , and is the Cartesian coordinates, as shown in Fig. 3 (A). The corresponding point on polar coordinate system is , as shown in Fig. 3 (C), where and are the radius and directional angle of the original point , respectively. The transform relation between the polar and Cartesian coordinates is as follow:

(4)

The height and width of transferred polar image are transformation radius and discretization . The polar transformation provides a pixel-wise representation of the original image in the polar coordinate system, which has the following properties:
1) Spatial Constraint: In the original fundus image, a useful geometry constraint is that the OC should be within the OD region, as shown in Fig. 3 (B). But this redial relationship is difficult to implement in the original Cartesian coordinate. By contrast, our polar transformation transfers this redial relationship to a spatial relationship, where the regions of cup, disc, and background appear the ordered layer structure, as shown in Fig. 3 (D). This layer-like spatial structure is convenient to use, especially some layer-based segmentation methods [36, 37] can be employed as the post-processing.
2) Equivalent Augmentation: Since the polar transformation is a pixel-wise mapping, the data augmentation on original fundus image is equivalent to that on polar coordinate. For example, moving the expansion center is equivalent to the drift cropping transformations on polar coordinate. Using different transformation radius is same as augmenting with the various scaling factor. Thus the data augmentation for deep learning can be done during the polar transformation with various parameters.
3) Balancing Cup Proportion: In the original fundus image, the distribution of OC/background pixels is heavily biased. Even in the cropped ROI, the cup region still accounts for a low proportion. Using Fig. 3 (B) as an example, the cup region only occupies about

. This extremely imbalance proportion easily leads the bias and overfitting in training the deep model. Our polar transformation flat the image based on OD center, that could enlarge the cup region by using interpolation and increase the OC proportion. As shown in Fig. 

3 (D), the ratio of cup region increases to over the ROI, which is more balanced than that in original fundus image. The balanced regions help avoid the overfitting during the model training and improve the segmentation performance further.

Note that the method in [38]

also utilizes the polar transformation to detect the cup outline based on the depth map estimated from stereo retinal fundus image pairs. Our work has significant difference compared to 

[38]. 1) Motivations are different. The polar transformation used in [38] aims at finding strongest depth edge in the radial direction as initial candidate points of cup border. In our work, we use polar transformation to obtain spatial constraint and augments the cup/disc region proportion. 2) Methods are different. The method in [38] detects the OD and OC boundaries sequentially, and polar transformation only use for OC segmentation. Our method segments OD and OC regions jointly, and considers their mutual relation under polar coordinate.

Iv Experiments

Iv-a Implementation

Our M-Net is implemented with Python based on Keras with Tensorflow backend. During training, we employ stochastic gradient descent (SGD) for optimizing the deep model. We use a gradually decreasing learning rate starting from

and a momentum of . red The transformation radius is set to , and the directional angles are draw into distinct bins, thus the size of transferred polar image is . The output of our M-Net is a 2-channel posterior probability map for OD and OC, where each pixel value represents the probability. A fixed threshold is employed to get a binary mask from the probability map. As same as that in previous works [10, 39], the largest connected region in OD/OC mask is selected and the ellipse fitting is utilized to generate the final segmentation result.

Iv-B Segmentation Experiments

We first evaluate the OD and OC segmentation performance. We employ the ORIGA dataset [40] containing 650 fundus images with 168 glaucomatous eyes and 482 normal eyes. The 650 images with manual ground truth boundaries are divided into 325 training images (including 73 glaucoma cases) and 325 testing images (including 95 glaucoma) as same as that in [26, 41]. To evaluate the segmentation performance, we use the overlapping error () and balanced accuracy (

) as the evaluation metrics for OD, OC, and rim regions:

(5)

with

(6)

where and denote the segmented mask and the manual ground truth, respectively. and denote the number of true positives and true negatives, respectively, and and denote the number of false positives and false negatives, respectively. Moreover, we follow the clinical convention to compute the vertical cup to disc ratio (CDR) which is an important indicator for glaucoma screening. When CDR is greater than a threshold, it is glaucomatous, otherwise, healthy. Thus an evaluation metric, named absolute CDR error , is defined as: where denotes the manual CDR from trained clinician, and is the CDR calculated by the segmented result.

(a) ORIGA dataset
(b) SCES dataset
Fig. 4: The ROC curves with AUC scores for glaucoma screening based on the vertical cup to disc ratio (CDR) on (A) ORIGA and (B) SCES datasets.
Method
R-Bend [9] 0.129 - 0.395 - - - 0.154
ASM [42] 0.148 - 0.313 - - - 0.107
Superpixel [10] 0.102 0.964 0.264 0.918 0.299 0.905 0.077
LRR [26] - - 0.244 - - - 0.078
QDSVM [39] 0.110 - - - - - -
U-Net [31] 0.115 0.959 0.287 0.901 0.303 0.921 0.102
Joint U-Net 0.108 0.961 0.285 0.913 0.325 0.903 0.083
Our M-Net 0.083 0.972 0.256 0.914 0.265 0.921 0.078
Joint U-Net + PT 0.072 0.979 0.251 0.914 0.250 0.935 0.074
Our M-Net + PT 0.071 0.983 0.230 0.930 0.233 0.941 0.071

TABLE I: Performance comparisons (%) of the different methods on ORIGA Dataset. (PT: Polar Transformation.)

We compare our M-Net with the several state-of-the-art OD/OC segmentation approaches: relevant-vessel bends (R-Bend) method in [9], active shape model (ASM) method in [42], superpixel-based classification (Superpixel) method in [10], quadratic divergence regularized SVM (QDSVM) method in [39], and low-rank superpixel representation (LRR) method in [26]. Additional, we compare with the deep learning method U-Net [31]. We report two results of U-Net, one is the original U-Net for segmenting OC and OD separately, and the other is U-Net utilized our multi-label loss function (Joint U-Net) for segmenting OC and OD jointly. We also provide segmentation results with/without the polar transformation (PT). The performances are shown in Table I.

R-Bend [9] provides a parameterization technique based on vessel bends, and ASM [42] employs the circular Hough transform initialization to segment the OD and OC regions. These two bottom-up methods extract the OD and OC regions separately, and do not perform well on the ORIGA dataset. Superpixel method in [10] utilizes superpixel classification to detect the OD and OC boundaries. It obtains a better performance than the other two bottom-up methods [9, 42]. The methods in LRR [26] and QDSVM [39] obtain good results. However, they only focus on either OD or OC segmentation, and can not calculate the CDR for glaucoma screening. Joint U-Net with our multi-label loss utilizes the mutual relation of OD and OC, and obtains a better performance than that in traditional U-Net [31]. Our M-Net with multi-scale input and side-output layers achieves a higher score than single-scale network and superpixel method [10]. It demonstrates that the multi-scale input and side-output layers are useful to guide the early layer training.

The polar transformation as one contribution of our work augments the proportion of cup region. One major advantage is that the polar transformation augments the proportion of cup region, and makes the area of the disc/cup and background more balance. The balanced regions help avoid the overfitting during the model training and improve the segmentation performance further. From Table I, polar transformation reduces about in Joint U-Net and in M-Net on scores. Note that the performance of Joint U-Net with PT is slight better than that in M-Net without PT. It shows that the gains of the polar transformation may be higher than that using multi-scale input and side-output layers. Finally, our M-Net with PT achieves the best performance, and outperforms other state-of-the-art methods.

Fig. 5: The visual examples of optic disc and cup segmentation, where the yellow and red region denote the cup and disc segmentations, respectively. From the left to right: fundus image, ground truth (GT), Joint U-Net, our M-Net and M-Net with polar transformation (PT). The last row shows the failed case.

Fig. 5 shows the visual examples of the segmentation results, where the first two rows are normal eyes and the rest rows are glaucoma cases. For the superpixel method [10], the segmented OC is smaller than ground truth in glaucoma cases, which may cause an under-estimated CDR. The deep learning methods (e.g., Joint U-Net and M-Net) obtain more accurate cup boundary, but it easily generates a larger OD. By contrast, our M-Net with PT can effectively and accurately segment OD and OC regions. The last row in Fig. 5 shows a challenging case for segmentation, where the image is blurred and has low-contrast for identifying the OC boundary. For this case, all the methods fail to produce accurate OC segmentation. This issue could potentially be addressed in future work through the use of more powerful network or additional image enhancement pre-processing.

Iv-C Glaucoma Screening

We also evaluate the proposed method on glaucoma screening by using the calculated CDR value. Two datasets are used, one is ORIGA dataset, the second is Singapore Chinese Eye Study (SCES) dataset. For ORIGA dataset, we employ 325 image for training and rest for testing, as same as used in segmentation experiment. For the SCES dataset, it consists of 1676 images, of which 46 () are glaucoma cases. Since the SCES dataset provides only clinical diagnoses, it will be used only to assess the diagnostic performance of our system. We use all the 650 images in ORIGA dataset for training and the whole 1676 images of SCES for testing. We report the Receiver Operating Characteristic (ROC) curve and area under ROC curve (AUC) as the overall measure of the diagnostic strength. The performances for glaucoma screening based on CDR are shown in Fig. 4.

From the glaucoma screening results, we have the following observations: 1) The non-deep learning method, superpixel [10], produces a competitive performance () on ORIGA dataset, which is better than M-Net (). But its performance is lower than others on SCES dataset. 2) The Joint U-Net with PT obtains higher scores than superpixel [10] and U-Net on both ORIGA and SCES datasets. 3) Our M-Net with PT achieves the best performances on ORIGA dataset () and SCES dataset (). Especially, our M-Net with PT has more than improvement on AUC than M-Net without PT, which demonstrates the effectiveness of polar transformation on glaucoma screening. 4) Our method also outperforms other deep learning based diagnostic method. For example, the deep learning method in [43] provides a glaucoma screening system by using deep visual feature, which obtained and on ORIGA and SCES dataset. However, it can not provide the CDR value as a clinical explanation. Our result of M-Net with PT is comparable to that of deep system [43]. 5) Final, all the deep learning based methods have better performances in SCES dataset than those in ORIGA dataset. One possible reason is that the size of training set on ORIGA is only 325. The more training data promotes the representation capability of deep learning.

Iv-D Discussion

Iv-D1 Running Time

The entire training phase of our method takes about 5 hours on a single NVIDIA Titan X GPU (100 iterations). However, the training phase could be done offline. In online testing, it costs only to generate the final segmentation map for one fundus image, which is faster than the existing methods, e.g., superpixel method [10] takes , ASM method [42] takes , R-Bend method [9] takes , sequential segmentation of OD and OC by using original U-Net [31] takes .

Iv-D2 Repeatability Experiment

Data Glaucoma (n=39) Normal (n=1481) All (n=1520)
Coefficient 0.8833 0.8262 0.8357
-value

TABLE II: Correlation Coefficients of Repeatability Experiment.
Fig. 6: Scatter plot of the CDR correspondence on the repeatability dataset.

In this experiment, we evaluate the repeatability of proposed method. We collect a repeatability dataset with two corresponding sets (A and B) consisting of 1520 fundus image pairs. For each pair, one image is selected from SCES dataset, and other one is the different image from the same visit. We run our proposed method on these two sets, and calculate the correlation coefficients of CDR values. Table II reports the repeatability test result, and the scatter plot of the CDR correspondence is shown in Fig. 6. As it can be seen, our method gets of -value and appears a good repeatability.

Iv-D3 Clinical Measurement

Our M-Net method segments the whole OD and OC regions, which could be used to calculate other clinical measurements. In this experiment, we evaluate the rim to disc area ratio (RDAR) [44] defined as: . The comparison of CDR and RDAR is shown in Table III, where our M-Net with PT obtains the best screening performance based on RDAR value on both datasets, which is consistent with the experiment based on CDR. Moreover, CDR measurement shows a better screening performance than RDAR. The possible reasons may be that the rim is calculated by subtracting of disc and cup regions, which contains the errors of both disc and cup segmentations. Thus the rim error is larger than cup error. This also observes from Table I, where the rim error () is larger than cup error () based on our M-Net with PT. Moreover, since the central retinal vessel trunk usually locates in the nasal optic disc sector [45], it renders difficult the automatic delineation of the optic disc region in horizontal. Thus, the vertical disc and cup may obtain the higher accuracy than that in horizontal.

ORIGA Dataset SCES Dataset
Method
Our M-Net 0.8019 0.7981 0.8397 0.8290
Joint U-Net + PT 0.8152 0.7921 0.8612 0.8003
Our M-Net + PT 0.8508 0.8425 0.8998 0.8488

TABLE III: AUC performance of the different Clinical Measurements. (CDR: vertical cup to disc ratio. RDAR: rim to disc area ratio)

V Conclusion

In this paper, we have developed a deep learning architecture, named M-Net, which solves the OD and OC segmentation jointly into a one-stage multi-label framework. The proposed M-Net employed the U-shape convolutional network as the body structure. And the multi-scale layer constructed an image pyramid to fed a multi-level inputs, while the side-output Layers acted as a early classifier to produce the companion local prediction maps for early scale layers. A multi-label loss function has been proposed to guarantee the final output for segmenting OD and OC together. For improving the segmentation result further, we also have introduced a polar transformation to transfer the original fundus image to the polar coordinate system. We have demonstrated that our system produces state-of-the-art segmentation results on ORIGA dataset. Simultaneously, the proposed method also obtained the satisfactory glaucoma screening performances by using the calculated CDR on both ORIGA and SCES datasets. The work implementation details are available at http://hzfu.github.io/proj_glaucoma_fundus.html.

References

  • [1] Y.-C. Tham, X. Li, T. Y. Wong, H. A. Quigley, T. Aung, and C.-Y. Cheng, “Global prevalence of glaucoma and projections of glaucoma burden through 2040: A systematic review and meta-analysis,” Ophthalmology, vol. 121, no. 11, pp. 2081–2090, 2014.
  • [2] D. F. Garway-Heath and R. A. Hitchings, “Quantitative evaluation of the optic nerve head in early glaucoma,” Br. J. Ophthalmol., vol. 82, no. 4, pp. 352–361, 1998.
  • [3] J. B. Jonas, A. Bergua, P. Schmitz–Valckenberg, K. I. Papastathopoulos, and W. M. Budde, “Ranking of optic disc variables for detection of glaucomatous optic nerve damage,” Invest. Ophthalmol. Vis. Sci., vol. 41, no. 7, pp. 1764–1773, 2000.
  • [4] M. D. HancoxO.D., “Optic disc size, an important consideration in the glaucoma evaluation,” Clinical Eye and Vision Care, vol. 11, no. 2, pp. 59 – 62, 1999.
  • [5] K. Lee, M. Niemeijer, M. K. Garvin, Y. H. Kwon, M. Sonka, and M. D. Abràmoff, “Segmentation of the optic disc in 3-D OCT scans of the optic nerve head,” IEEE Trans. Med. Imag., vol. 29, no. 1, pp. 159–168, 2010.
  • [6] M. Wu, T. Leng, L. de Sisternes, D. L. Rubin, and Q. Chen, “Automated segmentation of optic disc in sd-oct images and cup-to-disc ratios quantification by patch searching-based neural canal opening detection,” Opt. Express, vol. 23, no. 24, pp. 31 216–31 229, Nov 2015.
  • [7] H. Fu, D. Xu, S. Lin, D. W. K. Wong, and J. Liu, “Automatic Optic Disc Detection in OCT Slices via Low-Rank Reconstruction,” IEEE Trans. Biomed. Eng., vol. 62, no. 4, pp. 1151–1158, 2015.
  • [8] H. Fu, Y. Xu, S. Lin, X. Zhang, D. Wong, J. Liu, and A. Frangi, “Segmentation and Quantification for Angle-Closure Glaucoma Assessment in Anterior Segment OCT,” IEEE Trans. Med. Imag., vol. 36, no. 9, pp. 1930–1938, 2017.
  • [9] G. D. Joshi, J. Sivaswamy, and S. R. Krishnadas, “Optic Disk and Cup Segmentation from Monocular Colour Retinal Images for Glaucoma Assessment,” IEEE Trans. Med. Imag., vol. 30, no. 6, pp. 1192–1205, 2011.
  • [10] J. Cheng, J. Liu, Y. Xu, F. Yin, D. Wong, N. Tan, D. Tao, C.-Y. Cheng, T. Aung, and T. Wong, “Superpixel classification based optic disc and optic cup segmentation for glaucoma screening,” IEEE Trans. Med. Imag., vol. 32, no. 6, pp. 1019–1032, 2013.
  • [11] Y. Zheng, D. Stambolian, J. O’Brien, and J. Gee, “Optic Disc and Cup Segmentation from Color Fundus Photograph Using Graph Cut with Priors,” in Proc. MICCAI, 2013, pp. 75–82.
  • [12] A. Almazroa, R. Burman, K. Raahemifar, and V. Lakshminarayanan, “Optic Disc and Optic Cup Segmentation Methodologies for Glaucoma Image Detection: A Survey,” Journal of Ophthalmology, vol. 2015, 2015.
  • [13]

    A. Krizhevsky, I. Sutskever, and G. Hinton, “Imagenet classification with deep convolutional neural networks,” in

    Proc. NIPS, 2012, pp. 1097–1105.
  • [14] E. Shelhamer, J. Long, and T. Darrell, “Fully convolutional networks for semantic segmentation,” IEEE Trans. Pattern Anal. Mach. Intell., vol. 39, no. 4, pp. 640–651, 2017.
  • [15] V. Gulshan, L. Peng, M. Coram, and et.al, “Development and Validation of a Deep Learning Algorithm for Detection of Diabetic Retinopathy in Retinal Fundus Photographs,” Journal of the American Medical Association, vol. 304, no. 6, pp. 649–656, 2016.
  • [16] H. Fu, Y. Xu, D. Wong, and J. Liu, “Retinal Vessel Segmentation via Deep Learning Network and Fully-connected Conditional Random Fields,” in Proc. ISBI, 2016, pp. 698–701.
  • [17] H. Fu, Y. Xu, S. Lin, D. W. K. Wong, and J. Liu, “DeepVessel: Retinal Vessel Segmentation via Deep Learning and Conditional Random Field,” in Proc. MICCAI, 2016, pp. 132–139.
  • [18] K. Maninis, J. Pont-Tuset, P. Arbeláez, and L. V. Gool, “Deep Retinal Image Understanding,” in Proc. MICCAI, 2016, pp. 140–148.
  • [19] J. Lowell, A. Hunter, D. Steel, A. Basu, R. Ryder, E. Fletcher, and L. Kennedy, “Optic nerve head segmentation,” IEEE Trans. Med. Imag., vol. 23, no. 2, pp. 256–264, 2004.
  • [20]

    A. Aquino, M. E. Gegundez-Arias, and D. Marin, “Detecting the optic disc boundary in digital fundus images using morphological, edge detection, and feature extraction techniques,”

    IEEE Trans. Med. Imag., vol. 29, no. 11, pp. 1860–1869, 2010.
  • [21] S. Lu, “Accurate and efficient optic disc detection and segmentation by a circular transformation,” IEEE Trans. Med. Imag., vol. 30, no. 12, pp. 2126–2133, dec 2011.
  • [22] M. D. Abràmoff, W. L. M. Alward, E. C. Greenlee, L. Shuba, C. Y. Kim, J. H. Fingert, and Y. H. Kwon, “Automated segmentation of the optic disc from stereo color photographs using physiologically plausible features,” Invest. Ophthalmol. Vis. Sci., vol. 48, no. 4, p. 1665, 2007.
  • [23] D. W. K. Wong, J. Liu, J. H. Lim, X. Jia, F. Yin, H. Li, and T. Y. Wong, “Level-set based automatic cup-to-disc ratio determination using retinal fundus images in argali,” in Proc. EMBC, 2008, pp. 2266–2269.
  • [24] D. W. K. Wong, J. Liu, J. H. Lim, H. Li, and T. Y. Wong, “Automated detection of kinks from blood vessels for optic cup segmentation in retinal images,” in Proc. SPIE, vol. 7260, 2009, pp. 7260 – 7268.
  • [25] Y. Xu, D. Xu, S. Lin, J. Liu, J. Cheng, C. Y. Cheung, T. Aung, and T. Y. Wong, “Sliding window and regression based cup detection in digital fundus images for glaucoma diagnosis,” in Proc. MICCAI, 2011.
  • [26] Y. Xu, L. Duan, S. Lin, X. Chen, D. Wong, T. Wong, and J. Liu, “Optic Cup Segmentation for Glaucoma Detection Using Low-Rank Superpixel Representation,” in Proc. MICCAI, 2014.
  • [27] Y. Xu, J. Liu, S. Lin, D. Xu, C. Y. Cheung, T. Aung, and T. Y. Wong, “Efficient optic cup detection from intra-image learning with retinal structure priors,” in Proc. MICCAI, vol. 15, 2012, pp. 58–65.
  • [28] A. Sevastopolsky, “Optic disc and cup segmentation methods for glaucoma detection with modification of U-Net convolutional neural network,” Pattern Recognition and Image Analysis, vol. 27, no. 3, pp. 618–624, 2017.
  • [29] J. Zilly, J. M. Buhmann, and D. Mahapatra, “Glaucoma detection using entropy sampling and ensemble learning for automatic optic cup and disc segmentation,” Comput. Med. Imaging Graph., vol. 55, pp. 28–41, 2017.
  • [30] J. Xu, O. Chutatape, E. Sung, C. Zheng, and P. C. T. Kuan, “Optic disk feature extraction via modified deformable model technique for glaucoma analysis,” Pattern Recognit., vol. 40, no. 7, pp. 2063–2076, 2007.
  • [31] O. Ronneberger, P. Fischer, and T. Brox, “U-Net: Convolutional Networks for Biomedical Image Segmentation,” in Proc. MICCAI, 2015, pp. 234–241.
  • [32] G. Li and Y. Yu, “Visual saliency detection based on multiscale deep cnn features,” IEEE Trans. Image Process., vol. 25, no. 11, pp. 5012–5024, 2016.
  • [33] Y. Liu, M.-M. Cheng, X. Hu, K. Wang, and X. Bai, “Richer Convolutional Features for Edge Detection,” in Proc. CVPR, 2017.
  • [34] C. Lee, S. Xie, P. Gallagher, Z. Zhang, and Z. Tu, “Deeply-supervised nets,” in

    International Conference on Artificial Intelligence and Statistics

    , 2015.
  • [35] W. R. Crum, O. Camara, and D. L. G. Hill, “Generalized Overlap Measures for Evaluation and Validation in Medical Image Analysis,” IEEE Trans. Med. Imag., vol. 25, no. 11, pp. 1451–1461, 2006.
  • [36] A. Lang, A. Carass, M. Hauser, E. S. Sotirchos, P. A. Calabresi, H. S. Ying, and J. L. Prince, “Retinal layer segmentation of macular oct images using boundary classification,” Biomed. Opt. Express, vol. 4, no. 7, pp. 1133–1152, Jul 2013.
  • [37] P. a. Dufour, L. Ceklic, H. Abdillahi, S. Schröder, S. De Dzanet, U. Wolf-Schnurrbusch, and J. Kowal, “Graph-Based Multi-Surface Segmentation of OCT Data Using Trained Hard and Soft Constraints,” IEEE Trans. Med. Imag., vol. 32, no. 3, pp. 531–543, mar 2013.
  • [38] C. Muramatsu, T. Nakagawa, A. Sawada, Y. Hatanaka, T. Hara, T. Yamamoto, and H. Fujita, “Determination of cup-to-disc ratio of optical nerve head for diagnosis of glaucoma on stereo retinal fundus image pairs,” in Proc. SPIE, vol. 7260, 2009, p. 72603L.
  • [39] J. Cheng, D. Tao, D. W. K. Wong, and J. Liu, “Quadratic divergence regularized SVM for optic disc segmentation,” Biomed. Opt. Express, vol. 8, no. 5, pp. 2687–2696, 2017.
  • [40] Z. Zhang, F. Yin, J. Liu, W. Wong, N. Tan, B. Lee, J. Cheng, and T. Wong, “ORIGA(-light): an online retinal fundus image database for glaucoma analysis and research.” in Proc. EMBC, 2010, pp. 3065–3068.
  • [41] J. Cheng, Z. Zhang, D. Tao, D. Wong, J. Liu, M. Baskaran, T. Aung, and T. Wong, “Similarity regularized sparse group lasso for cup to disc ratio computation,” Biomed. Opt. Express, vol. 8, no. 8, pp. 1192–1205, 2017.
  • [42] F. Yin, J. Liu, S. H. Ong, Y. Sun, D. W. K. Wong, N. M. Tan, C. Cheung, M. Baskaran, T. Aung, and T. Y. Wong, “Model-based optic nerve head segmentation on retinal fundus images,” in Proc. EMBC, 2011, pp. 2626–2629.
  • [43] X. Chen, Y. Xu, S. Yan, D. Wing, T. Wong, and J. Liu, “Automatic Feature Learning for Glaucoma Detection Based on Deep Learning,” in Proc. MICCAI, 2015, pp. 669–677.
  • [44] J. B. Jonas, C. Y. Mardin, and A. E. Gründler, “Comparison of measurements of neuroretinal rim area between confocal laser scanning tomography and planimetry of photographs,” Br. J. Ophthalmol., vol. 82, no. 4, pp. 362–366, 1998.
  • [45] J. B. Jonas and W. M. Budde, “Is the nasal optic disc sector important for morphometric glaucoma diagnosis?” Br. J. Ophthalmol., vol. 86, no. 11, pp. 1232–1235, 2002.