Iterative Schemes to Solve Low-Dimensional Calibration Equations in Parallel MR Image Reconstruction with GRAPPA

GRAPPA (Generalized Autocalibrating Partially Parallel Acquisition) is a widely used parallel MRI reconstruction technique. The processing of data from multichannel receiver coils may increase the storage and computational requirements of GRAPPA reconstruction. Random projection on GRAPPA (RP-GRAPPA) uses random projection (RP) method to overcome the computational overheads of solving large linear equations in the calibration phase of GRAPPA, saving reconstruction time. However, RP-GRAPPA compromises the reconstruction accuracy in case of large reductions in the dimensions of calibration equations. In this paper, we present the implementation of GRAPPA reconstruction method using potential iterative solvers to estimate the reconstruction coefficients from the randomly projected calibration equations. Experimental results show that the proposed methods withstand the reconstruction accuracy (visually and quantitatively) against large reductions in the dimension of linear equations, when compared with RP-GRAPPA reconstruction. Particularly, the proposed method using conjugate gradient for least squares (CGLS) demonstrates more savings in the computational time of GRAPPA, without significant loss in the reconstruction accuracy, when compared with RP-GRAPPA. It is also demonstrated that the proposed method using CGLS complements the channel compression method for reducing the computational complexities associated with higher channel count, thereby resulting in additional memory savings and speedup.

[1]  Hammad Omer,et al.  A new parameterized architectural design for SENSE reconstruction , 2014 .

[2]  Allen Klinger Approximate pseudoinverse solutions to ill-conditioned linear systems , 1968 .

[3]  Heng Tao Shen,et al.  Principal Component Analysis , 2009, Encyclopedia of Biometrics.

[4]  Zhouping Xin,et al.  Step-sizes for the gradient method , 2008 .

[5]  Enrico Magli,et al.  Analysis of One-Time Random Projections for Privacy Preserving Compressed Sensing , 2016, IEEE Transactions on Information Forensics and Security.

[6]  M. Lustig,et al.  SPIRiT: Iterative self‐consistent parallel imaging reconstruction from arbitrary k‐space , 2010, Magnetic resonance in medicine.

[7]  Christina Triantafyllou,et al.  A 128‐channel receive‐only cardiac coil for highly accelerated cardiac MRI at 3 Tesla , 2008, Magnetic resonance in medicine.

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

[9]  Donglai Huo,et al.  Robust GRAPPA reconstruction , 2006, 3rd IEEE International Symposium on Biomedical Imaging: Nano to Macro, 2006..

[10]  Kenneth Ward Church,et al.  Margin-constrained random projections and very sparse random projections , 2006 .

[11]  Dimitris Achlioptas,et al.  Database-friendly random projections: Johnson-Lindenstrauss with binary coins , 2003, J. Comput. Syst. Sci..

[12]  Jianghong Wei,et al.  Fuzzy -Means and Cluster Ensemble with Random Projection for Big Data Clustering , 2016 .

[13]  Michael Lustig,et al.  Coil compression for accelerated imaging with Cartesian sampling , 2013, Magnetic resonance in medicine.

[14]  Xuelong Li,et al.  L1-Norm-Based 2DPCA , 2010, IEEE Transactions on Systems, Man, and Cybernetics, Part B (Cybernetics).

[15]  Scott B King,et al.  Optimum SNR data compression in hardware using an Eigencoil array , 2010, Magnetic resonance in medicine.

[16]  F.L. Fontaine,et al.  A Reconfigurable FPGA-based 16-Channel Front-End for MRI , 2006, 2006 Fortieth Asilomar Conference on Signals, Systems and Computers.

[17]  Michael A. Saunders,et al.  LSQR: An Algorithm for Sparse Linear Equations and Sparse Least Squares , 1982, TOMS.

[18]  Jyh-Shing Roger Jang,et al.  ANFIS: adaptive-network-based fuzzy inference system , 1993, IEEE Trans. Syst. Man Cybern..

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

[20]  Asoke K. Nandi,et al.  Blind adaptive carrier phase recovery for QAM systems , 2016, Digit. Signal Process..

[21]  Michael Schacht Hansen,et al.  Gadgetron: An open source framework for medical image reconstruction , 2013, Magnetic resonance in medicine.

[22]  Alberto D. Pascual-Montano,et al.  A survey of dimensionality reduction techniques , 2014, ArXiv.

[23]  Peter Boesiger,et al.  Array compression for MRI with large coil arrays , 2007, Magnetic resonance in medicine.

[24]  Santiago Aja-Fernández,et al.  Acquisition and Reconstruction of Magnetic Resonance Imaging , 2016 .

[25]  H. B. Curry The method of steepest descent for non-linear minimization problems , 1944 .

[26]  Heikki Mannila,et al.  Random projection in dimensionality reduction: applications to image and text data , 2001, KDD '01.

[27]  Wei Lin,et al.  A hybrid method for more efficient channel‐by‐channel reconstruction with many channels , 2012, Magnetic resonance in medicine.

[28]  Jim Ji,et al.  Channel reduction in massive array parallel MRI , 2009, 2009 Annual International Conference of the IEEE Engineering in Medicine and Biology Society.

[29]  L. Ying,et al.  Nonlinear GRAPPA: A kernel approach to parallel MRI reconstruction , 2012, Magnetic resonance in medicine.

[30]  Pascal Vincent,et al.  Representation Learning: A Review and New Perspectives , 2012, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[31]  Hammad Omer,et al.  Parallel MRI Reconstruction Algorithm Implementation on GPU , 2016 .

[32]  Kenneth Ward Church,et al.  Very sparse random projections , 2006, KDD '06.

[33]  Martin Blaimer,et al.  General formulation for quantitative G‐factor calculation in GRAPPA reconstructions , 2009, Magnetic resonance in medicine.

[34]  Lawrence L Wald,et al.  Massively parallel MRI detector arrays. , 2013, Journal of magnetic resonance.

[35]  M. Hestenes,et al.  Methods of conjugate gradients for solving linear systems , 1952 .

[36]  J. Shewchuk An Introduction to the Conjugate Gradient Method Without the Agonizing Pain , 1994 .

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

[38]  L. Ying,et al.  Fast GRAPPA reconstruction with random projection , 2015, Magnetic resonance in medicine.

[39]  Gianluca Bontempi,et al.  New Routes from Minimal Approximation Error to Principal Components , 2008, Neural Processing Letters.

[40]  Yu Li,et al.  A software channel compression technique for faster reconstruction with many channels. , 2008, Magnetic resonance imaging.

[41]  Michael E. Tipping,et al.  Probabilistic Principal Component Analysis , 1999 .

[42]  Nicole Seiberlich,et al.  Parallel MR imaging , 2012, Journal of magnetic resonance imaging : JMRI.

[43]  Yuchou Chang,et al.  Efficient GRAPPA reconstruction using random projection , 2013, 2013 IEEE 10th International Symposium on Biomedical Imaging.

[44]  Nicole Seiberlich,et al.  Improved radial GRAPPA calibration for real‐time free‐breathing cardiac imaging , 2011, Magnetic resonance in medicine.

[45]  Martin Blaimer,et al.  Reconstruction of undersampled non‐Cartesian data sets using pseudo‐Cartesian GRAPPA in conjunction with GROG , 2008, Magnetic resonance in medicine.

[46]  Stefan Skare,et al.  Comparison of reconstruction accuracy and efficiency among autocalibrating data‐driven parallel imaging methods , 2008, Magnetic resonance in medicine.