Publicly available framework for simulating and experimentally validating clinical PET systems.

Abstract

Monte Carlo (MC) simulations are a powerful tool to model medical imaging systems. However, before simulations can be considered the ground truth, they have to be validated with experiments. To provide a pipeline that models a clinical positron emission tomography (PET)/CT system using MC simulations after extensively validating the results against experimental measurements. A clinical four-ring PET imaging system was modeled using Geant4 application for tomographic emission (v. 9.0). To validate the simulations, PET images were acquired of a cylindrical phantom, point source, and image quality phantom with the modeled system and the simulations of the experimental procedures. For the purpose of validating the quantification capabilities and image quality provided by the simulation pipeline, the simulations were compared against the measurements in terms of their count rates and sensitivity as well as their image uniformity, resolution, recovery coefficients (RCs), coefficients of variation, contrast, and background variability. When compared to the measured data, the number of true detections in the MC simulations was within 5%. The scatter fraction was found to be 30.0%±2.2% and 28.8%±1.7% in the measured and simulated scans, respectively. Analyzing the measured and simulated sinograms, the sensitivities were found to be 8.2 and 7.8cps/kBq, respectively. The fraction of random coincidences were 19% in the measured data and 25% in the simulation. When calculating the image uniformity within the axial slices, the measured image exhibited a uniformity of 0.015±0.005, whereas the simulated image had a uniformity of 0.029±0.011. In the axial direction, the uniformity was measured to be 0.024±0.006 and 0.040±0.015 for the measured and simulated data, respectively. Comparing the image resolution, an average percentage difference of 2.9% was found between the measurements and simulations. The RCs calculated in both the measured and simulated images were found to be within the EARL ranges, except for that of the simulation of the smallest sphere. The coefficients of variation for the measured and simulated images were found to be 12% and 13%, respectively. Lastly, the background variability was consistent between the measurements and simulations, whereas the average percentage difference in the sphere contrasts was found to be 8.8%. The clinical PET/CT system was modeled and validated to provide a simulation pipeline for the community. The pipeline and the validation procedures have been made available (https://github.com/teaghan/PET_MonteCarlo).