A Stochastic Model for Automatic Extraction of 3D Neuronal Morphology

Tubular structures are frequently encountered in bio-medical images. The center-lines of these tubules provide an accurate representation of the topology of the structures. We introduce a stochastic Marked Point Process framework for fully automatic extraction of tubular structures requiring no user interaction or seed points for initialization. Our Marked Point Process model enables unsupervised network extraction by fitting a configuration of objects with globally optimal associated energy to the centreline of the arbors. For this purpose we propose special configurations of marked objects and an energy function well adapted for detection of 3D tubular branches. The optimization of the energy function is achieved by a stochastic, discrete-time multiple birth and death dynamics. Our method finds the centreline, local width and orientation of neuronal arbors and identifies critical nodes like bifurcations and terminals. The proposed model is tested on 3D light microscopy images from the DIADEM data set with promising results.

[1]  Pascal Fua,et al.  Reconstructing Geometrically Consistent Tree Structures from Noisy Images , 2010, MICCAI.

[2]  A. Baddeley,et al.  Stochastic geometry models in high-level vision , 1993 .

[3]  Tianxu Zhang,et al.  Marked Point Process for Vascular Tree Extraction on Angiogram , 2007, EMMCVPR.

[4]  Francis K. H. Quek,et al.  A review of vessel extraction techniques and algorithms , 2004, CSUR.

[5]  Erik Meijering,et al.  Neuron tracing in perspective , 2010, Cytometry. Part A : the journal of the International Society for Analytical Cytology.

[6]  Nicholas Ayache,et al.  Model-based multiscale detection of 3D vessels , 1998, Proceedings. 1998 IEEE Computer Society Conference on Computer Vision and Pattern Recognition (Cat. No.98CB36231).

[7]  Xavier Descombes,et al.  Author manuscript, published in "Journal of Mathematical Imaging and Vision (2009)" Object , 2008 .

[8]  Isabelle E. Magnin,et al.  Coronary tree extraction from X-ray angiograms using marked point processes , 2006, 3rd IEEE International Symposium on Biomedical Imaging: Nano to Macro, 2006..

[9]  Nassir Navab,et al.  Medical Image Computing and Computer-Assisted Intervention - MICCAI 2010, 13th International Conference, Beijing, China, September 20-24, 2010, Proceedings, Part III , 2010, MICCAI.

[10]  William M. Wells,et al.  Medical Image Computing and Computer-Assisted Intervention — MICCAI’98 , 1998, Lecture Notes in Computer Science.

[11]  Stephen R. Aylward,et al.  Initialization, noise, singularities, and scale in height ridge traversal for tubular object centerline extraction , 2002, IEEE Transactions on Medical Imaging.

[12]  P. Fua,et al.  Learning rotational features for filament detection , 2009, 2009 IEEE Conference on Computer Vision and Pattern Recognition.

[13]  Ju Lu,et al.  The DIADEM Data Sets: Representative Light Microscopy Images of Neuronal Morphology to Advance Automation of Digital Reconstructions , 2011, Neuroinformatics.

[14]  Xavier Descombes,et al.  A Point Process for Fully Automatic Road Network Detection in Satellite and Aerial Images , 2010 .

[15]  Horst Bischof,et al.  Multiscale Medialness for Robust Segmentation of 3 D Tubular Structures ∗ , 2005 .

[16]  Alejandro F. Frangi,et al.  Muliscale Vessel Enhancement Filtering , 1998, MICCAI.