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

UMC-PET: a fast and flexible Monte Carlo PET simulator

Objective.The GPU-based Ultra-fast Monte Carlo positron emission tomography simulator (UMC-PET) incorporates the physics of the emission, transport and detection of radiation in PET scanners. It includes positron range, non-colinearity, scatter and attenuation, as well as detector response. The objective of this work is to present and validate UMC-PET as a a multi-purpose, accurate, fast and flexible PET simulator.Approach.We compared UMC-PET against PeneloPET, a well-validated MC PET simulator, both in preclinical and clinical scenarios. Different phantoms for scatter fraction (SF) assessment following NEMA protocols were simulated in a 6R-SuperArgus and a Biograph mMR scanner, comparing energy histograms, NEMA SF, and sensitivity for different energy windows. A comparison with real data reported in the literature on the Biograph scanner is also shown.Main results.NEMA SF and sensitivity estimated by UMC-PET where within few percent of PeneloPET predictions. The discrepancies can be attributed to small differences in the physics modeling. Running in a 11 GB GeForce RTX 2080 Ti GPU, UMC-PET is ∼1500 to ∼2000 times faster than PeneloPET executing in a single core Intel(R) Xeon(R) CPU W-2155 @ 3.30 GHz.Significance.UMC-PET employs a voxelized scheme for the scanner, patient adjacent objects (such as shieldings or the patient bed), and the activity distribution. This makes UMC-PET extremely flexible. Its high simulation speed allows applications such as MC scatter correction, faster SRM estimation for complex scanners, or even MC iterative image reconstruction.

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.

Synthetic PET via Domain Translation of 3-D MRI

Historically, patient datasets have been used to develop and validate various reconstruction algorithms for PET/MRI and PET/CT. To enable such algorithm development, without the need for acquiring hundreds of patient exams, in this article we demonstrate a deep learning technique to generate synthetic but realistic whole-body PET sinograms from abundantly available whole-body MRI. Specifically, we use a dataset of 56 18F-FDG-PET/MRI exams to train a 3-D residual UNet to predict physiologic PET uptake from whole-body T1-weighted MRI. In training, we implemented a balanced loss function to generate realistic uptake across a large dynamic range and computed losses along tomographic lines of response to mimic the PET acquisition. The predicted PET images are forward projected to produce synthetic PET (sPET) time-of-flight (ToF) sinograms that can be used with vendor-provided PET reconstruction algorithms, including using CT-based attenuation correction (CTAC) and MR-based attenuation correction (MRAC). The resulting synthetic data recapitulates physiologic 18F-FDG uptake, e.g., high uptake localized to the brain and bladder, as well as uptake in liver, kidneys, heart, and muscle. To simulate abnormalities with high uptake, we also insert synthetic lesions. We demonstrate that this sPET data can be used interchangeably with real PET data for the PET quantification task of comparing CTAC and MRAC methods, achieving ≤ 7.6% error in mean-SUV compared to using real data. These results together show that the proposed sPET data pipeline can be reasonably used for development, evaluation, and validation of PET/MRI reconstruction methods.

β-amyloid PET harmonisation across longitudinal studies: Application to AIBL, ADNI and OASIS3

INTRODUCTION

The Centiloid scale was developed to harmonise the quantification of β-amyloid (Aβ) PET images across tracers, scanners, and processing pipelines. However, several groups have reported differences across tracers and scanners even after centiloid conversion. In this study, we aim to evaluate the impact of different pre and post-processing harmonisation steps on the robustness of longitudinal Centiloid data across three large international cohort studies.

METHODS

All Aβ PET data in AIBL (N = 3315), ADNI (N = 3442) and OASIS3 (N = 1398) were quantified using the MRI-based Centiloid standard SPM pipeline and the PET-only pipeline CapAIBL. SUVR were converted into Centiloids using each tracer's respective transform. Global Aβ burden from pre-defined target cortical regions in Centiloid units were quantified for both raw PET scans and PET scans smoothed to a uniform 8 mm full width half maximum (FWHM) effective smoothness. For Florbetapir, we assessed the performance of using both the standard Whole Cerebellum (WCb) and a composite white matter (WM)+WCb reference region. Additionally, our recently proposed quantification based on Non-negative Matrix Factorisation (NMF) was applied to all spatially and SUVR normalised images. Correlation with clinical severity measured by the Mini-Mental State Examination (MMSE) and effect size, as well as tracer agreement in 11C-PiB-18F-Florbetapir pairs and longitudinal consistency were evaluated.

RESULTS

The smoothing to a uniform resolution partially reduced longitudinal variability, but did not improve inter-tracer agreement, effect size or correlation with MMSE. Using a Composite reference region for 18F-Florbetapir improved inter-tracer agreement, effect size, correlation with MMSE, and longitudinal consistency. The best results were however obtained when using the NMF method which outperformed all other quantification approaches in all metrics used.

CONCLUSIONS

FWHM smoothing has limited impact on longitudinal consistency or outliers. A Composite reference region including subcortical WM should be used for computing both cross-sectional and longitudinal Florbetapir Centiloid. NMF improves Centiloid quantification on all metrics examined.

Impact of spill-in counts from off-target regions on [18F]Flortaucipir PET quantification

BACKGROUND

[¹⁸F]Flortaucipir (FTP) PET quantification is usually hindered by spill-in counts from off-target binding (OFF) regions. The present work aims to provide an in-depth analysis of the impact of OFF in FTP PET quantification, as well as to identify optimal partial volume correction (PVC) strategies to minimize this problem.

METHODS

309 amyloid-beta (Aβ) negative cognitively normal subjects were included in the study. Additionally, 510 realistic FTP images with different levels of OFF were generated using Monte Carlo simulation (MC). Images were corrected for PVC using both a simple two-compartment and a multi-region method including OFF regions. FTP standardized uptake value ratio (SUVR) was quantified in Braak Areas (BA), the hippocampus (which was not included in Braak I/II) and different OFF regions (caudate, putamen, pallidum, thalamus, choroid plexus (ChPlex), cerebellar white matter (cerebWM), hemispheric white matter (hemisWM) and cerebrospinal fluid (CSF)) using the lower portion of the cerebellum as a reference region. The correlations between OFF and cortical SUVRs were studied both in real and in simulated PET images, with and without PVC.

RESULTS

In-vivo, we found correlations between all OFF and target regions, especially strong for the hemisWM (slope>0.63, R²>0.4). All the correlations were attenuated but remained significant after applying PVC, except for the ChPlex. In MC simulations, the hemisWM and CSF were the main contributors to PVE in all BA (slopes 0.15-0.26 and 0.13-0.21 respectively). The hemisWM (slope=0.2), as well as the ChPlex (slope=0.02), influenced SUVRs in the hippocampus. The CerebWM was negatively correlated with all target regions (slope<-0.02, R²>0.8). While no other correlations between OFF and target regions were found, hemisWM was correlated with all OFF regions but the cerebWM (slopes 0.06-0.33). HemisWM correlations attenuated (slopes<0.06) when applying two-compartment PVC, but the hippocampus-ChPlex and the cerebWM correlations required more complex PVC with dedicated compartments for these regions. In-vivo, PVC removed a notably higher fraction of the correlation between OFF regions found to be affected by PVE in the simulation studies and BA (≈50%) than for OFF regions not affected by PVE (16%).

CONCLUSION

HemisWM is the main driver of spill-in effects in FTP PET, affecting both target regions and the rest of OFF regions. PVC successfully reduces PVE, even when using a simple two-compartment method. Despite PVC, non-zero correlations were still observed between target and OFF regions in vivo, which suggests the existence of biological or tracer-related contributions to these correlations.

Deep learning-based image reconstruction and post-processing methods in positron emission tomography for low-dose imaging and resolution enhancement

Image processing plays a crucial role in maximising diagnostic quality of positron emission tomography (PET) images. Recently, deep learning methods developed across many fields have shown tremendous potential when applied to medical image enhancement, resulting in a rich and rapidly advancing literature surrounding this subject. This review encapsulates methods for integrating deep learning into PET image reconstruction and post-processing for low-dose imaging and resolution enhancement. A brief introduction to conventional image processing techniques in PET is firstly presented. We then review methods which integrate deep learning into the image reconstruction framework as either deep learning-based regularisation or as a fully data-driven mapping from measured signal to images. Deep learning-based post-processing methods for low-dose imaging, temporal resolution enhancement and spatial resolution enhancement are also reviewed. Finally, the challenges associated with applying deep learning to enhance PET images in the clinical setting are discussed and future research directions to address these challenges are presented.

Musiré: multimodal simulation and reconstruction framework for the radiological imaging sciences

A software-based workflow is proposed for managing the execution of simulation and image reconstruction for SPECT, PET, CBCT, MRI, BLI and FMI packages in single and multimodal biomedical imaging applications. The workflow is composed of a Bash script, the purpose of which is to provide an interface to the user, and to organize data flow between dedicated programs for simulation and reconstruction. The currently incorporated simulation programs comprise GATE for Monte Carlo simulation of SPECT, PET and CBCT, SpinScenario for simulating MRI, and Lipros for Monte Carlo simulation of BLI and FMI. Currently incorporated image reconstruction programs include CASToR for SPECT and PET as well as RTK for CBCT. MetaImage (mhd) standard is used for voxelized phantom and image data format. Meshlab project (mlp) containers incorporating polygon meshes and point clouds defined by the Stanford triangle format (ply) are employed to represent anatomical structures for optical simulation, and to represent tumour cell inserts. A number of auxiliary programs have been developed for data transformation and adaptive parameter assignment. The software workflow uses fully automatic distribution to, and consolidation from, any number of Linux workstations and CPU cores. Example data are presented for clinical SPECT, PET and MRI systems using the Mida head phantom and for preclinical X-ray, PET and BLI systems employing the Digimouse phantom. The presented method unifies and simplifies multimodal simulation setup and image reconstruction management and might be of value for synergistic image research. This article is part of the theme issue 'Synergistic tomographic image reconstruction: part 2'.