Automatic segmentation of the choroid in enhanced depth imaging optical coherence tomography images

Jing Tian, Pina Marziliano, Mani Baskaran, Tin Aung Tun, Tin Aung, Jing Tian, Pina Marziliano, Mani Baskaran, Tin Aung Tun, Tin Aung

Abstract

Enhanced Depth Imaging (EDI) optical coherence tomography (OCT) provides high-definition cross-sectional images of the choroid in vivo, and hence is used in many clinical studies. However, the quantification of the choroid depends on the manual labelings of two boundaries, Bruch's membrane and the choroidal-scleral interface. This labeling process is tedious and subjective of inter-observer differences, hence, automatic segmentation of the choroid layer is highly desirable. In this paper, we present a fast and accurate algorithm that could segment the choroid automatically. Bruch's membrane is detected by searching the pixel with the biggest gradient value above the retinal pigment epithelium (RPE) and the choroidal-scleral interface is delineated by finding the shortest path of the graph formed by valley pixels using Dijkstra's algorithm. The experiments comparing automatic segmentation results with the manual labelings are conducted on 45 EDI-OCT images and the average of Dice's Coefficient is 90.5%, which shows good consistency of the algorithm with the manual labelings. The processing time for each image is about 1.25 seconds.

Keywords: (100.0100) Image processing; (100.2960) Image analysis; (110.4500) Optical coherence tomography; (170.4470) Ophthalmology.

Figures

Fig. 1
Fig. 1
Anatomy of the human eye: the choroid is the vascular layer between the sclera and the retina. Image courtesy of the National Eye Institute, National Institute of Health.
Fig. 2
Fig. 2
(a) EDI-OCT image in the macular region of the human eye. (b) The Bruch’s Membrane and the choroidal-scleral interface are manually labeled by the ophthalmologists.
Fig. 3
Fig. 3
The overview of the automatic segmentation of the choroid in EDI-OCT images.
Fig. 4
Fig. 4
BM detection algorithm. (a) Part of the EDI-OCT image, (b) The detected RPE (CRPE in red dotted line) and Bruch’s membrane (CBM in blue solid line). (c) The new image Îs has the flat Bruch’s membrane.
Fig. 5
Fig. 5
The valleys (local minimums) of the A-scans are used as the feature to detect the choroidal-scleral interface. However, there are also valleys caused by the speckle noise and the blood vessels in the choroid region.
Fig. 6
Fig. 6
Construction of graph G, each valley pixel is a vertex of the graph and is connected to the vertices in the next Nnh columns.
Fig. 7
Fig. 7
(a) A straighted EDI-OCT image. (b) The choroidal-scleral interface is sometimes invisible and the missing boundary is refereed as gaps (circles), where the valley pixels (cross) do not appear near the manual labeled choroidal-scleral interface (dashed lines). (c) The delineation result of the proposed algorithm by setting Nnh = 10 and Nnh = 20. When Nnh = 20, the value 3Nnh (60 pixels) is greater than the horizontal distance of the largest gap (57 pixels). Therefore, vs and ve is connected and the result agrees well with the manual labeling in (b);
Fig. 8
Fig. 8
(a) Construction of graph G, each valley pixel is a vertex of the graph and is connected to the vertex in the next Nnh columns; (b) Hard and soft thresholding methods for assigning penalty weights.
Fig. 9
Fig. 9
The Dice’s coefficient of the choroid segmentation between the output of the proposed algorithm and the manual labelings.
Fig. 10
Fig. 10
The ratio between the error and the manual labelings in estimating the choroidal thickness of 45 EDI-OCT images.
Fig. 11
Fig. 11
(a) A part of non-inverted image I43. (b) The comparison between the results of the proposed algorithm(indicated by the dash lines) and the ground truth labeled by the ophthalmologists (drawn in solid lines). (c) The enlarged major discrepancy region.
Fig. 12
Fig. 12
(a) A part of the inverted image I6. (b) The comparison between the results of the proposed algorithm(indicated by the dash lines) and the ground truth labeled by the ophthalmologists (drawn in solid lines). (c)The enlarged major discrepancy region.

References

    1. Margolis R., Spaide R. F., “A pilot study of enhanced depth imaging optical coherence tomography of the choroid in normal eyes,” Am. J. Ophthalmol. 147, 811–815 (2009).10.1016/j.ajo.2008.12.008
    1. Quigley H. A., “What’s the choroid got to do with angle closure?” Arch. Ophthalmol. 127(5), 693–4 (2009).10.1001/archophthalmol.2009.80
    1. Ikuno Y., Kawaguchi K., Nouchi T., Yasuno Y., “Choroidal thickness in healthy Japanese subjects,” Invest. Ophthalmol. Vis. Sci. 51, 2173–2176 (2010).10.1167/iovs.09-4383
    1. Ding X., Li J., Zeng J., Ma W., Liu R., Li T., Yu S., Tang S., “Choroidal thickness in healthy Chinese subjects,” Invest. Ophthalmol. Vis. Sci. 52, 9555–9560 (2011).10.1167/iovs.11-8076
    1. Manjunath V., Taha M., Fujimoto J. G., Duker J. S., “Choroidal thickness in normal eyes measured using Cirrus HD optical coherence tomography,” Am. J. Ophthalmol. 150, 325–329 (2010).10.1016/j.ajo.2010.04.018
    1. Nickla D. L., Wildsoet C., Wallman J., “The circadian rhythm in intraocular pressure and its relation to diurnal ocular growth changes in chicks,” Exp. Eye. Res. 66, 183–193 (1998).10.1006/exer.1997.0425
    1. Huang D., Swanson E. A., Lin C. P., Schuman J. S., Stinson W. G., Chang W., Hee M. R., Flotte T., Gregory K., Puliafito C. A., “Optical coherence tomography,” Science 254, 1178–1181 (1991).10.1126/science.1957169
    1. Wang R. K., “Signal degradation by multiple scattering in optical coherence tomography of dense tissue: a monte carlo study towards optical clearing of biotissues,” Phys. Med. Biol. 47, 2281–2299 (2002).10.1088/0031-9155/47/13/307
    1. Spaide R. F., Koizumi H., Pozzoni M. C., Pozonni M. C., “Enhanced depth imaging spectral-domain optical coherence tomography,” Am. J. Ophthalmol. 146, 496–500 (2008).10.1016/j.ajo.2008.05.032
    1. Chung S. E., Kang S. W., Lee J. H., Kim Y. T., “Choroidal thickness in polypoidal choroidal vasculopathy and exudative age-related macular degeneration,” Ophthalmology 118, 840–845 (2011).10.1016/j.ophtha.2010.09.012
    1. Fong A. H. C., Li K. K. W., Wong D., “Choroidal evaluation using enhanced depth imaging spectral-domain optical coherence tomography in Vogt-Koyanagi-Harada disease,” Retina 31, 502–509 (2011).10.1097/IAE.0b013e3182083beb
    1. Maruko I., Iida T., Sugano Y., Ojima A., Sekiryu T., “Subfoveal choroidal thickness in fellow eyes of patients with central serous chorioretinopathy,” Retina 31, 1603–1608 (2011).10.1097/IAE.0b013e31820f4b39
    1. Mwanza J.-C., Hochberg J. T., Banitt M. R., Feuer W. J., Budenz D. L., “Lack of association between glaucoma and macular choroidal thickness measured with enhanced depth-imaging optical coherence tomography,” Invest. Ophthalmol. Vis Sci. 52, 3430–3435 (2011).10.1167/iovs.10-6600
    1. Maruko I., Iida T., Sugano Y., Ojima A., Ogasawara M., Spaide R. F., “Subfoveal choroidal thickness after treatment of central serous chorioretinopathy,” Ophthalmology 117, 1792–1799 (2010).10.1016/j.ophtha.2010.01.023
    1. Maruko I., Iida T., Sugano Y., Oyamada H., Sekiryu T., Fujiwara T., Spaide R. F., “Subfoveal choroidal thickness after treatment of Vogt-Koyanagi-Harada disease,” Retina 31, 510–517 (2011).10.1097/IAE.0b013e3181eef053
    1. Yazdanpanah A., Hamar G., “Segmentation of intra-retinal layers from optical coherence tomgraphy images using an active contour approach,” IEEE Trans. Med. Imaging 30, 484–496 (2011).10.1109/TMI.2010.2087390
    1. Garvin M. K., Abramoff M. D., Kardon R., Russell S. R., Wu X., Sonka M., “Intraretinal layer segmentation of macular optical coherence tomography images using optimal 3-D graph search,” IEEE Trans. Med. Imaging 27, 1495–1505 (2008).10.1109/TMI.2008.923966
    1. Chiu S. J., Li X. T., Nicholas P., Toth C. A., Izatt J. A., Farsiu S., “Automatic segmentation of seven retinal layers in SD-OCT images congruent with expert manual segmentation,” Opt. Express 18, 19413–19428 (2010).10.1364/OE.18.019413
    1. Yang Q., Reisman C. A., Wang Z., Fukuma Y., Hangai M., Yoshimura N., Tomidokoro A., Araie M., Raza A. S., Hood D. C., Chan K., “Automated layer segmentation of macular OCT images using dual-scale gradient information.” Opt. Express 18, 21293–21307 (2010).10.1364/OE.18.021293
    1. Koozekanani D., Boyer K., Roberts C., “Retinal thickness measurements from optical coherence tomgraphy using a Markov boundary model,” IEEE Trans. Med. Imaging 20, 906–916 (2001).10.1109/42.952728
    1. Zhang L., Lee K., Niemeijer M., Mullins R. F., Sonka M., Abramoff M. D., “Automated segmentation of the choroid from clinical SD-OCT,” Invest. Ophthalmol. Vis. Sci. 53, 7510–7519 (2012).10.1167/iovs.12-10311
    1. Kajic V., Esmaeelpour M., Povazay B., Marshall D., Rosin P. L., Drexler W., “Automated choroidal segmentation of 1060 nm OCT in healthy and pathologic eyes using a statistical mode,” Biomed. Opt. Express 3, 86–103 (2012).10.1364/BOE.3.000086
    1. Torzicky T., Pircher M., Zotter S., Bonesi M., Gotzinger E., Hitzenberger C. K., “Automated measurement of choroidal thickness in the human eye by polarization sensitive optical coherence tomography,” Opt. Express 20, 7564–7574 (2012).10.1364/OE.20.007564
    1. Duan L., Yamanari M., Yasuno Y., “Automated phase retardation oriented segmentation of chorio-scleral interface by polarization sensitive optical coherence tomography,” Opt. Express 20, 3353–3366 (2012).10.1364/OE.20.003353
    1. Nickla D., Wallman J., “The multifunctional choroid,” Prog. Retinal Res. 29, 144–168 (2010).10.1016/j.preteyeres.2009.12.002
    1. Dijkstra E. W., “A note on two problems in connexion with graphs,” Numerische Math. 1, 269–271 (1959).10.1007/BF01386390
    1. Tian J., Marziliano P., “Location-based graph search algorithm for boundary detection in oct images,” to be submitted to IEEE Trans. Med. Imaging.

Source: PubMed

3
Subscribe