Resolution modeling in projection space using a factorized multi-block detector response function for PET image reconstruction

Evaluation of applying space-variant resolution modeling to attenuation correction in PET

Attenuation correction aims to recover the underestimated tracer uptake and improve the image contrast recovery in positron emission tomography (PET). However, traditional ray-tracing-based projection of attenuation maps is inaccurate as some physical effects are not considered, such as finite crystal size, inter-crystal penetration and inter-crystal scatter. In this study, we evaluated the effects of applying resolution modeling (RM) to attenuation correction by implementing space-variant RM to complement physical effects which are usually omitted in the traditional projection model. We verified this method on a brain PET scanner developed by our group, in both Monte Carlo simulation and real-world data, in comparison with space-invariant Gaussian RM, average-depth-of-interaction, and multi-ray tracing methods. The results indicate that the space-variant RM is superior in terms of artifacts reduction and contrast recovery.

High-throughput, accurate Monte Carlo simulation on CPU hardware for PET applications

Monte Carlo simulations (MCS) represent a fundamental approach to modelling the photon interactions in Positron Emission Tomography (PET). A variety of PET-dedicated MCS tools are available to assist and improve PET imaging applications. Of these, GATE has evolved into one of the most popular software for PET MCS because of its accuracy and flexibility. However, simulations are extremely time-consuming. The use of graphics processing units (GPU) has been proposed as a solution to this, with reported acceleration factors about 400-800. These factors refer to GATE benchmarks performed on a single CPU core. Consequently, CPU-based MCS can also be easily accelerated by one order of magnitude or beyond when exploiting multi-threading on powerful CPUs. Thus, CPU-based implementations become competitive when further optimisations can be achieved. In this context, we have developed a novel, CPU-based software called the PET Physics Simulator (PPS), which combines several efficient methods to significantly boost the performance. PPS flexibly applies GEANT4 cross-sections as a pre-calculated database, thus obtaining results equivalent to GATE. This is demonstrated for an elaborated PET scanner with 3-layer block detectors. All code optimisations yield an acceleration factor of  20 (single core). Multi-threading on a high-end CPU workstation (96 cores) further accelerates the PPS by a factor of 80. This results in a total speed-up factor of 1600, which outperforms comparable GPU-based MCS by a factor of  2. Optionally, the proposed method of coincidence multiplexing can further enhance the throughput by an additonal factor of 15. The combination of all optimisations corresponds to an acceleration factor of  24000. In this way, the PPS can simulate complex PET detector systems with an effective throughput of photon pairs in less than 10 milliseconds.

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.

SimPET-An open online platform for the Monte Carlo simulation of realistic brain PET data. Validation for 18 F-FDG scans

Purpose

SimPET (www.sim‐pet.org) is a free cloud‐based platform for the generation of realistic brain positron emission tomography (PET) data. In this work, we introduce the key features of the platform. In addition, we validate the platform by performing a comparison between simulated healthy brain FDG‐PET images and real healthy subject data for three commercial scanners (GE Advance NXi, GE Discovery ST, and Siemens Biograph mCT).

Methods

The platform provides a graphical user interface to a set of automatic scripts taking care of the code execution for the phantom generation, simulation (SimSET), and tomographic image reconstruction (STIR). We characterize the performance using activity and attenuation maps derived from PET/CT and MRI data of 25 healthy subjects acquired with a GE Discovery ST. We then use the created maps to generate synthetic data for the GE Discovery ST, the GE Advance NXi, and the Siemens Biograph mCT. The validation was carried out by evaluating Bland‐Altman differences between real and simulated images for each scanner. In addition, SPM voxel‐wise comparison was performed to highlight regional differences. Examples for amyloid PET and for the generation of ground‐truth pathological patients are included.

Results

The platform can be efficiently used for generating realistic simulated FDG‐PET images in a reasonable amount of time. The validation showed small differences between SimPET and acquired FDG‐PET images, with errors below 10% for 98.09% (GE Discovery ST), 95.09% (GE Advance NXi), and 91.35% (Siemens Biograph mCT) of the voxels. Nevertheless, our SPM analysis showed significant regional differences between the simulated images and real healthy patients, and thus, the use of the platform for converting control subject databases between different scanners requires further investigation.

Conclusions

The presented platform can potentially allow scientists in clinical and research settings to perform MC simulation experiments without the need for high‐end hardware or advanced computing knowledge and in a reasonable amount of time.

Spatially-variant image-based modeling of PSF deformations with application to a limited angle geometry from a dual-panel breast-PET imager

Dual-panel PET system configuration can lead to spatially variable point-spread functions (PSF) of considerable deformations due to depth-of-interaction effects and limited angular coverage. If not modelled properly, these effects result in decreased and inconsistent recovery of lesion activity across the field-of-view (FOV), as well as mispositioning of lesions in the reconstructed image caused by strong PSF asymmetries. We implemented and evaluated models of such PSF deformations with spatially-variant image-based resolution modeling (IRM) within reconstruction (varRM) using the Direct Image REConstruction for Time-of-flight (DIRECT) method and within post-reconstruction deconvolution methods. In addition, DIRECT reconstruction was performed with a spatially-invariant IRM (invRM) and without resolution modeling (noRM) for comparison. The methods were evaluated using simulated data for a realistic breast model with a set of 5 mm lesions located throughout the FOV of a dual-panel Breast-PET scanner. We simulated high-count data to focus on the ability of each method to correctly recover the PSF deformations, and a clinically realistic count level to assess the impact of low count data on the quantitative performance of the evaluated techniques. Performance of the methods evaluated herein was assessed by comparing lesion activity recovery (%BIAS), consistency (%SD) across the FOV, overall error (%RMSE), and recovery of each lesion location. As expected, all techniques using IRM provide considerable improvement over the noRM reconstruction. For the high-count cases, the overall quantitative performance of all IRM techniques, whether within reconstruction or within post-reconstruction, is similar if the lesion location misplacements are ignored. However, invRM provides less consistent performance on activity across lesions and is not able to recover accurate lesion locations. For a clinically realistic count level, varRM reconstruction consistently outperforms all compared approaches, while the post-reconstruction IRM approaches exhibit higher %SD and %RMSE values due to being more affected by the data noise than the within-reconstruction IRM approaches.