title: Neural tractography using an unscented Kalman filter.
abstract: 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.
authors.affiliation: Psychiatry Neuroimaging Laboratory, Harvard Medical School, Boston, MA, USA. malcolm@bwh.harvard.edu
authors.forename: James G
authors.initials: JG
authors.lastname: Malcolm
authors.forename: Martha E
authors.initials: ME
authors.lastname: Shenton
authors.forename: Yogesh
authors.initials: Y
authors.lastname: Rathi
mesh_terms: 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
references.citation: Med Image Comput Comput Assist Interv. 2008;11(Pt 1):305-12
references.pmid: 18979761
references.citation: Inf Process Med Imaging. 2003 Jul;18:672-83
references.pmid: 15344497
references.citation: Med Image Anal. 2009 Jun;13(3):432-44
references.pmid: 19269242
references.citation: Magn Reson Med. 2005 Nov;54(5):1216-25
references.pmid: 16200554
references.citation: Med Image Comput Comput Assist Interv. 2008;11(Pt 1):1-8
references.pmid: 18979725
references.citation: Inf Process Med Imaging. 2009;21:139-50
references.pmid: 19694259
references.citation: Magn Reson Med. 2001 May;45(5):770-80
references.pmid: 11323803
references.citation: Magn Reson Med. 2000 Oct;44(4):625-32
references.pmid: 11025519
references.citation: Neuroimage. 2002 Jun;16(2):378-88
references.pmid: 12030823
references.citation: IEEE Trans Vis Comput Graph. 2008 Nov-Dec;14(6):1635-42
references.pmid: 18989020
references.citation: Magn Reson Imaging. 2006 Nov;24(9):1263-70
references.pmid: 17071347
references.citation: Magn Reson Med. 2007 Sep;58(3):497-510
references.pmid: 17763358
references.citation: Neuroimage. 2004 Nov;23(3):1176-85
references.pmid: 15528117
references.citation: Neuroimage. 2006 Jul 1;31(3):1086-103
references.pmid: 16546404
references.citation: IEEE Trans Image Process. 2010 Feb;19(2):461-77
references.pmid: 19887312
references.citation: Philos Trans R Soc Lond B Biol Sci. 2005 May 29;360(1457):893-902
references.pmid: 16087434
references.citation: Neuroimage. 2009 Aug;47 Suppl 2:T98-106
references.pmid: 18657622
references.citation: Ann N Y Acad Sci. 2005 Dec;1064:113-33
references.pmid: 16394152
references.citation: Magn Reson Med. 2005 Nov;54(5):1194-206
references.pmid: 16161109
references.citation: Magn Reson Med. 2002 Oct;48(4):577-82
references.pmid: 12353272
references.citation: Magn Reson Med. 2002 Aug;48(2):331-40
references.pmid: 12210942
references.citation: Magn Reson Med. 2004 Dec;52(6):1358-72
references.pmid: 15562495
references.citation: Magn Reson Med. 2003 Nov;50(5):1077-88
references.pmid: 14587019
references.citation: Neuroimage. 2007 Jan 1;34(1):144-55
references.pmid: 17070705
references.citation: Magn Reson Med. 2002 Jun;47(6):1083-99
references.pmid: 12111955
references.citation: IEEE Trans Med Imaging. 2007 Nov;26(11):1464-71
references.pmid: 18041262
references.citation: Magn Reson Med. 2005 Dec;54(6):1480-9
references.pmid: 16265642
references.citation: Neuroimage. 2005 Oct 1;27(4):725-36
references.pmid: 16111897
references.citation: NMR Biomed. 2002 Nov-Dec;15(7-8):456-67
references.pmid: 12489095
references.citation: Proc IEEE Int Symp Biomed Imaging. 2006;2006:65-68
references.pmid: 19759891
references.citation: Neuroimage. 2007 Aug 15;37(2):474-88
references.pmid: 17596967
references.citation: Med Image Comput Comput Assist Interv. 2007;10(Pt 2):144-52
references.pmid: 18044563
references.citation: J Magn Reson. 2006 Dec;183(2):193-202
references.pmid: 16963296
references.citation: Magn Reson Med. 2006 Jul;56(1):104-17
references.pmid: 16755539
references.citation: Inf Process Med Imaging. 2003 Jul;18:684-95
references.pmid: 15344498
references.citation: Med Image Anal. 2008 Oct;12(5):527-34
references.pmid: 18664412
references.citation: IEEE Trans Med Imaging. 2006 Aug;25(8):965-78
references.pmid: 16894991