Evaluating sparsity penalty functions for combined compressed sensing and parallel MRI

The combination of compressed sensing (CS) and parallel magnetic resonance (MR) imaging enables further scan acceleration via undersampling than previously feasible. While many of these methods incorporate similar styles of CS, there remains significant variation in the particular choice of function used to promote sparsity. Having developed SpRING, a framework for combining CS and GRAPPA, a parallel MR image reconstruction method, we view the choice of penalty function as a design choice rather than a defining feature of the algorithm. For both simulated and real data, we compare different sparsity penalty functions to the empirical distribution of the reference images. Then, we perform reconstructions on uniformly undersampled data using a variety of penalty functions to illustrate the impact appropriately choosing the penalty function has on the performance of SpRING. These experiments demonstrate the importance of choosing an appropriate penalty function and how such a choice may differ between simulated data and real data.

[1]  P. Roemer,et al.  The NMR phased array , 1990, Magnetic resonance in medicine.

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

[3]  D. Donoho,et al.  SPIR-iT : Autocalibrating Parallel Imaging Compressed Sensing , 2008 .

[4]  Robin M Heidemann,et al.  Generalized autocalibrating partially parallel acquisitions (GRAPPA) , 2002, Magnetic resonance in medicine.

[5]  Vivek K. Goyal,et al.  Simultaneously Sparse Solutions to Linear Inverse Problems with Multiple System Matrices and a Single Observation Vector , 2009, SIAM J. Sci. Comput..

[6]  Emmanuel J. Candès,et al.  Quantitative Robust Uncertainty Principles and Optimally Sparse Decompositions , 2004, Found. Comput. Math..

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

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

[9]  V. K. Goyal,et al.  Combining nonconvex compressed sensing and GRAPPA using the nullspace method , 2009 .

[10]  J. R. Polimeni,et al.  Nullspace compressed sensing for accelerated imaging , 2008 .

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

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