Accelerated noncartesian sense reconstruction using a majorize-minimize algorithm combining variable-splitting

Magnetic resonance imaging (MRI) provides great flexibility in the choice of k-space sampling trajectories. NonCartesian trajectories exhibit several advantages over Cartesian ones but are less amenable to FFT-based manipulation of k-space data. Thus, existing iterative reconstruction methods for nonCartesian trajectories require relatively more computation (interpolation/gridding in addition to FFTs) and can be slow, especially for (undersampled) parallel MRI. In this work, we focus on SENSE-based regularized image reconstruction for nonCartesian trajectories and propose a majorize-minimize approach where we first majorize the SENSE data-fidelity term with a quadratic form involving a symmetric positive definite circulant matrix. For the minimization step, we apply a suitable variable splitting (VS) strategy combined with the augmented Lagrangian framework and alternating minimization that together decouple the circulant matrix from coil sensitivities and the regularizer. The resulting iterative algorithm admits simple update steps, is amenable to FFT-based matrix inversions due in part to the circulant matrix in the majorizer and provides a natural framework for incorporating a two-step procedure for acceleration. Simulations indicate that the proposed algorithm converges faster than some state-of-the-art VSbased iterative image reconstruction methods for the same problem.

[1]  Matthew I. Grivich,et al.  The magnetic field of current-carrying polygons: An application of vector field rotations , 2000 .

[2]  Marc Teboulle,et al.  A Fast Iterative Shrinkage-Thresholding Algorithm for Linear Inverse Problems , 2009, SIAM J. Imaging Sci..

[3]  A. Evans,et al.  MRI simulation-based evaluation of image-processing and classification methods , 1999, IEEE Transactions on Medical Imaging.

[4]  Xavier Bresson,et al.  Bregmanized Nonlocal Regularization for Deconvolution and Sparse Reconstruction , 2010, SIAM J. Imaging Sci..

[5]  Jeffrey A. Fessler,et al.  Parallel MR Image Reconstruction Using Augmented Lagrangian Methods , 2011, IEEE Transactions on Medical Imaging.

[6]  Jeffrey A. Fessler,et al.  Nonuniform fast Fourier transforms using min-max interpolation , 2003, IEEE Trans. Signal Process..

[7]  Jeffrey A. Fessler,et al.  Accelerated Regularized Estimation of MR Coil Sensitivities Using Augmented Lagrangian Methods , 2013, IEEE Transactions on Medical Imaging.

[8]  Alan C. Evans,et al.  MRI Simulation Based Evaluation and Classifications Methods , 1999, IEEE Trans. Medical Imaging.

[9]  Bo Liu,et al.  A statistical approach to SENSE regularization with arbitrary k‐space trajectories , 2008, Magnetic resonance in medicine.

[10]  T. Chan An Optimal Circulant Preconditioner for Toeplitz Systems , 1988 .

[11]  A. Bos Complex gradient and Hessian , 1994 .

[12]  Tom Goldstein,et al.  The Split Bregman Method for L1-Regularized Problems , 2009, SIAM J. Imaging Sci..

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

[14]  Martin Blaimer,et al.  Non‐Cartesian data reconstruction using GRAPPA operator gridding (GROG) , 2007, Magnetic resonance in medicine.

[15]  Klaas Paul Pruessmann,et al.  Realistic Analytical Phantoms for Parallel Magnetic Resonance Imaging , 2012, IEEE Transactions on Medical Imaging.

[16]  Raymond H. Chan,et al.  Conjugate Gradient Methods for Toeplitz Systems , 1996, SIAM Rev..