Efficient computation of Hessian-based enhancement filters for tubular structures in 3D images

Abstract This work presents guidelines for a computationally efficient implementation of multiscale image filters based on eigenanalysis of the Hessian matrix, for the enhancement of tubular structures. Our focus is the application to 3D medical images of blood vessels. The method uses matrix trace, determinant and sign to discard voxels unlikely to belong to vessels, prior to the calculation of the Hessian eigenvalues. As example of time savings, we provide results obtained in four computed tomography datasets (300 × 300 × 300 voxels) containing coronary and pulmonary arteries. The test based on the Hessian trace avoided the computation of the eigenvalues in half of the voxels on average, while the test combining the Hessian determinant and sign eliminated up to 10% additional voxels. The actual time savings depend on the algorithm used to compute the eigenvalues for the remaining voxels. With a very fast algorithm using a closed-form solution, the computational time was reduced from 20.5 to 12.5 seconds per scale, but the time gained thanks to the more complex of the two tests was negligible. However, this fast algorithm is prone to numerical instabilities. Accurate computation of the eigenvalues requires the use of iterative or hybrid algorithms. In this case, both tests produce time savings and the computational time can be reduced by several minutes per scale.

[1]  Xavier Tizon,et al.  Algorithms for the analysis of 3D magnetic resonance angiography images , 2004 .

[2]  Anshuman Razdan,et al.  Case Study: Cellar Scaffold Extraction Using Crest Point for Volume Rendering , 2003, VisSym.

[3]  Sang Uk Lee,et al.  Fast multiscale vessel enhancement filtering , 2008, SPIE Medical Imaging.

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

[5]  Lihong V. Wang,et al.  In vivo volumetric imaging of subcutaneous microvasculature by photoacoustic microscopy. , 2006, Optics express.

[6]  Mp Mirjam Rubbens,et al.  Quantification of collagen orientation in 3D engineered tissue , 2007 .

[7]  Miguel Angel Aragón Calvo,et al.  Morphology and dynamics of the cosmic web , 2007 .

[8]  P. Saggau,et al.  Towards automatic reconstruction of dendrite morphology from live neurons , 2004, The 26th Annual International Conference of the IEEE Engineering in Medicine and Biology Society.

[9]  Oliver K. Smith,et al.  Eigenvalues of a symmetric 3 × 3 matrix , 1961, Commun. ACM.

[10]  Nicholas Ayache,et al.  Preprocessing : data selection , pseudo ECG III . 3 − D centerlines reconstruction , 2011 .

[11]  Nicholas Ayache,et al.  Model-Based Detection of Tubular Structures in 3D Images , 2000, Comput. Vis. Image Underst..

[12]  Jürgen Weese,et al.  Multi-scale line segmentation with automatic estimation of width, contrast and tangential direction in 2D and 3D medical images , 1997, CVRMed.

[13]  David,et al.  [Lecture Notes in Computer Science] Computer Vision â ECCV 2008 Volume 5305 || Three Dimensional Curvilinear Structure Detection Using Optimally Oriented Flux , 2008 .

[14]  OrkiszMaciej,et al.  Models, algorithms and applications in vascular image segmentation , 2008 .

[15]  Anil A. Bharath,et al.  Retinal Blood Vessel Segmentation by Means of Scale-Space Analysis and Region Growing , 1999, MICCAI.

[16]  Wiro J Niessen,et al.  Semiautomatic segmentation and stenosis quantification of 3D contrast‐enhanced MR angiograms of the internal carotid artery , 2004, Magnetic resonance in medicine.

[17]  Marcela Hernández Hoyos,et al.  Models, algorithms and applications in vascular image segmentation , 2008 .

[18]  Max W. K. Law,et al.  Three Dimensional Curvilinear Structure Detection Using Optimally Oriented Flux , 2008, ECCV.

[19]  J. Koenderink The structure of images , 2004, Biological Cybernetics.

[20]  Qiang Li,et al.  Selective enhancement filters for nodules, vessels, and airway walls in two- and three-dimensional CT scans. , 2003, Medical physics.

[21]  Joachim Kopp EFFICIENT NUMERICAL DIAGONALIZATION OF HERMITIAN 3 × 3 MATRICES , 2006 .

[22]  Milan Sonka,et al.  Automated segmentation of pulmonary vascular tree from 3D CT images , 2004, SPIE Medical Imaging.

[23]  Guido Gerig,et al.  Three-dimensional multi-scale line filter for segmentation and visualization of curvilinear structures in medical images , 1998, Medical Image Anal..

[24]  Maxime Descoteaux,et al.  Bone Enhancement Filtering: Application to Sinus Bone Segmentation and Simulation of Pituitary Surgery , 2005, MICCAI.

[25]  Sebastian P. M. Dries,et al.  Automated Planning of Scan Geometries in Spine MRI Scans , 2007, MICCAI.

[26]  Max A. Viergever,et al.  Scale and the differential structure of images , 1992, Image Vis. Comput..