Probabilistic Brain Fiber Tractography on GPUs

Diffusion Tensor Magnetic Resonance Imaging (DT-MRI) is an emerging technique that explores the structural connectivity of the human brain. The probabilistic fiber tractography based on DT-MRI data behaves more robustly than deterministic approaches in the presence of fiber crossings, but requires more prohibitive computational time. In this work we present a GPU-based probabilistic framework for brain fiber tractography. The framework includes two main steps: 1) Markov-Chain Monte-Carlo (MCMC) sampling, and 2) probabilistic streamlining fiber tracking. We implement the Metropolis-Hastings sampling for local parameter estimation on GPU. In the probabilistic streamlining fiber tracking, we find that fiber lengths are exponentially distributed, and propose a novel segmenting strategy to improve the load balance. On mid-range GPUs, we achieve performance gains up to 34x and 50x over CPUs for the two steps respectively.

[1]  Junji Yamato,et al.  Real-time estimation of human visual attention with dynamic Bayesian network and MCMC-based particle filter , 2009, 2009 IEEE International Conference on Multimedia and Expo.

[2]  Heinz-Otto Peitgen,et al.  A GPU-based fiber tracking framework using geometry shaders , 2009, Medical Imaging.

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

[4]  Teresa H. Y. Meng,et al.  High-throughput Bayesian network learning using heterogeneous multicore computers , 2010, ICS '10.

[5]  Mark W. Woolrich,et al.  Advances in functional and structural MR image analysis and implementation as FSL , 2004, NeuroImage.

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

[7]  Eros Comunello,et al.  Enabling Interactive Brain Fiber Tracking with the GPU , 2008 .

[8]  Rüdiger Westermann,et al.  The application of GPU particle tracing to diffusion tensor field visualization , 2005, VIS 05. IEEE Visualization, 2005..

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

[10]  Mark W. Woolrich,et al.  Bayesian analysis of neuroimaging data in FSL , 2009, NeuroImage.

[11]  Ola Friman,et al.  Qualitative and quantitative analysis of probabilistic and deterministic fiber tracking , 2010, Medical Imaging.

[12]  William J. Dally,et al.  The GPU Computing Era , 2010, IEEE Micro.

[13]  Anna Vilanova,et al.  Visualization of DTI fibers using hair-rendering techniques , 2006 .

[14]  Derek K. Jones,et al.  Isotropic resolution diffusion tensor imaging with whole brain acquisition in a clinically acceptable time , 2002, Human brain mapping.

[15]  Mark W. Woolrich,et al.  Probabilistic diffusion tractography with multiple fibre orientations: What can we gain? , 2007, NeuroImage.

[16]  Geoffrey J M Parker,et al.  A framework for a streamline‐based probabilistic index of connectivity (PICo) using a structural interpretation of MRI diffusion measurements , 2003, Journal of magnetic resonance imaging : JMRI.

[17]  Heidi Johansen-Berg,et al.  Using diffusion imaging to study human connectional anatomy. , 2009, Annual review of neuroscience.

[18]  Pierre L'Ecuyer,et al.  Maximally equidistributed combined Tausworthe generators , 1996, Math. Comput..

[19]  Thomas Ertl,et al.  GPU-Based Hyperstreamlines for Diffusion Tensor Imaging , 2006, EuroVis.

[20]  Hubert Nguyen,et al.  GPU Gems 3 , 2007 .

[21]  Ross T. Whitaker,et al.  Interactive Visualization of Volumetric White Matter Connectivity in DT-MRI Using a Parallel-Hardware Hamilton-Jacobi Solver , 2007, IEEE Transactions on Visualization and Computer Graphics.

[22]  Gerik Scheuermann,et al.  Fast and Memory Efficient GPU-Based Rendering of Tensor Data , 2011 .

[23]  S. Chib,et al.  Understanding the Metropolis-Hastings Algorithm , 1995 .

[24]  Eros Comunello,et al.  Diffusion tensor fiber tracking on graphics processing units , 2008, Comput. Medical Imaging Graph..

[25]  Vid Petrovic,et al.  Visualizing Whole-Brain DTI Tractography with GPU-based Tuboids and LoD Management , 2007, IEEE Transactions on Visualization and Computer Graphics.

[26]  Carl-Fredrik Westin,et al.  Uncertainty in White Matter Fiber Tractography , 2005, MICCAI.

[27]  Anna Vilanova,et al.  Infants with perinatal hypoxic ischemia: feasibility of fiber tracking at birth and 3 months. , 2006, Radiology.

[28]  Hua Zhou,et al.  Graphics Processing Units and High-Dimensional Optimization. , 2010, Statistical science : a review journal of the Institute of Mathematical Statistics.

[29]  Timothy Edward John Behrens,et al.  Characterization and propagation of uncertainty in diffusion‐weighted MR imaging , 2003, Magnetic resonance in medicine.

[30]  Cliburn Chan,et al.  Understanding GPU Programming for Statistical Computation: Studies in Massively Parallel Massive Mixtures , 2010, Journal of computational and graphical statistics : a joint publication of American Statistical Association, Institute of Mathematical Statistics, Interface Foundation of North America.

[31]  M. E. Muller,et al.  A Note on the Generation of Random Normal Deviates , 1958 .

[32]  E. V. Aart,et al.  Acceleration of a geodesic fiber-tracking algorithm for diffusion tensor imaging using CUDA , 2010 .

[33]  M. Raichle,et al.  Tracking neuronal fiber pathways in the living human brain. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[34]  Mario A. R. Dantas,et al.  Design and Implementation of Brain Fiber Tracking for GPUs and PC Clusters , 2009, 2009 21st International Symposium on Computer Architecture and High Performance Computing.

[35]  Jarad Niemi,et al.  Efficient Bayesian inference in stochastic chemical kinetic models using graphical processing units , 2011, 1101.4242.

[36]  Monte Carlo Integration.  Markov Chain Monte Carlo and Gibbs Sampling , 2002 .

[37]  Jean-Philippe Thiran,et al.  DTI mapping of human brain connectivity: statistical fibre tracking and virtual dissection , 2003, NeuroImage.

[38]  Carl-Fredrik Westin,et al.  A Bayesian approach for stochastic white matter tractography , 2006, IEEE Transactions on Medical Imaging.