Compressed sensing in MRI with a Markov random field prior for spatial clustering of subband coefficients

Recent work in compressed sensing of magnetic resonance images (CS-MRI) concentrates on encoding structured sparsity in acquisition or in the reconstruction stages. Subband coefficients of typical images obey a certain structure, which can be viewed in terms of fixed groups (like wavelet trees) or statistically (certain configurations are more likely than others). Approaches using wavelet tree-sparsity have already demonstrated excellent performance in MRI. However, the use of statistical models for spatial clustering of the subband coefficients has not been studied well in CS-MRI yet, although the potentials of such an approach have been indicated. In this paper, we design a practical reconstruction algorithm as a variant of the proximal splitting methods, making use of a Markov Random Field prior model for spatial clustering of subband coefficients. The results for different undersampling patterns demonstrate an improved reconstruction performance compared to both standard CS-MRI methods and methods based on wavelet tree sparsity.

[1]  Olgica Milenkovic,et al.  Subspace Pursuit for Compressive Sensing Signal Reconstruction , 2008, IEEE Transactions on Information Theory.

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

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

[4]  CONSTRUCTION OF REGULAR AND IRREGULAR SHEARLET FRAMES , 2007 .

[5]  Aleksandra Pizurica,et al.  A joint inter- and intrascale statistical model for Bayesian wavelet based image denoising , 2002, IEEE Trans. Image Process..

[6]  Junzhou Huang,et al.  Compressive Sensing MRI with Wavelet Tree Sparsity , 2012, NIPS.

[7]  A. Pižurica,et al.  Efficient design of a low redundant Discrete Shearlet Transform , 2009, 2009 International Workshop on Local and Non-Local Approximation in Image Processing.

[8]  Michael I. Jordan,et al.  Loopy Belief Propagation for Approximate Inference: An Empirical Study , 1999, UAI.

[9]  Fei Yang,et al.  Compressed magnetic resonance imaging based on wavelet sparsity and nonlocal total variation , 2012, 2012 9th IEEE International Symposium on Biomedical Imaging (ISBI).

[10]  M. Ramasubba Reddy,et al.  An iterative hard thresholding algorithm for CS MRI , 2012, Medical Imaging.

[11]  Deanna Needell,et al.  CoSaMP: Iterative signal recovery from incomplete and inaccurate samples , 2008, ArXiv.

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

[13]  Thomas Blumensath,et al.  Accelerated iterative hard thresholding , 2012, Signal Process..

[14]  Junzhou Huang,et al.  Exploiting the wavelet structure in compressed sensing MRI. , 2014, Magnetic resonance imaging.

[15]  J. Besag On the Statistical Analysis of Dirty Pictures , 1986 .

[16]  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.

[17]  24th European Signal Processing Conference, EUSIPCO 2016, Budapest, Hungary, August 29 - September 2, 2016 , 2016, European Signal Processing Conference.

[18]  Michael A. Saunders,et al.  Atomic Decomposition by Basis Pursuit , 1998, SIAM J. Sci. Comput..

[19]  Volkan Cevher,et al.  Structured Sparsity: Discrete and Convex approaches , 2015, ArXiv.

[20]  Hiêp Quang Luong,et al.  Augmented Lagrangian based reconstruction of non-uniformly sub-Nyquist sampled MRI data , 2011, Signal Process..

[21]  Wilfried Philips,et al.  On structured sparsity and selected applications in tomographic imaging , 2011, Optical Engineering + Applications.

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

[23]  Volkan Cevher,et al.  Sparse Signal Recovery Using Markov Random Fields , 2008, NIPS.

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

[25]  José M. Bioucas-Dias,et al.  An Augmented Lagrangian Approach to the Constrained Optimization Formulation of Imaging Inverse Problems , 2009, IEEE Transactions on Image Processing.

[26]  Stan Z. Li Markov Random Field Modeling in Image Analysis , 2009, Advances in Pattern Recognition.

[27]  Volkan Cevher,et al.  Low-Dimensional Models for Dimensionality Reduction and Signal Recovery: A Geometric Perspective , 2010, Proceedings of the IEEE.

[28]  José M. Bioucas-Dias,et al.  A fast algorithm for the constrained formulation of compressive image reconstruction and other linear inverse problems , 2009, 2010 IEEE International Conference on Acoustics, Speech and Signal Processing.

[29]  Ben Adcock,et al.  The quest for optimal sampling: Computationally efficient, structure-exploiting measurements for compressed sensing , 2014, ArXiv.