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.

Full field spatially-variant image-based resolution modelling reconstruction for the HRRT.

Accurate characterisation of the scanner's point spread function across the entire field of view (FOV) is crucial in order to account for spatially dependent factors that degrade the resolution of the reconstructed images. The HRRT users' community resolution modelling reconstruction software includes a shift-invariant resolution kernel, which leads to transaxially non-uniform resolution in the reconstructed images. Unlike previous work to date in this field, this work is the first to model the spatially variant resolution across the entire FOV of the HRRT, which is the highest resolution human brain PET scanner in the world. In this paper we developed a spatially variant image-based resolution modelling reconstruction dedicated to the HRRT, using an experimentally measured shift-variant resolution kernel. Previously, the system response was measured and characterised in detail across the entire FOV of the HRRT, using a printed point source array. The newly developed resolution modelling reconstruction was applied on measured phantom, as well as clinical data and was compared against the HRRT users' community resolution modelling reconstruction, which is currently in use. Results demonstrated improvements both in contrast and resolution recovery, particularly for regions close to the edges of the FOV, with almost uniform resolution recovery across the entire transverse FOV. In addition, because the newly measured resolution kernel is slightly broader with wider tails, compared to the deliberately conservative kernel employed in the HRRT users' community software, the reconstructed images appear to have not only improved contrast recovery (up to 20% for small regions), but also better noise characteristics.

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 [(18)F]FDG and [(15)O]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.

Evaluation of a direct 4D reconstruction method using generalised linear least squares for estimating nonlinear micro-parametric maps.

OBJECTIVE: Estimation of nonlinear micro-parameters is a computationally demanding and fairly challenging process, since it involves the use of rather slow iterative nonlinear fitting algorithms and it often results in very noisy voxel-wise parametric maps. Direct reconstruction algorithms can provide parametric maps with reduced variance, but usually the overall reconstruction is impractically time consuming with common nonlinear fitting algorithms. METHODS: In this work we employed a recently proposed direct parametric image reconstruction algorithm to estimate the parametric maps of all micro-parameters of a two-tissue compartment model, used to describe the kinetics of [[Formula: see text]F]FDG. The algorithm decouples the tomographic and the kinetic modelling problems, allowing the use of previously developed post-reconstruction methods, such as the generalised linear least squares (GLLS) algorithm. RESULTS: Results on both clinical and simulated data showed that the proposed direct reconstruction method provides considerable quantitative and qualitative improvements for all micro-parameters compared to the conventional post-reconstruction fitting method. Additionally, region-wise comparison of all parametric maps against the well-established filtered back projection followed by post-reconstruction non-linear fitting, as well as the direct Patlak method, showed substantial quantitative agreement in all regions. CONCLUSIONS: The proposed direct parametric reconstruction algorithm is a promising approach towards the estimation of all individual microparameters of any compartment model. In addition, due to the linearised nature of the GLLS algorithm, the fitting step can be very efficiently implemented and, therefore, it does not considerably affect the overall reconstruction time.

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

PURPOSE: 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. METHODS: Using an array of 12 × 7 printed point sources, along with a custom made phantom, and with the MR magnet on, the system's spatially variant image-based PSF was characterized in detail. Moreover, basis function parameters were systematically varied during reconstruction (list-mode TF OSEM) to evaluate their impact on the reconstructed resolution and the image space PSF. Following the spatial resolution optimization, phantom, and clinical studies were subsequently reconstructed using representative basis function parameters. RESULTS: Based on the analysis and under standard basis function parameters, the axial and tangential components of the PSF were found to be almost invariant under spatial transformations (~4 mm) while the radial component varied modestly from 4 to 6.7 mm. Using a systematic investigation into the basis function parameter space, the spatial resolution was found to degrade for basis functions with a large radius and small shape parameter. However, it was found that optimizing the spatial resolution in the reconstructed PET images, while having a good basis function superposition and keeping the image representation error to a minimum, is feasible, with the parameter combination range depending upon the scanner's intrinsic resolution characteristics. CONCLUSIONS: Using the printed point source array as a MR compatible methodology for experimentally measuring the scanner's PSF, the system's spatially variant resolution properties were successfully evaluated in image space. Overall the PET subsystem exhibits excellent resolution characteristics mainly due to the fact that the raw data are not under-sampled/rebinned, enabling the spatial resolution to be dictated by the scanner's intrinsic resolution and the image reconstruction parameters. Due to the impact of these parameters on the resolution properties of the reconstructed images, the image space PSF varies both under spatial transformations and due to basis function parameter selection. Nonetheless, for a range of basis function parameters, the image space PSF remains unaffected, with the range depending on the scanner's intrinsic resolution properties.

Isotope specific resolution recovery image reconstruction in high resolution PET imaging.

PURPOSE: Measuring and incorporating a scanner-specific point spread function (PSF) within image reconstruction has been shown to improve spatial resolution in PET. However, due to the short half-life of clinically used isotopes, other long-lived isotopes not used in clinical practice are used to perform the PSF measurements. As such, non-optimal PSF models that do not correspond to those needed for the data to be reconstructed are used within resolution modeling (RM) image reconstruction, usually underestimating the true PSF owing to the difference in positron range. In high resolution brain and preclinical imaging, this effect is of particular importance since the PSFs become more positron range limited and isotope-specific PSFs can help maximize the performance benefit from using resolution recovery image reconstruction algorithms. METHODS: In this work, the authors used a printing technique to simultaneously measure multiple point sources on the High Resolution Research Tomograph (HRRT), and the authors demonstrated the feasibility of deriving isotope-dependent system matrices from fluorine-18 and carbon-11 point sources. Furthermore, the authors evaluated the impact of incorporating them within RM image reconstruction, using carbon-11 phantom and clinical datasets on the HRRT. RESULTS: The results obtained using these two isotopes illustrate that even small differences in positron range can result in different PSF maps, leading to further improvements in contrast recovery when used in image reconstruction. The difference is more pronounced in the centre of the field-of-view where the full width at half maximum (FWHM) from the positron range has a larger contribution to the overall FWHM compared to the edge where the parallax error dominates the overall FWHM. CONCLUSIONS: Based on the proposed methodology, measured isotope-specific and spatially variant PSFs can be reliably derived and used for improved spatial resolution and variance performance in resolution recovery image reconstruction. The benefits are expected to be more substantial for more energetic positron emitting isotopes such as Oxygen-15 and Rubidium-82.