Quantitative validation of Monte Carlo SPECT simulation: application to a Mediso AnyScan GATE simulation.

BACKGROUND: Monte Carlo (MC) simulations are used in nuclear medicine imaging as they provide unparalleled insight into processes that are not directly experimentally measurable, such as scatter and attenuation in an acquisition. Whilst MC is often used to provide a 'ground-truth', this is only the case if the simulation is fully validated against experimental data. This work presents a quantitative validation for a MC simulation of a single-photon emission computed tomography (SPECT) system. METHODS: An MC simulation model of the Mediso AnyScan SCP SPECT system installed at the UK National Physical Laboratory was developed in the GATE (Geant4 Application for Tomographic Emission) toolkit. Components of the detector head and two collimator configurations were modelled according to technical specifications and physical measurements. Experimental detection efficiency measurements were collected for a range of energies, permitting an energy-dependent intrinsic camera efficiency correction function to be determined and applied to the simulation on an event-by-event basis. Experimental data were collected in a range of geometries with [Formula: see text]Tc for comparison to simulation. The procedure was then repeated with [Formula: see text]Lu to determine how the validation extended to another isotope and set of collimators. RESULTS: The simulation's spatial resolution, sensitivity, energy spectra and the projection images were compared with experimental measurements. The simulation and experimental uncertainties were determined and propagated to all calculations, permitting the quantitative agreement between simulated and experimental SPECT acquisitions to be determined. Statistical agreement was seen in sinograms and projection images of both [Formula: see text]Tc and [Formula: see text]Lu data. Average simulated and experimental sensitivity ratios of ([Formula: see text]) were seen for emission and scatter windows of [Formula: see text]Tc, and ([Formula: see text]) and ([Formula: see text]) for the 113 and 208 keV emissions of [Formula: see text]Lu, respectively. CONCLUSIONS: MC simulations will always be an approximation of a physical system and the level of agreement should be assessed. A validation method is presented to quantify the level of agreement between a simulation model and a physical SPECT system.

Positional dependence of activity determination in single photon emission computed tomography.

Accurate image quantification requires accurate calibration of the detector and is vital if dosimetry is to be performed in molecular radiotherapy. A dependence on the position of calibration has been observed in single photon emission computed tomography images when attenuation correction (AC) and scatter correction are applied. This work investigates the origin of this dependence in single photon emission computed tomography scans of phantom inserts filled with Lu solution. A 113 ml sphere and inserts representing a mathematical model of a spleen and an anatomical model of a patient spleen were imaged at the centre and edge of elliptical phantoms. For these inserts, the difference in calibration factor between the positions was around 10% for images reconstructed with AC and triple energy window scatter correction. A combination of experimental imaging and Monte Carlo simulation was used to isolate possible causes due to imaging or reconstruction in turn. Inconsistent application of AC between different reconstruction systems was identified as the origin of the positional dependence.

Improving molecular radiotherapy dosimetry using anthropomorphic calibration.

The optimised delivery of Molecular Radiotherapy requires individualised calculation of absorbed dose to both targeted lesions and neighbouring healthy tissue. To achieve this, accurate quantification of the activity distribution in the patient by external detection is vital. METHODS: This work extends specific anatomy-related calibration to true organ shapes. A set of patient-specific 3D printed organ inserts based on a diagnostic CT scan was produced, comprising the liver, spleen and both kidneys. The inserts were used to calculate patient-specific calibration factors for 177Lu. These calibration factors were compared with previously reported calibration factors for corresponding organ models based on the Cristy and Eckerman phantom series and for a comparably sized sphere. Monte Carlo calculations of the patient-specific radiation dose were performed for comparison with current clinical dosimetry methods for these data. RESULTS: Patient-specific calibration factors are shown to be dependent on the volume, shape and position of the organ containing activity with a corresponding impact on the calculation of the dose to the patient. The impact of organ morphology on calculated dose is reduced when the dominant contributor to dose is beta particles. This is due to the small range of beta particles in tissue. Overestimations of recovered activity and hence dose of up to 135% are observed. CONCLUSION: For accurate quantification to be performed calibration factors accounting for organ size, shape and position must be used. Such quantification is vital if accurate, patient-specific dosimetry is to be achieved.