Automated segmentation and area estimation of neural foramina with boundary regression model

Abstract Accurate segmentation and area estimation of neural foramina from both CT and MR images are essential to clinical diagnosis of neural foramina stenosis. Existing clinical routine, relying on physician's purely manual segmentation, becomes very tedious, laborious, and inefficient. Automated segmentation is highly desirable but faces big challenges from diverse boundary, local weak/no boundary, and intra/inter-modality intensity inhomogeneity. In this paper, a novel boundary regression segmentation framework is proposed for fully automated and multi-modal segmentation of neural foramina. It creatively formulates the segmentation task as a boundary regression problem which models a highly nonlinear mapping function from substantially diverse neural foramina images directly to desired object boundaries. By leveraging a seamless combination of multiple output support vector regression (MSVR) and multiple kernel learning (MKL), the proposed framework enables the domain knowledge learning in a holistic fashion which successfully handles the extreme diversity posing a tremendous challenge to conventional segmentation methods. The performance evaluation was conducted on a dataset including 912 MR images and 306 CT images collected from 152 subjects. Experimental results show that the proposed automated segmentation framework is highly consistent with physician with average DSI (dice similarity index) as high as 0.9005 (CT), 0.8984 (MR), 0.8935 (MR+CT) and BD (boundary distance) as low as 0.6393 mm (CT), 0.6586 mm (MR), 0.6881 mm (MR+CT). Based on this accurate automated segmentation, the estimated areas, highly correlated to their independent ground truth, have been achieved with correlation coefficient: 0.9154 (CT) and 0.8789 (MR). Hence, the proposed approach enables an efficient, accurate and convenient tool for clinical diagnosis of neural foramina stenosis.

[1]  Shuo Li,et al.  Graph Cuts with Invariant Object-Interaction Priors: Application to Intervertebral Disc Segmentation , 2011, IPMI.

[3]  Olga Veksler,et al.  Fast Approximate Energy Minimization via Graph Cuts , 2001, IEEE Trans. Pattern Anal. Mach. Intell..

[4]  Franco Postacchini,et al.  Stenosis of Lumbar Intervertebral Foramen: Anatomic Study on Predisposing Factors , 2002, Spine.

[5]  Cordelia Schmid,et al.  Evaluation of GIST descriptors for web-scale image search , 2009, CIVR '09.

[6]  Manik Varma,et al.  More generality in efficient multiple kernel learning , 2009, ICML '09.

[7]  Caroline Petitjean,et al.  A review of segmentation methods in short axis cardiac MR images , 2011, Medical Image Anal..

[8]  J. Evans,et al.  Lumbar Intervertebral Foramens: An in Vitro Study of Their Shape in Relation to Intervertebral Disc Pathology , 1991, Spine.

[9]  G. G. Stokes "J." , 1890, The New Yale Book of Quotations.

[10]  Fernando Pérez-Cruz,et al.  SVM multiregression for nonlinear channel estimation in multiple-input multiple-output systems , 2004, IEEE Transactions on Signal Processing.

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

[12]  Benoit M. Dawant,et al.  Morphometric analysis of white matter lesions in MR images: method and validation , 1994, IEEE Trans. Medical Imaging.

[13]  Constantin Schizas,et al.  Influence of anatomical variations on lumbar foraminal stenosis pathogenesis , 2015, European Spine Journal.

[14]  John A. Hipp,et al.  Assessment of Magnetic Resonance Imaging in the Diagnosis of Lumbar Spine Foraminal Stenosis—A Surgeon's Perspective , 2006, Journal of spinal disorders & techniques.

[15]  V. Haughton,et al.  Anatomic Changes of the Spinal Canal and Intervertebral Foramen Associated With Flexion‐Extension Movement , 1996, Spine.

[16]  Hakan Tuna,et al.  Morphometric analysis of the roots and neural foramina of the lumbar vertebrae. , 2006, Surgical neurology.

[17]  L. Costaridou,et al.  Combining 2D wavelet edge highlighting and 3D thresholding for lung segmentation in thin-slice CT. , 2007, The British journal of radiology.

[18]  Aaron Fenster,et al.  Spine Image Fusion Via Graph Cuts , 2013, IEEE Transactions on Biomedical Engineering.

[19]  R. Delamarter,et al.  Spinal Fusion in the United States: Analysis of Trends From 1998 to 2008 , 2012, Spine.

[20]  Takeo Kanade,et al.  Real-time topometric localization , 2012, 2012 IEEE International Conference on Robotics and Automation.

[21]  Suketaka Momoshima,et al.  Morphometric analysis of the lumbar intervertebral foramen in patients with degenerative lumbar scoliosis by multidetector-row computed tomography , 2012, European Spine Journal.

[22]  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.

[23]  Philip Sedgwick,et al.  Pearson’s correlation coefficient , 2012, BMJ : British Medical Journal.

[24]  W. Bartynski,et al.  Lumbar root compression in the lateral recess: MR imaging, conventional myelography, and CT myelography comparison with surgical confirmation. , 2003, AJNR. American journal of neuroradiology.

[25]  Khoa N. Le A mathematical approach to edge detection in hyperbolic-distributed and Gaussian-distributed pixel-intensity images using hyperbolic and Gaussian masks , 2011, Digit. Signal Process..

[26]  Zhongyi Hu,et al.  Multiple-output support vector regression with a firefly algorithm for interval-valued stock price index forecasting , 2014, Knowl. Based Syst..

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

[28]  Kurt Konolige,et al.  CenSurE: Center Surround Extremas for Realtime Feature Detection and Matching , 2008, ECCV.

[29]  Luis Alonso,et al.  Multioutput Support Vector Regression for Remote Sensing Biophysical Parameter Estimation , 2011, IEEE Geoscience and Remote Sensing Letters.

[30]  M. Adams,et al.  THE BIOMECHANICS OF BACK PAIN , 2003 .

[31]  Antonio Torralba,et al.  Modeling the Shape of the Scene: A Holistic Representation of the Spatial Envelope , 2001, International Journal of Computer Vision.

[32]  Gunnar Rätsch,et al.  Large Scale Multiple Kernel Learning , 2006, J. Mach. Learn. Res..

[33]  Bill Triggs,et al.  Histograms of oriented gradients for human detection , 2005, 2005 IEEE Computer Society Conference on Computer Vision and Pattern Recognition (CVPR'05).

[34]  Hakan Tuna,et al.  The lumbar roots and pedicles: a morphometric analysis and anatomical features , 2008, Journal of Clinical Neuroscience.

[35]  Manohar M. Panjabi,et al.  Dynamic Intervertebral Foramen Narrowing During Simulated Rear Impact , 2006, Spine.

[36]  P. Fieguth,et al.  Adaboost and Support Vector Machines for White Matter Lesion Segmentation in MR Images , 2005, 2005 IEEE Engineering in Medicine and Biology 27th Annual Conference.

[37]  Yi Wang,et al.  Automatic Left Ventricle Segmentation Using Iterative Thresholding and an Active Contour Model With Adaptation on Short-Axis Cardiac MRI , 2010, IEEE Transactions on Biomedical Engineering.

[38]  Fernando Pérez-Cruz,et al.  Multi-dimensional Function Approximation and Regression Estimation , 2002, ICANN.

[39]  Heung Sik Kang,et al.  A practical MRI grading system for lumbar foraminal stenosis. , 2010, AJR. American journal of roentgenology.

[40]  Olivier Chapelle,et al.  Training a Support Vector Machine in the Primal , 2007, Neural Computation.