Plane-dependent ML scatter scaling: 3D extension of the 2D simulated single scatter (SSS) estimate

Improved Positron Emission Tomography Quantification: Evaluation of a Maximum-Likelihood Scatter Scaling Algorithm

Incorrect scatter scaling of positron emission tomography (PET) images can lead to halo artifacts, quantitative bias, or reconstruction failure. Tail-fitted scatter scaling (TFSS) possesses performance limitations in multiple cases. This study aims to investigate a novel method for scatter scaling: maximum-likelihood scatter scaling (MLSS) in scenarios where TFSS tends to induce artifacts or are observed to cause reconstruction abortion. [68Ga]Ga-RGD PET scans of nine patients were included in cohort 1 in the scope of investigating the reduction of halo artifacts relative to the scatter estimation method. PET scans of 30 patients administrated with [68Ga]Ga-uPAR were included in cohort 2, used for an evaluation of the robustness of MLSS in cases where TFSS-integrated reconstructions are observed to fail. A visual inspection of MLSS-corrected images scored higher than TFSS-corrected reconstructions of cohort 1. The quantitative investigation near the bladder showed a relative difference in tracer uptake of up to 94.7%. A reconstruction of scans included in cohort 2 resulted in failure in 23 cases when TFSS was used. The lesion uptake values of cohort 2 showed no significant difference. MLSS is suggested as an alternative scatter-scaling method relative to TFSS with the aim of reducing halo artifacts and a robust reconstruction process.

MCGPU-PET: An Open-Source Real-Time Monte Carlo PET Simulator

Monte Carlo (MC) simulations are commonly used to model the emission, transmission, and/or detection of radiation in Positron Emission Tomography (PET). In this work, we introduce a new open-source MC software for PET simulation, MCGPU-PET, which has been designed to fully exploit the computing capabilities of modern GPUs to simulate the acquisition of more than 100 million coincidences per second from voxelized sources and material distributions. The new simulator is an extension of the PENELOPE-based MCGPU code previously used in cone-beam CT and mammography applications. We validated the accuracy of the accelerated code by comparing it to GATE and PeneloPET simulations achieving an agreement within 10 percent approximately. As an example application of the code for fast estimation of PET coincidences, a scan of the NEMA IQ phantom was simulated. A fully 3D sinogram with 6382 million true coincidences and 731 million scatter coincidences was generated in 54 seconds in one GPU. MCGPU-PET provides an estimation of true and scatter coincidences and spurious background (for positron-gamma emitters such as 124I) at a rate 3 orders of magnitude faster than CPU-based MC simulators. This significant speed-up enables the use of the code for accurate scatter and prompt-gamma background estimations within an iterative image reconstruction process.

Comparison of maximum likelihood and conventional PET scatter scaling methods for 18 F-FDG and 68 Ga-DOTATATE PET/CT

PURPOSE

We aim to quantify differences between a new maximum likelihood (ML) background scaling (MLBS) algorithm and two conventional scatter scaling methods for clinical PET/CT. A common source of reduced image quantification with conventional scatter corrections is attributed to erroneous scaling of the initial scatter estimate to match acquired scattered events in the sinogram. MLBS may have performance advantages over conventional methods by using all available data intersecting the subject.

METHODS

A retrospective analysis was performed on subjects injected with ¹⁸ F-FDG (N = 71) and ⁶⁸ Ga-DOTATATE (N = 11) and imaged using time-of-flight (TOF) PET/CT. The scatter distribution was estimated with single scatter simulation approaches. Conventional scaling algorithms included (a) tail fitted background scaling (TFBS), which scales the scatter to "tails" outside the emission support, and (b) absolute scatter correction (ABS), which utilizes the simulated scatter distribution with no scaling applied. MLBS consisted of an alternating iterative reconstruction with a TOF-based ML activity image update allowing negative values (NEG-ML) and nested loop ML scatter scaling estimation. Scatter corrections were compared using reconstructed images as follows: (a) normalized relative difference images were generated and used for voxel-wise analysis, (b) liver and suspected lesion ROIs were drawn to compute mean SUVs, and (c) a qualitative analysis of overall diagnostic image quality, impact of artifacts, and lesion conspicuity was performed. Absolute quantification and normalized relative differences were also assessed with an ¹⁸ F-FDG phantom study.

RESULTS

For human subjects ¹⁸ F-FDG data, Bland-Altman plots demonstrated that the largest normalized voxel-wise differences were observed close to the lower limit (SUV = 1.0). MLBS reconstructions trended towards higher scatter fractions compared to TFBS and ABS images, with median voxel differences across all subjects for TFBS-MLBS measured at 1.7% and 7.6% for ¹⁸ F-FDG and ⁶⁸ Ga-DOTATATE, respectively. For mean SUV analysis, there was a high degree of correlation between the scatter corrections. For ¹⁸ F-FDG, ABS scatter correction reconstructions trended towards higher liver mean SUVs relative to MLBS. The qualitative image analysis revealed no significant differences between TFBS and MLBS image reconstructions. For a uniformly filled relatively large 37 cm diameter phantom, MLBS produced the lowest bias in absolute quantification, while normalized voxel-wise differences showed a trend in scatter correction performance consistent with the human subjects study.

CONCLUSIONS

For ¹⁸ F-FDG, MLBS is at least a valid substitute to TFBS, providing reconstructed image performance comparable to TFBS in most subjects but exhibiting quantitative differences in cases where TFBS is typically prone to inaccuracies (e.g., due to patient motion and CT-based attenuation map truncation). Particularly for low contrast regions, quantification differs for ABS compared to MLBS and TFBS, and caution should be taken when utilizing ABS for decision-making based on quantitative metrics.

Time efficient scatter correction for time-of-flight PET: the immediate scatter approximation

Utilization of time-of-flight (TOF) information allows us to improve image quality and convergence rate in iterative PET image reconstruction. In order to obtain quantitatively correct images accurate scatter correction (SC) is required that accounts for the non-uniform distribution of scatter events over the TOF bins. However, existing simplified TOF-SC algorithms frequently exhibit limited accuracy while the currently accepted reference method-the TOF extension of the single scatter simulation approach (TOF-SSS)-is computationally demanding and can substantially slow down the reconstruction. In this paper we propose and evaluate a new accelerated TOF-SC algorithm in order to improve this situation. The key idea of the algorithm is the use of an immediate scatter approximation (ISA) for scatter time distribution calculation which speeds up estimation of the required TOF scatter by a factor of up to five in comparison to TOF-SSS. The proposed approach was evaluated in dedicated phantom measurements providing challenging high activity contrast conditions as well as in representative clinical patient data sets. Our results show that ISA is a viable alternative to TOF-SSS. The reconstructed images are in excellent quantitative agreement with those obtained with TOF-SSS while overall reconstruction time can be reduced by a factor of two in whole-body studies, even when using a listmode reconstruction not optimized for speed.

Joint Reconstruction of Activity and Attenuation in Time-of-Flight PET: A Quantitative Analysis

Methods for joint activity reconstruction and attenuation reconstruction of time-of-flight (TOF) PET data provide an effective solution to attenuation correction when no (or incomplete or inaccurate) information on attenuation is available. One of the main barriers limiting use of these methods in clinical practice is their lack of validation in a relatively large patient database. In this contribution, we aim to validate reconstruction performed with maximum-likelihood activity reconstruction and attenuation registration (MLRR) in a whole-body patient dataset. Furthermore, a partial validation (because the scale problem of the algorithm is avoided for now) of reconstruction performed with maximum-likelihood activity and attenuation (MLAA) is also provided. We present a quantitative comparison between these 2 methods of joint reconstruction and the current clinical gold standard, maximum-likelihood expectation maximization (MLEM) with CT-based attenuation correction. Methods: The whole-body TOF PET emission data of each patient dataset were processed as a whole to reconstruct an activity volume covering all the acquired bed positions, helping reduce the problem of a scale per bed position in MLAA to a global scale for the entire activity volume. Three reconstruction algorithms were used: MLEM, MLRR, and MLAA. A maximum-likelihood scaling of the single-scatter simulation estimate to the emission data was used for scatter correction. The reconstruction results for various regions of interest were then analyzed. Results: The joint reconstructions of the whole-body patient dataset provided better quantification than the gold standard in cases of PET and CT misalignment caused by patient or organ motion. Our quantitative analysis showed a difference of -4.2% ± 2.3% between MLRR and MLEM and a difference of -7.5% ± 4.6% between MLAA and MLEM, averaged over all regions of interest. Conclusion: Joint reconstruction of activity and attenuation provides a useful means to estimate tracer distribution when CT-based-attenuation images are subject to misalignment or are not available. With an accurate estimate of the scatter contribution in the emission measurements, the joint reconstructions of TOF PET data are within clinically acceptable accuracy.

Clinical Evaluation of 68Ga-PSMA-II and 68Ga-RM2 PET Images Reconstructed With an Improved Scatter Correction Algorithm

OBJECTIVE

Gallium-68-labeled radiopharmaceuticals pose a challenge for scatter estimation because their targeted nature can produce high contrast in these regions of the kidneys and bladder. Even small errors in the scatter estimate can result in washout artifacts. Administration of diuretics can reduce these artifacts, but they may result in adverse events. Here, we investigated the ability of algorithmic modifications to mitigate washout artifacts and eliminate the need for diuretics or other interventions.

MATERIALS AND METHODS

The model-based scatter algorithm was modified to account for PET/MRI scanner geometry and challenges of non-FDG tracers. Fifty-three clinical 68Ga-RM2 and 68Ga-PSMA-11 whole-body images were reconstructed using the baseline scatter algorithm. For comparison, reconstruction was also processed with modified sampling in the single-scatter estimation and with an offset in the scatter tail-scaling process. None of the patients received furosemide to attempt to decrease the accumulation of radiopharmaceuticals in the bladder. The images were scored independently by three blinded reviewers using the 5-point Likert scale.

RESULTS

The scatter algorithm improvements significantly decreased or completely eliminated the washout artifacts. When comparing the baseline and most improved algorithm, the image quality increased and image artifacts were reduced for both 68Ga-RM2 and for 68Ga-PSMA-11 in the kidneys and bladder regions.

CONCLUSION

Image reconstruction with the improved scatter correction algorithm mitigated washout artifacts and recovered diagnostic image quality in 68Ga PET, indicating that the use of diuretics may be avoided.

Scatter Artifact with Ga-68-PSMA-11 PET: Severity Reduced With Furosemide Diuresis and Improved Scatter Correction

PURPOSE

To assess the utility of furosemide diuresis and the role of an improved scatter correction algorithm in reducing scatter artifact severity on Ga-68- Prostate-specific membrane antigen (PSMA)-11 positron emission tomography (PET).

MATERIALS AND METHODS

A total of 139 patients underwent Ga-68-PSMA-11 PET imaging for prostate cancer: 47 non-time-of-flight (non-TOF) PET/computed tomography, 51 PET/magnetic resonance imaging (MRI) using the standard TOF scatter correction algorithm, and 41 PET/MRI using an improved TOF scatter correction algorithm. Whole-body PET acquisitions were subdivided into 3 regions: around kidneys; between kidneys and bladder; and around bladder. The images were reviewed, and scatter artifact severity was rated using a Likert-type scale.

RESULTS

The worst scatter occurred when using non-TOF scatter correction without furosemide, where 42.1% of patients demonstrated severe scatter artifacts in 1 or more regions. Improved TOF scatter correction resulted in the smallest percentage of studies with severe scatter (6.5%). Scatter ratings by region were lowest using improved TOF scatter correction. Furosemide reduced mean scatter severity when using non-TOF and standard TOF.

CONCLUSIONS

Both furosemide and scatter correction algorithm play a role in reducing scatter in PSMA PET. Improved TOF scatter correction resulted in the lowest scatter severity.

Simultaneous whole-body 18F-PSMA-1007-PET/MRI with integrated high-resolution multiparametric imaging of the prostatic fossa for comprehensive oncological staging of patients with prostate cancer: a pilot study

INTRODUCTION

The aim of the present study was to explore the clinical feasibility and reproducibility of a comprehensive whole-body ¹⁸F-PSMA-1007-PET/MRI protocol for imaging prostate cancer (PC) patients.

METHODS

Eight patients with high-risk biopsy-proven PC underwent a whole-body PET/MRI (3 h p.i.) including a multi-parametric prostate MRI after ¹⁸F-PSMA-1007-PET/CT (1 h p.i.) which served as reference. Seven patients presented with non-treated PC, whereas one patient presented with biochemical recurrence. SUV_mean-quantification was performed using a 3D-isocontour volume-of-interest. Imaging data was consulted for TNM-staging and compared with histopathology. PC was confirmed in 4/7 patients additionally by histopathology after surgery. PET-artifacts, co-registration of pelvic PET/MRI and MRI-data were assessed (PI-RADS 2.0).

RESULTS

The examinations were well accepted by patients and comprised 1 h. SUV_mean-values between PET/CT (1 h p.i.) and PET/MRI (3 h p.i.) were significantly correlated (p < 0.0001, respectively) and similar to literature of ¹⁸F-PSMA-1007-PET/CT 1 h vs 3 h p.i. The dominant intraprostatic lesion could be detected in all seven patients in both PET and MRI. T2c, T3a, T3b and T4 features were detected complimentarily by PET and MRI in five patients. PET/MRI demonstrated moderate photopenic PET-artifacts surrounding liver and kidneys representing high-contrast areas, no PET-artifacts were observed for PET/CT. Simultaneous PET-readout during prostate MRI achieved optimal co-registration results.

CONCLUSIONS

The presented ¹⁸F-PSMA-1007-PET/MRI protocol combines efficient whole-body assessment with high-resolution co-registered PET/MRI of the prostatic fossa for comprehensive oncological staging of patients with PC.

Investigation of the halo-artifact in 68Ga-PSMA-11-PET/MRI

OBJECTIVES

Combined positron emission tomography (PET) and magnetic resonance imaging (MRI) targeting the prostate-specific membrane antigen (PSMA) with a 68Ga-labelled PSMA-analog (68Ga-PSMA-11) is discussed as a promising diagnostic method for patients with suspicion or history of prostate cancer. One potential drawback of this method are severe photopenic (halo-) artifacts surrounding the bladder and the kidneys in the scatter-corrected PET images, which have been reported to occur frequently in clinical practice. The goal of this work was to investigate the occurrence and impact of these artifacts and, secondly, to evaluate variants of the standard scatter correction method with regard to halo-artifact suppression.

METHODS

Experiments using a dedicated pelvis phantom were conducted to investigate whether the halo-artifact is modality-, tracer-, and/or concentration-dependent. Furthermore, 31 patients with history of prostate cancer were selected from an ongoing 68Ga-PSMA-11-PET/MRI study. For each patient, PET raw data were reconstructed employing six different variants of PET scatter correction: absolute scatter scaling, relative scatter scaling, and relative scatter scaling combined with prompt gamma correction, each of which was combined with a maximum scatter fraction (MaxSF) of MaxSF = 75% or MaxSF = 40%. Evaluation of the reconstructed images with regard to halo-artifact suppression was performed both quantitatively using statistical analysis and qualitatively by two independent readers.

RESULTS

The phantom experiments did not reveal any modality-dependency (PET/MRI vs. PET/CT) or tracer-dependency (68Ga vs. 18F-FDG). Patient- and phantom-based data indicated that halo-artifacts derive from high organ-to-background activity ratios (OBR) between bladder/kidneys and surrounding soft tissue, with a positive correlation between OBR and halo size. Comparing different variants of scatter correction, reducing the maximum scatter fraction from the default value MaxSF = 75% to MaxSF = 40% was found to efficiently suppress halo-artifacts in both phantom and patient data. In 1 of 31 patients, reducing the maximum scatter fraction provided new PET-based information changing the patient's diagnosis.

CONCLUSION

Halo-artifacts are particularly observed for 68Ga-PSMA-11-PET/MRI due to 1) the biodistribution of the PSMA-11-tracer resulting in large OBRs for bladder and kidneys and 2) inaccurate scatter correction methods currently used in clinical routine, which tend to overestimate the scatter contribution. If not compensated for, 68Ga-PSMA-11 uptake pathologies may be masked by halo-artifacts leading to false-negative diagnoses. Reducing the maximum scatter fraction was found to efficiently suppress halo-artifacts.