IGT logo

National Center for Image Guided Therapy

The Publication Database hosted by SPL

All Publications | Upload | Advanced Search | Gallery View | Download Statistics | Help | Import | Log in

Tensor Kernels for Simultaneous Fiber Model Estimation and Tractography

Institution:
1Psychiatry Neuroimaging Laboratory, Brigham and Women's Hospital, Harvard Medical School, Boston, MA, USA. yogesh@bwh.harvard.edu
2Department of Electrical Engineering, University of Waterloo, Canada
3Laboratory for Mathematical Imaging, Brigham and Women’s Hospital, Harvard Medical School, Boston, MA, USA
4VA Boston Healthcare System, Brockton Division, Brockton, MA, USA
Publication Date:
Jul-2010
Journal:
Magn Reson Med
Volume Number:
64
Issue Number:
1
Pages:
138-148
Citation:
Magn Reson Med. 2010 Jul;64(1):138-48.
PubMed ID:
20572129
PMCID:
PMC3043656
Appears in Collections:
PNL, LMI, NA-MIC, NAC
Sponsors:
P41 RR13218-14/RR/NCRR NIH HHS/United States
R01 MH050740-15/MH/NIMH NIH HHS/United States
R01 MH074794-05/MH/NIMH NIH HHS/United States
R01 MH082918-03/MH/NIMH NIH HHS/United States
U54 EB005149-07/EB/NIBIB NIH HHS/United States
Generated Citation:
Rathi Y., Malcolm J.G., Michailovich O., Westin C-F., Shenton M.E., Bouix S. Tensor Kernels for Simultaneous Fiber Model Estimation and Tractography. Magn Reson Med. 2010 Jul;64(1):138-48. PMID: 20572129. PMCID: PMC3043656.
Downloaded: 636 times. [view map]
Paper: Download, View online
Export citation:

This paper proposes a novel framework for joint orientation distribution function estimation and tractography based on a new class of tensor kernels. Existing techniques estimate the local fiber orientation at each voxel independently so there is no running knowledge of confidence in the measured signal or estimated fiber orientation. In this work, fiber tracking is formulated as recursive estimation: at each step of tracing the fiber, the current estimate of the orientation distribution function is guided by the previous. To do this, second-and higher-order tensor-based kernels are employed. A weighted mixture of these tensor kernels is used for representing crossing and branching fiber structures. While tracing a fiber, the parameters of the mixture model are estimated based on the orientation distribution function at that location and a smoothness term that penalizes deviation from the previous estimate along the fiber direction. This ensures smooth estimation along the direction of propagation of the fiber. In synthetic experiments, using a mixture of two and three components it is shown that this approach improves the angular resolution at crossings. In vivo experiments using two and three components examine the corpus callosum and corticospinal tract and confirm the ability to trace through regions known to contain such crossing and branching.

Additional Material
1 File (34.35kB)
Rathi-MagResMed2010.jpg (34.35kB)