Divide et impera: Acceleration of DTI tractography using multi‐GPU parallel processing

Diffusion tensor imaging (DTI) tractography technique represents neural fiber pathways by using local tensor information based on water diffusion anisotropy in brain white matter. However, DTI tractography is often unable to reconstruct crossing, kissing, and branching fiber trajectories due to intrinsic limitations of DTI. Increasingly complex tractography algorithms provided reliable and visually pleasing results, yet at an increasing computational cost in comparison with simple tractography algorithms. To shorten the computation time, we developed multi‐GPU (graphics processing unit)‐based parallelized versions of deterministic and probabilistic tractography algorithms to investigate their utility for near‐real time tractography. We were able to dramatically reduce the computation time using multiple GPUs (three NVIDIA TESLA C1060s) in comparison to the central processing unit (CPU) sequential processing. Deterministic tractography could accelerate 101 times faster, and probabilistic tractography could accelerate 63 times faster. The results showed that parallel tractography algorithm is well suited with GPU which has fundamentally parallelized architecture. © 2013 Wiley Periodicals, Inc. Int J Imaging Syst Technol, 23, 256–264, 2013

[1]  Jie Cheng,et al.  Programming Massively Parallel Processors. A Hands-on Approach , 2010, Scalable Comput. Pract. Exp..

[2]  Timothy Edward John Behrens,et al.  Non-invasive mapping of connections between human thalamus and cortex using diffusion imaging , 2003, Nature Neuroscience.

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

[4]  David G. Norris,et al.  An Investigation of Functional and Anatomical Connectivity Using Magnetic Resonance Imaging , 2002, NeuroImage.

[5]  Alan Connelly,et al.  Diffusion-weighted magnetic resonance imaging fibre tracking using a front evolution algorithm , 2003, NeuroImage.

[6]  Zeynep Firat,et al.  High b-value diffusion-weighted MR imaging of normal brain at 3T. , 2009, European journal of radiology.

[7]  Yonggang Lu,et al.  Improved fiber tractography with Bayesian tensor regularization , 2006, NeuroImage.

[8]  Rachid Deriche,et al.  Deterministic and Probabilistic Tractography Based on Complex Fibre Orientation Distributions , 2009, IEEE Transactions on Medical Imaging.

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

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

[11]  P. Kingsley,et al.  Introduction to diffusion tensor imaging mathematics: Part I. Tensors, rotations, and eigenvectors , 2006 .

[12]  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.

[13]  C. Beaulieu,et al.  The basis of anisotropic water diffusion in the nervous system – a technical review , 2002, NMR in biomedicine.

[14]  Jie Cheng,et al.  CUDA by Example: An Introduction to General-Purpose GPU Programming , 2010, Scalable Comput. Pract. Exp..

[15]  Kunihiko Fukuda,et al.  Usefulness of high-b-value diffusion-weighted imaging in acute cerebral infarction , 2007, European Radiology.

[16]  Daniel C. Alexander,et al.  Probabilistic Monte Carlo Based Mapping of Cerebral Connections Utilising Whole-Brain Crossing Fibre Information , 2003, IPMI.

[17]  Michael J. Flynn,et al.  Some Computer Organizations and Their Effectiveness , 1972, IEEE Transactions on Computers.