Systematic Evaluation of the Impact of Involuntary Motion in Whole Body Dynamic PET

Involuntary patient motion can happen in dynamic whole body (DWB) PET due to long scanning times, which may cause inaccurate quantification of tissue parameters. To quantify the impact on Patlak parameters, we simulated dynamic data using patient-derived motion fields, systematically introducing the motion at different passes of the dynamic scan, both inter and intra-frame. Estimated parameters are compared against the ground truth. Results show that errors can be large, even for small motion. Caution is advised when quantitatively evaluating DWB-PET images, if any motion has been detected

NEMA NU 2-2018 performance evaluation of a new generation 30-cm axial field-of-view Discovery MI PET/CT.

PURPOSE The DMI PET/CT is a modular silicon photomultiplier-based scanner with an axial field-of-view (FOV) between 15 and 25 cm depending on ring configuration (3, 4, or 5 rings). A new generation of the system includes a reengineered detector module, featuring improved electronics and an additional 6th ring, extending the axial FOV to 30 cm. We report on the performance evaluation of the 6-ring upgraded Generation 2 (Gen2) system while values are also reported for the 5-ring configuration of the very same system prior to the upgrade. METHODS PET performance was evaluated using the NEMA NU 2-2018 standard for spatial resolution, sensitivity, image quality, count rate performance, timing resolution, and image co-registration accuracy. Patient images were used to assess image quality. RESULTS The average system sensitivity was measured at 32.76 cps/kBq (~ 47% increase to 5 rings at 22.29 cps/kBq) while noise equivalent count rate peaked at 434.3 kcps corresponding to 23.6 kBq/mL (~ 60% increase to Generation 1 (Gen1) and 39% to Gen2 5 rings). Contrast recovery ranged between 54.5 and 85.8% similar to 5 rings, while the 6 rings provided lower background variability (2.3-8.5% for 5 rings vs 1.9-6.8% for 6 rings) and lower lung error (4.0% for the 5 rings and 3.16% for the 6 rings). Transverse/axial full width at half-maximum (FWHM) at 1 cm (3.79/4.26 mm) and 10 cm (4.29/4.55 mm), scatter fraction (40.2%), energy resolution (9.63%), and time-of-flight (TOF) resolution (389.6 ps at 0 kBq/mL) were in line to previously reported values measured across different system configurations. Improved patient image quality is obtained with the 6 rings compared to the 5 rings, while image quality is retained even at reduced scan times, enabling WB dynamic acquisitions. CONCLUSIONS The higher sensitivity of the 6-ring DMI compared to the 5-ring configuration may lead to improved image quality of clinical images at reduced scan time. Additionally, it could equally be used to allow improved temporal sampling and/or reduced overall scan time in dynamic acquisitions. Conversely, temporal sampling and scan time could be traded per application to further drive injected dose at lower levels

An Investigation of Stochastic Variance Reduction Algorithms for Relative Difference Penalised 3D PET Image Reconstruction

Penalised PET image reconstruction algorithms are often accelerated during early iterations with the use of subsets. However, these methods may exhibit limit cycle behaviour at later iterations due to variations between subsets. Desirable converged images can be achieved for a subclass of these algorithms via the implementation of a relaxed step size sequence, but the heuristic selection of parameters will impact the quality of the image sequence and algorithm convergence rates. In this work, we demonstrate the adaption and application of a class of stochastic variance reduction gradient algorithms for PET image reconstruction using the relative difference penalty and numerically compare convergence performance to BSREM. The two investigated algorithms are: SAGA and SVRG. These algorithms require the retention in memory of recently computed subset gradients, which are utilised in subsequent updates. We present several numerical studies based on Monte Carlo simulated data and a patient data set for fully 3D PET acquisitions. The impact of the number of subsets, different preconditioners and step size methods on the convergence of regions of interest values within the reconstructed images is explored. We observe that when using constant preconditioning, SAGA and SVRG demonstrate reduced variations in voxel values between subsequent updates and are less reliant on step size hyper-parameter selection than BSREM reconstructions. Furthermore, SAGA and SVRG can converge significantly faster to the penalised maximum likelihood solution than BSREM, particularly in low count data

An Investigation of Stochastic Variance Reduction Algorithms for Relative Difference Penalized 3D PET Image Reconstruction

Penalised PET image reconstruction algorithms are often accelerated during early iterations with the use of subsets. However, these methods may exhibit limit cycle behaviour at later iterations due to variations between subsets. Desirable converged images can be achieved for a subclass of these algorithms via the implementation of a relaxed step size sequence, but the heuristic selection of parameters will impact the quality of the image sequence and algorithm convergence rates. In this work, we demonstrate the adaption and application of a class of stochastic variance reduction gradient algorithms for PET image reconstruction using the relative difference penalty and numerically compare convergence performance to BSREM. The two investigated algorithms are: SAGA and SVRG. These algorithms require the retention in memory of recently computed subset gradients, which are utilised in subsequent updates. We present several numerical studies based on Monte Carlo simulated data and a patient data set for fully 3D PET acquisitions. The impact of the number of subsets, different preconditioners and step size methods on the convergence of regions of interest values within the reconstructed images is explored. We observe that when using constant preconditioning, SAGA and SVRG demonstrate reduced variations in voxel values between subsequent updates and are less reliant on step size hyper-parameter selection than BSREM reconstructions. Furthermore, SAGA and SVRG can converge significantly faster to the penalised maximum likelihood solution than BSREM, particularly in low count data

Impact of different image reconstructions on PET quantification in non-small cell lung cancer: a comparison of adenocarcinoma and squamous cell carcinoma

OBJECTIVE: Positron emission tomography (PET) using 18F-fluordeoxyglucose (F-FDG) is an established imaging modality for tumor staging in patients with non-small cell lung cancer (NSCLC). There is a growing interest in using F-FDG PET for therapy response assessment in NSCLC which relies on quantitative PET parameters such as standardized uptake values (SUV). Different reconstruction algorithms in PET may affect SUV. We sought to determine the variation of SUV in patients with NSCLC when using ordered subset expectation maximization (OSEM) and block sequential regularized expectation maximization (BSREM) in latest-generation digital PET/CT, including a subanalysis for adenocarcinoma and squamous cell carcinoma. METHODS: A total of 58 patients (34 = adenocarcinoma, 24 = squamous cell carcinoma) that underwent a clinically indicated F-FDG PET/CT for staging were reviewed. PET images were reconstructed with OSEM and BSREM reconstruction with noise penalty strength β-levels of 350, 450, 600, 800 and 1200. Lung tumors maximum standardized uptake value (SUV) were compared. RESULTS: Lung tumors SUV were significantly lower in adenocarcinomas compared to squamous cell carcinomas in all reconstructions evaluated (all p 0.05). There was a statistically significant difference of the relative increase of SUV in adenocarcinoma (mean + 34.8%) and squamous cell carcinoma (mean 23.4%), when using BSREM instead of OSEM (p < 0.05). CONCLUSIONS: In NSCLC the relative change of SUV when using BSREM instead of OSEM is significantly higher in adenocarcinoma as compared to squamous cell carcinoma. ADVANCES IN KNOWLEDGE: The impact of BSREM on SUV may vary in different histological subtypes of NSCLC. This highlights the importance for careful standardisation of β-value used for serial F-FDG PET scans when following-up NSCLC patients

A 5D computational phantom for pharmacokinetic simulation studies in dynamic emission tomography

Introduction: Dynamic image acquisition protocols are increasingly used in emission tomography for drug development and clinical research. As such, there is a need for computational phantoms to accurately describe both the spatial and temporal distribution of radiotracers, also accounting for periodic and non-periodic physiological processes occurring during data acquisition. Methods: A new 5D anthropomorphic digital phantom was developed based on a generic simulation platform, for accurate parametric imaging simulation studies in emission tomography. The phantom is based on high spatial and temporal information derived from real 4D MR data and a detailed multi-compartmental pharmacokinetic modelling simulator. Results: The proposed phantom is comprised of three spatial and two temporal dimensions, including periodic physiological processes due to respiratory motion and non-periodic functional processes due to tracer kinetics. Example applications are shown in parametric [18F]FDG and [15O]H2O PET imaging, successfully generating realistic macro- and micro-parametric maps. Conclusions: The envisaged applications of this digital phantom include the development and evaluation of motion correction and 4D image reconstruction algorithms in PET and SPECT, development of protocols and methods for tracer and drug development as well as new pharmacokinetic parameter estimation algorithms, amongst others. Although the simulation platform is primarily developed for generating dynamic phantoms for emission tomography studies, it can easily be extended to accommodate dynamic MR and CT imaging simulation protocols

Experimental evaluation and basis function optimization of the spatially variant image-space PSF on the Ingenuity PET/MR scanner

The Ingenuity time-of-flight (TF) PET/MR is a recently developed hybrid scanner combining the molecular imaging capabilities of PET with the excellent soft tissue contrast of MRI. It is becoming common practice to characterize the system's point spread function (PSF) and understand its variation under spatial transformations to guide clinical studies and potentially use it within resolution recovery image reconstruction algorithms. Furthermore, due to the system's utilization of overlapping and spherical symmetric Kaiser-Bessel basis functions during image reconstruction, its image space PSF and reconstructed spatial resolution could be affected by the selection of the basis function parameters. Hence, a detailed investigation into the multidimensional basis function parameter space is needed to evaluate the impact of these parameters on spatial resolution