Accelerating nonlinear diffusion tensor estimation for medical image processing using high performance GPU clusters

Diffusion Tensor Imaging (DTI) is a non-invasive magnetic resonance technique that produces in vivo images of biological tissues with local microstructural characteristics such as water diffusion. It can be used, for example, to localize white matter lesions, or in neuro-navigation surgery of brain tumors. Diffusion tensor maps are usually computed on a voxel-by-voxel basis by fitting the signal intensities of diffusion weighted images as a function of their corresponding data acquisition parameters. This processing is highly computation-intensive and can be time-consuming which constraints the clinical use of DTI. This study presents the application of using high performance GPU clusters in diffusion tensor estimation by accelerating the multivariate non-linear regression. The results are tested in simulated DTI brain datasets and show significant performance gain in tensor fitting in addition to favorable scalability characteristics. The proposed GPU implementation framework can further promote the clinical use of DTI, and can be used to accelerate statistical analysis of DTI where Monte Carlo simulations are employed, or readily applied to quantitative assessment of DTI using bootstrap analysis.

[1]  D. Marquardt An Algorithm for Least-Squares Estimation of Nonlinear Parameters , 1963 .

[2]  Kenneth Levenberg A METHOD FOR THE SOLUTION OF CERTAIN NON – LINEAR PROBLEMS IN LEAST SQUARES , 1944 .

[3]  Inanc Senocak,et al.  An MPI-CUDA Implementation for Massively Parallel Incompressible Flow Computations on Multi-GPU Clusters , 2010 .

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

[5]  Erkki Oja,et al.  Solving Large Regression Problems using an Ensemble of GPU-accelerated ELMs , 2010, ESANN.

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

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

[8]  Carlo Pierpaoli,et al.  Variance of estimated DTI‐derived parameters via first‐order perturbation methods , 2007, Magnetic resonance in medicine.

[9]  Sinisa Pajevic,et al.  Color schemes to represent the orientation of anisotropic tissues from diffusion tensor data: Application to white matter fiber tract mapping in the human brain , 1999, Magnetic resonance in medicine.

[10]  P. Basser,et al.  Microstructural and physiological features of tissues elucidated by quantitative-diffusion-tensor MRI. , 1996, Journal of magnetic resonance. Series B.

[11]  Christopher Nimsky,et al.  Hybrid Visualization for White Matter Tracts using Triangle Strips and Point Sprites , 2006, IEEE Transactions on Visualization and Computer Graphics.

[12]  P. Basser,et al.  The b matrix in diffusion tensor echo‐planar imaging , 1997, Magnetic resonance in medicine.

[13]  Cheng Guan Koay,et al.  A simple scheme for generating nearly uniform distribution of antipodally symmetric points on the unit sphere , 2011, J. Comput. Sci..

[14]  I. Peña,et al.  Matrix computations on graphics processors and clusters of gpus , 2011 .

[15]  P. Basser,et al.  Diffusion tensor MR imaging of the human brain. , 1996, Radiology.

[16]  P. Basser,et al.  Estimation of the effective self-diffusion tensor from the NMR spin echo. , 1994, Journal of magnetic resonance. Series B.