A Maximum Likelihood Approach to Parallel Imaging With Coil Sensitivity Noise

Parallel imaging is a powerful technique to speed up magnetic resonance (MR) image acquisition via multiple coils. Both the received signal of each coil and its sensitivity map, which describes its spatial response, are needed during reconstruction. Widely used schemes such as SENSE assume that sensitivity maps of the coils are noiseless while the only errors are in coil outputs. In practice, however, sensitivity maps are subject to a wide variety of errors. At first glance, sensitivity noise appears to result in an errors-in-variables problem of the kind that is typically solved using total least squares (TLSs). However, existing TLS algorithms are in general inappropriate for the specific type of block structure that arises in parallel imaging. In this paper, we take a maximum likelihood approach to the problem of parallel imaging in the presence of independent Gaussian sensitivity noise. This results in a quasi-quadratic objective function, which can be efficiently minimized. Experimental evidence suggests substantial gains over conventional SENSE, especially in nonideal imaging conditions like low signal-to-noise ratio (SNR), high g-factors and large acceleration, using sensitivity maps suffering from misalignment, ringing, and random noise.

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

[2]  Sabine Van Huffel,et al.  On the equivalence of constrained total least squares and structured total least squares , 1996, IEEE Trans. Signal Process..

[3]  Nikolas P. Galatsanos,et al.  Regularized constrained total least-squares image restoration , 1994, Other Conferences.

[4]  Sabine Van Huffel,et al.  Total Least Squares and Errors-in-variables Modeling , 2007, Comput. Stat. Data Anal..

[5]  Yi Wang,et al.  Description of parallel imaging in MRI using multiple coils , 2000, Magnetic resonance in medicine.

[6]  N. Pelc,et al.  Making Better SENSE : Wavelet Denoising , Tikhonov Regularization , and Total Least Squares , 2001 .

[7]  C A McKenzie,et al.  Coil‐by‐coil image reconstruction with SMASH , 2001, Magnetic resonance in medicine.

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

[9]  Raja Muthupillai,et al.  Implications of SENSE MR in routine clinical practice. , 2003, European journal of radiology.

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

[11]  Peter Boesiger,et al.  2D sense for faster 3D MRI , 2002, Magnetic Resonance Materials in Physics, Biology and Medicine.

[12]  Robin M Heidemann,et al.  SMASH, SENSE, PILS, GRAPPA: How to Choose the Optimal Method , 2004, Topics in magnetic resonance imaging : TMRI.

[13]  R. Fletcher Practical Methods of Optimization , 1988 .

[14]  William H. Press,et al.  Numerical Recipes in C, 2nd Edition , 1992 .

[15]  Roland Bammer,et al.  Diffusion tensor imaging using single‐shot SENSE‐EPI , 2002, Magnetic resonance in medicine.

[16]  Nikolas P. Galatsanos,et al.  Iterative maximum a posteriori (MAP) restoration from partially-known blur for tomographic reconstruction , 1995, Proceedings., International Conference on Image Processing.

[17]  William H. Press,et al.  Numerical recipes in C , 2002 .

[18]  S. Schoenberg,et al.  Single breath-hold real-time cine MR imaging: improved temporal resolution using generalized autocalibrating partially parallel acquisition (GRAPPA) algorithm , 2003, European Radiology.

[19]  Demetri Terzopoulos,et al.  Regularization of Inverse Visual Problems Involving Discontinuities , 1986, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[20]  Daniel K. Sodickson,et al.  Recent advances in image reconstruction, coil sensitivity calibration, and coil array design for SMASH and generalized parallel MRI , 2001, Magnetic Resonance Materials in Physics, Biology and Medicine.

[21]  P. Boesiger,et al.  Electrodynamics and ultimate SNR in parallel MR imaging , 2004, Magnetic resonance in medicine.

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

[23]  M. Peruggia Total Least Squares and Errors-in-Variables Modeling: Analysis, Algorithms and Applications , 2003 .

[24]  Jerry M. Mendel,et al.  The constrained total least squares technique and its applications to harmonic superresolution , 1991, IEEE Trans. Signal Process..

[25]  Jacques A. de Guise,et al.  A method for modeling noise in medical images , 2004, IEEE Transactions on Medical Imaging.

[26]  Gene H. Golub,et al.  Matrix computations , 1983 .

[27]  K. Kwong,et al.  Parallel imaging reconstruction using automatic regularization , 2004, Magnetic resonance in medicine.

[28]  S. Huffel,et al.  Total Least Squares and Errors-in-Variables Modeling : Analysis, Algorithms and Applications , 2002 .