The effect of regularization in motion compensated PET image reconstruction: a realistic numerical 4D simulation study

Following continuous improvement in PET spatial resolution, respiratory motion correction has become an important task. Two of the most common approaches that utilize all detected PET events to motion-correct PET data are the reconstruct-transform-average method (RTA) and motion-compensated image reconstruction (MCIR). In RTA, separate images are reconstructed for each respiratory frame, subsequently transformed to one reference frame and finally averaged to produce a motion-corrected image. In MCIR, the projection data from all frames are reconstructed by including motion information in the system matrix so that a motion-corrected image is reconstructed directly. Previous theoretical analyses have explained why MCIR is expected to outperform RTA. It has been suggested that MCIR creates less noise than RTA because the images for each separate respiratory frame will be severely affected by noise. However, recent investigations have shown that in the unregularized case RTA images can have fewer noise artefacts, while MCIR images are more quantitatively accurate but have the common salt-and-pepper noise. In this paper, we perform a realistic numerical 4D simulation study to compare the advantages gained by including regularization within reconstruction for RTA and MCIR, in particular using the median-root-prior incorporated in the ordered subsets maximum a posteriori one-step-late algorithm. In this investigation we have demonstrated that MCIR with proper regularization parameters reconstructs lesions with less bias and root mean square error and similar CNR and standard deviation to regularized RTA. This finding is reproducible for a variety of noise levels (25, 50, 100 million counts), lesion sizes (8 mm, 14 mm diameter) and iterations. Nevertheless, regularized RTA can also be a practical solution for motion compensation as a proper level of regularization reduces both bias and mean square error.

[1]  C. J. Thompson,et al.  Motion correction of PET images using multiple acquisition frames , 1997, IEEE Transactions on Medical Imaging.

[2]  Jeffrey A. Fessler Penalized weighted least-squares image reconstruction for positron emission tomography , 1994, IEEE Trans. Medical Imaging.

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

[4]  Ronald H. Huesman,et al.  List mode reconstruction for PET with motion compensation: a simulation study , 2002, Proceedings IEEE International Symposium on Biomedical Imaging.

[5]  Guobao Wang,et al.  Penalized Likelihood PET Image Reconstruction Using Patch-Based Edge-Preserving Regularization , 2012, IEEE Transactions on Medical Imaging.

[6]  Anastasios Gaitanis,et al.  Modeling and Simulation of 4D PET-CT and PET-MR Images. , 2013, PET clinics.

[7]  Xiaodong Tao,et al.  Lesion detectability in motion compensated image reconstruction of respiratory gated PET/CT , 2006, 2006 IEEE Nuclear Science Symposium Conference Record.

[8]  S. Nehmeh,et al.  Respiratory Motion Correction Strategies in Thoracic PET-CT Imaging. , 2013, PET clinics.

[9]  A J Reader,et al.  List-mode-based reconstruction for respiratory motion correction in PET using non-rigid body transformations , 2007, Physics in medicine and biology.

[10]  C Tsoumpas,et al.  Analysis and comparison of two methods for motion correction in PET imaging. , 2012, Medical physics.

[11]  F. Fazio,et al.  Implementation and evaluation of a 3D one-step late reconstruction algorithm for 3D positron emission tomography brain studies using median root prior , 2002, European Journal of Nuclear Medicine and Molecular Imaging.

[12]  T Schaeffter,et al.  Fast generation of 4D PET-MR data from real dynamic MR acquisitions , 2011, Physics in medicine and biology.

[13]  Jeffrey A. Fessler,et al.  Noise Properties of Motion-Compensated Tomographic Image Reconstruction Methods , 2013, IEEE Transactions on Medical Imaging.

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

[15]  John W. Clark,et al.  A motion-incorporated reconstruction method for gated PET studies , 2006, Physics in medicine and biology.

[16]  Donald W. Wilson,et al.  Noise properties of the EM algorithm. I. Theory , 1994 .

[17]  L. Furenlid,et al.  SPECT detectors: the Anger Camera and beyond , 2011, Physics in medicine and biology.

[18]  Ravindra Manjeshwar,et al.  Theoretical Comparison of Motion Correction Techniques for PET Image Reconstruction , 2006, 2006 IEEE Nuclear Science Symposium Conference Record.

[19]  Jeffrey A. Fessler,et al.  Spatial Resolution Properties of Motion-Compensated Tomographic Image Reconstruction Methods , 2012, IEEE Transactions on Medical Imaging.

[20]  Julian C. Matthews,et al.  Bias in iterative reconstruction of low-statistics PET data: benefits of a resolution model , 2009, 2009 IEEE Nuclear Science Symposium Conference Record (NSS/MIC).

[21]  P. Green On Use of the EM Algorithm for Penalized Likelihood Estimation , 1990 .

[22]  Tobias Schaeffter,et al.  Hierarchical adaptive local affine registration for fast and robust respiratory motion estimation , 2011, Medical Image Anal..

[23]  B. Tsui,et al.  Noise properties of the EM algorithm: II. Monte Carlo simulations. , 1994, Physics in medicine and biology.

[24]  S. Alenius,et al.  Bayesian image reconstruction for emission tomography based on median root prior , 1997, European Journal of Nuclear Medicine.

[25]  Florian Büther,et al.  Optimal number of respiratory gates in positron emission tomography: a cardiac patient study. , 2009, Medical physics.

[26]  Jeffrey A. Fessler,et al.  A penalized-likelihood image reconstruction method for emission tomography, compared to postsmoothed maximum-likelihood with matched spatial resolution , 2003, IEEE Transactions on Medical Imaging.

[27]  Ulla Ruotsalainen,et al.  Generalization of median root prior reconstruction , 2002, IEEE Transactions on Medical Imaging.

[28]  Habib Zaidi,et al.  Four-Dimensional Image Reconstruction Strategies in Cardiac-Gated and Respiratory-Gated PET Imaging. , 2013, PET clinics.

[29]  Paul Kinahan,et al.  Noise and signal properties in PSF-based fully 3D PET image reconstruction: an experimental evaluation , 2010, Physics in medicine and biology.

[30]  Nikolaos Dikaios,et al.  Improved motion-compensated image reconstruction for PET using sensitivity correction per respiratory gate and an approximate tube-of-response backprojector. , 2011, Medical physics.

[31]  C. Tsoumpas,et al.  STIR: software for tomographic image reconstruction release 2 , 2012, 2006 IEEE Nuclear Science Symposium Conference Record.

[32]  Klaus P Schäfers,et al.  The promise of nuclear medicine technology: status and future perspective of high-resolution whole-body PET. , 2008, Physica medica : PM : an international journal devoted to the applications of physics to medicine and biology : official journal of the Italian Association of Biomedical Physics.

[33]  P. K. Marsden,et al.  Investigation of MR-Based Attenuation Correction and Motion Compensation for Hybrid PET/MR , 2012, IEEE Transactions on Nuclear Science.

[34]  L. Eriksson,et al.  Comparison of TOF and non-TOF iterative reconstruction at low statistics , 2012, Nuclear Science Symposium and Medical Imaging Conference.