Tensor Splines for Interpolation and Approximation of DT-MRI With Applications to Segmentation of Isolated Rat Hippocampi

In this paper, we present novel algorithms for statistically robust interpolation and approximation of diffusion tensors-which are symmetric positive definite (SPD) matrices-and use them in developing a significant extension to an existing probabilistic algorithm for scalar field segmentation, in order to segment diffusion tensor magnetic resonance imaging (DT-MRI) datasets. Using the Riemannian metric on the space of SPD matrices, we present a novel and robust higher order (cubic) continuous tensor product of -splines algorithm to approximate the SPD diffusion tensor fields. The resulting approximations are appropriately dubbed tensor splines. Next, we segment the diffusion tensor field by jointly estimating the label (assigned to each voxel) field, which is modeled by a Gauss Markov measure field (GMMF) and the parameters of each smooth tensor spline model representing the labeled regions. Results of interpolation, approximation, and segmentation are presented for synthetic data and real diffusion tensor fields from an isolated rat hippocampus, along with validation. We also present comparisons of our algorithms with existing methods and show significantly improved results in the presence of noise as well as outliers.

[1]  O. Faugeras,et al.  Statistics on Multivariate Normal Distributions: A Geometric Approach and its Application to Diffusion Tensor MRI , 2004 .

[2]  James C. Gee,et al.  Spatial transformations of diffusion tensor magnetic resonance images , 2001, IEEE Transactions on Medical Imaging.

[3]  Carl-Fredrik Westin,et al.  Riemannian Mean Curvature Flow , 2005, ISVC.

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

[5]  J. E. Tanner,et al.  Spin diffusion measurements : spin echoes in the presence of a time-dependent field gradient , 1965 .

[6]  Geoffrey E. Hinton,et al.  A View of the Em Algorithm that Justifies Incremental, Sparse, and other Variants , 1998, Learning in Graphical Models.

[7]  R. Clark,et al.  The medial temporal lobe. , 2004, Annual review of neuroscience.

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

[9]  S. Helgason.  LIE GROUPS AND SYMMETRIC SPACES. , 1968 .

[10]  Mariano Rivera,et al.  Entropy Controlled Gauss-Markov Random Measure Field Models for Early Vision , 2005, VLSM.

[11]  P. Thomas Fletcher,et al.  Principal Geodesic Analysis on Symmetric Spaces: Statistics of Diffusion Tensors , 2004, ECCV Workshops CVAMIA and MMBIA.

[12]  Ricardo Insausti,et al.  CHAPTER 23 – Hippocampal Formation , 2004 .

[13]  S. Helgason.  Differential Geometry, Lie Groups, and Symmetric Spaces , 1978 .

[14]  Guillermo Sapiro,et al.  Robust anisotropic diffusion , 1998, IEEE Trans. Image Process..

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

[16]  Joachim Weickert,et al.  Tensor Field Interpolation with PDEs , 2006, Visualization and Processing of Tensor Fields.

[17]  C. D. Boor,et al.  On Calculating B-splines , 1972 .

[18]  A. Terras.  Harmonic Analysis on Symmetric Spaces and Applications I , 1985 .

[19]  David E. Breen,et al.  Level set modeling and segmentation of diffusion tensor magnetic resonance imaging brain data , 2003, J. Electronic Imaging.

[20]  Mariano Rivera,et al.  Gauss-Markov Measure Field Models for Low-Level Vision , 2001, IEEE Trans. Pattern Anal. Mach. Intell..

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

[22]  Zhizhou Wang,et al.  An affine invariant tensor dissimilarity measure and its applications to tensor-valued image segmentation , 2004, Proceedings of the 2004 IEEE Computer Society Conference on Computer Vision and Pattern Recognition, 2004. CVPR 2004..

[23]  D. Amaral,et al.  The Hippocampal Formation , 2009 .

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

[25]  David E. Breen,et al.  Level Set Segmentation and Modeling of DT-MRI human brain data , 2003 .

[26]  Zhizhou Wang,et al.  A Constrained Variational Principle for Direct Estimation and Smoothing of the Diffusion Tensor Field from DWI , 2003, IPMI.

[27]  Jun Zhang.  The mean field theory in EM procedures for Markov random fields , 1992, IEEE Trans. Signal Process..

[28]  Nicholas Ayache,et al.  Fast and Simple Calculus on Tensors in the Log-Euclidean Framework , 2005, MICCAI.

[29]  Juan Ruiz-Alzola,et al.  Tensor Field Regularization using Normalized Convolution and Markov Random Fields in a Bayesian Framework , 2006, Visualization and Processing of Tensor Fields.

[30]  P. Basser,et al.  A continuous tensor field approximation of discrete DT-MRI data for extracting microstructural and architectural features of tissue. , 2002, Journal of magnetic resonance.

[31]  P. Basser,et al.  Estimation of the effective self-diffusion tensor from the NMR spin echo. , 1994, Journal of magnetic resonance. Series B.

[32]  Michael A. King,et al.  Structural insights from high-resolution diffusion tensor imaging and tractography of the isolated rat hippocampus , 2006, NeuroImage.

[33]  Maher Moakher,et al.  A rigorous framework for diffusion tensor calculus , 2005, Magnetic resonance in medicine.

[34]  John F. Hughes,et al.  Smooth interpolation of orientations with angular velocity constraints using quaternions , 1992, SIGGRAPH.