Efficient Recursive Algorithms for Computing the Mean Diffusion Tensor and Applications to DTI Segmentation

Computation of the mean of a collection of symmetric positive definite (SPD) matrices is a fundamental ingredient of many algorithms in diffusion tensor image (DTI) processing. For instance, in DTI segmentation, clustering, etc. In this paper, we present novel recursive algorithms for computing the mean of a set of diffusion tensors using several distance/divergence measures commonly used in DTI segmentation and clustering such as the Riemannian distance and symmetrized Kullback-Leibler divergence. To the best of our knowledge, to date, there are no recursive algorithms for computing the mean using these measures in literature. Recursive algorithms lead to a gain in computation time of several orders in magnitude over existing non-recursive algorithms. The key contributions of this paper are: (i) we present novel theoretical results on a recursive estimator for Karcher expectation in the space of SPD matrices, which in effect is a proof of the law of large numbers (with some restrictions) for the manifold of SPD matrices. (ii) We also present a recursive version of the symmetrized KL-divergence for computing the mean of a collection of SPD matrices. (iii) We present comparative timing results for computing the mean of a group of SPD matrices (diffusion tensors) depicting the gains in compute time using the proposed recursive algorithms over existing non-recursive counter parts. Finally, we also show results on gains in compute times obtained by applying these recursive algorithms to the task of DTI segmentation.

[1]  Zhizhou Wang,et al.  Tensor Field Segmentation Using Region Based Active Contour Model , 2004, ECCV.

[2]  Lasse Riis Østergaard,et al.  Active Surface Approach for Extraction of the Human Cerebral Cortex from MRI , 2006, MICCAI.

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

[4]  René Vidal,et al.  Segmenting Fiber Bundles in Diffusion Tensor Images , 2008, ECCV.

[5]  Hanqing Lu,et al.  Real-Time Visual Tracking via Incremental Covariance Model Update on Log-Euclidean Riemannian Manifold , 2009, 2009 Chinese Conference on Pattern Recognition.

[6]  Andrew J. Davison,et al.  Active Matching , 2008, ECCV.

[7]  Baba C. Vemuri,et al.  A unified framework for estimating diffusion tensors of any order with symmetric positive-definite constraints , 2010, 2010 IEEE International Symposium on Biomedical Imaging: From Nano to Macro.

[8]  Maher Moakher,et al.  Symmetric Positive-Definite Matrices: From Geometry to Applications and Visualization , 2006, Visualization and Processing of Tensor Fields.

[9]  Jiří Matas,et al.  Computer Vision - ECCV 2004 , 2004, Lecture Notes in Computer Science.

[10]  Carl-Fredrik Westin,et al.  Segmentation of Thalamic Nuclei from DTI Using Spectral Clustering , 2006, MICCAI.

[11]  Rachid Deriche,et al.  DTI segmentation by statistical surface evolution , 2006, IEEE Transactions on Medical Imaging.

[12]  Xavier Pennec,et al.  A Riemannian Framework for Tensor Computing , 2005, International Journal of Computer Vision.

[13]  Guillermo Sapiro,et al.  Geodesic Active Contours , 1995, International Journal of Computer Vision.

[14]  Bengt Jönsson,et al.  Restricted Diffusion in Cylindrical Geometry , 1995 .

[15]  Zhizhou Wang,et al.  DTI segmentation using an information theoretic tensor dissimilarity measure , 2005, IEEE Transactions on Medical Imaging.

[16]  N. Ayache,et al.  Log‐Euclidean metrics for fast and simple calculus on diffusion tensors , 2006, Magnetic resonance in medicine.

[17]  H. Karcher.  Riemannian center of mass and mollifier smoothing , 1977 .

[18]  Baba C. Vemuri,et al.  Tensor Splines for Interpolation and Approximation of DT-MRI With Applications to Segmentation of Isolated Rat Hippocampi , 2007, IEEE Transactions on Medical Imaging.

[19]  W. Ballmann.  Manifolds of non positive curvature , 1985 .

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

[21]  Ghassan Hamarneh,et al.  DT-MRI segmentation using graph cuts , 2007, SPIE Medical Imaging.

[22]  J. Weickert,et al.  Level-Set Methods for Tensor-Valued Images , 2003 .

[23]  Anthony J. Yezzi,et al.  Curve evolution implementation of the Mumford-Shah functional for image segmentation, denoising, interpolation, and magnification , 2001, IEEE Trans. Image Process..

[24]  Frank Nielsen,et al.  This article has been accepted for publication in a future issue of this journal, but has not been fully edited. Content may change prior to final publication. 1 Total Bregman Divergence and its Applications to DTI Analysis , 2022 .

[25]  Hans Hagen,et al.  Visualization and Processing of Tensor Fields , 2014 .