Real-Time Reconstruction of Sensitivity Encoded Radial Magnetic Resonance Imaging Using a Graphics Processing Unit

A barrier to the adoption of non-Cartesian parallel magnetic resonance imaging for real-time applications has been the times required for the image reconstructions. These times have exceeded the underlying acquisition time thus preventing real-time display of the acquired images. We present a reconstruction algorithm for commodity graphics hardware (GPUs) to enable real time reconstruction of sensitivity encoded radial imaging (radial SENSE). We demonstrate that a radial profile order based on the golden ratio facilitates reconstruction from an arbitrary number of profiles. This allows the temporal resolution to be adjusted on the fly. A user adaptable regularization term is also included and, particularly for highly undersampled data, used to interactively improve the reconstruction quality. Each reconstruction is fully self-contained from the profile stream, i.e., the required coil sensitivity profiles, sampling density compensation weights, regularization terms, and noise estimates are computed in real-time from the acquisition data itself. The reconstruction implementation is verified using a steady state free precession (SSFP) pulse sequence and quantitatively evaluated. Three applications are demonstrated; real-time imaging with real-time SENSE 1) or k-t SENSE 2) reconstructions, and 3) offline reconstruction with interactive adjustment of reconstruction settings.

[1]  K. Kwong,et al.  Parallel imaging reconstruction using automatic regularization , 2004, Magnetic resonance in medicine.

[2]  Christof Baltes,et al.  k‐t BLAST reconstruction from non‐Cartesian k‐t space sampling , 2006, Magnetic resonance in medicine.

[3]  Dwight G. Nishimura,et al.  Rapid gridding reconstruction with a minimal oversampling ratio , 2005, IEEE Transactions on Medical Imaging.

[4]  David Atkinson,et al.  Cartesian SENSE and k‐t SENSE reconstruction using commodity graphics hardware , 2008, Magnetic resonance in medicine.

[5]  Sébastien Roujol,et al.  Online real‐time reconstruction of adaptive TSENSE with commodity CPU/GPU hardware , 2009, Magnetic resonance in medicine.

[6]  F H Epstein,et al.  Adaptive sensitivity encoding incorporating temporal filtering (TSENSE) † , 2001, Magnetic resonance in medicine.

[7]  Peter Boesiger,et al.  k‐t BLAST and k‐t SENSE: Dynamic MRI with high frame rate exploiting spatiotemporal correlations , 2003, Magnetic resonance in medicine.

[8]  Olaf Dössel,et al.  An Optimal Radial Profile Order Based on the Golden Ratio for Time-Resolved MRI , 2007, IEEE Transactions on Medical Imaging.

[9]  Daniel K Sodickson,et al.  Self‐calibrating parallel imaging with automatic coil sensitivity extraction , 2002, Magnetic resonance in medicine.

[10]  Satoshi Matsuoka,et al.  Bandwidth intensive 3-D FFT kernel for GPUs using CUDA , 2008, 2008 SC - International Conference for High Performance Computing, Networking, Storage and Analysis.

[11]  Stefan Kunis,et al.  Using NFFT 3---A Software Library for Various Nonequispaced Fast Fourier Transforms , 2009, TOMS.

[12]  Jeffrey L Duerk,et al.  Image‐guided and ‐monitored renal artery stenting using only MRI , 2006, Journal of magnetic resonance imaging : JMRI.

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

[14]  P. Hansen Rank-Deficient and Discrete Ill-Posed Problems: Numerical Aspects of Linear Inversion , 1987 .

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

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

[17]  J. Pipe,et al.  Sampling density compensation in MRI: Rationale and an iterative numerical solution , 1999, Magnetic resonance in medicine.

[18]  P. Boesiger,et al.  SENSE: Sensitivity encoding for fast MRI , 1999, Magnetic resonance in medicine.

[19]  P. Boesiger,et al.  Advances in sensitivity encoding with arbitrary k‐space trajectories , 2001, Magnetic resonance in medicine.

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

[21]  Warren J Manning,et al.  Magnetic Resonance–Guided Coronary Artery Stent Placement in a Swine Model , 2002, Circulation.

[22]  H Eggers,et al.  Projection reconstruction balanced fast field echo for interactive real‐time cardiac imaging , 2001, Magnetic resonance in medicine.