Trajectory optimized NUFFT: Faster non‐Cartesian MRI reconstruction through prior knowledge and parallel architectures
Authors:
Open Access
- 17 October 2018
- journal article
- research article
- Published by Wiley in Magnetic Resonance in Medicine
- Vol. 81 (3) , 2064-2071
Abstract
Purpose: The non‐uniform fast Fourier transform (NUFFT) involves interpolation of non‐uniformly sampled Fourier data onto a Cartesian grid, an interpolation that is slowed by complex, non‐local data access patterns. A faster NUFFT would increase the clinical relevance of the plethora of advanced non‐Cartesian acquisition methods. Methods: Here we customize the NUFFT procedure for a radial trajectory and GPU architecture to eliminate the bottlenecks encountered when allowing for arbitrary trajectories and hardware. We call the result TRON, for TR ajectory O ptimized N UFFT. We benchmark the speed and accuracy TRON on a Shepp‐Logan phantom and on whole‐body continuous golden‐angle radial MRI. Results: TRON was 6–30× faster than the closest competitor, depending on test data set, and was the most accurate code tested. Conclusions: Specialization of the NUFFT algorithm for a particular trajectory yielded significant speed gains. TRON can be easily extended to other trajectories, such as spiral and PROPELLER. TRON can be downloaded at http://github.com/davidssmith/TRON .Keywords
Funding Information
- National Cancer Institute (K25 CA176219)
- National Institute of Diabetes and Digestive and Kidney Diseases (R21 DK096282)
- National Center for Advancing Translational Sciences (UL1 TR000445)
This publication has 16 references indexed in Scilit:
- XD‐GRASP: Golden‐angle radial MRI with reconstruction of extra motion‐state dimensions using compressed sensingMagnetic Resonance in Medicine, 2015
- Continuously moving table MRI with golden angle radial samplingMagnetic Resonance in Medicine, 2014
- More IMPATIENT: A gridding-accelerated Toeplitz-based strategy for non-Cartesian high-resolution 3D MRI on GPUsJournal of Parallel and Distributed Computing, 2013
- Self-gated Radial MRI for Respiratory Motion Compensation on Hybrid PET/MR Systems2013
- The nonequispaced FFT on graphics processing unitsProceedings in Applied Mathematics and Mechanics, 2012
- A CUDA-based reverse gridding algorithm for MR reconstructionMagnetic Resonance Imaging, 2012
- Accelerating the Nonequispaced Fast Fourier Transform on Commodity Graphics HardwareIEEE Transactions on Medical Imaging, 2008
- Image quality assessment: from error visibility to structural similarityIEEE Transactions on Image Processing, 2004
- Nonuniform fast fourier transforms using min-max interpolationIEEE Transactions on Signal Processing, 2003
- Selection of a convolution function for Fourier inversion using gridding (computerised tomography application)IEEE Transactions on Medical Imaging, 1991