A CUDA-based reverse gridding algorithm for MR reconstruction.

MR raw data collected using non-Cartesian method can be transformed on Cartesian grids by traditional gridding algorithm (GA) and reconstructed by Fourier transform. However, its runtime complexity is O(K×N(2)), where resolution of raw data is N×N and size of convolution window (CW) is K. And it involves a large number of matrix calculation including modulus, addition, multiplication and convolution. Therefore, a Compute Unified Device Architecture (CUDA)-based algorithm is proposed to improve the reconstruction efficiency of PROPELLER (a globally recognized non-Cartesian sampling method). Experiment shows a write-write conflict among multiple CUDA threads. This induces an inconsistent result when synchronously convoluting multiple k-space data onto the same grid. To overcome this problem, a reverse gridding algorithm (RGA) was developed. Different from the method of generating a grid window for each trajectory as in traditional GA, RGA calculates a trajectory window for each grid. This is what "reverse" means. For each k-space point in the CW, contribution is cumulated to this grid. Although this algorithm can be easily extended to reconstruct other non-Cartesian sampled raw data, we only implement it based on PROPELLER. Experiment illustrates that this CUDA-based RGA has successfully solved the write-write conflict and its reconstruction speed is 7.5 times higher than that of traditional GA.

[1]  Hossein Sedarat,et al.  On the optimality of the gridding reconstruction algorithm , 2000, IEEE Transactions on Medical Imaging.

[2]  Qing Huo Liu,et al.  An Efficient MR Image Reconstruction Method for Arbitrary K-space Trajectories Without Density Compensation , 2006, 2006 International Conference of the IEEE Engineering in Medicine and Biology Society.

[3]  C. Ahn,et al.  High-Speed Spiral-Scan Echo Planar NMR Imaging-I , 1986, IEEE Transactions on Medical Imaging.

[4]  Hairong Zheng,et al.  A Novel Gridding Algorithm Using NUFFT with Applications to Ultrasound Diffraction Tomography , 2008, 2008 2nd International Conference on Bioinformatics and Biomedical Engineering.

[5]  R. Turner,et al.  Echo-planar imaging: magnetic resonance imaging in a fraction of a second. , 1991, Science.

[6]  Jeffrey A. Fessler,et al.  Model-Based Image Reconstruction for MRI , 2010, IEEE Signal Processing Magazine.

[7]  Klaus Mueller,et al.  Methods for efficient, high quality volume resampling in the frequency domain , 2004, IEEE Visualization 2004.

[8]  S J Riederer,et al.  Current technical development of magnetic resonance imaging. , 2000, IEEE engineering in medicine and biology magazine : the quarterly magazine of the Engineering in Medicine & Biology Society.

[9]  H. Levine.  Medical Imaging , 2010, Annals of Biomedical Engineering.

[10]  Adam B Kerr,et al.  Spiral imaging artifact reduction: A comparison of two k‐trajectory measurement methods , 2009, Journal of magnetic resonance imaging : JMRI.

[11]  R. R. Ernst,et al.  NMR Fourier zeugmatography. , 2011, Journal of magnetic resonance.

[12]  Martin Cadík,et al.  FFT and Convolution Performance in Image Filtering on GPU , 2006, Tenth International Conference on Information Visualisation (IV'06).

[13]  David P. Luebke,et al.  CUDA: Scalable parallel programming for high-performance scientific computing , 2008, 2008 5th IEEE International Symposium on Biomedical Imaging: From Nano to Macro.

[14]  James G Pipe,et al.  Brain imaging in the unsedated pediatric patient: comparison of periodically rotated overlapping parallel lines with enhanced reconstruction and single-shot fast spin-echo sequences. , 2003, AJNR. American journal of neuroradiology.

[15]  A. A. Tamhane,et al.  Accelerating PROPELLER-MRI by Means of Under-Sampling and Iterative Image Reconstruction using the Non-Uniform Fast Fourier Transform , 2007 .

[16]  J G Pipe,et al.  PROPELLER MRI: Clinical testing of a novel technique for quantification and compensation of head motion , 2001, Journal of magnetic resonance imaging : JMRI.

[17]  Hongyu Guo,et al.  GPU Acceleration of PROPELLER MRI Using CUDA , 2009, 2009 3rd International Conference on Bioinformatics and Biomedical Engineering.

[18]  J. Pipe.  Motion correction with PROPELLER MRI: Application to head motion and free‐breathing cardiac imaging , 1999, Magnetic resonance in medicine.

[19]  Christopher Kumar Anand,et al.  Random Volumetric MRI Trajectories via Genetic Algorithms , 2008, Int. J. Biomed. Imaging.

[20]  P. Mansfield.  Multi-planar image formation using NMR spin echoes , 1977 .

[21]  Stuart Crozier,et al.  A new method in accelerating PROPELLER MRI , 2008, 2008 30th Annual International Conference of the IEEE Engineering in Medicine and Biology Society.

[22]  Dirk van Ormondt,et al.  Dynamic MR-Imaging with Radial Scanning, a Post-Acquisition Keyhole Approach , 2003, EURASIP J. Adv. Signal Process..

[23]  Justin P. Haldar,et al.  Accelerating advanced MRI reconstructions on GPUs , 2008, J. Parallel Distributed Comput..

[24]  James R. MacFall,et al.  Correction of MR k-space data corrupted by spike noise , 2000, IEEE Transactions on Medical Imaging.

[25]  H. D. Morris,et al.  K‐space in the clinic , 2004, Journal of magnetic resonance imaging : JMRI.

[26]  Steven C.R. Williams,et al.  Proc. Intl. Soc. Mag. Reson. Med. , 2012 .

[27]  Volker Rasche,et al.  Resampling of data between arbitrary grids using convolution interpolation , 1999, IEEE Transactions on Medical Imaging.

[28]  Paul C. Lauterbur.  Medical Imaging by Nuclear Magnetic Resonance Zeugmatography , 1979, IEEE Transactions on Nuclear Science.

[29]  Xiao Han,et al.  GPU-accelerated, gradient-free MI deformable registration for atlas-based MR brain image segmentation , 2009, 2009 IEEE Computer Society Conference on Computer Vision and Pattern Recognition Workshops.

[30]  Justin P. Haldar,et al.  Accelerating iterative field-compensated MR image reconstruction on GPUs , 2010, 2010 IEEE International Symposium on Biomedical Imaging: From Nano to Macro.

[31]  Qasim Zaidi,et al.  Three-dimensional shape from non-homogeneous textures: carved and stretched surfaces. , 2004, Journal of vision.

[32]  C. McKinstry,et al.  Nuclear magnetic resonance imaging in medicine , 1986, The Ulster medical journal.

[33]  Kyehyun Kim,et al.  Fast 2D-3D registration using GPU-based preprocessing , 2005, Proceedings of 7th International Workshop on Enterprise networking and Computing in Healthcare Industry, 2005. HEALTHCOM 2005..

[34]  Qing Huo Liu,et al.  An improved MRI reconstruction method based on table-lookup gridding , 2009, 2009 International Conference on Image Analysis and Signal Processing.

[35]  John D. Owens,et al.  GPU Computing , 2008, Proceedings of the IEEE.

[36]  John L. Prince,et al.  Discussion and Future Work , 1994 .

[37]  P. Jezzard,et al.  Sources of distortion in functional MRI data , 1999, Human brain mapping.

[38]  P. C. Lauterbur,et al.  NMR imaging in medicine: history, principles and prospects , 1983 .

[39]  Jan Timmer,et al.  The gridding method for image reconstruction by Fourier transformation , 1995, IEEE Trans. Medical Imaging.

[40]  Jeffrey A. Fessler,et al.  3D Forward and Back-Projection for X-Ray CT Using Separable Footprints , 2010, IEEE Transactions on Medical Imaging.

[41]  Hongyu Guo,et al.  Fast iterative reconstruction method for PROPELLER MRI , 2009, International Symposium on Multispectral Image Processing and Pattern Recognition.

[42]  Bob S. Hu,et al.  Fast Spiral Coronary Artery Imaging , 1992, Magnetic resonance in medicine.

[43]  A. Filler.  The History, Development and Impact of Computed Imaging in Neurological Diagnosis and Neurosurgery: CT, MRI, and DTI , 2009 .

[44]  R Mezrich,et al.  A perspective on K-space. , 1995, Radiology.

[45]  Tarek S. Abdelrahman,et al.  hiCUDA: a high-level directive-based language for GPU programming , 2009, GPGPU-2.

[46]  Jens Frahm,et al.  Spiral imaging: A critical appraisal , 2005, Journal of magnetic resonance imaging : JMRI.

[47]  Martin K. Purvis,et al.  GPU as a General Purpose Computing Resource , 2008, 2008 Ninth International Conference on Parallel and Distributed Computing, Applications and Technologies.

[48]  Mark J. Harris.  Many-core GPU computing with NVIDIA CUDA , 2008, ICS '08.

[49]  James A. Patton,et al.  Nuclear Magnetic Resonance (NMR) Imaging: Potential Impact On Medical Diagnosis And Digital Data Communications , 1982, Other Conferences.

[50]  Tobias Schaeffter,et al.  Accelerating the Nonequispaced Fast Fourier Transform on Commodity Graphics Hardware , 2008, IEEE Transactions on Medical Imaging.

[51]  A. Macovski,et al.  Selection of a convolution function for Fourier inversion using gridding [computerised tomography application]. , 1991, IEEE transactions on medical imaging.

[52]  P. Lauterbur,et al.  Image Formation by Induced Local Interactions: Examples Employing Nuclear Magnetic Resonance , 1973, Nature.

[53]  J. D. O'Sullivan,et al.  A Fast Sinc Function Gridding Algorithm for Fourier Inversion in Computer Tomography , 1985, IEEE Transactions on Medical Imaging.

[54]  Justin P. Haldar,et al.  Accelerating MR image reconstruction on GPUs , 2009, 2009 IEEE International Symposium on Biomedical Imaging: From Nano to Macro.

[55]  Rik Van de Walle,et al.  Reconstruction of MR images from data acquired on a general nonregular grid by pseudoinverse calculation , 2000, IEEE Transactions on Medical Imaging.