|
|
||||||||
1From the UPMC Eye Center, Department of Ophthalmology, University of Pittsburgh School of Medicine, Pittsburgh, Pennsylvania; the 2New England Eye Center, Tufts-New England Medical Center, Tufts University School of Medicine, Boston, Massachusetts; and the 3Massachusetts Institute of Technology, Cambridge, Massachusetts.
| Abstract |
|---|
|
|
|---|
METHODS. Four layer structures within the retina were defined: the macular nerve fiber layer (mNFL), the inner retinal complex (IRC; retinal ganglion cell [RGC] layer+inner plexiform and nuclear layers), outer plexiform layer (OPL), and outer retinal complex (ORC; outer nuclear layer+photoreceptor layer). Normal and glaucomatous eyes underwent fast macular map and fast NFL OCT scans. Linear macular images were analyzed using the developed algorithm, and the results were compared with the cpNFL thickness measurement.
RESULTS. Forty-seven subjects (23 normal and 24 with glaucoma) were analyzed. mNFL, cpNFL, IRC, and the total retinal thicknesses were significantly greater in normal than in glaucomatous eyes (P
0.0002; Wilcoxon), whereas OPL thickness did not show a significant difference (P = 0.46). ORC thickness was significantly greater in glaucomatous than normal eyes (P = 0.035). Areas under the receiver operator characteristic curve (AROCs) for discriminating normal from glaucomatous eyes were highest with mNFL+IRC (0.97) and lowest with OPL (0.56). AROCs for OPL and ORC were significantly smaller than those for mNFL, IRC, mNFL+IRC, and cpNFL (P
0.01). AROCs for IRC, mNFL+IRC, and cpNFL were significantly larger than for retinal thickness (P
0.049). Among the best-performing parameters (mNFL, IRC, mNFL+IRC, and cpNFL) there was no significant difference in AROCs (P
0.15).
CONCLUSIONS. The newly developed macular segmentation algorithm described herein demonstrated its ability to quantify objectively the glaucomatous damage to RGCs and NFL and to discriminate between glaucomatous and normal eyes. Further algorithm refinement and improvements in resolution and image quality may yield a more powerful methodology for clinical glaucoma evaluation.
Since Zeimer et al.1 suggested the concept that macular thickness could be used as a measure of glaucoma, several studies have explored this idea by using emerging ocular imaging modalities, such as the retinal thickness analyzer (RTA)2 3 and OCT.4 5 6 7 8 The theory behind this concept is that assessment of macular thickness may be superior to circumpapillary nerve fiber layer (cpNFL) measurements, because retinal ganglion cell (RGC) bodies are 10 to 20 times the diameter of their axons and because the RGC layer is more than one cell thick in the macula.
Although the studies mentioned showed that visual function correlates significantly with cpNFL thickness and macular thickness, the cpNFL thickness measurement outperform the macular thickness in terms of magnitude of association with visual function.4 7 This may be because the measured macular thickness represents the total retinal thickness. In this case retinal structures other than the RGC layer may confound the data so that the macular thickness measurement was not as sensitive as cpNFL thickness measurement.
The currently available commercial OCT unit (StratusOCT; Carl Zeiss Meditec, Inc., Dublin, CA) is the third generation of OCT technology, with a resolution of 8 to 10 µm and a scanning speed of approximately 400 Hz. With this high-resolution OCT, retinal layer structures can be distinguished in great detail. The purpose of this study was to develop a software algorithm to perform automated segmentation of retinal layer structures on linear macular OCT scan images and to test its performance in discriminating normal from glaucomatous eyes, in comparison with conventional cpNFL thickness measurement.
| Materials and Methods |
|---|
|
|
|---|
The summarized flow of the algorithm follows.
Preprocessing.
Before segmentation, images were processed to yield better algorithm performance.
|
|
|
Learning Set.
A small group of subjects (5 normal eyes and 11 glaucomatous eyes) were used to subjectively assess the performance of the algorithm in proper border detection. These subjects were not used for the testing of the algorithm described in the next section.
Testing of the Algorithm
Subjects.
Sixty-four subjects (27 normal and 37 with glaucoma) were scanned with a StratusOCT unit using the fast macular mapping and fast NFL scan modes (a set of three consecutive peripapillary circular scans 3.4 mm in diameter). Institutional Review Board (IRB) and Ethics Committee approval was obtained for the study, and all participants gave their informed consent to participate study. The study adhered to the principles of the Declaration of Helsinki.
All the participants had comprehensive ocular examination, reliable Swedish Interactive Thresholding Algorithm (SITA) standard 24-2 perimetry results (Carl Zeiss Meditec, Inc.), and OCT scanning of the macular and peripapillary regions at the same visit. All eyes had best corrected visual acuity of 20/40 or better, refractive error between 6.00 and +3.00 D, and no media opacities. Patients with diabetes, any medical condition that might affect visual field (VF) other than glaucoma, or treatment with medications that might affect retinal thickness were excluded from the study.
Normal subjects were healthy volunteers with normal ocular examination findings and VF glaucoma hemifield test (GHT) results within normal limits. All subjects with glaucoma had VF GHT results outside normal limits and a VF defect defined as one with fewer than 30% fixation losses, false-positive, or false-negative responses. They also had more than one of the following conditions: optic disc rim notching, cup asymmetry, large cupping (>0.7 vertical cup/disc ratio), NFL defect, or intraocular pressure >32 mm Hg. These glaucoma subjects had no history of ocular trauma or surgery other than glaucoma interventions or uncomplicated cataract extraction and no ocular diseases not attributed to glaucoma. One eye was randomly selected if both eyes qualified for the study.
Image Quality Control/Assessment.
During acquisition, OCT scans had to demonstrate signal-to-noise ratio (SNR)
35 dB. Because SNR does not take into account all aspects of image quality (e.g., signal intensity distribution), we also used our previously presented quality index (QI) to further differentiate acceptable/unacceptable scans (Ishikawa H, et al. IOVS 2003;44:ARVO E-Abstract 3358; Ishikawa H, et al. IOVS 2004;45:ARVO E-Abstract 3317). QI was defined as a product of intensity ratio (IR) and tissue signal ratio (TSR), where IR was equivalent to SNR calculated from the logged raw data, and TSR was a ratio of the number of the pixels within a highly reflective area versus low reflective area. The cutoff of QI for poor-quality images on the OCT macular scan was 17.5 based on our unpublished data. Any subjects with a mean QI rating <17.5 for the six macular scans was excluded from the study.
Data and Statistical Analysis.
All linear macular OCT images were analyzed by using the developed algorithm. cpNFL thickness was calculated by the OCT system software version 3.0 (Carl Zeiss Meditec, Inc.). Both algorithm outputs were subjectively evaluated by two experts (HI, DS) independently to detect algorithm errors. Criteria of algorithm error were obvious disruption of the detected border, and/or border wandering (detected border jumping to and from different anatomic structures) for >5% consecutive (i.e., an uninterrupted error) or 20% cumulative (i.e., adding up all errors amounts to 20% of the image width) of the entire image.
Thicknesses of four different retinal segments (mNFL, IRC, OPL, and ORC) were calculated on each linear macular image, and thickness maps were plotted on each set of fast macular scans. Overall thicknesses of each segment were then calculated, incorporating interpolated point values.
Differences in thicknesses between normal and glaucomatous eyes were evaluated using the Wilcoxon test. The capability of discriminating normal from glaucomatous eyes was assessed according to the area under the receiver operator characteristic curve (AROC).
| Results |
|---|
|
|
|---|
Normal subjects were younger than subjects with glaucoma (52.1 ± 14.8 normal vs. 62.4 ± 14.5 glaucoma; P = 0.02, t-test).
Algorithm failure occurred more often with poor-quality (55.6%) than with good-quality (10.2%) images (Table 1) . Expert assessment of algorithm failure showed high agreement between observers (
= 0.89). The ILM border was the most reliable border, whereas the IRC border was the least reliable one (Table 2) .
|
|
0.0002, Wilcoxon), whereas OPL thickness did not show significant difference (P = 0.46; Table 3 , Fig. 4 ). ORC thickness was significantly thicker in glaucomatous than in normal eyes (P = 0.035). The power was 0.998 to detect a difference in mNFL, IRC, or mNFL+IRC measurements between glaucomatous and normal eyes with
= 0.05. Age showed statistically significant linear correlation with all thickness parameters (Pearson correlation, all P < 0.01) except for ORC and OPL (both P > 0.2). QI also showed statistically significant linear correlation with all thickness parameters (all P < 0.02) except for OPL and total retinal thickness (both P > 0.6).
|
|
0.01). AROCs for IRC, mNFL+IRC, and cpNFL were significantly larger than that for total retinal thickness (P
0.049). Among the top-performing parameters (mNFL, IRC, mNFL+IRC, and cpNFL) there was no significant difference in AROCs (P
0.15).
|
|
| Discussion |
|---|
|
|
|---|
The use of the fast macular mapping mode, which limits the number of sampling lines (A-scans) per image to 128, may be controversial. The improved visualization of the retinal layers with StratusOCT in comparison with older-generation OCT units is not only due to its higher axial resolution (which is negligible) but also to its denser transverse sampling (up to 512 sampling lines per image), feasible because of the StratusOCTs 400 A-scan/s scan rate. Hence higher-density scans may theoretically have better retinal layer separation. However, the price for the higher-density scan, which is decreased scanning speed, should not be overlooked. Six images must be acquired for macular mapping, which makes the overall scanning time approximately 7.5 seconds or longer. StratusOCT does not provide a scan mode that performs six high-density scan in a rapid succession without stopping between each image acquisition. Therefore, scan registration (centering at the fovea) is less accurate than the fast macular mapping scan, which completes the scanning of six 128-sampling images in just under 2 seconds. We also found that, for segmentation purposes, the smoothed images (Fig. 1b) using 128, 256, and 512 A scans/image showed no significant difference in appearance or in segmentation algorithm performance (Ishikawa H, unpublished data, 2003). Therefore, we decided to use the most commonly used fast macular mapping scanning mode for this study.
To improve the stability of the algorithm, we combined several retinal layers so that detected borders showed minimum variability. We found that the outer borders of the mNFL, inner nuclear layer, and OPL were the most stable of all possible borders. Thus, four different segments (mNFL, IRC, OPL, and ORC) were detected by the algorithm.
Even after this simplification of macular segmentation, we experienced algorithm failure of at least one detected border in approximately 10% of the good-quality images. Although a higher failure rate in poor-quality images was expected, we were surprised to find that the algorithm failed on more than half of them. When considering the most common areas for failure, we found that the outer borders of the IRC and OPL were least reliable. This may have been because the IRC and OPL had less pixel intensity and blurrier borders than did the better-defined layers such as ILM, NFL, and RPE. Further refinement of the algorithm is required.
In the present population, 88.9% of the OCT images from normal and 64.9% from glaucoma eyes were successfully enrolled in the study as good-quality images. This implies that the diseased eyes were harder to get scanned in good quality.
All parameters but OPL showed significant differences between normal and glaucomatous eyesa finding that agreed with previous histologic studies.11 12 As expected, the total retinal thickness was greater in normal than in glaucomatous eyes, also in agreement with previous studies.3 5 6 7 8 Because the normal and glaucoma groups were not age matched, these discrepancies could be due to age differences. We found statistically significant correlations between age and thickness parameters, except for OPL and ORC. However, based on the correlation analysis, the measured differences far exceeded the expected differences due to age (e.g., the measured difference in mNFL+IRC was approximately 30 µm, whereas the expected age-related difference was approximately 6.3 µm).
To our surprise, although the glaucoma discrimination power of the ORC was the second worst, the ORC was thicker in glaucomatous eyes than in normal ones. The ORC is a combination of the outer nuclear layer and photoreceptor layer. According to the current understanding of pathophysiology, these layers are thought not to be affected by glaucoma. A potential explanation of this finding is that our definition of ORC on the OCT A-scan profile was inaccurate, perhaps encompassing more structural elements than expected. Further investigation is needed, especially in the correlation of histologic structures and OCT images.
Contrary to expectations, mNFL+IRC thickness was the best-performing macular segmentation parameter in the present study. Our expectation was that IRC thickness would be the best performer, because it was the closest to the RGC layers thickness.
The biggest sources of instability were speckle noise and uneven tissue reflectivity across sampling lines within an image. There are two approaches to minimize these factors; higher resolution, and improved signal quality (signal-to-noise ratio). Of course, further sophistication of the algorithm by using better methods and procedures for preprocessing filtering and border detection may be another key to improving the overall performance of macular segmentation analysis.
Data from our group published by Guedes et al.7 showed that the AROC of macular thickness (total retinal thickness) was 0.80 for discriminating normal from those with advanced glaucoma (mean MD, 12.0 dB). Although these data were based on the first-generation OCT unit, which provided lower resolution images than the StratusOCT we used in the present study, our results agree with these nearly perfectly. In other words, even with the higher-resolution images, total macular thickness showed quite similar discriminating power. A major reason for this is probably that macular thickness measurement did not take advantage of the improved resolution since the ILM and RPE were the most clearly identifiable layers, even on the lower resolution OCT images. One can argue this also in another way. As mentioned earlier, glaucomatous eyes may tend to have worse image quality than normal eyes, and all thickness parameters except for OPL and total retinal thickness showed significant correlations with QI. The combination of these facts may have yielded artifactiously higher sensitivities for mNFL+IRC and cpNFL, as well as a lower sensitivity for total retinal thickness.
Guedes et al.7 and Wollstein et al.4 reported that cpNFL was the best parameter among many different ocular imaging parameters in terms of glaucoma discrimination. In the present study, we found that mNFL, IRC, and mNFL+IRC performed at least as well as cpNFL and better than the whole retinal thickness. This was an unexpected finding, as the cpNFL had the following theoretical advantages over the mNFL: First, the OCT scanning pattern in the macula employs a spoke pattern configuration, whereas continuous circular scanning is used in the circumpapillary region. The spoke pattern may miss information between each linear scan. Whereas at the scan patterns center the points are overlapping, at the farthest extreme, there is a 1.6-mm space between spokes. Second, since all retinal nerve fibers radiate toward the optic nerve head, circular scanning around the optic nerve head is capable of detecting damage in areas that are not covered in the macular scan (i.e., the representation of the 50% of the eyes ganglion cells that do not reside in the macula). Finally, injury to the RGC may not affect the detectable structure of the cell bodies early on, but could result in thinning of their axons. Further investigation is needed.
In conclusion, our newly developed automated algorithm for macular segmentation showed potential usefulness for objective quantification of glaucomatous macular RGC and NFL damage. Improved algorithm sophistication and in scan resolution and image quality may yield a powerful, novel method for clinical glaucoma evaluation.
| Footnotes |
|---|
Submitted for publication March 24, 2004; revised June 23 and August 10, 2004, and January 4 and February 18, 2005; accepted February 28, 2005.
Disclosure: H. Ishikawa, None; D.M. Stein, None; G. Wollstein, None; S. Beaton, None; J.G. Fujimoto, Carl Zeiss Meditec (F); J.S. Schuman, Carl Zeiss Meditec (F)
The publication costs of this article were defrayed in part by page charge payment. This article must therefore be marked "advertisement" in accordance with 18 U.S.C.
1734 solely to indicate this fact.
Corresponding author: Hiroshi Ishikawa, UPMC Eye Center, Department of Ophthalmology, University of Pittsburgh School of Medicine, 203 Lothrop Street, Eye and Ear Institute, Room 841, Pittsburgh, PA 15213; ishikawah{at}upmc.edu.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
H. W. van Dijk, P. H. B. Kok, M. Garvin, M. Sonka, J. H. DeVries, R. P. J. Michels, M. E. J. van Velthoven, R. O. Schlingemann, F. D. Verbraak, and M. D. Abramoff Selective Loss of Inner Retinal Layer Thickness in Type 1 Diabetic Patients with Minimal Diabetic Retinopathy Invest. Ophthalmol. Vis. Sci., July 1, 2009; 50(7): 3404 - 3409. [Abstract] [Full Text] [PDF] |
||||
![]() |
I Gorczynska, V J Srinivasan, L N Vuong, R W S Chen, J J Liu, E Reichel, M Wojtkowski, J S Schuman, J S Duker, and J G Fujimoto Projection OCT fundus imaging for visualising outer retinal pathology in non-exudative age-related macular degeneration Br. J. Ophthalmol., May 1, 2009; 93(5): 603 - 609. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Ishikawa, J. Kim, T. R. Friberg, G. Wollstein, L. Kagemann, M. L. Gabriele, K. A. Townsend, K. R. Sung, J. S. Duker, J. G. Fujimoto, et al. Three-Dimensional Optical Coherence Tomography (3D-OCT) Image Enhancement with Segmentation-Free Contour Modeling C-Mode Invest. Ophthalmol. Vis. Sci., March 1, 2009; 50(3): 1344 - 1349. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. M. Nolan, J. M. Stringham, S. Beatty, and D. M. Snodderly Spatial Profile of Macular Pigment and Its Relationship to Foveal Architecture Invest. Ophthalmol. Vis. Sci., May 1, 2008; 49(5): 2134 - 2142. [Abstract] [Full Text] [PDF] |
||||
![]() |
V. J. Srinivasan, B. K. Monson, M. Wojtkowski, R. A. Bilonick, I. Gorczynska, R. Chen, J. S. Duker, J. S. Schuman, and J. G. Fujimoto Characterization of Outer Retinal Morphology with High-Speed, Ultrahigh-Resolution Optical Coherence Tomography Invest. Ophthalmol. Vis. Sci., April 1, 2008; 49(4): 1571 - 1579. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Tappeiner, D. Barthelmes, M. H. Abegg, S. Wolf, and J. C. Fleischhauer Impact of Optic Media Opacities and Image Compression on Quantitative Analysis of Optical Coherence Tomography Invest. Ophthalmol. Vis. Sci., April 1, 2008; 49(4): 1609 - 1614. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. R. Sadda, S. Joeres, Z. Wu, P. Updike, P. Romano, A. T. Collins, and A. C. Walsh Error Correction and Quantitative Subanalysis of Optical Coherence Tomography Data Using Computer-Assisted Grading Invest. Ophthalmol. Vis. Sci., February 1, 2007; 48(2): 839 - 848. [Abstract] [Full Text] [PDF] |
||||
![]() |
V. Senatorov, I. Malyukova, R. Fariss, E. F. Wawrousek, S. Swaminathan, S. K. Sharan, and S. Tomarev Expression of Mutated Mouse Myocilin Induces Open-Angle Glaucoma in Transgenic Mice. J. Neurosci., November 15, 2006; 26(46): 11903 - 11914. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Ishikawa, M. L. Gabriele, G. Wollstein, R. D. Ferguson, D. X. Hammer, L. A. Paunescu, S. A. Beaton, and J. S. Schuman Retinal nerve fiber layer assessment using optical coherence tomography with active optic nerve head tracking. Invest. Ophthalmol. Vis. Sci., March 1, 2006; 47(3): 964 - 967. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. Barthelmes, F. K. Sutter, M. M. Kurz-Levin, M. M. Bosch, H. Helbig, G. Niemeyer, and J. C. Fleischhauer Quantitative Analysis of OCT Characteristics in Patients with Achromatopsia and Blue-Cone Monochromatism. Invest. Ophthalmol. Vis. Sci., March 1, 2006; 47(3): 1161 - 1166. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |