Estimating Crossing Fibers: A Tensor Decomposition Approach

Diffusion weighted magnetic resonance imaging is a unique tool for non-invasive investigation of major nerve fiber tracts. Since the popular diffusion tensor (DT-MRI) model is limited to voxels with a single fiber direction, a number of high angular resolution techniques have been proposed to provide information about more diverse fiber distributions. Two such approaches are Q-Ball imaging and spherical deconvolution, which produce orientation distribution functions (ODFs) on the sphere. For analysis and visualization, the maxima of these functions have been used as principal directions, even though the results are known to be biased in case of crossing fiber tracts. In this paper, we present a more reliable technique for extracting discrete orientations from continuous ODFs, which is based on decomposing their higher-order tensor representation into an isotropic component, several rank-1 terms, and a small residual. Comparing to ground truth in synthetic data shows that the novel method reduces bias and reliably reconstructs crossing fibers which are not resolved as individual maxima in the ODF We present results on both Q-Ball and spherical deconvolution data and demonstrate that the estimated directions allow for plausible fiber tracking in a real data set.

[1]  Alan Connelly,et al.  Robust determination of the fibre orientation distribution in diffusion MRI: Non-negativity constrained super-resolved spherical deconvolution , 2007, NeuroImage.

[2]  Carl-Fredrik Westin,et al.  3D Histological Reconstruction of Fiber Tracts and Direct Comparison with Diffusion Tensor MRI Tractography , 2006, MICCAI.

[3]  Phillip A. Regalia,et al.  On the Best Rank-1 Approximation of Higher-Order Supersymmetric Tensors , 2001, SIAM J. Matrix Anal. Appl..

[4]  Pierre Comon,et al.  Decomposition of quantics in sums of powers of linear forms , 1996, Signal Process..

[5]  L. Lathauwer,et al.  On the Best Rank-1 and Rank-( , 2004 .

[6]  Philip A. Cook,et al.  Exploiting peak anisotropy for tracking through complex structures , 2007, 2007 IEEE 11th International Conference on Computer Vision.

[7]  R. Deriche,et al.  Regularized, fast, and robust analytical Q‐ball imaging , 2007, Magnetic resonance in medicine.

[8]  Alan Connelly,et al.  Direct estimation of the fiber orientation density function from diffusion-weighted MRI data using spherical deconvolution , 2004, NeuroImage.

[9]  S. Arridge,et al.  Detection and modeling of non‐Gaussian apparent diffusion coefficient profiles in human brain data , 2002, Magnetic resonance in medicine.

[10]  Gerik Scheuermann,et al.  HOT-lines: tracking lines in higher order tensor fields , 2005, VIS 05. IEEE Visualization, 2005..

[11]  Peter J. Basser,et al.  Spectral decomposition of a 4th-order covariance tensor: Applications to diffusion tensor MRI , 2007, Signal Process..

[12]  T. Mareci,et al.  Generalized diffusion tensor imaging and analytical relationships between diffusion tensor imaging and high angular resolution diffusion imaging , 2003, Magnetic resonance in medicine.

[13]  Thomas R. Knösche,et al.  Parametric spherical deconvolution: Inferring anatomical connectivity using diffusion MR imaging , 2007, NeuroImage.

[14]  D. Parker,et al.  Analysis of partial volume effects in diffusion‐tensor MRI , 2001, Magnetic resonance in medicine.

[15]  P. Callaghan,et al.  RAPID COMMUNICATION: NMR microscopy of dynamic displacements: k-space and q-space imaging , 1988 .

[16]  Bernd Hamann,et al.  Tensor Lines in Tensor Fields of Arbitrary Order , 2007, ISVC.

[17]  B. Vemuri,et al.  Generalized scalar measures for diffusion MRI using trace, variance, and entropy , 2005, Magnetic resonance in medicine.

[18]  Yihong Yang,et al.  How accurately can the diffusion profiles indicate multiple fiber orientations? A study on general fiber crossings in diffusion MRI. , 2006, Journal of magnetic resonance.

[19]  N. Makris,et al.  High angular resolution diffusion imaging reveals intravoxel white matter fiber heterogeneity , 2002, Magnetic resonance in medicine.

[20]  D. Tuch.  Q‐ball imaging , 2004, Magnetic resonance in medicine.

[21]  Ching-Po Lin,et al.  Validation of Diffusion Tensor Magnetic Resonance Axonal Fiber Imaging with Registered Manganese-Enhanced Optic Tracts , 2001, NeuroImage.

[22]  Mark W. Woolrich,et al.  Probabilistic diffusion tractography with multiple fibre orientations: What can we gain? , 2007, NeuroImage.

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

[24]  Rachid Deriche,et al.  Probabilistic Q-ball tractography solves crossings of callosal fibers , 2007 .

[25]  Gene H. Golub,et al.  Symmetric Tensors and Symmetric Tensor Rank , 2008, SIAM J. Matrix Anal. Appl..

[26]  Liqun Qi,et al.  On the successive supersymmetric rank‐1 decomposition of higher‐order supersymmetric tensors , 2007, Numer. Linear Algebra Appl..

[27]  P. Basser,et al.  Microstructural and physiological features of tissues elucidated by quantitative-diffusion-tensor MRI. , 1996, Journal of magnetic resonance. Series B.

[28]  Duan Xu,et al.  Q‐ball reconstruction of multimodal fiber orientations using the spherical harmonic basis , 2006, Magnetic resonance in medicine.

[29]  Carl-Fredrik Westin,et al.  Geometrically constrained two-tensor model for crossing tracts in DWI. , 2006, Magnetic resonance imaging.

[30]  Carl-Fredrik Westin,et al.  Diffusion Tensor Analysis With Invariant Gradients and Rotation Tangents , 2007, IEEE Transactions on Medical Imaging.

[31]  S. Wakana,et al.  Fiber tract-based atlas of human white matter anatomy. , 2004, Radiology.

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

[33]  若菜 勢津.  Fiber tract-based atlas of human white matter anatomy , 2006 .

[34]  V. Wedeen,et al.  Mapping fiber orientation spectra in cerebral white matter with Fourier-transform diffusion MRI , 2000 .

[35]  Vin de Silva,et al.  Tensor rank and the ill-posedness of the best low-rank approximation problem , 2006, math/0607647.

[36]  Carl-Fredrik Westin,et al.  Diffusion k-tensor Estimation from Q-ball Imaging Using Discretized Principal Axes , 2006, MICCAI.

[37]  J. Thiran,et al.  Diffusion Spectrum Imaging tractography in complex cerebral white matter: an investigation of the centrum semiovale , 2004 .

[38]  P. Grenier,et al.  MR imaging of intravoxel incoherent motions: application to diffusion and perfusion in neurologic disorders. , 1986, Radiology.

[39]  F. L. Hitchcock.  Multiple Invariants and Generalized Rank of a P‐Way Matrix or Tensor , 1928 .

[40]  Gene H. Golub,et al.  Genericity And Rank Deficiency Of High Order Symmetric Tensors , 2006, 2006 IEEE International Conference on Acoustics Speech and Signal Processing Proceedings.

[41]  A. Anderson.  Measurement of fiber orientation distributions using high angular resolution diffusion imaging , 2005, Magnetic resonance in medicine.

[42]  B W Kreher,et al.  Multitensor approach for analysis and tracking of complex fiber configurations , 2005, Magnetic resonance in medicine.

[43]  L. Frank.  Characterization of anisotropy in high angular resolution diffusion‐weighted MRI , 2002, Magnetic resonance in medicine.

[44]  M. Horsfield,et al.  Optimal strategies for measuring diffusion in anisotropic systems by magnetic resonance imaging , 1999, Magnetic resonance in medicine.