Compartment modeling of dynamic brain PET--the impact of scatter corrections on parameter errors.

PURPOSE The aim of this study was to investigate the effect of scatter and its correction on kinetic parameters in dynamic brain positron emission tomography (PET) tumor imaging. The 2-tissue compartment model was used, and two different reconstruction methods and two scatter correction (SC) schemes were investigated. METHODS The gate Monte Carlo (MC) software was used to perform 2 × 15 full PET scan simulations of a voxelized head phantom with inserted tumor regions. The two sets of kinetic parameters of all tissues were chosen to represent the 2-tissue compartment model for the tracer 3'-deoxy-3'-((18)F)fluorothymidine (FLT), and were denoted FLT1 and FLT2. PET data were reconstructed with both 3D filtered back-projection with reprojection (3DRP) and 3D ordered-subset expectation maximization (OSEM). Images including true coincidences with attenuation correction (AC) and true+scattered coincidences with AC and with and without one of two applied SC schemes were reconstructed. Kinetic parameters were estimated by weighted nonlinear least squares fitting of image derived time-activity curves. Calculated parameters were compared to the true input to the MC simulations. RESULTS The relative parameter biases for scatter-eliminated data were 15%, 16%, 4%, 30%, 9%, and 7% (FLT1) and 13%, 6%, 1%, 46%, 12%, and 8% (FLT2) for K1, k2, k3, k4, Va, and Ki, respectively. As expected, SC was essential for most parameters since omitting it increased biases by 10 percentage points on average. SC was not found necessary for the estimation of Ki and k3, however. There was no significant difference in parameter biases between the two investigated SC schemes or from parameter biases from scatter-eliminated PET data. Furthermore, neither 3DRP nor OSEM yielded the smallest parameter biases consistently although there was a slight favor for 3DRP which produced less biased k3 and Ki estimates while OSEM resulted in a less biased Va. The uncertainty in OSEM parameters was about 26% (FLT1) and 12% (FLT2) larger than for 3DRP although identical postfilters were applied. CONCLUSIONS SC was important for good parameter estimations. Both investigated SC schemes performed equally well on average and properly corrected for the scattered radiation, without introducing further bias. Furthermore, 3DRP was slightly favorable over OSEM in terms of kinetic parameter biases and SDs.

[1]  D. Newport,et al.  Evaluation of simulation-based scatter correction for 3-D PET cardiac imaging , 1995 .

[2]  R. Boellaard,et al.  Experimental and clinical evaluation of iterative reconstruction (OSEM) in dynamic PET: quantitative characteristics and effects on kinetic modeling. , 2001, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[3]  L. Schwartz,et al.  Promise and pitfalls of quantitative imaging in oncology clinical trials. , 2012, Magnetic resonance imaging.

[4]  T. Cloughesy,et al.  18 F-fluorothymidine kinetics of malignant brain tumors , 2007 .

[5]  B. Gulyás,et al.  New PET radiopharmaceuticals beyond FDG for brain tumor imaging. , 2012, The quarterly journal of nuclear medicine and molecular imaging : official publication of the Italian Association of Nuclear Medicine (AIMN) [and] the International Association of Radiopharmacology (IAR), [and] Section of the Society of....

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

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

[8]  J. Humm,et al.  Impact of attenuation and scatter correction in estimating tumor hypoxia-related kinetic parameters for FMISO dynamic animal-PET imaging , 2008, 2008 IEEE Nuclear Science Symposium Conference Record.

[9]  Michael E. Phelps,et al.  18F-FLT    and 18F-FDOPA PET kinetics in recurrent brain tumors , 2014, European Journal of Nuclear Medicine and Molecular Imaging.

[10]  D. Louis Collins,et al.  Design and construction of a realistic digital brain phantom , 1998, IEEE Transactions on Medical Imaging.

[11]  Janet F. Eary,et al.  NCI-Sponsored Trial for the Evaluation of Safety and Preliminary Efficacy of 3′-Deoxy-3′-[18F]fluorothymidine (FLT) as a Marker of Proliferation in Patients with Recurrent Gliomas: Preliminary Efficacy Studies , 2009, Molecular Imaging and Biology.

[12]  Jan Axelsson,et al.  The influence of time sampling scheme on kinetic parameters obtained from compartmental modeling of a dynamic PET study - A Monte Carlo study , 2012, 2012 IEEE Nuclear Science Symposium and Medical Imaging Conference Record (NSS/MIC).

[13]  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.

[14]  T. Cloughesy,et al.  3′-Deoxy-3′-18F-Fluorothymidine PET and MRI for Early Survival Predictions in Patients with Recurrent Malignant Glioma Treated with Bevacizumab , 2012, The Journal of Nuclear Medicine.

[15]  Otto Muzik,et al.  Imaging proliferation in vivo with [F-18]FLT and positron emission tomography , 1998, Nature Medicine.

[16]  Paul Kinahan,et al.  Analytic 3D image reconstruction using all detected events , 1989 .

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

[18]  R. Buchert,et al.  Evaluation of non-uniform weighting in non-linear regression for pharmacokinetic neuroreceptor modelling , 2008, Nuclear medicine communications.

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

[20]  Eric R. Ziegel,et al.  Engineering Statistics , 2004, Technometrics.

[21]  Sung-Cheng Huang,et al.  18F-fluorothymidine kinetics of malignant brain tumors , 2007, European Journal of Nuclear Medicine and Molecular Imaging.

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

[23]  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.

[24]  D. Newport,et al.  A Single Scatter Simulation Technique for Scatter Correction in 3D PET , 1996 .

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

[26]  Mark Muzi,et al.  Kinetic analysis of 3'-deoxy-3'-18F-fluorothymidine in patients with gliomas. , 2006, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[27]  T. Cloughesy,et al.  Kinetics of 3′-Deoxy-3′-18F-Fluorothymidine During Treatment Monitoring of Recurrent High-Grade Glioma , 2010, Journal of Nuclear Medicine.

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

[29]  I. Kawrakow,et al.  Influence of photon energy cuts on PET Monte Carlo simulation results. , 2012, Medical physics.

[30]  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.

[31]  C. Tsoumpas,et al.  Evaluation of the single scatter simulation algorithm implemented in the STIR library , 2004, IEEE Symposium Conference Record Nuclear Science 2004..

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

[33]  Marvin Bergsneider,et al.  Imaging proliferation in brain tumors with 18F-FLT PET: comparison with 18F-FDG. , 2005, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

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