ML-reconstruction for TOF-PET with simultaneous estimation of the attenuation factors

In positron emission tomography (PET), attenuation correction is typically done based on information obtained from transmission tomography. Recently, it has been shown that stable maximum-likelihood reconstruction of both the attenuation and the activity from time-of-flight (TOF) PET emission data is possible. Mathematical analysis revealed that the TOF-PET data determine the attenuation correction factors uniquely except for a scale factor. Here, we propose a maximum likelihood algorithm (called MLACF) that jointly estimates the image of the activity distribution and the sinogram with the attenuation factors. This method avoids the reconstruction of the attenuation image. If additive contributions (such as scatter and randoms) can be ignored, the algorithm even does not require storage of the attenuation correction factors. However, in contrast, this algorithm does not impose the consistency of the attenuation sinogram, which may result in increased noise propagation. This paper presents the derivation of the algorithm, an (incomplete) theoretical analysis of the corresponding likelihood function, and first results on 2D and 3D simulations.

[1]  R. D. Badawi,et al.  Self-normalization of emission data in 31) PET , 1999 .

[2]  Alvaro R. De Pierro,et al.  Activity and attenuation recovery from activity data only in emission computed tomography , 2006 .

[3]  Frank Natterer,et al.  Determination of tissue attenuation in emission tomography of optically dense media , 1993 .

[4]  Johan Nuyts,et al.  ML-Reconstruction for TOF-PET With Simultaneous Estimation of the Attenuation Factors , 2014, IEEE Transactions on Medical Imaging.

[5]  V. Sossi,et al.  Effect of normalization method on image uniformity and binding potential estimates on microPET , 2004, IEEE Symposium Conference Record Nuclear Science 2004..

[6]  G.L. Zeng,et al.  A method of attenuation map and emission activity reconstructions from emission data , 1999, 1999 IEEE Nuclear Science Symposium. Conference Record. 1999 Nuclear Science Symposium and Medical Imaging Conference (Cat. No.99CH37019).

[7]  Til Aach,et al.  Simultaneous Reconstruction of Activity and Attenuation for PET/MR , 2011, IEEE Transactions on Medical Imaging.

[8]  R. Badawi,et al.  Self-normalisation of emission data in 3D-PET , 1998, 1998 IEEE Nuclear Science Symposium Conference Record. 1998 IEEE Nuclear Science Symposium and Medical Imaging Conference (Cat. No.98CH36255).

[9]  Donald Sashin,et al.  Efficiency normalization techniques for 3D PET data , 1995, 1995 IEEE Nuclear Science Symposium and Medical Imaging Conference Record.

[10]  Grant T. Gullberg,et al.  Toward accurate attenuation correction in SPECT without transmission measurements , 1997, IEEE Transactions on Medical Imaging.

[11]  C. Watson.  Extension of Single Scatter Simulation to Scatter Correction of Time-of-Flight PET , 2007, IEEE Transactions on Nuclear Science.

[12]  F. Crepaldi,et al.  Activity and Attenuation Reconstruction for Positron Emission Tomography Using Emission Data Only Via Maximum Likelihood and Iterative Data Refinement , 2007, IEEE Transactions on Nuclear Science.

[13]  A Geissbuhler,et al.  A normalization technique for 3D PET data. , 1991, Physics in medicine and biology.

[14]  M. Defrise,et al.  Time-of-flight PET data determine the attenuation sinogram up to a constant , 2012, Physics in medicine and biology.

[15]  Maurizio Conti,et al.  Simultaneous Reconstruction of Activity and Attenuation in Time-of-Flight PET , 2012, IEEE Transactions on Medical Imaging.

[16]  D. Townsend,et al.  Physical and clinical performance of the mCT time-of-flight PET/CT scanner , 2011, Physics in medicine and biology.

[17]  Frederic H Fahey,et al.  Data acquisition in PET imaging. , 2002, Journal of nuclear medicine technology.

[18]  A. V. Bronnikov,et al.  Reconstruction of attenuation map using discrete consistency conditions , 2000, IEEE Transactions on Medical Imaging.

[19]  Yuji Nakamoto,et al.  Respiratory motion artifacts on PET emission images obtained using CT attenuation correction on PET-CT , 2003, European Journal of Nuclear Medicine and Molecular Imaging.

[20]  Y. Censor,et al.  A New Approach to the Emission Computerized Tomography Problem: Simultaneous Calculation of Attenuation and Activity Coefficients , 1979, IEEE Transactions on Nuclear Science.

[21]  Patrick Dupont,et al.  Simultaneous maximum a posteriori reconstruction of attenuation and activity distributions from emission sinograms , 1999, IEEE Transactions on Medical Imaging.

[22]  Paul Kinahan,et al.  Attenuation correction for a combined 3D PET/CT scanner. , 1998, Medical physics.

[23]  René M. Botnar,et al.  A Self-Normalization Reconstruction Technique for PET Scans Using the Positron Emission Data , 2012, IEEE Transactions on Medical Imaging.

[24]  H. Kudo,et al.  A new approach to SPECT attenuation correction without transmission measurements , 2000, 2000 IEEE Nuclear Science Symposium. Conference Record (Cat. No.00CH37149).

[25]  Vladimir Y. Panin,et al.  Reconstruction of uniform sensitivity emission image with partially known axial attenuation information in PET-CT scanners , 2012, 2012 IEEE Nuclear Science Symposium and Medical Imaging Conference Record (NSS/MIC).

[26]  Vladimir Y. Panin,et al.  LSO background radiation as a transmission source using time of flight information , 2013 .

[27]  Johan Nuyts,et al.  Transmission-less attenuation correction in time-of-flight PET: analysis of a discrete iterative algorithm , 2014, Physics in medicine and biology.

[28]  James E. Bowsher,et al.  An EM algorithm for estimating SPECT emission and transmission parameters from emission data only , 2001, IEEE Transactions on Medical Imaging.

[29]  B. H. Hasegawa,et al.  Attenuation Correction Strategies in Emission Tomography , 2006 .

[30]  Alvaro R. De Pierro,et al.  A modified expectation maximization algorithm for penalized likelihood estimation in emission tomography , 1995, IEEE Trans. Medical Imaging.