An iterative reconstruction method of complex images using expectation maximization for radial parallel MRI.

In MRI (magnetic resonance imaging), signal sampling along a radial k-space trajectory is preferred in certain applications due to its distinct advantages such as robustness to motion, and the radial sampling can be beneficial for reconstruction algorithms such as parallel MRI (pMRI) due to the incoherency. For radial MRI, the image is usually reconstructed from projection data using analytic methods such as filtered back-projection or Fourier reconstruction after gridding. However, the quality of the reconstructed image from these analytic methods can be degraded when the number of acquired projection views is insufficient. In this paper, we propose a novel reconstruction method based on the expectation maximization (EM) method, where the EM algorithm is remodeled for MRI so that complex images can be reconstructed. Then, to optimize the proposed method for radial pMRI, a reconstruction method that uses coil sensitivity information of multichannel RF coils is formulated. Experiment results from synthetic and in vivo data show that the proposed method introduces better reconstructed images than the analytic methods, even from highly subsampled data, and provides monotonic convergence properties compared to the conjugate gradient based reconstruction method.

[1]  B. F. Logan,et al.  The Fourier reconstruction of a head section , 1974 .

[2]  B. Logan,et al.  Reconstructing Interior Head Tissue from X-Ray Transmissions , 1974 .

[3]  D. Rubin,et al.  Maximum likelihood from incomplete data via the EM - algorithm plus discussions on the paper , 1977 .

[4]  L. Shepp,et al.  Maximum Likelihood Reconstruction for Emission Tomography , 1983, IEEE Transactions on Medical Imaging.

[5]  K. Lange,et al.  EM reconstruction algorithms for emission and transmission tomography. , 1984, Journal of computer assisted tomography.

[6]  P. Green Bayesian reconstructions from emission tomography data using a modified EM algorithm. , 1990, IEEE transactions on medical imaging.

[7]  H. Malcolm Hudson,et al.  Accelerated image reconstruction using ordered subsets of projection data , 1994, IEEE Trans. Medical Imaging.

[8]  C. Byrne Iterative algorithms for deblurring and deconvolution with constraints , 1998 .

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

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

[11]  Walter F Block,et al.  Time‐resolved contrast‐enhanced imaging with isotropic resolution and broad coverage using an undersampled 3D projection trajectory , 2002, Magnetic resonance in medicine.

[12]  Peter M. Jakob,et al.  Direct parallel imaging reconstruction of radially sampled data using GRAPPA with relative shifts , 2003 .

[13]  Bida Zhang,et al.  Convergence behavior of iterative SENSE reconstruction with non‐Cartesian trajectories , 2005, Magnetic resonance in medicine.

[14]  J Velikina,et al.  Highly constrained backprojection for time‐resolved MRI , 2006, Magnetic resonance in medicine.

[15]  Bida Zhang,et al.  An improved iterative SENSE reconstruction method , 2007 .

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

[17]  J. C. Ye,et al.  Projection reconstruction MR imaging using FOCUSS , 2007, Magnetic resonance in medicine.

[18]  HyunWook Park,et al.  Improvement of the diffusion‐weighted images acquired with radial trajectories using projection data regeneration , 2007, Journal of magnetic resonance imaging : JMRI.

[19]  P. Qu,et al.  Interative Back-Projection Reconstruction for Radial SENSE , 2007 .

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

[21]  Sean B Fain,et al.  Iterative projection reconstruction of time‐resolved images using highly‐constrained back‐projection (HYPR) , 2008, Magnetic resonance in medicine.

[22]  Kevin M Johnson,et al.  Ultrashort TE spectroscopic imaging (UTESI) using complex highly‐constrained backprojection with local reconstruction (HYPR LR) , 2009, Magnetic resonance in medicine.

[23]  E. Kholmovski,et al.  Acquisition and reconstruction of undersampled radial data for myocardial perfusion magnetic resonance imaging , 2009, Journal of magnetic resonance imaging : JMRI.

[24]  Andrew J Reader,et al.  AB-OSEM reconstruction for improved Patlak kinetic parameter estimation: a simulation study , 2010, Physics in medicine and biology.

[25]  Wei Lin,et al.  A rapid and robust numerical algorithm for sensitivity encoding with sparsity constraints: Self‐feeding sparse SENSE , 2010, Magnetic resonance in medicine.

[26]  R. Buxton The physics of functional magnetic resonance imaging (fMRI) , 2013, Reports on progress in physics. Physical Society.