Neural tractography using an unscented Kalman filter.
We describe a technique to simultaneously estimate a local neural fiber model and trace out its path. Existing techniques estimate the local fiber orientation at each voxel independently so there is no running knowledge of confidence in the estimated fiber model. We formulate fiber tracking as recursive estimation: at each step of tracing the fiber, the current estimate is guided by the previous. To do this we model the signal as a mixture of Gaussian tensors and perform tractography within a filter framework. Starting from a seed point, each fiber is traced to its termination using an unscented Kalman filter to simultaneously fit the local model and propagate in the most consistent direction. Despite the presence of noise and uncertainty, this provides a causal estimate of the local structure at each point along the fiber. Synthetic experiments demonstrate that this approach reduces signal reconstruction error and significantly improves the angular resolution at crossings and branchings. In vivo experiments confirm the ability to trace out fibers in areas known to contain such crossing and branching while providing inherent path regularization.
Malcolm
James G
JG
Psychiatry Neuroimaging Laboratory, Harvard Medical School, Boston, MA, USA. malcolm@bwh.harvard.edu
Shenton
Martha E
ME
Rathi
Yogesh
Y
D000465:Algorithms; D001185:Artificial Intelligence; D001921:Brain / Q000033:anatomy & histology*; D003198:Computer Simulation; D038524:Diffusion Magnetic Resonance Imaging / Q000379:methods*; D006801:Humans; D007089:Image Enhancement / Q000379:methods; D007090:Image Interpretation, Computer-Assisted / Q000379:methods*; D021621:Imaging, Three-Dimensional / Q000379:methods*; D008959:Models, Neurological; D015233:Models, Statistical; D009413:Nerve Fibers, Myelinated / Q000648:ultrastructure*; D010363:Pattern Recognition, Automated / Q000379:methods*; D015203:Reproducibility of Results; D012680:Sensitivity and Specificity; D012815:Signal Processing, Computer-Assisted
Med Image Comput Comput Assist Interv. 2008;11(Pt 1):305-12
18979761
Inf Process Med Imaging. 2003 Jul;18:672-83
15344497
Med Image Anal. 2009 Jun;13(3):432-44
19269242
Magn Reson Med. 2005 Nov;54(5):1216-25
16200554
Med Image Comput Comput Assist Interv. 2008;11(Pt 1):1-8
18979725
Inf Process Med Imaging. 2009;21:139-50
19694259
Magn Reson Med. 2001 May;45(5):770-80
11323803
Magn Reson Med. 2000 Oct;44(4):625-32
11025519
Neuroimage. 2002 Jun;16(2):378-88
12030823
IEEE Trans Vis Comput Graph. 2008 Nov-Dec;14(6):1635-42
18989020
Magn Reson Imaging. 2006 Nov;24(9):1263-70
17071347
Magn Reson Med. 2007 Sep;58(3):497-510
17763358
Neuroimage. 2004 Nov;23(3):1176-85
15528117
Neuroimage. 2006 Jul 1;31(3):1086-103
16546404
IEEE Trans Image Process. 2010 Feb;19(2):461-77
19887312
Philos Trans R Soc Lond B Biol Sci. 2005 May 29;360(1457):893-902
16087434
Neuroimage. 2009 Aug;47 Suppl 2:T98-106
18657622
Ann N Y Acad Sci. 2005 Dec;1064:113-33
16394152
Magn Reson Med. 2005 Nov;54(5):1194-206
16161109
Magn Reson Med. 2002 Oct;48(4):577-82
12353272
Magn Reson Med. 2002 Aug;48(2):331-40
12210942
Magn Reson Med. 2004 Dec;52(6):1358-72
15562495
Magn Reson Med. 2003 Nov;50(5):1077-88
14587019
Neuroimage. 2007 Jan 1;34(1):144-55
17070705
Magn Reson Med. 2002 Jun;47(6):1083-99
12111955
IEEE Trans Med Imaging. 2007 Nov;26(11):1464-71
18041262
Magn Reson Med. 2005 Dec;54(6):1480-9
16265642
Neuroimage. 2005 Oct 1;27(4):725-36
16111897
NMR Biomed. 2002 Nov-Dec;15(7-8):456-67
12489095
Proc IEEE Int Symp Biomed Imaging. 2006;2006:65-68
19759891
Neuroimage. 2007 Aug 15;37(2):474-88
17596967
Med Image Comput Comput Assist Interv. 2007;10(Pt 2):144-52
18044563
J Magn Reson. 2006 Dec;183(2):193-202
16963296
Magn Reson Med. 2006 Jul;56(1):104-17
16755539
Inf Process Med Imaging. 2003 Jul;18:684-95
15344498
Med Image Anal. 2008 Oct;12(5):527-34
18664412
IEEE Trans Med Imaging. 2006 Aug;25(8):965-78
16894991