Compressed Sensing With Wavelet Domain Dependencies for Coronary MRI: A Retrospective Study

Coronary magnetic resonance imaging (MRI) is a noninvasive imaging modality for diagnosis of coronary artery disease. One of the limitations of coronary MRI is its long acquisition time due to the need of imaging with high spatial resolution and constraints on respiratory and cardiac motions. Compressed sensing (CS) has been recently utilized to accelerate image acquisition in MRI. In this paper, we develop an improved CS reconstruction method, Bayesian least squares-Gaussian scale mixture (BLS-GSM), that uses dependencies of wavelet domain coefficients to reduce the observed blurring and reconstruction artifacts in coronary MRI using traditional l1 regularization. Images of left and right coronary MRI was acquired in 7 healthy subjects with fully-sampled k-space data. The data was retrospectively undersampled using acceleration rates of 2, 4, 6, and 8 and reconstructed using l1 thresholding, l1 minimization and BLS-GSM thresholding. Reconstructed right and left coronary images were compared with fully-sampled reconstructions in vessel sharpness and subjective image quality (1-4 for poor-excellent). Mean square error (MSE) was also calculated for each reconstruction. There were no significant differences between the fully sampled image score versus rate 2, 4, or 6 for BLS-GSM for both right and left coronaries (=N.S.). However, for l1 thresholding significant differences (p <; 0.05) were observed for rates higher than 2 and 4 for right and left coronaries respectively. l1 minimization also yields images with lower scores compared to the reference for rates higher than 4 for both coronaries. These results were consistent with the quantitative vessel sharpness readings. BLS-GSM allows acceleration of coronary MRI with acceleration rates beyond what can be achieved with l1 regularization.

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

[2]  D. Mozaffarian,et al.  Executive summary: heart disease and stroke statistics--2010 update: a report from the American Heart Association. , 2010, Circulation.

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

[4]  Emmanuel J. Candès,et al.  The curvelet transform for image denoising , 2002, IEEE Trans. Image Process..

[5]  F H Epstein,et al.  Adaptive sensitivity encoding incorporating temporal filtering (TSENSE) † , 2001, Magnetic resonance in medicine.

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

[7]  K. T. Block,et al.  Undersampled radial MRI with multiple coils. Iterative image reconstruction using a total variation constraint , 2007, Magnetic resonance in medicine.

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

[9]  Edward H. Adelson,et al.  Shiftable multiscale transforms , 1992, IEEE Trans. Inf. Theory.

[10]  René M. Botnar,et al.  “Soap‐Bubble” visualization and quantitative analysis of 3D coronary magnetic resonance angiograms , 2002, Magnetic resonance in medicine.

[11]  David P. Wipf,et al.  Iterative Reweighted 1 and 2 Methods for Finding Sparse Solutions , 2010, IEEE J. Sel. Top. Signal Process..

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

[13]  Leslie Ying,et al.  Accelerating SENSE using distributed compressed sensing , 2009 .

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

[15]  Himanshu Bhat,et al.  Contrast-Enhanced Whole-Heart Coronary Magnetic Resonance Angiography at 3 T Using Interleaved Echo Planar Imaging , 2010, Investigative radiology.

[16]  W. Manning,et al.  Simultaneous acquisition of spatial harmonics (SMASH): Fast imaging with radiofrequency coil arrays , 1997, Magnetic resonance in medicine.

[17]  René M. Botnar,et al.  Double-oblique free-breathing high resolution three-dimensional coronary magnetic resonance angiography. , 1999, Journal of the American College of Cardiology.

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

[19]  C. Sidney Burrus,et al.  Nonlinear processing of a shift-invariant discrete wavelet transform (DWT) for noise reduction , 1995, Defense, Security, and Sensing.

[20]  Leon Axel,et al.  Combination of Compressed Sensing and Parallel Imaging for Highly-Accelerated 3 D First-Pass Cardiac Perfusion MRI , 2009 .

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

[22]  I. Johnstone,et al.  Ideal spatial adaptation by wavelet shrinkage , 1994 .

[23]  Thoralf Niendorf,et al.  Short breath‐hold, volumetric coronary MR angiography employing steady‐state free precession in conjunction with parallel imaging , 2005, Magnetic resonance in medicine.

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

[25]  Thoralf Niendorf,et al.  Toward single breath‐hold whole‐heart coverage coronary MRA using highly accelerated parallel imaging with a 32‐channel MR system , 2006, Magnetic resonance in medicine.

[26]  C. Robert The Bayesian choice : a decision-theoretic motivation , 1996 .

[27]  Michael P. Friedlander,et al.  Probing the Pareto Frontier for Basis Pursuit Solutions , 2008, SIAM J. Sci. Comput..

[28]  Peter M. Jakob,et al.  A Combination of Nonconvex Compressed Sensing and GRAPPA (CS-GRAPPA) , 2009 .

[29]  Albert Macovski,et al.  Coronary Angiography with Magnetization‐Prepared T2 Contrast , 1995, Magnetic resonance in medicine.

[30]  X. Rabasseda,et al.  A report from the American Heart Association Scientific Sessions 2011 (November 12-16, 2011, Orlando, Florida, USA). , 2012, Drugs of today.

[31]  Richard Baraniuk,et al.  The Dual-tree Complex Wavelet Transform , 2007 .

[32]  Volkan Cevher,et al.  Model-Based Compressive Sensing , 2008, IEEE Transactions on Information Theory.

[33]  Martin J. Wainwright,et al.  Image denoising using scale mixtures of Gaussians in the wavelet domain , 2003, IEEE Trans. Image Process..

[34]  Debiao Li,et al.  Whole‐heart coronary magnetic resonance angiography at 3 Tesla in 5 minutes with slow infusion of Gd‐BOPTA, a high‐relaxivity clinical contrast agent , 2007, Magnetic resonance in medicine.

[35]  Shrikanth S. Narayanan,et al.  Accelerated three‐dimensional upper airway MRI using compressed sensing , 2009, Magnetic resonance in medicine.

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

[37]  M. Lustig,et al.  Improving non‐contrast‐enhanced steady‐state free precession angiography with compressed sensing , 2009, Magnetic resonance in medicine.

[38]  M. West On scale mixtures of normal distributions , 1987 .

[39]  Bob S. Hu,et al.  Real‐time interactive coronary MRA , 2001, Magnetic resonance in medicine.

[40]  Rachid Deriche,et al.  Fast algorithms for low-level vision , 1988, [1988 Proceedings] 9th International Conference on Pattern Recognition.

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

[42]  Jerome M. Shapiro,et al.  Embedded image coding using zerotrees of wavelet coefficients , 1993, IEEE Trans. Signal Process..

[43]  Peter Boesiger,et al.  Sensitivity‐encoded coronary MRA at 3T , 2004, Magnetic resonance in medicine.

[44]  René M. Botnar,et al.  Coronary magnetic resonance angiography for the detection of coronary stenoses. , 2001, The New England journal of medicine.