Exploiting rank deficiency and transform domain sparsity for MR image reconstruction.

The reconstruction of magnetic resonance (MR) images from the partial samples of their k-space data using compressed sensing (CS)-based methods has generated a lot of interest in recent years. To reconstruct the MR images, these techniques exploit the sparsity of the image in a transform domain (wavelets, total variation, etc.). In a recent work, it has been shown that it is also possible to reconstruct MR images by exploiting their rank deficiency. In this work, it will be shown that, instead of exploiting the sparsity of the image or rank deficiency alone, better reconstruction results can be achieved by combining transform domain sparsity with rank deficiency. To reconstruct an MR image using its transform domain sparsity and its rank deficiency, this work proposes a combined l(1)-norm (of the transform coefficients) and nuclear norm (of the MR image matrix) minimization problem. Since such an optimization problem has not been encountered before, this work proposes and derives a first-order algorithm to solve it. The reconstruction results show that the proposed approach yields significant improvements, in terms of both visual quality as well as the signal to noise ratio, over previous works that reconstruct MR images either by exploiting rank deficiency or by the standard CS-based technique popularly known as the 'Sparse MRI.'

[1]  A. Majumdar,et al.  An algorithm for sparse MRI reconstruction by Schatten p-norm minimization. , 2011, Magnetic resonance imaging.

[2]  G. Sapiro,et al.  A collaborative framework for 3D alignment and classification of heterogeneous subvolumes in cryo-electron tomography. , 2013, Journal of structural biology.

[3]  Xiaodong Li,et al.  Dense error correction for low-rank matrices via Principal Component Pursuit , 2010, 2010 IEEE International Symposium on Information Theory.

[4]  Mathews Jacob,et al.  Accelerated Dynamic MRI Exploiting Sparsity and Low-Rank Structure: k-t SLR , 2011, IEEE Transactions on Medical Imaging.

[5]  Rick Chartrand,et al.  Exact Reconstruction of Sparse Signals via Nonconvex Minimization , 2007, IEEE Signal Processing Letters.

[6]  Rick Chartrand,et al.  Fast algorithms for nonconvex compressive sensing: MRI reconstruction from very few data , 2009, 2009 IEEE International Symposium on Biomedical Imaging: From Nano to Macro.

[7]  Krishnamurthy Dvijotham,et al.  A nullspace analysis of the nuclear norm heuristic for rank minimization , 2010, 2010 IEEE International Conference on Acoustics, Speech and Signal Processing.

[8]  Paulo J. G. Lisboa,et al.  An adaptive hybrid image compression method and its application to medical images , 2008, 2008 5th IEEE International Symposium on Biomedical Imaging: From Nano to Macro.

[9]  Pablo A. Parrilo,et al.  Rank-Sparsity Incoherence for Matrix Decomposition , 2009, SIAM J. Optim..

[10]  Weiyu Xu,et al.  Null space conditions and thresholds for rank minimization , 2011, Math. Program..

[11]  Javier Villanueva-Oller,et al.  SVD lossy adaptive encoding of 3D digital images for ROI progressive transmission , 2010, Image Vis. Comput..

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

[13]  E. Candès,et al.  Stable signal recovery from incomplete and inaccurate measurements , 2005, math/0503066.

[14]  Wei Lin,et al.  A rapid and robust numerical algorithm for sensitivity encoding with sparsity constraints: Self‐feeding sparse SENSE , 2010, Magnetic resonance in medicine.

[15]  Armando Manduca,et al.  Highly Undersampled Magnetic Resonance Image Reconstruction via Homotopic $\ell_{0}$ -Minimization , 2009, IEEE Transactions on Medical Imaging.

[16]  Panos Nasiopoulos,et al.  This article has been accepted for publication in a future issue of this journal, but has not been fully edited. Content may change prior to final publication. 1 Symmetry-Based Scalable Lossless Compression of 3D Medical Image Data , 2022 .

[17]  D. Donoho For most large underdetermined systems of linear equations the minimal 𝓁1‐norm solution is also the sparsest solution , 2006 .

[18]  A. Manduca,et al.  Wavelet compression of medical images. , 1998, Radiology.

[19]  L. Ying,et al.  Accelerating SENSE using compressed sensing , 2009, Magnetic resonance in medicine.

[20]  E. DiBella,et al.  Reconstruction of dynamic contrast enhanced magnetic resonance imaging of the breast with temporal constraints. , 2010, Magnetic resonance imaging.

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

[22]  Justin P. Haldar,et al.  Spatiotemporal imaging with partially separable functions: A matrix recovery approach , 2010, 2010 IEEE International Symposium on Biomedical Imaging: From Nano to Macro.

[23]  Emmanuel J. Candès,et al.  A Singular Value Thresholding Algorithm for Matrix Completion , 2008, SIAM J. Optim..

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

[25]  Mário A. T. Figueiredo,et al.  Signal restoration with overcomplete wavelet transforms: comparison of analysis and synthesis priors , 2009, Optical Engineering + Applications.

[26]  R. Frayne,et al.  MR image reconstruction of sparsely sampled 3D k-space data by projection-onto-convex sets. , 2006, Magnetic resonance imaging.

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

[28]  Felix J. Herrmann,et al.  Non-parametric seismic data recovery with curvelet frames , 2008 .

[29]  George Carayannis,et al.  MRI scan time reduction through non-uniform sampling and SVD-based estimation , 1996, Signal Process..

[30]  Peter Boesiger,et al.  Compressed sensing in dynamic MRI , 2008, Magnetic resonance in medicine.

[31]  Huai Li,et al.  Optimization of wavelet decomposition for image compression and feature preservation , 2003, IEEE Transactions on Medical Imaging.

[32]  Weiyu Xu,et al.  Necessary and sufficient conditions for success of the nuclear norm heuristic for rank minimization , 2008, 2008 47th IEEE Conference on Decision and Control.

[33]  Yi Ma,et al.  Robust principal component analysis? , 2009, JACM.

[34]  M Braun,et al.  Compression and reconstruction of MRI images using 2D DCT. , 1992, Magnetic resonance imaging.

[35]  Amel Benazza-Benyahia,et al.  A wavelet-based regularized reconstruction algorithm for SENSE parallel MRI with applications to neuroimaging , 2011, Medical Image Anal..

[36]  Shiqian Ma,et al.  An efficient algorithm for compressed MR imaging using total variation and wavelets , 2008, 2008 IEEE Conference on Computer Vision and Pattern Recognition.