Cloud-based serverless computing enables accelerated monte carlo simulations for nuclear medicine imaging.

Objective.This study investigates the potential of cloud-based serverless computing to accelerate Monte Carlo (MC) simulations for nuclear medicine imaging tasks. MC simulations can pose a high computational burden-even when executed on modern multi-core computing servers. Cloud computing allows simulation tasks to be highly parallelized and considerably accelerated.Approach.We investigate the computational performance of a cloud-based serverless MC simulation of radioactive decays for positron emission tomography imaging using Amazon Web Service (AWS) Lambda serverless computing platform for the first time in scientific literature. We provide a comparison of the computational performance of AWS to a modern on-premises multi-thread reconstruction server by measuring the execution times of the processes using between105and2·1010simulated decays. We deployed two popular MC simulation frameworks-SimSET and GATE-within the AWS computing environment. Containerized application images were used as a basis for an AWS Lambda function, and local (non-cloud) scripts were used to orchestrate the deployment of simulations. The task was broken down into smaller parallel runs, and launched on concurrently running AWS Lambda instances, and the results were postprocessed and downloaded via the Simple Storage Service.Main results.Our implementation of cloud-based MC simulations with SimSET outperforms local server-based computations by more than an order of magnitude. However, the GATE implementation creates more and larger output file sizes and reveals that the internet connection speed can become the primary bottleneck for data transfers. Simulating 109decays using SimSET is possible within 5 min and accrues computation costs of about $10 on AWS, whereas GATE would have to run in batches for more than 100 min at considerably higher costs.Significance.Adopting cloud-based serverless computing architecture in medical imaging research facilities can considerably improve processing times and overall workflow efficiency, with future research exploring additional enhancements through optimized configurations and computational methods.

Super-resolution reconstruction of γ-ray CT images for PET-enabled dual-energy CT imaging.

Dual-energy computed tomography (DECT) enables material decomposition for tissues and produces additional information for PET/CT imaging to potentially improve the characterization of diseases. PET-enabled DECT (PDECT) allows the generation of PET and DECT images simultaneously with a conventional PET/CT scanner without the need for a second x-ray CT scan. In PDECT, high-energy γ-ray CT (GCT) images at 511 keV are obtained from time-of-flight (TOF) PET data and are combined with the existing x-ray CT images to form DECT imaging. We have developed a kernel-based maximum-likelihood attenuation and activity (MLAA) method that uses x-ray CT images as a priori information for noise suppression. However, our previous studies focused on GCT image reconstruction at the PET image resolution which is coarser than the image resolution of the x-ray CT. In this work, we explored the feasibility of generating super-resolution GCT images at the corresponding CT resolution. The study was conducted using both phantom and patient scans acquired with the uEXPLORER total-body PET/CT system. GCT images at the PET resolution with a pixel size of 4.0 mm × 4.0 mm and at the CT resolution with a pixel size of 1.2 mm × 1.2 mm were reconstructed using both the standard MLAA and kernel MLAA methods. The results indicated that the GCT images at the CT resolution had sharper edges and revealed more structural details compared to the images reconstructed at the PET resolution. Furthermore, images from the kernel MLAA method showed substantially improved image quality compared to those obtained with the standard MLAA method.

Feasibility of PET-enabled dual-energy CT imaging: First physical phantom and patient results.

X-ray computed tomography (CT) in PET/CT is commonly operated with a single energy, resulting in a limitation of lacking tissue composition information. Dual-energy (DE) spectral CT enables material decomposition by using two different x-ray energies and may be combined with PET for improved multimodality imaging, but would either require hardware upgrade or increase radiation dose due to the added second x-ray CT scan. Recently proposed PET-enabled DECT method allows dual-energy spectral imaging using a conventional PET/CT scanner without the need for a second x-ray CT scan. A gamma-ray CT (gCT) image at 511 keV can be generated from the existing time-of-flight PET data with the maximum-likelihood attenuation and activity (MLAA) approach and is then combined with the low-energy x-ray CT image to form dual-energy spectral imaging. To improve the image quality of gCT, a kernel MLAA method was further proposed by incorporating x-ray CT as a priori information. The concept of this PET-enabled DECT has been validated using simulation studies, but not yet with 3D real data. In this work, we developed a general open-source implementation for gCT reconstruction from PET data and use this implementation for the first real data validation with both a physical phantom study and a human subject study on a uEXPLORER total-body PET/CT system. These results have demonstrated the feasibility of this method for spectral imaging and material decomposition.

Development of a Monte Carlo-based scatter correction method for total-body PET using the uEXPLORER PET/CT scanner.

Objective.This study presents and evaluates a robust Monte Carlo-based scatter correction (SC) method for long axial field of view (FOV) and total-body positron emission tomography (PET) using the uEXPLORER total-body PET/CT scanner.Approach.Our algorithm utilizes the Monte Carlo (MC) tool SimSET to compute SC factors in between individual image reconstruction iterations within our in-house list-mode and time-of-flight-based image reconstruction framework. We also introduced a unique scatter scaling technique at the detector block-level for optimal estimation of the scatter contribution in each line of response. First image evaluations were derived from phantom data spanning the entire axial FOV along with image data from a human subject with a large body mass index. Data was evaluated based on qualitative inspections, and contrast recovery, background variability, residual scatter removal from cold regions, biases and axial uniformity were quantified and compared to non-scatter-corrected images.Main results.All reconstructed images demonstrated qualitative and quantitative improvements compared to non-scatter-corrected images: contrast recovery coefficients improved by up to 17.2% and background variability was reduced by up to 34.3%, and the residual lung error was between 1.26% and 2.08%. Low biases throughout the axial FOV indicate high quantitative accuracy and axial uniformity of the corrections. Up to 99% of residual activity in cold areas in the human subject was removed, and the reliability of the method was demonstrated in challenging body regions like in the proximity of a highly attenuating knee prosthesis.Significance.The MC SC method employed was demonstrated to be accurate and robust in TB-PET. The results of this study can serve as a benchmark for optimizing the quantitative performance of future SC techniques.

Numerical investigation reveals challenges in measuring the contrast recovery coefficients in PET.

Objective.Contrast recovery coefficient (CRC) is essential for image quality (IQ) assessment in positron emission tomography (PET), typically measured according to the National Electrical Manufacturers Association (NEMA) NU 2 standard. This study quantifies systematic uncertainties of the CRC measurement by a numerical investigation of the effects from scanner-independent parameters like voxel size, region-of-interest (ROI) misplacement, and sphere position on the underlying image grid.Approach.CRC measurements with 2D and 3D ROIs were performed on computer-generated images of a NEMA IQ-like phantom, using voxel sizes of 1-4 mm for sphere diameters of 5-40 mm-first in absence of noise and blurring, then with simulated spatial resolution and image noise with varying noise levels. The systematic uncertainties of the CRC measurement were quantified from above variations of scanner-independent parameters. Subsampled experimental images of a NEMA IQ phantom were additionally used to investigate the impact of ROI misplacement at different noise levels.Main results.In absence of noise and blurring, systematic uncertainties were up to 28.8% and 31.0% with 2D and 3D ROIs, respectively, for the 10 mm sphere, with the highest impact from ROI misplacement. In all cases, smaller spheres showed higher uncertainties with larger voxels. Contrary to prior assumptions, the use of 3D ROIs did not exhibit less susceptibility for parameter changes. Experimental and computer-generated images both demonstrated considerable variations on individual CRC measurements when background coefficient-of-variation exceeded 20%, despite negligible effects on mean CRC.Significance.This study underscores the effect of scanner-independent parameters on reliability, reproducibility, and comparability of CRC measurements. Our findings highlight the trade-off between the benefits of smaller voxel sizes and noise-induced CRC fluctuations, which is not considered in the current version of the NEMA IQ standards. The results furthermore warrant adjustments to the standard to accommodate the advances in sensitivity and spatial resolution of current-generation PET scanners.

Relating18F-FDG image signal-to-noise ratio to time-of-flight noise-equivalent count rate in total-body PET using the uEXPLORER scanner.

Objective.This work assessed the relationship between image signal-to-noise ratio (SNR) and total-body noise-equivalent count rate (NECR)-for both non-time-of-flight (TOF) NECR and TOF-NECR-in a long uniform water cylinder and 14 healthy human subjects using the uEXPLORER total-body PET/CT scanner.Approach.A TOF-NEC expression was modified for list-mode PET data, and both the non-TOF NECR and TOF-NECR were compared using datasets from a long uniform water cylinder and 14 human subjects scanned up to 12 h after radiotracer injection.Main results.The TOF-NECR for the uniform water cylinder was found to be linearly proportional to the TOF-reconstructed image SNR2in the range of radioactivity concentrations studied, but not for non-TOF NECR as indicated by the reducedR2value. The results suggest that the use of TOF-NECR to estimate the count rate performance of TOF-enabled PET systems may be more appropriate for predicting the SNR of TOF-reconstructed images.Significance.Image quality in PET is commonly characterized by image SNR and, correspondingly, the NECR. While the use of NECR for predicting image quality in conventional PET systems is well-studied, the relationship between SNR and NECR has not been examined in detail in long axial field-of-view total-body PET systems, especially for human subjects. Furthermore, the current NEMA NU 2-2018 standard does not account for count rate performance gains due to TOF in the NECR evaluation. The relationship between image SNR and total-body NECR in long axial FOV PET was assessed for the first time using the uEXPLORER total-body PET/CT scanner.