Monte Carlo simulation of digital photon counting PET

Design and proof of concept of a double-panel TOF-PET system

OBJECTIVE

Positron Emission Tomography (PET) is a well-known imaging technology for the diagnosis, treatment, and monitoring of several diseases. Most PET scanners use a Ring-Shaped Detector Configuration (RSDC), which helps obtain homogeneous image quality but are restricted to an invariable Field-of-View (FOV), scarce spatial resolution, and low sensitivity. Alternatively, few PET systems use Open Detector Configurations (ODC) to permit an accessible FOV adaptable to different target sizes, thus optimizing sensitivity. Yet, to compensate the lack of angular coverage in ODC-PET, developing a detector with high-timing performance is mandatory to enable Time-of-Flight (TOF) techniques during reconstruction. The main goal of this work is to provide a proof of concept PET scanner appropriate for constructing the new generation of ODC-PET suitable for biopsy guidance and clinical intervention during acquisition. The designed detector has to be compact and robust, and its requirements in terms of performance are spatial and time resolutions < 2 mm and < 200 ps, respectively.

METHODS

The present work includes a simulation study of an ODC-PET based on 2-panels with variable distance. The image quality (IQ) and Derenzo phantoms have been simulated and evaluated. The phantom simulations have also been performed using a ring-shaped PET for comparison purposes of the ODC approach with conventional systems. Then, an experimental evaluation of a prototype detector that has been designed following the simulation results is presented. This study focused on tuning the ASIC parameters and evaluating the scintillator surface treatment (ESR and TiO2), and configuration that yields the best Coincidence Time Resolution (CTR). Moreover, the scalability of the prototype to a module of 64 × 64mm2 and its preliminary evaluation regarding pixel identification are provided.

RESULTS

The simulation results reported sensitivity (%) values at the center of the FOV of 1.96, 1.63, and 1.18 for panel distances of 200, 250, and 300 mm, respectively. The IQ reconstructed image reported good uniformity (87%) and optimal CRC values, and the Derenzo phantom reconstruction suggests a system resolution of 1.6-2 mm. The experimental results demonstrate that using TiO2 coating yielded better detector performance than ESR. Acquired data was filtered by applying an energy window of ± 30% at the photopeak level. After filtering, best CTR of 230 ± 2 ps was achieved for an 8 × 8 LYSO pixel block with 2 × 2 × 12mm3 each. The detector performance remained constant after scaling-up the prototype to a module of 64 × 64mm2, and the flood map demonstrates the module's capabilities to distinguish the small pixels; thus, a spatial resolution < 2 mm (pixel size) is achieved.

CONCLUSIONS

The simulated results of this biplanar scanner show high performance in terms of image quality and sensitivity. These results are comparable to state-of-the-art PET technology and, demonstrate that including TOF information minimizes the image artifacts due to the lack of angular projections. The experimental results concluded that using TiO2 coating provide the best performance. The results suggest that this scanner may be suitable for organ study, breast, prostate, or cardiac applications, with good uniformity and CRC.

PET digitization chain for Monte Carlo simulation in GATE

Objective. We introduce a versatile methodology for the accurate modelling of PET imaging systems via Monte Carlo simulations, using the Geant4 application for tomographic emission (GATE) platform. Accurate Monte Carlo modelling involves the incorporation of a complete analytical signal processing chain, called the digitizer in GATE, to emulate the different count rates encountered in actual positron emission tomography (PET) systems.Approach. The proposed approach consists of two steps: (1) modelling the digitizer to replicate the detection chain of real systems, covering all available parameters, whether publicly accessible or supplied by manufacturers; (2) estimating the remaining parameters, i.e. background noise level, detection efficiency, and pile-up, using optimisation techniques based on experimental single and prompt event rates. We show that this two-step optimisation reproduces the other experimental count rates (true, scatter, and random), without the need for additional adjustments. This method has been applied and validated with experimental data derived from the NEMA count losses test for three state-of-the-art SiPM-based time-of-flight (TOF)-PET systems: Philips Vereos, Siemens Biograph Vision 600 and GE Discovery MI 4-ring.Main results. The results show good agreement between experiments and simulations for the three PET systems, with absolute relative discrepancies below 3%, 6%, 6%, 7% and 12% for prompt, random, true, scatter and noise equivalent count rates, respectively, within the 0-10 kBq·ml-1activity concentration range typically observed in whole-body18F scans.Significance. Overall, the proposed digitizer optimisation method was shown to be effective in reproducing count rates and NECR for three of the latest generation SiPM-based TOF-PET imaging systems. The proposed methodology could be applied to other PET scanners.

Validation of a discovery MI 4-ring model according to the NEMA NU 2-2018 standards: from Monte Carlo simulations to clinical-like reconstructions

BACKGROUND

We propose a comprehensive evaluation of a Discovery MI 4-ring (DMI) model, using a Monte Carlo simulator (GATE) and a clinical reconstruction software package (PET toolbox). The following performance characteristics were compared with actual measurements according to NEMA NU 2-2018 guidelines: system sensitivity, count losses and scatter fraction (SF), coincidence time resolution (CTR), spatial resolution (SR), and image quality (IQ). For SR and IQ tests, reconstruction of time-of-flight (TOF) simulated data was performed using the manufacturer's reconstruction software.

RESULTS

Simulated prompt, random, true, scatter and noise equivalent count rates closely matched the experimental rates with maximum relative differences of 1.6%, 5.3%, 7.8%, 6.6%, and 16.5%, respectively, in a clinical range of less than 10 kBq/mL. A 3.6% maximum relative difference was found between experimental and simulated sensitivities. The simulated spatial resolution was better than the experimental one. Simulated image quality metrics were relatively close to the experimental results.

CONCLUSIONS

The current model is able to reproduce the behaviour of the DMI count rates in the clinical range and generate clinical-like images with a reasonable match in terms of contrast and noise.

TOPAS-imaging: extensions to the TOPAS simulation toolkit for medical imaging systems

Objective. The TOol for PArticle Simulation (TOPAS) is a Geant4-based Monte Carlo software application that has been used for both research and clinical studies in medical physics. So far, most users of TOPAS have focused on radiotherapy-related studies, such as modeling radiation therapy delivery systems or patient dose calculation. Here, we present the first set of TOPAS extensions to make it easier for TOPAS users to model medical imaging systems.Approach. We used the extension system of TOPAS to implement pre-built, user-configurable geometry components such as detectors (e.g. flat-panel and multi-planar detectors) for various imaging modalities and pre-built, user-configurable scorers for medical imaging systems (e.g. digitizer chain).Main results. We developed a flexible set of extensions that can be adapted to solve research questions for a variety of imaging modalities. We then utilized these extensions to model specific examples of cone-beam CT (CBCT), positron emission tomography (PET), and prompt gamma (PG) systems. The first of these new geometry components, the FlatImager, was used to model example CBCT and PG systems. Detected signals were accumulated in each detector pixel to obtain the intensity of x-rays penetrating objects or prompt gammas from proton-nuclear interaction. The second of these new geometry components, the RingImager, was used to model an example PET system. Positron-electron annihilation signals were recorded in crystals of the RingImager and coincidences were detected. The simulated data were processed using corresponding post-processing algorithms for each modality and obtained results in good agreement with the expected true signals or experimental measurement.Significance. The newly developed extension is a first step to making it easier for TOPAS users to build and simulate medical imaging systems. Together with existing TOPAS tools, this extension can help integrate medical imaging systems with radiotherapy simulations for image-guided radiotherapy.

Design and performance simulation studies of a breast PET insert integrable into a clinical whole-body PET/MRI scanner

Objective. Three different breast positron emission tomography (PET) insert geometries are proposed for integration into an existing magnetic resonance imaging (MRI) breast coil (Breast Biopsy Coil, NORAS MRI products) to be used inside a whole-body PET/MRI scanner (Biograph mMR, Siemens Healthineers) to enhance the sensitivity and spatial resolution of imaging inside the breast.Approach. Monte Carlo simulations were performed to predict and compare the performance characteristics of the three geometries in terms of the sensitivity, spatial resolution, scatter fraction, and noise equivalent count rate (NECR). In addition, the background single count rate due to organ uptake in a clinical scan scenario was predicted using a realistic anthropomorphic phantom.Main results. In the center of the field of view (cFOV), absolute sensitivities of 3.1%, 2.7%, and 2.2% were found for Geometry A (detectors arranged in two cylinders), Geometry B (detectors arranged in two partial cylinders), and Geometry C (detectors arranged in two half cylinders combined with two plates), respectively. The full width at half maximum spatial resolution was determined to be 1.7 mm (Geometry A), 1.8 mm (Geometry B) and 2.0 mm (Geometry C) at 5 mm from the cFOV. Designs with multiple scintillation-crystal layers capable of determining the depth of interaction (DOI) strongly improved the spatial resolution at larger distances from the transaxial cFOV. The system scatter fractions were 33.1% (Geometries A and B) and 32.3% (Geometry C). The peak NECRs occurred at source activities of 300 MBq (Geometry A), 310 MBq (Geometry B) and 340 MBq (Geometry C). The background single-event count rates were 17.1 × 10⁶cps (Geometry A), 15.3 × 10⁶cps (Geometry B) and 14.8 × 10⁶cps (Geometry C). Geometry A in the three-layer DOI variant exhibited the best PET performance characteristics but could be challenging to manufacture. Geometry C had the lowest impact on the spatial resolution and the lowest sensitivity among the investigated geometries.Significance. Geometry B in the two-layer DOI variant represented an effective compromise between the PET performance and manufacturing difficulty and was found to be a promising candidate for the future breast PET insert.

Autonomous digitizer calibration of a Monte Carlo detector model through evolutionary simulation

Simulating the response of a radiation detector is a modelling challenge due to the stochastic nature of radiation, often complex geometries, and multi-stage signal processing. While sophisticated tools for Monte Carlo simulation have been developed for radiation transport, emulating signal processing and data loss must be accomplished using a simplified model of the electronics called the digitizer. Due to a large number of free parameters, calibrating a digitizer quickly becomes an optimisation problem. To address this, we propose a novel technique by which evolutionary algorithms calibrate a digitizer autonomously. We demonstrate this by calibrating six free parameters in a digitizer model for the ADAC Forte. The accuracy of solutions is quantified via a cost function measuring the absolute percent difference between simulated and experimental coincidence count rates across a robust characterisation data set, including three detector configurations and a range of source activities. Ultimately, this calibration produces a count rate response with 5.8% mean difference to the experiment, improving from 18.3% difference when manually calibrated. Using evolutionary algorithms for model calibration is a notable advancement because this method is novel, autonomous, fault-tolerant, and achieved through a direct comparison of simulation to reality. The software used in this work has been made freely available through a GitHub repository.

Towards efficient Monte Carlo N-Particle simulation of a positron emission tomography (PET) via source volume definition

Monte Carlo N-Particle (MCNP) simulation has been extensively proven in nuclear medicine imaging systems, most notably in designing and optimizing new medical imaging tools. It enables more complicated geometries and the simulation of particles passing through and interacting with materials. However, a relatively long simulation time is a drawback of Monte Carlo simulation, mainly when complex geometry exists. The current study presents an alternative variance reduction technique for a modeled positron emission tomography (PET) camera by reducing the height of the source volume definition while maintaining the geometry of the simulated model. The National Electrical Manufacturers Association (NEMA) of the International Electrotechnical Commission (IEC) PET's phantom was used with a 1 cm diameter and 7 cm height of line source placed in the middle. The first geometry was fully filled the line source with 0.50 mCi radioactivity. In contrast, the second geometry decreased the source definition to 2.4 cm in height, covering 1 cm above and below the sub-block detector level. The source volume definition approach led to a 71% reduction in the total photons to be simulated. Results showed that the proposed variance reduction strategy could produce spatial resolution as precise as fully filled geometry and sped up the simulation time by approximately 65%. Hence, this strategy can be utilized for further PET optimizing simulation studies.

Validation of a computational chain from PET Monte Carlo simulations to reconstructed images

The study aimed to create a pipeline from Monte Carlo simulated projections of a Gate PET system to reconstructed images. The PET system was modelled after the GE Discovery MI (DMI) PET/CT, and the simulated projections were reconstructed with the stand-alone reconstruction software CASToR. Attenuation correction, normalisation calibration, random estimation, and scatter estimation for the simulations were computed with in-house programs. The pipeline was compared in both projection and image space with data acquired on a clinical DMI and reconstructed with GE's off-line PET reconstruction software (PET Toolbox) and CASToR. The simulated and measured data were compared for the number of prompt coincidences, scatter fraction, contrast recovery coefficient (CRC), signal-to-noise ratio (SNR), background variability, residual lung error, and image profiles. A slight discrepancy was noted in the projection space, but good agreements were generally achieved in image space between simulated and measured data. The CRC was found to be 81 % for Gate – CASToR, 84 % for GE – CASToR, and 84 % for GE - PET Toolbox for the largest sphere of the NEMA image quality (IQ) phantom, and the SNR was found to be 98 for Gate – CASToR, 91 for GE – CASToR, and 93 for GE – PET Toolbox. Profiles drawn over the spheres for the NEMA IQ phantom and the Data Spectrum (DS) phantom show a good match between measurement and simulation. The results indicate feasibility to utilise the pipeline as a tool for off-line simulation-based studies. A complete pipeline introduces possibilities to study the impact of single parameters in the whole chain from simulation to reconstructed images.

Shining Damaged Hearts: Immunotherapy-Related Cardiotoxicity in the Spotlight of Nuclear Cardiology

The emerging use of immunotherapies in cancer treatment increases the risk of immunotherapy-related cardiotoxicity. In contrast to conventional chemotherapy, these novel therapies have expanded the forms and presentations of cardiovascular damage to a broad spectrum from asymptomatic changes to fulminant short- and long-term complications in terms of cardiomyopathy, arrythmia, and vascular disease. In cancer patients and, particularly, cancer patients undergoing (immune-)therapy, cardio-oncological monitoring is a complex interplay between pretherapeutic risk assessment, identification of impending cardiotoxicity, and post-therapeutic surveillance. For these purposes, the cardio-oncologist can revert to a broad spectrum of nuclear cardiological diagnostic workup. The most promising commonly used nuclear medicine imaging techniques in relation to immunotherapy will be discussed in this review article with a special focus on the continuous development of highly specific molecular markers and steadily improving methods of image generation. The review closes with an outlook on possible new developments of molecular imaging and advanced image evaluation techniques in this exciting and increasingly growing field of immunotherapy-related cardiotoxicity.

Monte Carlo evaluation of hypothetical long axial field-of-view PET scanner using GE discovery MI PET front-end architecture

PURPOSE

The development of total-body PET scanners is of growing interest in the PET community. Investigation into the imaging properties of a hypothetical extended axial field-of-view (AFOV) GE Healthcare SiPM-based Discovery MI (DMI) system architecture has not yet been performed. In this work, we assessed its potential as a whole-body scanner using Monte Carlo simulations. The aim of this work was to (1) develop and validate a Monte Carlo model of a 4-ring scanner and (2) extend its AFOV up to 2 m to evaluate performance gain through NEMA-based evaluation.

METHODS

The DMI 4-ring geometry and its pulse digitization scheme were modeled within the GATE Monte Carlo platform using published literature. The GATE scanner model was validated by comparing results against published NEMA performance measurements. Following the validation of the 4-ring model, the model was extended to simulate 8, 20, 30, and 40-ring systems. Spatial resolution, sensitivity, NECR, and scatter fraction were characterized with modified NEMA NU-2 2018 standards; however, the image quality measurements were not acquired due to computational limitations. Spatial resolutions were simulated for all scanner ring configurations using point sources to examine the effects of parallax errors. NEMA count rates were estimated using a standard 70 cm scatter phantom and an extended version of scatter phantom of length 200 cm with (1-800) MBq of ¹⁸ F for all scanners. Sensitivity was evaluated using NEMA methods with a 70 cm standard and a 200 cm long line source.

RESULTS

The average FWHM of the radial/tangential/axial spatial resolution reconstructed with filtered back-projection at 1 and 10 cm from the scanner center were 3.94/4.10/4.41 mm and 5.29/4.89/5.90 mm for the 4-ring scanner. Sensitivity was determined to be 14.86 cps/kBq at the center of the FOV for the 4-ring scanner using a 70 cm line source. Sensitivity enhancement up to 21-fold and 60-fold were observed for 1 m and 2 m AFOV scanners compared to 4-ring scanner using a 200 cm long line source. Spatial resolution simulations in a 2 m AFOV scanner suggest a maximum degradation of ∼23.8% in the axial resolution compared to the 4-ring scanner. However, the transverse resolution was found to be relatively constant when increasing the axial acceptance angle up to ±70°. The peak NECR was 212.92 kcps at 22.70 kBq/mL with a scatter fraction of 38.9% for a 4-ring scanner with a 70 cm scatter phantom. Comparison of peak NECR using the 200 cm long scatter phantom relative to the 4-ring scanner resulted in a NECR gain of 15 for the 20-ring and 28 for the 40-ring geometry. Spatial resolution, sensitivity, and scatter fraction showed an agreement within ∼7% compared with published measured values.

CONCLUSIONS

The 4-ring DMI scanner simulation was successfully validated against published NEMA measurements. Sensitivity and NECR performance of extended 1 and 2 meters AFOV scanners based upon the DMI architecture were subsequently simulated. Increases in sensitivity and count-rate performance are consistent with prior simulation studies utilizing extensions of the Siemens mCT architecture and published NEMA measurements with the uEXPLORER system. This article is protected by copyright. All rights reserved.

Evaluation of simulators for x-ray speckle-based phase contrast imaging

X-ray phase contrast imaging (PCI) denotes a group of highly sensitive imaging techniques that permits imaging at scales ranging from nanoscopic to the medical. Recently introduced, speckle-based imaging has seen a rapid development because of its experimental simplicity and its capability to retrieve the refraction, the scattering and the absorption of a sample using a conventional x-ray set-up. Precise simulation would permit to optimise the imaging setups for different applications, but until now works on simulation of x-ray speckle-based PCI have been very few. In this work we evaluate different simulation codes, based on Monte-Carlo, analytical ray-tracing and wave-optics Fresnel propagation. The simulation results are compared to both synchrotron and conventional imaging experiments to permits their validation. We obtain a strong similarity between simulated and experimental data. We discuss the validity and applicability of each approach.

Yttrium-90 quantitative phantom study using digital photon counting PET

BACKGROUND

PET imaging of 90Y-microsphere distribution following radioembolisation is challenging due to the count-starved statistics from the low branching ratio of e+/e- pair production during 90Y decay. PET systems using silicon photo-multipliers have shown better 90Y image quality compared to conventional photo-multiplier tubes. The main goal of the present study was to evaluate reconstruction parameters for different phantom configurations and varying listmode acquisition lengths to improve quantitative accuracy in 90Y dosimetry, using digital photon counting PET/CT.

METHODS

Quantitative PET and dosimetry accuracy were evaluated using two uniform cylindrical phantoms specific for PET calibration validation. A third body phantom with a 9:1 hot sphere-to-background ratio was scanned at different activity concentrations of 90Y. Reconstructions were performed using OSEM algorithm with varying parameters. Time-of-flight and point-spread function modellings were included in all reconstructions. Absorbed dose calculations were carried out using voxel S-values convolution and were compared to reference Monte Carlo simulations. Dose-volume histograms and root-mean-square deviations were used to evaluate reconstruction parameter sets. Using listmode data, phantom and patient datasets were rebinned into various lengths of time to assess the influence of count statistics on the calculation of absorbed dose. Comparisons between the local energy deposition method and the absorbed dose calculations were performed.

RESULTS

Using a 2-mm full width at half maximum post-reconstruction Gaussian filter, the dosimetric accuracy was found to be similar to that found with no filter applied but also reduced noise. Larger filter sizes should not be used. An acquisition length of more than 10 min/bed reduces image noise but has no significant impact in the quantification of phantom or patient data for the digital photon counting PET. 3 iterations with 10 subsets were found suitable for large spheres whereas 1 iteration with 30 subsets could improve dosimetry for smaller spheres.

CONCLUSION

The best choice of the combination of iterations and subsets depends on the size of the spheres. However, one should be careful on this choice, depending on the imaging conditions and setup. This study can be useful in this choice for future studies for more accurate 90Y post-dosimetry using a digital photon counting PET/CT.

Advanced Monte Carlo simulations of emission tomography imaging systems with GATE

Built on top of the Geant4 toolkit, GATE is collaboratively developed for more than 15 years to design Monte Carlo simulations of nuclear-based imaging systems. It is, in particular, used by researchers and industrials to design, optimize, understand and create innovative emission tomography systems. In this paper, we reviewed the recent developments that have been proposed to simulate modern detectors and provide a comprehensive report on imaging systems that have been simulated and evaluated in GATE. Additionally, some methodological developments that are not specific for imaging but that can improve detector modeling and provide computation time gains, such as Variance Reduction Techniques and Artificial Intelligence integration, are described and discussed.

Fully digital PET is unaffected by any deterioration in TOF resolution and TOF image quality in the wide range of routine PET count rates

PURPOSE

Digital PET involving silicon photomultipliers (SiPM) provides an enhanced time-of-flight (TOF) resolution as compared with photomultiplier (PMT)-based PET, but also a better prevention of the count-related rises in dead time and pile-up effects mainly due to smaller trigger domains (i.e., the detection surfaces associated with each trigger circuit). This study aimed to determine whether this latter property could help prevent against deteriorations in TOF resolution and TOF image quality in the wide range of PET count rates documented in clinical routine.

METHODS

Variations, according to count rates, in timing resolution and in TOF-related enhancement of the quality of phantom images were compared between the first fully digital PET (Vereos) and a PMT-based PET (Ingenuity). Single-count rate values were additionally extracted from the list-mode data of routine analog- and digital-PET exams at each 500-ms interval, in order to determine the ranges of routine PET count rates.

RESULTS

Routine PET count rates were lower for the Vereos than for the Ingenuity. For Ingenuity, the upper limits were estimated at approximately 21.7 and 33.2 Mcps after injection of respectively 3 and 5 MBq.kg^-1 of current ¹⁸F-labeled tracers. At 5.8 Mcps, corresponding to the lower limit of the routine count rates documented with the Ingenuity, timing resolutions provided by the scatter phantom were 326 and 621 ps for Vereos and Ingenuity, respectively. At higher count rates, timing resolution was remarkably stable for Vereos but exhibited a progressive deterioration for Ingenuity, respectively reaching 732 and 847 ps at the upper limits of 21.7 and 33.2 Mcps. The averaged TOF-related gain in signal/noise ratio was stable at approximately 2 for Vereos but decreased from 1.36 at 5.8 Mcps to 1.14 and 1.00 at respectively 21.7 and 33.2 Mcps for Ingenuity.

CONCLUSION

Contrary to the Ingenuity PMT-based PET, the Vereos fully digital PET is unaffected by any deterioration in TOF resolution and consequently, in the quality of TOF images, in the wide range of routine PET count rates. This advantage is even more striking with higher count-rates for which the preferential use of digital PET should be further recommended (i.e., dynamic PET recording, higher injected activities).