This Article 
   
 Share 
   
 Bibliographic References 
   
 Add to: 
 
Digg
Furl
Spurl
Blink
Simpy
Google
Del.icio.us
Y!MyWeb
 
 Search 
   
Estimating Crossing Fibers: A Tensor Decomposition Approach
November/December 2008 (vol. 14 no. 6)
pp. 1635-1642
Thomas Schultz, MPI Informatik
Hans-Peter Seidel, MPI Informatik
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] A. L. Alexander, K. M. Hasan, M. Lazar, J. S. Tsuruda, and D. L. Parker, Analysis of partial volume effects in diffusion-tensor MRI. Magnetic Resonance in Medicine 45: 770–780, 2001.
[2] D. C. Alexander, G. J. Barker, and S. R. Arridge, Detection and modeling of non-gaussian apparent diffusion coefficient profiles in human brain data. Magnetic Resonance in Medicine, 48: 331–340, 2002.
[3] A. W. Anderson, Measurement of fiber orientation distributions using high angular resolution diffusion imaging. Magnetic Resonance in Medicine, 54 (5): 1194–1206, 2005.
[4] A. Anwander, M. Descoteaux, and R. Deriche, Probabilistic Q-Ball tractography solves crossings of callosal fibers. In Proc. Human Brain Mapping, page 342, 2007.
[5] P. J. Basser, J. Mattiello, and D. L. Bihan, Estimation of the effective self-diffusion tensor from the NMR spin echo. Journal of Magnetic Resonance, B(103): 247–254, 1994.
[6] P. J. Basser and S. Pajevic, Spectral decomposition of a 4th-order covariance tensor: Applications to diffusion tensor MRI. Signal Processing, 87: 220–236, 2007.
[7] P. J. Basser and C. Pierpaoli, Microstructural and physiological features of tissues elucidated by quantitative-diffusion-tensor MRI. Journal of Magnetic Resonance, B(111): 209–219, 1996.
[8] T. Behrens, H. Johansen-Berg, S. Jbabdi, M. Rushworth, and M. Woolrich, Probabilistic diffusion tractography with multiple fibre orientations: What can we gain? NeuroImage, 34: 144–155, 2007.
[9] Ø. Bergmann, G. Kindlmann, A. Lundervold, and C.-F. Westin, Diffusion k-tensor estimation from Q-Ball imaging using discretized principal axes. In Proc. Medical Image Computing and Computer-Assisted Intervention (MICCAI'06), volume 4191 of LNCS, pages 268–275, 2006.
[10] P. T. Callaghan, C. D. Eccles, and Y. Xia, NMR microscopy of dynamic displacements: k-space and q-space imaging. Journal of Physics E, 21 (8): 820–822, 1988.
[11] P. Comon, G. Golub, L.-H. Lim, and B. Mourrain, Symmetric tensors and symmetric tensor rank. Technical Report SCCM-06-02, Stanford Scientific Computing and Computational Mathematics (SCCM), 2006. To appear in: SIAM Journal on Matrix Analysis and Applications.
[12] P. Comon and B. Mourrain, Decomposition of quantics in sums of powers of linear forms. Signal Processing, 53 (2): 96–107, September 1996.
[13] J. Dauguet, S. Peled, V. Berezovskii, T. Delzescaux, S. K. Warfield, R. Born, and C.-F. Westin, 3D histological reconstruction of fiber tracts and direct comparison with diffusion tensor MRI tractography. In R. Larsen, M. Nielsen, and J. Sporring, editors, Medical Image Computing and Computer-Assisted Intervention (MICCAI'06), volume 4190 of LNCS, pages 109–116. Springer, 2006.
[14] L. De Lathauwer, B. De Moor, and J. Vandewalle, On the best rank-1 and rank-(r1, r2, , rn) approximation of higher-order tensors. SIAM Journal on Matrix Analysis and Applications, 21 (4): 1324–1342, 2000.
[15] V. de Silva and L.-H. Lim, Tensor rank and the ill-posedness of the best low-rank approximation problem. Technical Report SCCM-06-06, Stanford Scientific Computing and Computational Mathematics (SCCM), 2006. To appear in: SIAM Journal on Matrix Analysis and Applications.
[16] M. Descoteaux, E. Angelino, S. Fitzgibbons, and R. Deriche, Regularized, fast, and robust analytical Q-Ball imaging. Magnetic Resonance in Medicine, 58: 497–510, 2007.
[17] L. R. Frank, Characterization of anisotropy in high angular resolution diffusion-weighted MRI. Magnetic Resonance in Medicine, 47: 1083–1099, 2002.
[18] P. Hagmann, T. G. Reese, W.-Y. I. Tseng, R. Meuli, J.-P. Thiran, and V. J. Wedeen, Diffusion spectrum imaging tractography in complex cerebral white matter: an investigation of the centrum semiovale. In Proc. International Society of Magnetic Resonance in Medicine (ISMRM'04), page 623, 2004.
[19] C. P. Hess, P. Mukherjee, E. T. Han, D. Xu, and D. B. Vigneron, Q-Ball reconstruction of multimodal fiber orientations using the spherical harmonic basis. Magnetic Resonance in Medicine, 56: 104–117, 2006.
[20] F. L. Hitchcock, Multiple invariants and generalized rank of a p-way matrix or tensor. Journal of Mathematics and Physics, 7 (1): 39–79, 1927.
[21] M. Hlawitschka and G. Scheuermann, HOT-lines: Tracking lines in higher order tensor fields. In C. Silva, E. Gröller, and H. Rushmeier, editors, Proceedings of IEEE Visualization 2005, pages 27–34, 2005.
[22] M. Hlawitschka, G. Scheuermann, A. Anwander, T. Knösche, M. Tittgemeyer, and B. Hamann, Tensor lines in tensor fields of arbitrary order. In G. Bebis et al., editor, Advances in Visual Computing (Proc. ISVC'07), volume 4841 of LNCS, pages 341–350. Springer, 2007.
[23] D. K. Jones, M. A. Horsfield, and A. Simmons, Optimal strategies for measuring diffusion in anisotropic systems by magnetic resonance imaging. Magnetic Resonance in Medicine, 42: 515–525, 1999.
[24] E. Kaden, T. R. Knösche, and A. Anwander, Parametric spherical deconvolution: Inferring anatomical connectivity using diffusion MR imaging. NeuroImage, 37: 474–488, 2007.
[25] G. Kindlmann, D. Ennis, R. Whitaker, and C.-F. Westin, Diffusion tensor analysis with invariant gradients and rotation tangents. IEEE Transactions on Medical Imaging, 26 (11): 1483–1499, 2007.
[26] E. Kofidis and P. A. Regalia, On the best rank-1 approximation of higher-order supersymmetric tensors. SIAM Journal on Matrix Analysis and Applications, 23 (3): 863–884, 2002.
[27] B. Kreher, J. Schneider, I. Mader, E. Martin, J. Hennig, and K. Il'yasov, Multitensor approach for analysis and tracking of complex fiber configurations. Magnetic Resonance in Medicine, 54: 1216–1225, 2005.
[28] D. Le Bihan, E. Breton, D. Lallemand, P. Grenier, E. Cabanis, and M. Laval-Jeantet, MR imaging of intravoxel incoherent motions: Application to diffusion and perfusion in neurologic disorders. Radiology, 161 (2): 401–407, 1986.
[29] C.-P. Lin, W.-Y. I. Tseng, H.-C. Cheng, and J.-H. Chen, Validation of diffusion tensor magnetic resonance axonal fiber imaging with registered manganese-enhanced optic tracts. NeuroImage, 14 (5): 1035–1047, 2001.
[30] E. Özarslan and T. Mareci, Generalized diffusion tensor imaging and analytical relationships between diffusion tensor imaging and high angular resolution diffusion imaging. Magnetic Resonance in Medicine, 50: 955–965, 2003.
[31] E. Özarslan, B. C. Vemuri, and T. H. Mareci, Generalized scalar measures for diffusion MRI using trace, variance, and entropy. Magnetic Resonance in Medicine, 53: 866–876, 2005.
[32] S. Peled, O. Friman, F. Jolesz, and C.-F. Westin, Geometrically constrained two-tensor model for crossing tracts in DWI. Magnetic Resonance Imaging, 24 (9): 1263–1270, 2006.
[33] K. K. Seunarine, P. A. Cook, M. G. Hall, K. V. Embleton, G. J. M. Parker, and D. C. Alexander, Exploiting peak anisotropy for tracking through complex structures. In Proc. IEEE International Conference on Computer Vision: Workshop Mathematical Methods in Biomedical Image Analysis (MMBIA), pages 1–8, 2007.
[34] E. Stejskal and J. Tanner, Spin diffusion measurements: Spin echoes in the presence of a time-dependent field gradient. Journal of Chemical Physics, 42: 288–292, 1965.
[35] J.-D. Tournier, F. Calamante, and A. Connelly, Robust determination of the fibre orientation distribution in diffusion MRI: Non-negativity constrained super-resolved spherical deconvolution. NeuroImage, 35: 1459–1472, 2007.
[36] J.-D. Tournier, F. Calamante, D. G. Gadian, and A. Connelly, Direct estimation of the fiber orientation density function from diffusion-weighted MRI data using spherical deconvolution. NeuroImage, 23: 1176–1185, 2004.
[37] D. S. Tuch, Q-Ball imaging. Magnetic Resonance in Medicine, 52: 1358–1372, 2004.
[38] D. S. Tuch, T. G. Reese, M. R. Wiegell, N. Makris, J. W. Belliveau, and V. J. Wedeen, High angular resolution diffusion imaging reveals intravoxel white matter fiber heterogeneity. Magnetic Resonance in Medicine, 48: 577–582, 2002.
[39] S. Wakana, H. Jiang, L. M. Nagae-Poetscher, P. C. M. van Zijl, and S. Mori, Fiber tract-based atlas of human white matter anatomy. Radiology, 230: 77–87, 2004.
[40] Y. Wang and L. Qi, On the successive supersymmetric rank-1 decomposition of higher-order supersymmetric tensors. Numerical Linear Algebra with Applications, 14: 503–519, 2007.
[41] V. J. Wedeen, T. G. Reese, D. S. Tuch, M. R. Weigel, J.-G. Dou, R. M. Weisskoff, and D. Chessler, Mapping fiber orientation spectra in cerebral white matter with fourier-transform diffusion MRI. In Proc. Intl. Soc. Magn. Res. Med., volume 8, page 82, 2000.
[42] W. Zhan and Y. Yang, How accurately can the diffusion profiles indicate multiple fiber orientations? a study on general fiber crossings in diffusion MRI. Journal of Magnetic Resonance, 183: 193–202, 2006.

Index Terms:
Index Terms—DW-MRI, Q-Ball, spherical deconvolution, fiber tracking, higher-order tensor, tensor decomposition.
Citation:
Thomas Schultz, Hans-Peter Seidel, "Estimating Crossing Fibers: A Tensor Decomposition Approach," IEEE Transactions on Visualization and Computer Graphics, vol. 14, no. 6, pp. 1635-1642, Nov.-Dec. 2008, doi:10.1109/TVCG.2008.128
Usage of this product signifies your acceptance of the Terms of Use.