Development of a Semi-Automatic Segmentation Method for Retinal OCT Images Tested in Patients with Diabetic Macular Edema

Transcription

1 Development of a Semi-Automatic Segmentation Method for Retinal OCT Images Tested in Patients with Diabetic Macular Edema Yijun Huang 1,2 *, Ronald P. Danis 1,2, Jeong W. Pak 1, Shiyu Luo 1, James White 1, Xian Zhang 3, Ashwini Narkar 1, Amitha Domalpally 1 1 Fundus Photograph Reading Center, Department of Ophthalmology and Visual Sciences, University of Wisconsin-Madison, Madison, Wisconsin, United States of America, 2 McPherson Eye Research Institute, University of Wisconsin-Madison, Madison, Wisconsin, United States of America, 3 Department of Psychiatry, Yale University, New Haven, Connecticut, United States of America Abstract Purpose: To develop EdgeSelect, a semi-automatic method for the segmentation of retinal layers in spectral domain optical coherence tomography images, and to compare the segmentation results with a manual method. Methods: SD-OCT (Heidelberg Spectralis) scans of 28 eyes (24 patients with diabetic macular edema and 4 normal subjects) were imported into a customized MATLAB application, and were manually segmented by three graders at the layers corresponding to the inner limiting membrane (ILM), the inner segment/ellipsoid interface (ISe), the retinal/retinal pigment epithelium interface (RPE), and the Bruch s membrane (BM). The scans were then segmented independently by the same graders using EdgeSelect, a semi-automated method allowing the graders to guide/correct the layer segmentation interactively. The inter-grader reproducibility and agreement in locating the layer positions between the manual and EdgeSelect methods were assessed and compared using the Wilcoxon signed rank test. Results: The inter-grader reproducibility using the EdgeSelect method for retinal layers varied from 0.15 to 1.21 mm, smaller than those using the manual method ( mm). The Wilcoxon test indicated the EdgeSelect method had significantly better reproducibility than the manual method. The agreement between the manual and EdgeSelect methods in locating retinal layers ranged from 0.08 to 1.32 mm. There were small differences between the two methods in locating the ILM (p = 0.012) and BM layers (p,0.001), but these were statistically indistinguishable in locating the ISe (p = 0.896) and RPE layers (p = 0.771). Conclusions: The EdgeSelect method resulted in better reproducibility and good agreement with a manual method in a set of eyes of normal subjects and with retinal disease, suggesting that this approach is feasible for OCT image analysis in clinical trials. Citation: Huang Y, Danis RP, Pak JW, Luo S, White J, et al. (2013) Development of a Semi-Automatic Segmentation Method for Retinal OCT Images Tested in Patients with Diabetic Macular Edema. PLoS ONE 8(12): e doi: /journal.pone Editor: Orhan Aktas, University of Düsseldorf, Germany Received April 21, 2013; Accepted October 29, 2013; Published December 26, 2013 Copyright: ß 2013 Huang et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Funding: This research was supported by an unrestricted fund provided by Research to Prevent Blindness Foundation to Department of Ophthalmology, University of Wisconsin - Madison; the funder had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. The authors were grateful that the data was collected and provided by DRCR.net for analysis; DRCR.net has no role in the study design, data analysis, decision to publish, or preparation of the manuscript. Competing Interests: The authors have declared that no competing interests exist. * yhuang@rc.ophth.wisc.edu Introduction Optical Coherence Tomography (OCT) - determined retinal layer thickness measurements have used in clinical trials as quantitative, morphologic endpoints for the diagnosis and classification of retinal diseases and for monitoring treatment effects [1,2,3,4,5,6], including studies in patients with diabetic macular edema. In order for the retinal layer thickness be measured properly, accuracy of the applied layer segmentation method becomes an important determinant especially when structural-altering retinal lesion is present. With the advent of spectral domain OCT, a number of research groups and manufacturers have engaged in developing layer segmentation strategies [7,8,9,10,11,12,13]. To date, published and commercially available segmentation methods fall into two groups, either fully automated or manual. With a fully-automated method, computer algorithm determines the desired layers with no human supervision. While it is very convenient and practical in the traditional clinical practice settings, the segmentation results using automatic algorithms are prone to layer misidentification errors, especially in eyes with intermediate to severe retinal lesions where the error rate ranges from 30% to 80% [14,15,16,17,18]. The layer misidentification, while confounded by other OCT imaging and operating errors such as weak signal quality, eye movement, and decentration, are mostly caused by the complexity of morphological configuration and reflectivity changes of retinal PLOS ONE 1 December 2013 Volume 8 Issue 12 e82922

2 lesions that are beyond the reach of traditional image segmentation techniques. Manual segmentation by human graders has been considered the gold standard in many previous reports [12,13,19]. Manual segmentation by human graders usually requires the grader to identify the layers either by free-hand drawing [12], or placing seed points and the computer interpolating the layers via pointfitting algorithms [13]. While the gross errors of layer misidentification frequented in the automated methods are avoided, the manual methods are usually time- and labor-intensive, and generally yield higher inter-grader variability. Taking into consideration the advantages and drawbacks of both approaches, we developed a semi-automatic, interactive segmentation method (called herein EdgeSelect). We tested the performance of the EdgeSelect method and compared it against a manual method for detection of the inner and outer boundaries of the retina. Methods Study dataset and preparation Retina volume scans from 28 eyes consisting of 4 normal subjects and 24 patients with diabetic macular edema were obtained using Heidelberg Spectralis OCT devices (Heidelberg Engineering Inc, Heidelberg, Germany). The patient scans were provided by the Diabetic Retinopathy Clinical Research Network (DRCR.net) from participants enrolled in a study comparing the measurements from spectral domain and time domain OCTs. The volume scans consist of horizontal raster lines covering a mm area centered at the fovea (Figure 1A), of which each B scan was captured with 9-frame averaging (ART = 9). Scans were obtained with certified photographers to minimize the OCT data acquisition artifacts [15,20]. The data samples were saved in the Heidelberg proprietary.e2e format. They were exported from a Heidelberg Heyex review software (version 5.1) in.vol format and converted to the DICOM (Digital Imaging and Communication in Medicine) [21] OPT (ophthalmic tomography) format using a custom application built in MATLAB (MATLAB R2011b, The MathWorks, Natick, MA). The standardized OCT images were then segmented by three graders independently for the layers of the inner limiting membrane (ILM), the inner segment/ellipsoid interface (ISe), the retina/retinal pigment epithelium interface (RPE), and Bruch s membrane (BM), using both manual and EdgeSelect segmentation methods (Figure 1B). The study was conducted in accordance with Health Insurance Portability and Accountability Act (HIPAA) requirements and the tenets of the Declaration of Helsinki. All images were obtained under informed written consent and de-identified through the DRCR.net clinical trial protocol, and the study protocol was approved by the institutional review board of the University of Wisconsin-Madison. Figure 1. A representative OCT image with layers segmented using EdgeSelect. (A) OCT images were captured with horizontal raster lines covering mm area centered at the fovea. (B) A representative B scan is shown by the red line in (A). The inner limiting membrane (ILM), the inner-segment/ellipsoid interface (ISe), the retinal/ retinal pigment epithelium interface (RPE), and the Bruch s membrane (BM), were segmented. doi: /journal.pone g001 Segmentation Procedures The manual method. The manual segmentation method used in this study was similar to the method described by Hood, et. al. [13] Graders used the computer mouse to place the seed points at a desired OCT layer, and the computer automatically interpolated the layer based on a spline fitting algorithm. The graders continued to add/delete seed points until the resultant spline line adequately demarcated the layer. The graders segmented each of the four layers for all B-scans in the volumetric data before continuing on the next layer. The EdgeSelect method. Several weeks after the graders performed the manual segmentation of the scans, they resegmented them for the same four layers using the EdgeSelect method. The general steps of the EdgeSelect method to generate the retinal layers are illustrated and described in Figure 2 by segmenting the ILM layer in a representative B-scan. It is of note that the steps leading to interactive edge selection and layer generation are based on automatic image processing routines and require no grader intervention. 1. To reduce speckle noise, OCT images were smoothed by applying a 363 pixel boxcar averaging filter. 2. The layer locations in OCT images were defined as the transition between two regions with difference in reflectance intensity, which were identified in EdgeSelect as the zerocrossing pixels when a Laplacian of Gaussian filter was applied [22]. The zero-crossing pixels correspond to the local maximum gradient in intensity. Figure 2B shows the zerocrossing binary image derived from the original OCT image. 3. Each of the zero-crossing pixels were assigned values with a weighting scheme based on the intensity and gradient of the neighboring pixels. The weighting scheme was designed based on the characteristics of the desired layer. For ILM, a gradient filter enhancing the transition from low reflectivity to high reflectivity was used. The weighted zero-crossing map was plotted in Figure 2C; many spurious edges were suppressed by the weighting scheme. 4. The weighted zero-crossing map was then applied with a Canny-edge detection scheme [23]: only the pixels connected with pixels with strong weight were considered as edge. PLOS ONE 2 December 2013 Volume 8 Issue 12 e82922

3 Figure 2D illustrates the edge candidates superimposed on the original OCT B-scan. 5. To generate the layer location, EdgeSelect used the Dijkstra s algorithm [24] interactively based on the edges selected by the user. The program automatically selected one edge candidate from the leftmost of the image as the source, and one from the rightmost of the image as the destination, and generated the shortest path from the source to the destination. The choices of the source and destination edges depends on the relative position of the desired layer within the retinal tissue. For example, if ILM layer is segmented, the program will choose the edges near the vitreal high reflectivity band; if ISe, RPE, or BM layers are segmented, the program will choose the edges near the posterior high reflectivity band. If the generated layer was satisfactory to demarcate the desired layer, no further human intervention was necessary. However, in diseased retinas where the layered structure was altered, the grader could pick additional or remove edge candidates from the initial set for the Dijkstra s algorithm to re-route the shortest path until the resultant layer was correct. This process is illustrated in Figure 2E and 2F. Figure 2E shows the resultant ILM layer from initial automatic edge selection. Because of the presence of a large cyst in the OCT image, the initial shortestpath failed to demarcate correctly the layer on the center portion of the image. The grader intervened by selecting the correct edge on the top of the cyst, and the ILM layer was recalculated automatically with the inclusion of the new edge candidate (Figure 2F). Data Analysis Comparison of the manual and EdgeSelect methods was performed by evaluating the inter-grader reproducibility for each method, and the agreement of the layer locations between the two methods. The inter-grader reproducibility of the mean layer location was calculated with concordance correlation coefficients. In addition, we also evaluated the variability at the pixel level. For a particular layer, the absolute difference in boundary location (DBL) at each pixel was calculated between any grader pairs, and the inter-grader reproducibility was defined as the value of the DBL averaged across the entire 3D volumetric data and among the three grader pairs. Scatter plots of the inter-grader reproducibility of each data sample was used to compare the manual and EdgeSelect method visually, and Wilcoxon signed rank test [25] was used to determine if the inter-grader reproducibility between the two methods differed, assuming the data sample distribution was not normal. We were also interested in determining whether the layers were accurately segmented using the EdgeSelect method, i.e., if the layer locations identified by the EdgeSelect method were in agreement with those by the manual method. At each pixel, we computed the layer location averaged among the three graders for the EdgeSelect and the manual methods independently. The agreement was measured by computing the absolute difference of the mean layer locations between the manual and the EdgeSelect methods at each pixel, and then averaged across the 3D data. Similarly, scattered plot of agreement of each data sample was used for visual inspection, and Wilcoxon signed rank test was used to determine the statistical difference between the two measurements. Results Inter-grader reproducibility Both the Manual and EdgeSelect methods exhibited high concordance correlation among graders in identifying the mean layer locations of the ILM, ISe, RPE and BM (Table 1). At the pixel level, reproducibility for the EdgeSelect method varied from 0.15 to 1.22 mm; in contrast, reproducibility for the manual method ranged from 3.36 to 6.43 mm. To further illustrate the difference in reproducibility between these two methods, the scatter plots of reproducibility for the four layers are shown in Figure 3. In each of the 28 data samples, reproducibility of the EdgeSelect method pixel-wise was better than that of the manual method. Wilcoxon test indicated the improvement was statistically significant (p,0.001) for each layer. Figure 2. Graphic representation of segmenting the inner limiting membrane (ILM) interface in a representative B scan using the EdgeSelect method. (A) A representative OCT B scan. (B) Detection of image contrast change using zero-crossing of a Laplacian-of-Gaussian (LoG) filter. (C) The detected edges are assigned different weights based on intensity/gradient characteristics. Higher intensity represents larger weight, corresponding to strong edges; lower intensity signals are assigned lesser weight corresponding to weak edges. (D) An edge candidate map (blue lines) is generated using a Canny-like filtering scheme, and is superimposed on the original OCT image. (E) The start and end edge candidates are initiated, and the path of the shortest distance via Dijkstra s algorithm (red line) is computed. (F) Human grader intervention adds additional seed edges, and the program regenerates automatically the updated shortest path until the proper layer segmentation is reached. doi: /journal.pone g002 PLOS ONE 3 December 2013 Volume 8 Issue 12 e82922

4 Figure 3. Comparison of the inter-grader reproducibility between the manual and the EdgeSelect methods in all 28 data samples. The filled symbols are data points from patients and the open symbols are from normal subjects. ILM: inner limiting membrane; ISe: the innersegment/ellipsoid interface; RPE: the retinal/retinal pigment epithelium interface; BM: the Bruch s membrane. doi: /journal.pone g003 It was also noted that the inter-grader reproducibility reported here was similar in magnitude to the MEAN (mean(dlbl)) reported by Hood et. al. [13], in which it was reported to be mm. Since the MEAN (mean(dlbl)) by Hood et. al. was defined as the absolute difference between grader and the mean location of all graders, and in our report the inter-grader reproducibility was defined as the absolute difference between grader pairs, the value of the inter-grader reproducibility is twice Table 1. Mean, standard deviation, and concordance correlation of the inter-grader reproducibility of the EdgeSelect and the Manual methods. Layer Concordance Mean (mm) Std. Dev. (mm) Corr. EdgeSelect ILM ISe RPE BM Manual ILM ISe RPE BM Layers are: inner limiting membrane (ILM), inner segment/ellipsoid interface (ISe), retina/retinal pigment epithelium interface (RPE), and Bruch s membrane (BM). doi: /journal.pone t001 as large as the MEAN (mean(dlbl)). Not surprisingly, the range of the inter-grader reproducibility for the manual method, mm, was congruent with those reported by Hood et. al., as the manual methods employed in these two reports were virtually identical. Agreement Figure 4 shows the agreement in identifying layer locations between the manual and EdgeSelect methods. The mean difference for the ILM, ISe, RPE, and BM layers were mm (mean 6 1 std. dev.), mm, mm, and mm, respectively. There was a small difference between the manual and EdgeSelect methods in identifying the ILM (p = 0.012, Wilcoxon test) and BM layers (p,0.001), but the ISe and RPE layers were statistically indistinguishable between the two methods (p = and p = 0.771). The small difference between the two methods in locating the ILM and BM layers was presumably due to the disparity that EdgeSelect method objectively found the location of the local maximum gradient in reflectance, while the results of the manual method were influenced by the human perception of an edge, in particular where there was an asymmetrical intensity profile of the neighboring pixels. Additional analysis is done for comparing the inter-grader variability and agreement in patient and normal subject groups independently. The results is provided as Table S1. PLOS ONE 4 December 2013 Volume 8 Issue 12 e82922

5 Figure 4. Agreement in segmented layer locations between the manual and EdgeSelect methods. Solid lines indicate the mean difference of the agreement between the two methods and the dashed lines indicate mean 6 2 standard deviations. The filled symbols are data points from patients and the open symbols are from normal subjects. ILM: inner limiting membrane; ISe: the inner-segment/ellipsoid interface; RPE: the retinal/retinal pigment epithelium interface; BM: the Bruch s membrane. doi: /journal.pone g004 Discussion Accurate segmentation of OCT retinal images has been a challenge for OCT device manufacturers and research groups. An ideal algorithm or strategy incorporates elements to ensure high accuracy in layer identification in various clinical conditions as well as low inter-session variability. Among the existing methods, automatic segmentation algorithms have the computational superiority and usually produce perfect inter-session reliability when applied to the same OCT image. However, in diseased eyes where retinal lesions are heterogeneous and complex, automated methods frequently fail to identify retinal layers properly, and the human observer remains the best decision-maker for definition and control of the desired segmented results. In contrast, manual methods usually have limited layer misidentification errors, but tend to have high inter-grader variability and are usually time and effort consuming. In this report, we described the EdgeSelect method, which allows graders to incorporate their clinical knowledge to guide the selection of the proper edge candidates of a particular layer, but delegates the calculation of the exact pixel location of the layer path to the computer via using the Dijkstra s algorithm interactively. When comparing with a manual method, the EdgeSelect method demonstrated better inter-grader reproducibility, while maintaining good agreement with the manual method. When compared to an automatic method, the EdgeSelect method should perform at the same level of efficiency in the retinas with normal or near-normal structure, as the initial automatic edge selection is likely accurate and hence no human intervention is needed. In diseased retinas, the advantage of the interactive nature of the EdgeSelect method becomes evident, especially where patchy or locally discontinuous layer signals are common. Using EdgeSelect, we anticipate increased efficiency in layer segmentation, which should allow segmentation of large 3D data sets to become feasible for graders. Lastly, the EdgeSelect method relies on local transition of the reflectance signal to determine the proper edges, which is less likely to be device dependent. Together with methods of data standardization and homogenization [26], EdgeSelect can be applied universally to images from different SD-OCT devices, which would allow all OCT images to be processed using the same software algorithm for harmonization of the measurements. Supporting Information Table S1 The inter-grader variability and agreement between EdgeSelect and manual measurement is analyzed independently for the patient and normal subjects groups. (DOCX) PLOS ONE 5 December 2013 Volume 8 Issue 12 e82922

6 Author Contributions Conceived and designed the experiments: YH RPD. Performed the experiments: YH RPD JWP SL JW AD. Analyzed the data: YH RPD XZ AN. Wrote the paper: YH. References 1. Drexler W, Fujimoto JG (2008) State-of-the-art retinal optical coherence tomography. Progress in retinal and eye research 27: Browning DJ, Glassman AR, Aiello LP, Bressler NM, Bressler SB, et al. (2008) Optical coherence tomography measurements and analysis methods in optical coherence tomography studies of diabetic macular edema. Ophthalmology 115: , 1371 e Keane PA, Patel PJ, Liakopoulos S, Heussen FM, Sadda SR, et al. (2012) Evaluation of age-related macular degeneration with optical coherence tomography. Survey of ophthalmology 57: Lim JI, Tan O, Fawzi AA, Hopkins JJ, Gil-Flamer JH, et al. (2008) A pilot study of Fourier-domain optical coherence tomography of retinal dystrophy patients. American journal of ophthalmology 146: Albrecht P, Ringelstein M, Müller A, Keser N, Dietlein T, et al. (2012) Degeneration of retinal layers in multiple sclerosis subtypes quantified by optical coherence tomography. Multiple Sclerosis Journal 18: Nolan R, Gelfand JM, Green AJ (2013) Fingolimod treatment in multiple sclerosis leads to increased macular volume. Neurology 80: Chiu SJ, Li XT, Nicholas P, Toth CA, Izatt JA, et al. (2010) Automatic segmentation of seven retinal layers in SDOCT images congruent with expert manual segmentation. Optics express 18: Garvin MK, Abramoff MD, Wu X, Russell SR, Burns TL, et al. (2009) Automated 3-D intraretinal layer segmentation of macular spectral-domain optical coherence tomography images. IEEE transactions on medical imaging 28: Vermeer KA, van der Schoot J, Lemij HG, de Boer JF (2011) Automated segmentation by pixel classification of retinal layers in ophthalmic OCT images. Biomedical optics express 2: Yang Q, Reisman CA, Chan K, Ramachandran R, Raza A, et al. (2011) Automated segmentation of outer retinal layers in macular OCT images of patients with retinitis pigmentosa. Biomedical optics express 2: Ahlers C, Simader C, Geitzenauer W, Stock G, Stetson P, et al. (2008) Automatic segmentation in three-dimensional analysis of fibrovascular pigmentepithelial detachment using high-definition optical coherence tomography. The British journal of ophthalmology 92: Sadda SR, Joeres S, Wu Z, Updike P, Romano P, et al. (2007) Error correction and quantitative subanalysis of optical coherence tomography data using computer-assisted grading. Investigative ophthalmology & visual science 48: Hood DC, Cho J, Raza AS, Dale EA, Wang M (2011) Reliability of a computeraided manual procedure for segmenting optical coherence tomography scans. Optometry and vision science : official publication of the American Academy of Optometry 88: Sadda SR, Wu Z, Walsh AC, Richine L, Dougall J, et al. (2006) Errors in retinal thickness measurements obtained by optical coherence tomography. Ophthalmology 113: Han IC, Jaffe GJ (2010) Evaluation of artifacts associated with macular spectraldomain optical coherence tomography. Ophthalmology 117: e Ho J, Sull AC, Vuong LN, Chen Y, Liu J, et al. (2009) Assessment of artifacts and reproducibility across spectral- and time-domain optical coherence tomography devices. Ophthalmology 116: Lammer J, Scholda C, Prunte C, Benesch T, Schmidt-Erfurth U, et al. (2011) Retinal thickness and volume measurements in diabetic macular edemfa: a comparison of four optical coherence tomography systems. Retina 31: Krebs I, Hagen S, Brannath W, Haas P, Womastek I, et al. (2010) Repeatability and reproducibility of retinal thickness measurements by optical coherence tomography in age-related macular degeneration. Ophthalmology 117: Seigo MA, Sotirchos ES, Newsome S, Babiarz A, Eckstein C, et al. (2012) In vivo assessment of retinal neuronal layers in multiple sclerosis with manual and automated optical coherence tomography segmentation techniques. Journal of neurology 259: Tewarie P, Balk L, Costello F, Green A, Martin R, et al. (2012) The OSCAR-IB consensus criteria for retinal OCT quality assessment. PLoS One 7: e NEMA (2013) Digital Imaging and Communication in Medicine (DICOM) Standard. 22. Rogowska J (2008) Overview and Fundamentals of Medical Image Segmentation. In: Bankman I, editor. Handbook of Medical Image Processing and Analysis: Academic press. pp Canny J (1986) A computational approach to edge detection. IEEE transactions on pattern analysis and machine intelligence 8: Sniedovich M (2006) Dijkstra s algorithm revisited: the dynamic programming connexion. Control and Cybernetics 35: Siegel S (1956) Non-parametric statistics for the behavioral sciences. New York: McGraw-Hill. 26. Huang Y, Gangaputra S, Lee KE, Narkar AR, Klein R, et al. (2012) Signal quality assessment of retinal optical coherence tomography images. Investigative ophthalmology & visual science 53: PLOS ONE 6 December 2013 Volume 8 Issue 12 e82922