Accelerating Global Tractography Using Parallel Markov Chain Monte Carlo.

Global tractography estimates brain connectivity by determining the optimal configuration of signal-generating fiber segments that best describes the measured diffusion-weighted data, promising better stability than local greedy methods with respect to imaging noise. However, global tractography is computationally very demanding and requires computation times that are often prohibitive for clinical applications. We present here a reformulation of the global tractography algorithm for fast parallel implementation amendable to acceleration using multi-core CPUs and general-purpose GPUs. Our method is motivated by the key observation that each fiber segment is affected by a limited spatial neighborhood. That is, a fiber segment is influenced only by the fiber segments that are (or can potentially be) connected to its both ends and also by the diffusion-weighted signal in its proximity. This observation makes it possible to parallelize the Markov chain Monte Carlo (MCMC) algorithm used in the global tractography algorithm so that updating of independent fiber segments can be done concurrently. The experiments show that the proposed algorithm can significantly speed up global tractography, while at the same time maintain or improve tractography performance.

[1]  W. K. Hastings,et al.  Monte Carlo Sampling Methods Using Markov Chains and Their Applications , 1970 .

[2]  Daoqiang Zhang,et al.  Identification of MCI individuals using structural and functional connectivity networks , 2012, NeuroImage.

[3]  P. V. van Zijl,et al.  Three‐dimensional tracking of axonal projections in the brain by magnetic resonance imaging , 1999, Annals of neurology.

[4]  Paul M. Thompson,et al.  Automatic clustering of white matter fibers in brain diffusion MRI with an application to genetics , 2014, NeuroImage.

[5]  Susanne Schnell,et al.  Global fiber reconstruction becomes practical , 2011, NeuroImage.

[6]  Dinggang Shen,et al.  Development Trends of White Matter Connectivity in the First Years of Life , 2011, PloS one.

[7]  Maxime Descoteaux,et al.  Quantitative evaluation of 10 tractography algorithms on a realistic diffusion MR phantom , 2011, NeuroImage.

[8]  Peter F. Neher,et al.  MITK global tractography , 2012, Medical Imaging.

[9]  Dinggang Shen,et al.  Human Brain Connectomics: Networks, Techniques, and Applications [Life Sciences] , 2010, IEEE Signal Processing Magazine.

[10]  Dinggang Shen,et al.  Longitudinal Tractography with Application to Neuronal Fiber Trajectory Reconstruction in Neonates , 2011, MICCAI.

[11]  Emile H. L. Aarts,et al.  Simulated annealing and Boltzmann machines - a stochastic approach to combinatorial optimization and neural computing , 1990, Wiley-Interscience series in discrete mathematics and optimization.

[12]  P. Basser,et al.  In vivo fiber tractography using DT‐MRI data , 2000, Magnetic resonance in medicine.

[13]  Wenbin Li,et al.  Enriched white matter connectivity networks for accurate identification of MCI patients , 2011, NeuroImage.

[14]  Randy L. Gollub,et al.  Reproducibility of quantitative tractography methods applied to cerebral white matter , 2007, NeuroImage.

[15]  Dinggang Shen,et al.  Altered Structural Connectivity in Neonates at Genetic Risk for Schizophrenia: a Combined Study Using Morphological and White Matter Networks , 2022 .

[16]  Guido Gerig,et al.  Probabilistic white matter fiber tracking using particle filtering and von Mises-Fisher sampling , 2009, Medical Image Anal..

[17]  Paul M. Thompson,et al.  Labeling white matter tracts in hardi by fusing multiple tract atlases with applications to genetics , 2013, 2013 IEEE 10th International Symposium on Biomedical Imaging.

[18]  Dinggang Shen,et al.  PopTract: Population-Based Tractography , 2011, IEEE Transactions on Medical Imaging.

[19]  Dinggang Shen,et al.  Spatial Transformation of DWI Data Using Non-Negative Sparse Representation , 2012, IEEE Transactions on Medical Imaging.

[20]  Tim McGraw,et al.  Stochastic DT-MRI Connectivity Mapping on the GPU , 2007, IEEE Transactions on Visualization and Computer Graphics.

[21]  V. Kiselev,et al.  Gibbs tracking: A novel approach for the reconstruction of neuronal pathways , 2008, Magnetic resonance in medicine.

[22]  Jungsoo Lee,et al.  Divide et impera: Acceleration of DTI tractography using multi‐GPU parallel processing , 2013, Int. J. Imaging Syst. Technol..

[23]  Chong Wang,et al.  Asymptotically Exact, Embarrassingly Parallel MCMC , 2013, UAI.