Weighted-average alternating minimization method for magnetic resonance image reconstruction based on compressive sensing

The problem of compressive-sensing (CS) L2-L1-TV reconstruction of magnetic resonance (MR) scans from undersampled $k$-space data has been addressed in numerous studies. However, the regularization parameters in models of CS L2-L1-TV reconstruction are rarely studied. Once the regularization parameters are given, the solution for an MR reconstruction model is fixed and is less effective in the case of strong noise. To overcome this shortcoming, we present a new alternating formulation to replace the standard L2-L1-TV reconstruction model. A weighted-average alternating minimization method is proposed based on this new formulation and a convergence analysis of the method is carried out. The advantages of and the motivation for the proposed alternating formulation are explained. Experimental results demonstrate that the proposed formulation yields better reconstruction results in the case of strong noise and can improve image reconstruction via flexible parameter selection.

[1]  Balas K. Natarajan,et al.  Sparse Approximate Solutions to Linear Systems , 1995, SIAM J. Comput..

[2]  Benar Fux Svaiter,et al.  General Projective Splitting Methods for Sums of Maximal Monotone Operators , 2009, SIAM J. Control. Optim..

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

[4]  Shiqian Ma,et al.  Fast Multiple-Splitting Algorithms for Convex Optimization , 2009, SIAM J. Optim..

[5]  Yonggui Zhu,et al.  A Fast Method for Reconstruction of Total-Variation MR Images With a Periodic Boundary Condition , 2013, IEEE Signal Processing Letters.

[6]  Patrick L. Combettes,et al.  Signal Recovery by Proximal Forward-Backward Splitting , 2005, Multiscale Model. Simul..

[7]  Z. Opial.  Weak convergence of the sequence of successive approximations for nonexpansive mappings , 1967 .

[8]  Paul Tseng,et al.  A coordinate gradient descent method for nonsmooth separable minimization , 2008, Math. Program..

[9]  Mário A. T. Figueiredo,et al.  Gradient Projection for Sparse Reconstruction: Application to Compressed Sensing and Other Inverse Problems , 2007, IEEE Journal of Selected Topics in Signal Processing.

[10]  P. L. Combettes,et al.  A proximal decomposition method for solving convex variational inverse problems , 2008, 0807.2617.

[11]  B. Mercier,et al.  A dual algorithm for the solution of nonlinear variational problems via finite element approximation , 1976 .

[12]  J. Moreau.  Proximité et dualité dans un espace hilbertien , 1965 .

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

[14]  Yunmei Chen,et al.  A novel method and fast algorithm for MR image reconstruction with significantly under-sampled data , 2010 .

[15]  Yonggui Zhu,et al.  Fast Alternating Minimization Method for Compressive Sensing MRI under Wavelet Sparsity and TV Sparsity , 2011, 2011 Sixth International Conference on Image and Graphics.

[16]  J. Spingarn.  Partial inverse of a monotone operator , 1983 .

[17]  Xiaoming Huo,et al.  Uncertainty principles and ideal atomic decomposition , 2001, IEEE Trans. Inf. Theory.

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

[19]  D. Gabay.  Applications of the method of multipliers to variational inequalities , 1983 .

[20]  Emmanuel J. Candès,et al.  Robust uncertainty principles: exact signal reconstruction from highly incomplete frequency information , 2004, IEEE Transactions on Information Theory.

[21]  Junfeng Yang,et al.  A New Alternating Minimization Algorithm for Total Variation Image Reconstruction , 2008, SIAM J. Imaging Sci..

[22]  E. Candès,et al.  Sparsity and incoherence in compressive sampling , 2006, math/0611957.

[23]  P. L. Combettes.  Iterative construction of the resolvent of a sum of maximal monotone operators , 2009 .

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

[25]  Junfeng Yang,et al.  A Fast Alternating Direction Method for TVL1-L2 Signal Reconstruction From Partial Fourier Data , 2010, IEEE Journal of Selected Topics in Signal Processing.

[26]  Junzhou Huang,et al.  Efficient MR Image Reconstruction for Compressed MR Imaging , 2010, MICCAI.

[27]  Wotao Yin,et al.  Bregman Iterative Algorithms for (cid:2) 1 -Minimization with Applications to Compressed Sensing ∗ , 2008 .

[28]  Xue-Cheng Tai,et al.  Augmented Lagrangian Method, Dual Methods and Split Bregman Iteration for ROF Model , 2009, SSVM.

[29]  R. Glowinski,et al.  Augmented Lagrangian and Operator-Splitting Methods in Nonlinear Mechanics , 1987 .

[30]  Jean-Jacques Fuchs,et al.  On sparse representations in arbitrary redundant bases , 2004, IEEE Transactions on Information Theory.

[31]  David L Donoho,et al.  Compressed sensing , 2006, IEEE Transactions on Information Theory.

[32]  L. Rudin,et al.  Nonlinear total variation based noise removal algorithms , 1992 .

[33]  Paul Tseng,et al.  A Modified Forward-backward Splitting Method for Maximal Monotone Mappings 1 , 1998 .