Accelerating Non-Cartesian MRI Reconstruction Convergence Using k-Space Preconditioning

We propose a k-space preconditioning formulation for accelerating the convergence of iterative Magnetic Resonance Imaging (MRI) reconstructions from non-uniformly sampled k-space data. Existing methods either use sampling density compensations which sacrifice reconstruction accuracy, or circulant preconditioners which increase per-iteration computation. Our approach overcomes both shortcomings. Concretely, we show that viewing the reconstruction problem in the dual formulation allows us to precondition in k-space using density-compensation-like operations. Using the primal-dual hybrid gradient method, the proposed preconditioning method does not have inner loops and are competitive in accelerating convergence compared to existing algorithms. We derive $\ell 2$ -optimized preconditioners, and demonstrate through experiments that the proposed method converges in about ten iterations in practice.

[1]  G H Glover,et al.  Projection Reconstruction Techniques for Reduction of Motion Effects in MRI , 1992, Magnetic resonance in medicine.

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

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

[4]  Stephen P. Boyd,et al.  Distributed Optimization and Statistical Learning via the Alternating Direction Method of Multipliers , 2011, Found. Trends Mach. Learn..

[5]  D. Nishimura,et al.  Reduced aliasing artifacts using variable‐density k‐space sampling trajectories , 2000, Magnetic resonance in medicine.

[6]  G. Zaharchuk,et al.  Recommended implementation of arterial spin-labeled perfusion MRI for clinical applications: A consensus of the ISMRM perfusion study group and the European consortium for ASL in dementia. , 2015, Magnetic resonance in medicine.

[7]  Alexandru Onose,et al.  An accelerated splitting algorithm for radio-interferometric imaging: when natural and uniform weighting meet , 2017, 1701.01748.

[8]  Jeffrey A. Fessler,et al.  Fast, iterative image reconstruction for MRI in the presence of field inhomogeneities , 2003, IEEE Transactions on Medical Imaging.

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

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

[11]  Antonin Chambolle,et al.  Diagonal preconditioning for first order primal-dual algorithms in convex optimization , 2011, 2011 International Conference on Computer Vision.

[12]  Dwight G Nishimura,et al.  Design and analysis of a practical 3D cones trajectory , 2006, Magnetic resonance in medicine.

[13]  M. Lustig,et al.  SPIRiT: Iterative self‐consistent parallel imaging reconstruction from arbitrary k‐space , 2010, Magnetic resonance in medicine.

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

[15]  Kevin M. Johnson,et al.  Optimized 3D ultrashort echo time pulmonary MRI , 2013, Magnetic resonance in medicine.

[16]  R. Hoge,et al.  Density compensation functions for spiral MRI , 1997, Magnetic resonance in medicine.

[17]  Michael Elad,et al.  ESPIRiT—an eigenvalue approach to autocalibrating parallel MRI: Where SENSE meets GRAPPA , 2014, Magnetic resonance in medicine.

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

[19]  Q. Liu,et al.  An accurate algorithm for nonuniform fast Fourier transforms (NUFFT's) , 1998 .

[20]  Gaël Varoquaux,et al.  The NumPy Array: A Structure for Efficient Numerical Computation , 2011, Computing in Science & Engineering.

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

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

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

[24]  I. M. Otivation Playing with Duality: An Overview of Recent Primal-Dual Approaches for Solving Large-Scale Optimization Problems , 2018 .

[25]  Peter Börnert,et al.  Accelerating compressed sensing in parallel imaging reconstructions using an efficient circulant preconditioner for cartesian trajectories , 2019, Magnetic resonance in medicine.

[26]  Andrew E. Yagle New fast preconditioners for Toeplitz-like linear systems , 2002, 2002 IEEE International Conference on Acoustics, Speech, and Signal Processing.

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

[28]  Jeffrey A. Fessler,et al.  Augmented Lagrangian with Variable Splitting for Faster Non-Cartesian ${\rm L}_{1}$-SPIRiT MR Image Reconstruction , 2014, IEEE Transactions on Medical Imaging.

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

[30]  Antonin Chambolle,et al.  Stochastic Primal-Dual Hybrid Gradient Algorithm with Arbitrary Sampling and Imaging Applications , 2017, SIAM J. Optim..

[31]  Antonin Chambolle,et al.  A First-Order Primal-Dual Algorithm for Convex Problems with Applications to Imaging , 2011, Journal of Mathematical Imaging and Vision.

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

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

[34]  Yunhong Shu,et al.  3616 A Preconditioned ADMM Strategy for Field-Corrected Non-Cartesian MRI Reconstruction , 2013 .

[35]  D. Donoho,et al.  Sparse MRI: The application of compressed sensing for rapid MR imaging , 2007, Magnetic resonance in medicine.

[36]  Frank Ong,et al.  Motion robust high resolution 3D free‐breathing pulmonary MRI using dynamic 3D image self‐navigator , 2018, Magnetic resonance in medicine.

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

[38]  S. Hido,et al.  CuPy : A NumPy-Compatible Library for NVIDIA GPU Calculations , 2017 .