2D and 3D multiphase active contours without edges based algorithms for simultaneous segmentation of retinal layers from OCT images

iv Acknowledgements I would like to thank my supervisor, Dr. Grady Rylander, for his guidance and support. I would also like to thank Bingqing Wang for his help throughout my research. Thanks as well to Dr. Mia Markey for serving as a reader for this thesis. A special thanks to my family and friends for their support throughout my college experience. Glaucoma is a common disease that is difficult to diagnosis early using only visual field tests. Current research indicates that determination of retinal nerve fiber layer thickness (RNFLT) can serve as an early indicator of glaucoma [1]. RNFLT is measured by segmenting non-invasive optical coherence tomography images. However, high speckle noise and presence of artifacts in the images cause traditional layer detection and segmentation methods to fail. Multiphase active contours segmentation methods utilize region intensity and shape terms to produce multiple continuous boundaries simultaneously in noisy environments. A 2D and 3D multiphase active contours based algorithm was created to segment synthetic and real human retina OCT images. The 2D multiphase algorithm segmented eight simultaneous layers with a 3.14% mean A-scan error rate per layer. The 3D approach performed qualitatively accurate segmentation of a 20 image stack vi simultaneously. In an artificial, high-noise image stack the incorporation of more pixels per layer allowed improved segmentation using the 3D algorithm over the 2D. These results indicate that 2D and 3D multiphase active contours algorithms can be used to accurately segment retina layers. With further development to reduce computation time and automate initialization, these algorithms could be used to provide close to real-time clinical retinal image segmentation.

[1]  Bernd Hamann,et al.  Segmentation of Three-dimensional Retinal Image Data , 2007, IEEE Transactions on Visualization and Computer Graphics.

[2]  Heinrich Niemann,et al.  Automated segmentation of the optic nerve head for diagnosis of glaucoma , 2005, Medical Image Anal..

[3]  Jitendra Malik,et al.  Scale-Space and Edge Detection Using Anisotropic Diffusion , 1990, IEEE Trans. Pattern Anal. Mach. Intell..

[4]  Antonin Chambolle,et al.  Implementation of an adaptive finite-element approximation of the Mumford-Shah functional , 2000, Numerische Mathematik.

[5]  Timothy F. Cootes,et al.  Active shape models , 1998 .

[6]  Rémi Ronfard,et al.  Region-based strategies for active contour models , 1994, International Journal of Computer Vision.

[7]  Tony F. Chan,et al.  Level set based shape prior segmentation , 2005, 2005 IEEE Computer Society Conference on Computer Vision and Pattern Recognition (CVPR'05).

[8]  Bingqing Wang,et al.  Retinal Nerve Fiber Layer Reflectance for Early Glaucoma Diagnosis , 2014, Journal of glaucoma.

[9]  Xiaodong Wu,et al.  Automated 3-D Intraretinal Layer Segmentation of Macular Spectral-Domain Optical Coherence Tomography Images , 2009, IEEE Transactions on Medical Imaging.

[10]  Milan Sonka,et al.  Segmentation of the Surfaces of the Retinal Layer from OCT Images , 2006, MICCAI.

[11]  Alan L. Yuille,et al.  Region Competition: Unifying Snakes, Region Growing, and Bayes/MDL for Multiband Image Segmentation , 1996, IEEE Trans. Pattern Anal. Mach. Intell..

[12]  Kaleem Siddiqi,et al.  Area and length minimizing flows for shape segmentation , 1998, IEEE Trans. Image Process..

[13]  Baba C. Vemuri,et al.  Shape Modeling with Front Propagation: A Level Set Approach , 1995, IEEE Trans. Pattern Anal. Mach. Intell..

[14]  John F. Canny,et al.  A Computational Approach to Edge Detection , 1986, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[15]  James A. Sethian,et al.  Level Set Methods and Fast Marching Methods: Evolving Interfaces in Computational Geometry, Fluid , 2012 .

[16]  M. Bashkansky,et al.  Statistics and reduction of speckle in optical coherence tomography. , 2000, Optics letters.

[17]  Jerry L. Prince,et al.  Snakes, shapes, and gradient vector flow , 1998, IEEE Trans. Image Process..

[18]  Bingqing Wang,et al.  Thickness, phase retardation, birefringence, and reflectance of the retinal nerve fiber layer in normal and glaucomatous non-human primates. , 2012, Investigative ophthalmology & visual science.

[19]  Guillermo Sapiro,et al.  3D active contours , 1996 .

[20]  T. Dawber,et al.  The Framingham Eye Study monograph: An ophthalmological and epidemiological study of cataract, glaucoma, diabetic retinopathy, macular degeneration, and visual acuity in a general population of 2631 adults, 1973-1975. , 1980, Survey of ophthalmology.

[21]  S. Osher,et al.  Level set methods: an overview and some recent results , 2001 .

[22]  Laurent D. Cohen,et al.  Avoiding local minima for deformable curves in image analysis , 1997 .

[23]  H. Quigley Number of people with glaucoma worldwide. , 1996, The British journal of ophthalmology.

[24]  Jayant Shah,et al.  A common framework for curve evolution, segmentation and anisotropic diffusion , 1996, Proceedings CVPR IEEE Computer Society Conference on Computer Vision and Pattern Recognition.

[25]  Guido Gerig,et al.  User-guided 3D active contour segmentation of anatomical structures: Significantly improved efficiency and reliability , 2006, NeuroImage.

[26]  J. Duker,et al.  Imaging of macular diseases with optical coherence tomography. , 1995, Ophthalmology.

[27]  D Marr,et al.  Theory of edge detection , 1979, Proceedings of the Royal Society of London. Series B. Biological Sciences.

[28]  Teresa C. Chen,et al.  Retinal nerve fiber layer thickness map determined from optical coherence tomography images. , 2005, Optics express.

[29]  A. Tannenbaum,et al.  Despeckling of medical ultrasound images , 2006, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[30]  Woo Kyung Moon,et al.  Segmentation of breast tumor in three-dimensional ultrasound images using three-dimensional discrete active contour model. , 2003, Ultrasound in medicine & biology.

[31]  J. Fujimoto,et al.  Optical coherence tomography of the human retina. , 1995, Archives of ophthalmology.

[32]  Tony F. Chan,et al.  Active contours without edges , 2001, IEEE Trans. Image Process..

[33]  J. Fujimoto,et al.  Optical coherence tomography: A new tool for glaucoma diagnosis , 1995, Current opinion in ophthalmology.

[34]  Mila Nikolova,et al.  Algorithms for Finding Global Minimizers of Image Segmentation and Denoising Models , 2006, SIAM J. Appl. Math..

[35]  Alexander Wong,et al.  Intra-retinal layer segmentation in optical coherence tomography images. , 2009, Optics express.

[36]  Tony F. Chan,et al.  A Multiphase Level Set Framework for Image Segmentation Using the Mumford and Shah Model , 2002, International Journal of Computer Vision.

[37]  Olivier D. Faugeras,et al.  Reconciling Distance Functions and Level Sets , 1999, Scale-Space.

[38]  Xavier Bresson,et al.  Fast Global Minimization of the Active Contour/Snake Model , 2007, Journal of Mathematical Imaging and Vision.

[39]  A. Govardhan,et al.  Active Contours and Image Segmentation: The Current State Of the Art , 2012 .

[40]  Michel Barlaud,et al.  Inward and outward curve evolution using level set method , 1999, Proceedings 1999 International Conference on Image Processing (Cat. 99CH36348).

[41]  S. Osher,et al.  Motion of multiple junctions: a level set approach , 1994 .

[42]  Eduardo Fernández,et al.  Webvision: The Organization of the Retina and Visual System , 1995 .

[43]  Demetri Terzopoulos,et al.  Snakes: Active contour models , 2004, International Journal of Computer Vision.

[44]  Ghassan Hamarneh,et al.  Intra-retinal Layer Segmentation in Optical Coherence Tomography Using an Active Contour Approach , 2009, MICCAI.

[45]  Bernd Hamann,et al.  Adaptation of a support vector machine algorithm for segmentation and visualization of retinal structures in volumetric optical coherence tomography data sets. , 2007, Journal of biomedical optics.

[46]  Delia Cabrera Fernandez,et al.  Delineating fluid-filled region boundaries in optical coherence tomography images of the retina , 2005, IEEE Transactions on Medical Imaging.

[47]  Ghassan Hamarneh,et al.  Adaptive Regularization for Image Segmentation Using Local Image Curvature Cues , 2010, ECCV.

[48]  Xiaodong Wu,et al.  Automated segmentation of intraretinal layers from macular optical coherence tomography images , 2007, SPIE Medical Imaging.

[49]  David A. Atchison,et al.  Optics of the Human Eye , 2023 .

[50]  G. Ripandelli,et al.  Optical coherence tomography. , 1998, Seminars in ophthalmology.

[51]  J. Hornegger,et al.  Retinal Nerve Fiber Layer Segmentation on FD-OCT Scans of Normal Subjects and Glaucoma Patients , 2010, Biomedical optics express.

[52]  Maciej Wojtkowski,et al.  High-definition and 3-dimensional imaging of macular pathologies with high-speed ultrahigh-resolution optical coherence tomography. , 2006, Ophthalmology.

[53]  Kim L. Boyer,et al.  Retinal thickness measurements from optical coherence tomography using a Markov boundary model , 2001, IEEE Transactions on Medical Imaging.

[54]  M. Brezinski Optical Coherence Tomography: Principles and Applications , 2006 .

[55]  D. Mumford,et al.  Optimal approximations by piecewise smooth functions and associated variational problems , 1989 .

[56]  Nikos Paragios,et al.  Shape Priors for Level Set Representations , 2002, ECCV.

[57]  Xiaodong Wu,et al.  Intraretinal Layer Segmentation of Macular Optical Coherence Tomography Images Using Optimal 3-D Graph Search , 2008, IEEE Transactions on Medical Imaging.

[58]  F. Medeiros,et al.  Detection of glaucoma progression with stratus OCT retinal nerve fiber layer, optic nerve head, and macular thickness measurements. , 2009, Investigative ophthalmology & visual science.

[59]  T. Chan,et al.  A Variational Level Set Approach to Multiphase Motion , 1996 .

[60]  Alex M. Andrew,et al.  Level Set Methods and Fast Marching Methods: Evolving Interfaces in Computational Geometry, Fluid Mechanics, Computer Vision, and Materials Science (2nd edition) , 2000 .

[61]  D. Zack,et al.  Retinal ganglion cell death in experimental glaucoma and after axotomy occurs by apoptosis. , 1995, Investigative ophthalmology & visual science.

[62]  Carmen A Puliafito,et al.  Automated detection of retinal layer structures on optical coherence tomography images. , 2005, Optics express.

[63]  Boris Hermann,et al.  Robust segmentation of intraretinal layers in the normal human fovea using a novel statistical model based on texture and shape analysis. , 2010, Optics express.

[64]  Carlos Vázquez,et al.  Multiregion competition: A level set extension of region competition to multiple region image partitioning , 2006, Comput. Vis. Image Underst..

[65]  D CohenLaurent On active contour models and balloons , 1991 .

[66]  G. Wollstein,et al.  Optical coherence tomography and histologic measurements of nerve fiber layer thickness in normal and glaucomatous monkey eyes. , 2007, Investigative ophthalmology & visual science.

[67]  R. Ansari,et al.  Thickness profiles of retinal layers by optical coherence tomography image segmentation. , 2008, American journal of ophthalmology.

[68]  H. Quigley,et al.  The number of people with glaucoma worldwide in 2010 and 2020 , 2006, British Journal of Ophthalmology.

[69]  Ronen Basri,et al.  Fast multiscale image segmentation , 2000, Proceedings IEEE Conference on Computer Vision and Pattern Recognition. CVPR 2000 (Cat. No.PR00662).

[70]  Thomas E. Milner,et al.  Automated method for RNFL segmentation in spectral domain OCT , 2008, SPIE BiOS.

[71]  Mark Sussman,et al.  An Efficient, Interface-Preserving Level Set Redistancing Algorithm and Its Application to Interfacial Incompressible Fluid Flow , 1999, SIAM J. Sci. Comput..

[72]  Joseph A Izatt,et al.  CORRELATION OF PATHOLOGIC FEATURES IN SPECTRAL DOMAIN OPTICAL COHERENCE TOMOGRAPHY WITH CONVENTIONAL RETINAL STUDIES , 2008, Retina.

[73]  J. Schmitt,et al.  Speckle in optical coherence tomography. , 1999, Journal of biomedical optics.

[74]  Tony F. Chan,et al.  Image segmentation using level sets and the piecewise-constant Mumford-Shah model , 2000 .

[75]  J. Sethian,et al.  FRONTS PROPAGATING WITH CURVATURE DEPENDENT SPEED: ALGORITHMS BASED ON HAMILTON-JACOB1 FORMULATIONS , 2003 .

[76]  Daniel Cremers,et al.  Shape priors in variational image segmentation: Convexity, Lipschitz continuity and globally optimal solutions , 2008, 2008 IEEE Conference on Computer Vision and Pattern Recognition.

[77]  Laurent D. Cohen,et al.  On active contour models and balloons , 1991, CVGIP Image Underst..

[78]  D Kucera,et al.  Segmentation of sequences of echocardiographic images using a simplified 3D active contour model with region-based external forces. , 1997, Computerized medical imaging and graphics : the official journal of the Computerized Medical Imaging Society.

[79]  L. Vese,et al.  A level set algorithm for minimizing the Mumford-Shah functional in image processing , 2001, Proceedings IEEE Workshop on Variational and Level Set Methods in Computer Vision.