Umeå universitets logga

Positron emission tomography (PET) is increasingly used for delineation of tumour tissue in, for example, radiotherapy treatment planning. The most common method used is to outline volumes with a certain per cent uptake over background in a static image. However, PET data can also be collected dynamically and analysed by kinetic models, which potentially represent the underlying biology better. In the present study, a three-parameter kinetic model was used for voxel-wise evaluation of (11)C-acetate data of head/neck tumours. These parameters which represent the tumour blood volume, the uptake rate and the clearance rate of the tissue were derived for each voxel using a linear regression method and used for segmentation of active tumour tissue. This feasibility study shows that it is possible to segment images based on derived model parameters. There is, however, room for improvements concerning the PET data acquisition, noise reduction and the kinetic modelling. In conclusion, this early study indicates a strong potential of the method even though no 'true' tumour volume was available for validation.

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) softwarewas 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′-(18F)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 therewas 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 parameterswas 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.

Compartmental modeling of dynamic PET data enables quantifi- cation of tracer kinetics in vivo, through the calculated model parameters. In this study, we aimed to investigate the effect of early frame sampling and reconstruction method on pharmacokinetic parameters obtained from a 2-tissue model, in terms of bias and uncertainty (SD). Methods: The GATE Monte Carlo software was used to simulate 2 · 15 dynamic 3′-deoxy-3′-18F-fluorothymidine (18F-FLT) brain PET studies, typical in terms of noise level and kinetic parameters. The data were reconstructed by both 3- dimensional (3D) filtered backprojection with reprojection (3DRP) and 3D ordered-subset expectation maximization (OSEM) into 6 dynamic image sets with different early frame durations of 1, 2, 4, 6, 10, and 15 s. Bias and SD were evaluated for fitted parameter estimates, calculated from regions of interest. Results: The 2-tissue-model parameter estimates K1, k2, and fraction of arterial blood in tissue depended on early frame sampling, and a sampling of 6–15 s generally minimized bias and SD. The shortest sampling of 1 s yielded a 25% and 42% larger bias than the other schemes, for 3DRP and OSEM, respectively, and a parameter uncertainty that was 10%–70% higher. The schemes from 4 to 15 s were generally not significantly different in regards to bias and SD. Typically, the reconstruction method 3DRP yielded less framesampling dependence and less uncertain results, compared with OSEM, but was on average more biased. Conclusion: Of the 6 sampling schemes investigated in this study, an early frame duration of 6–15 s generally kept both bias and uncertainty to a minimum, for both 3DRP and OSEM reconstructions. Veryshort frames of 1 s should be avoided because they typically resulted in the largest parameter bias and uncertainty. Furthermore, 3DRP may be p

Purpose: This work describes PETSTEP (PET Simulator of Tracers via Emission Projection): a faster and more accessible alternative to Monte Carlo (MC) simulation generating realistic PET images, for studies assessing image features and segmentation techniques.

Methods: PETSTEP was implemented within Matlab as open source software. It allows generating threedimensional PET images from PET/CT data or synthetic CT and PET maps, with user-drawn lesions and user-set acquisition and reconstruction parameters. PETSTEP was used to reproduce images of the NEMA body phantom acquired on a GE Discovery 690 PET/CT scanner, and simulated with MC for the GE Discovery LS scanner, and to generate realistic Head and Neck scans. Finally the sensitivity (S) and Positive Predictive Value (PPV) of three automatic segmentation methods were compared when applied to the scanner-acquired and PETSTEP-simulated NEMA images.

Results: PETSTEP produced 3D phantom and clinical images within 4 and 6 min respectively on a single core 2.7 GHz computer. PETSTEP images of the NEMA phantom had mean intensities within 2% of the scanner-acquired image for both background and largest insert, and 16% larger background Full Width at Half Maximum. Similar results were obtained when comparing PETSTEP images to MC simulated data. The S and PPV obtained with simulated phantom images were statistically significantly lower than for the original images, but led to the same conclusions with respect to the evaluated segmentation methods.

Conclusions: PETSTEP allows fast simulation of synthetic images reproducing scanner-acquired PET data and shows great promise for the evaluation of PET segmentation methods.