ABSTRACT
Monte-Carlo simulations of SPECT images are notoriously slow to converge due to the large ratio between the number of photons emitted and detected in the collimator. This work proposes a method to accelerate the simulations based on fixed forced detection (FFD) combined with an analytical response of the detector. FFD is based on a Monte-Carlo simulation but forces the detection of a photon in each detector pixel weighted by the probability of emission (or scattering) and transmission to this pixel. The method was evaluated with numerical phantoms and on patient images. We obtained differences with analog Monte Carlo lower than the statistical uncertainty. The overall computing time gain can reach up to five orders of magnitude. Source code and examples are available in the Gate V8.0 release.
Subject(s)
Computer Simulation , Monte Carlo Method , Phantoms, Imaging , Photons , Tomography, Emission-Computed, Single-Photon/methods , HumansABSTRACT
Simulations of planar whole body acquisitions in therapeutic procedures are often extensively time-consuming and therefore rarely used. However, optimising tools and variance reduction techniques can be employed to overcome this problem. In this paper, a variety of features available in GATE are explored and their capabilities to reduce simulation time are evaluated. For this purpose, the male XCAT phantom was used as a virtual patient with 177Lu-DOTATATE pharmacokinetic for whole body planar acquisition simulations in a Siemens Symbia T2 model. Activity distribution was divided into 8 compartments that were simulated separately. GATE optimization techniques included reducing the amount of time spent in both voxel and detector tracking. Some acceleration techniques led to a decrease of CPU-time by a factor of 167, while image statistics were kept constant. In that context, the simulation of therapeutic procedure imaging would still require 46days on a single CPU, but this could be reduced to hours on a dedicated cluster.
