Do scatter and random corrections affect the errors in kinetic parameters in dynamic PET? - a Monte Carlo study

Dynamic positron emission tomography (PET) data can be evaluated by compartmental models, yielding model specific kinetic parameters. For the parameters to be of quantitative use however, understanding and estimation of errors and uncertainties associated with them are crucial. The aim in this study was to investigate the effects of the inclusion of scattered and random counts and their respective corrections on kinetic parameter errors. The MC software GATE was used to simulate two dynamic PET scans of a phantom containing three regions; blood, tissue and a static background. The two sets of time-activity-curves (TACs) used were generated for a 2-tissue compartment model with preset parameter values (K1, k2, k3, k4 and Va). The PET data was reconstructed into 19 frames by both ordered-subset expectation maximization (OSEM) and 3D filtered backprojection with reprojection (3DFBPRP) with normalization and additional corrections (A=attenuation, R=random, S=scatter, C=correction): True counts (AC), true+random counts (ARC), true+scattered counts (ASC) and total counts (ARSC). The results show that parameter estimates from true counts (AC), true+random counts (ARC), true+scattered counts (ASC) and total counts (ARSC) were not significantly different, with the exception of Va where the bias increased with added corrections. Thus, the inclusion of and correction for scattered and random counts did not affect the bias in parameter estimates K1, k2, k3, k4 and Ki. Uncorrected total counts (only AC) resulted in biases of hundreds or even thousands of percent, emphasizing the need for proper corrections. Reconstructions with 3DFBPRP resulted in overall 20-40% less biased estimates compared to OSEM.

[1]  Yuichi Kimura,et al.  PET kinetic analysis: error consideration of quantitative analysis in dynamic studies , 2008, Annals of nuclear medicine.

[2]  Simon R. Cherry,et al.  A comparison of PET detector modules employing rectangular and round photomultiplier tubes , 1995 .

[3]  Mark Muzi,et al.  Kinetic modeling of 3'-deoxy-3'-fluorothymidine in somatic tumors: mathematical studies. , 2005, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[4]  S. Nehmeh,et al.  Validation of GATE Monte Carlo simulations of the GE Advance/Discovery LS PET scanners. , 2005, Medical physics.

[5]  Sadek Nehmeh,et al.  SU‐FF‐I‐54: A Monte Carlo Model of the Discovery ST PET Scanner , 2006 .

[6]  Hiroshi Watabe,et al.  Monte Carlo estimation of scatter effects on quantitative myocardial blood flow and perfusable tissue fraction using 3D-PET and 15O-water , 2012, Physics in medicine and biology.

[7]  Kris Thielemans,et al.  STIR: Software for Tomographic Image Reconstruction Release 2 , 2006 .

[8]  Simon R. Cherry,et al.  Effects of scatter on model parameter estimates in 3D PET studies of the human brain , 1995 .

[9]  Lula Rosso,et al.  Evaluation of limited blood sampling population input approaches for kinetic quantification of [18F]fluorothymidine PET data , 2012, EJNMMI Research.

[10]  S. Libutti,et al.  Parametric images of blood flow in oncology PET studies using [15O]water. , 2000, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[11]  C Lartizien,et al.  GATE: a simulation toolkit for PET and SPECT. , 2004, Physics in medicine and biology.