Monotone FISTA With Variable Acceleration for Compressed Sensing Magnetic Resonance Imaging

An improvement of the monotone fast iterative shrinkage-thresholding algorithm (MFISTA) for faster convergence is proposed in this paper. Our motivation is to reduce the reconstruction time of compressed sensing problems in magnetic resonance imaging. The proposed modification introduces an extra term, which is a multiple of the proximal-gradient step, into the so-called momentum formula used for the computation of the next iterate in MFISTA. In addition, the modified algorithm selects the next iterate as a possibly improved point obtained by any other procedure, such as an arbitrary shift, a line search, or other methods. As an example, an arbitrary-length shift in the direction from the previous iterate to the output of the proximal-gradient step is considered. The resulting algorithm accelerates MFISTA in a manner that varies with the iterative steps. Convergence analysis shows that the proposed modification provides improved theoretical convergence bounds, and that it has more flexibility in its parameters than the original MFISTA. Since such problems need to be studied in the context of functions of several complex variables, a careful extension of FISTA-like methods to complex variables is provided.

[1]  Sungheon Kim,et al.  Golden‐angle radial sparse parallel MRI: Combination of compressed sensing, parallel imaging, and golden‐angle radial sampling for fast and flexible dynamic volumetric MRI , 2014, Magnetic resonance in medicine.

[2]  Jeffrey A. Fessler,et al.  An optimized first-order method for image restoration , 2015, 2015 IEEE International Conference on Image Processing (ICIP).

[3]  J. Pauly,et al.  Accelerating parameter mapping with a locally low rank constraint , 2015, Magnetic resonance in medicine.

[4]  M. Lustig,et al.  Compressed Sensing MRI , 2008, IEEE Signal Processing Magazine.

[5]  Robert D. Nowak,et al.  An EM algorithm for wavelet-based image restoration , 2003, IEEE Trans. Image Process..

[6]  Daniel K Sodickson,et al.  Low‐rank plus sparse matrix decomposition for accelerated dynamic MRI with separation of background and dynamic components , 2015, Magnetic resonance in medicine.

[7]  I. Daubechies,et al.  An iterative thresholding algorithm for linear inverse problems with a sparsity constraint , 2003, math/0307152.

[8]  E.J. Candes,et al.  An Introduction To Compressive Sampling , 2008, IEEE Signal Processing Magazine.

[9]  Y. Nesterov A method for unconstrained convex minimization problem with the rate of convergence o(1/k^2) , 1983 .

[10]  H. Chandarana,et al.  Volumetric multicomponent T1ρ relaxation mapping of the human liver under free breathing at 3T , 2019, Magnetic resonance in medicine.

[11]  Donghwan Kim,et al.  Optimized first-order methods for smooth convex minimization , 2014, Mathematical Programming.

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

[13]  R. Otazo,et al.  Accelerated mono‐ and biexponential 3D‐T1ρ relaxation mapping of knee cartilage using golden angle radial acquisitions and compressed sensing , 2020, Magnetic resonance in medicine.

[14]  Emmanuel J. Candès,et al.  Adaptive Restart for Accelerated Gradient Schemes , 2012, Foundations of Computational Mathematics.

[15]  Paul C. Lauterbur,et al.  Principles of magnetic resonance imaging : a signal processing perspective , 1999 .

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

[17]  Jeffrey A Fessler,et al.  On NUFFT-based gridding for non-Cartesian MRI. , 2007, Journal of magnetic resonance.

[18]  Stephen P. Boyd,et al.  Proximal Algorithms , 2013, Found. Trends Optim..

[19]  Marc Teboulle,et al.  Fast Gradient-Based Algorithms for Constrained Total Variation Image Denoising and Deblurring Problems , 2009, IEEE Transactions on Image Processing.

[20]  S. Yun,et al.  An accelerated proximal gradient algorithm for nuclear norm regularized linear least squares problems , 2009 .

[21]  G. Chang,et al.  Biexponential T1ρ relaxation mapping of human knee cartilage in vivo at 3 T , 2017, NMR in biomedicine.

[22]  Adrien B. Taylor,et al.  Exact Worst-Case Performance of First-Order Methods for Composite Convex Optimization , 2015, SIAM J. Optim..

[23]  Jeffrey A. Fessler,et al.  ITERATIVE RECONSTRUCTION METHODS FOR NON-CARTESIAN MRI , 2007 .

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

[25]  Li Feng,et al.  Compressed sensing for body MRI , 2017, Journal of magnetic resonance imaging : JMRI.

[26]  Marc Teboulle,et al.  Performance of first-order methods for smooth convex minimization: a novel approach , 2012, Mathematical Programming.

[27]  Alvaro R. De Pierro,et al.  Fast and exact unidimensional L2-L1 optimization as an accelerator for iterative reconstruction algorithms , 2016, Digit. Signal Process..

[28]  R. Otazo,et al.  Accelerating 3D‐T1ρ mapping of cartilage using compressed sensing with different sparse and low rank models , 2018, Magnetic resonance in medicine.

[29]  Elias S. Helou,et al.  Accelerating Overrelaxed and Monotone Fast Iterative Shrinkage-Thresholding Algorithms With Line Search for Sparse Reconstructions , 2017, IEEE Transactions on Image Processing.

[30]  Michael Elad,et al.  L1-L2 Optimization in Signal and Image Processing , 2010, IEEE Signal Processing Magazine.

[31]  I. Yamada,et al.  Over-relaxation of the fast iterative shrinkage-thresholding algorithm with variable stepsize , 2011 .

[32]  Kyle J. Myers,et al.  Foundations of Image Science , 2003, J. Electronic Imaging.