A hybrid scatter correction for 3D PET based on an estimation of the distribution of unscattered coincidences: implementation on the ECAT EXACT HR+

Quantification of brain oxygen extraction and metabolism with [15O]-gas PET: A technical review in the era of PET/MRI

Oxygen extraction fraction (OEF) and the cerebral metabolic rate of oxygen (CMRO₂) are key cerebral physiological parameters to identify at-risk cerebrovascular patients and understand brain health and function. PET imaging with [¹⁵O]-oxygen tracers, either through continuous or bolus inhalation, provides non-invasive assessment of OEF and CMRO₂. Numerous tracer delivery, PET acquisition, and kinetic modeling approaches have been adopted to map brain oxygenation. The purpose of this technical review is to critically evaluate different methods for [¹⁵O]-gas PET and its impact on the accuracy and reproducibility of OEF and CMRO₂ measurements. We perform a meta-analysis of brain oxygenation PET studies in healthy volunteers and compare between continuous and bolus inhalation techniques. We also describe OEF metrics that have been used to detect hemodynamic impairment in cerebrovascular disease. For these patients, advanced techniques to accelerate the PET scans and potential synthesis with MRI to avoid arterial blood sampling would facilitate broader use of [¹⁵O]-oxygen PET for brain physiological assessment.

Validation of a simplified scatter correction method for 3D brain PET with 15O

Objective

Positron emission tomography (PET) enables quantitative measurements of various biological functions. Accuracy in data acquisition and processing schemes is a prerequisite for this. The correction of scatter is especially important when a 3D PET scanner is used. The aim of this study was to validate the use of a simplified calculation-based scatter correction method for ¹⁵O studies in the brain.

Methods

We applied two scatter correction methods to the same ¹⁵O PET data acquired from patients with cerebrovascular disease (n = 10): a hybrid dual-energy-window scatter correction (reference method), and a deconvolution scatter correction (simplified method). The PET study included three sequential scans for ¹⁵O-CO, ¹⁵O-O₂, and ¹⁵O-H₂O, from which the following quantitative parameters were calculated, cerebral blood flow, cerebral blood volume, cerebral metabolic rate of oxygen, and oxygen extraction fraction.

Results

Both scatter correction methods provided similar reconstruction images with almost identical image noise, although there were slightly greater differences in white-matter regions compared with gray matter regions. These differences were also greater for ¹⁵O-CO than for ¹⁵O-H₂O and ¹⁵O-O₂. Region of interest analysis of the quantitative parameters demonstrated that the differences were less than 10 % (except for cerebral blood volume in white-matter regions), and the agreement between the methods was excellent, with intraclass correlation coefficients above 0.95 for all the parameters.

Conclusions

The deconvolution scatter correction despite its simplified implementation provided similar results to the hybrid dual-energy-window scatter correction. We consider it suitable for application in a clinical ¹⁵O brain study using a 3D PET scanner.

Applicability of emission-based attenuation map for rapid CBF, OEF, and CMRO2 measurements using gaseous (15)O-labeled compounds

Background

Cerebral blood flow (CBF), oxygen extraction fraction (OEF), and cerebral metabolic rate of oxygen (CMRO₂) images have facilitated understanding of the pathophysiological basis of cerebrovascular disorders. Such parametric images can be rapidly, measured within around 15 min, using positron emission tomography (PET) with sequentially administered ¹⁵O-labeled oxygen and water. For further shortening, one option is to eliminate the transmission scan by applying an emission-based attenuation correction.

Methods

The validity of the present method was tested by comparing parametric values with emission-based attenuation correction to those with transmission-based correction. This was applied to 27 subjects who were diagnosed with or without cerebrovascular disorders. All subjects received the rapid CBF/OEF/CMRO₂ PET measurements. An emission-based attenuation map was generated by estimating the edge of the brain tissue contour on an obtained sinogram and by assuming the uniform tissue coefficient to be 0.1 cm⁻¹. Then images were reconstructed, and parametric images were computed.

Results

No difference was apparent between the emission- and transmission-based methods. Paired t-test showed no significant differences in CBF, OEF, or CMRO₂ values between the emission- and transmission-based methods, except in the parietal and occipital and cerebellum and occipital regions, and the differences were less than 10%. The regression analysis showed a close correlation of r = 0.89 to 0.99.

Conclusions

The present study revealed that the attenuation correction can be performed by the emission-based estimation method and clinical PET duration can be shortened for the CBF, OEF, and CMRO₂ gas study.

Comparative evaluation of scatter correction in 3D PET using different scatter-level approximations

OBJECTIVE

In 3D PET, scatter of the gamma photons is one of the most significant physical factors which degrades not only image quality but also quantification. The currently most used scatter estimation method is the analytic single scatter simulation (SSS) which usually accommodates for multiple scattering by scaling the single scatter estimation. However, it has not been clear yet how accurate this approximation is for cases where multiple scatter is significant, raising the question: "How important is correction for multiple scattered photons, and how accurately do we need to simulate all scattered events by appropriate scaling?" This study answers these questions and evaluates the accuracy of SSS implementation in the open-source library STIR.

METHODS

Different scatter orders approximations are evaluated including different levels of scattering and different scaling approaches using Monte Carlo (i.e. SimSET) data. SimSET simulations of a large anthropomorphic phantom were reconstructed with iterative reconstruction algorithms. Images reconstructed with 3D filtered back-projection reprojection algorithm have been compared quantitatively in order to clarify the errors due to different scatter order approximations.

RESULTS

Quantification in regions has improved by scatter correction. For example, in the heart the ideal value was 3, whereas before scatter correction the standard uptake value (SUV) was 4.0, after single scatter correction was 3.3 and after single and double scatter correction was 3.0. After correction by scaling single scatter with tail-fit, the SUV was 3.1, whereas with total-fit it was 3.0. Similarly, for the SSS correction methodology implemented in STIR using tail-fit the heart SUV was 3.1 whereas using total-fit it was 3.0.

CONCLUSIONS

The results demonstrate that correction for double scatter improves image contrast and therefore it is required for the accurate estimation of activity distribution in PET imaging. However, it has been also shown that scaling the single scatter distribution is a reasonable approximation to compensate for total scatter. Finally, scatter correction with STIR has shown excellent agreement with Monte Carlo simulations.

Novel scatter compensation of list-mode PET data using spatial and energy dependent corrections

With the widespread use of positron emission tomography (PET) crystals with greatly improved energy resolution (e.g., 11.5% with LYSO as compared to 20% with BGO) and of list-mode acquisitions, the use of the energy of individual events in scatter correction schemes becomes feasible. We propose a novel scatter approach that incorporates the energy of individual photons in the scatter correction and reconstruction of list-mode PET data in addition to the spatial information presently used in clinical scanners. First, we rewrite the Poisson likelihood function of list-mode PET data including the energy distributions of primary and scatter coincidences and show that this expression yields an MLEM reconstruction algorithm containing both energy and spatial dependent corrections. To estimate the spatial distribution of scatter coincidences we use the single scatter simulation (SSS). Next, we derive two new formulae which allow estimation of the 2-D (coincidences) energy probability density functions (E-PDF) of primary and scatter coincidences from the 1-D (photons) E-PDFs associated with each photon. We also describe an accurate and robust object-specific method for estimating these 1-D E-PDFs based on a decomposition of the total energy spectra detected across the scanner into primary and scattered components. Finally, we show that the energy information can be used to accurately normalize the scatter sinogram to the data. We compared the performance of this novel scatter correction incorporating both the position and energy of detected coincidences to that of the traditional approach modeling only the spatial distribution of scatter coincidences in 3-D Monte Carlo simulations of a medium cylindrical phantom and a large, nonuniform NCAT phantom. Incorporating the energy information in the scatter correction decreased bias in the activity distribution estimation by ~20% and ~40% in the cold regions of the large NCAT phantom at energy resolutions 11.5% and 20% at 511 keV, respectively, compared to when using the spatial information alone.

A scatter-compensated crystal interference factor in component-based normalization for high-resolution whole-body PET

On a positron emission tomography (PET) scanner consisting of block detectors, coincidence responses to scattered radiation may differ from those to true depending on the crystal pair position within a coincidence block pair. Furthermore, these differences are considered to vary according to the radial position of the coincidence block pair. These conditions create ringing artifacts in the reconstructed image due to the lack of scatter compensation in detector normalization. In component-based normalization, a scatter-compensated crystal interference factor is therefore required in addition to the scatter-compensated block profile and intrinsic crystal efficiencies. In this study, we propose a scatter-compensated component-based normalization scheme using an annulus phantom, which provides true and scattered radiations over a large transaxial field of view, and evaluates the quality of three different-sized phantom images with whole-body PET. The results showed that the proposed normalization method significantly reduces the ringing artifacts in reconstructed images with different scattered/true fractions. The proposed algorithm, which introduced the scatter-compensated crystal interference factor, worked well under different scattered/true ratio conditions and was considered to be a robust, practical normalization method in high-resolution whole-body PET.

Evaluation of dynamic row-action maximum likelihood algorithm reconstruction for quantitative 15O brain PET

OBJECTIVE

A modified version of row-action maximum likelihood algorithm (RAMLA) using a 'subset-dependent' relaxation parameter for noise suppression, or dynamic RAMLA (DRAMA), has been proposed. The aim of this study was to assess the capability of DRAMA reconstruction for quantitative (15)O brain positron emission tomography (PET).

METHODS

Seventeen healthy volunteers were studied using a 3D PET scanner. The PET study included 3 sequential PET scans for C(15)O, (15)O(2) and H (2) (15) O. First, the number of main iterations (N (it)) in DRAMA was optimized in relation to image convergence and statistical image noise. To estimate the statistical variance of reconstructed images on a pixel-by-pixel basis, a sinogram bootstrap method was applied using list-mode PET data. Once the optimal N (it) was determined, statistical image noise and quantitative parameters, i.e., cerebral blood flow (CBF), cerebral blood volume (CBV), cerebral metabolic rate of oxygen (CMRO(2)) and oxygen extraction fraction (OEF) were compared between DRAMA and conventional FBP. DRAMA images were post-filtered so that their spatial resolutions were matched with FBP images with a 6-mm FWHM Gaussian filter.

RESULTS

Based on the count recovery data, N (it) = 3 was determined as an optimal parameter for (15)O PET data. The sinogram bootstrap analysis revealed that DRAMA reconstruction resulted in less statistical noise, especially in a low-activity region compared to FBP. Agreement of quantitative values between FBP and DRAMA was excellent. For DRAMA images, average gray matter values of CBF, CBV, CMRO(2) and OEF were 46.1 +/- 4.5 (mL/100 mL/min), 3.35 +/- 0.40 (mL/100 mL), 3.42 +/- 0.35 (mL/100 mL/min) and 42.1 +/- 3.8 (%), respectively. These values were comparable to corresponding values with FBP images: 46.6 +/- 4.6 (mL/100 mL/min), 3.34 +/- 0.39 (mL/100 mL), 3.48 +/- 0.34 (mL/100 mL/min) and 42.4 +/- 3.8 (%), respectively.

CONCLUSION

DRAMA reconstruction is applicable to quantitative (15)O PET study and is superior to conventional FBP in terms of image quality.

Impact of Ge-68∕Ga-68-based versus CT-based attenuation correction on PET

Transmission (Tx) scans are used in PET for attenuation correction (AC). For standalone PET this is typically done using Ge-68/Ga-68 sources, for PET-CT using CT. Therefore, standalone PET suffers from emission contamination during Tx scans, PET-CT does not. Here, we studied the effects of AC across the two systems. With a cylindrical phantom (Jaszczak Phantom, Data Spectrum Corp.) with hollow spheres (diameter 10-60 mm) two studies were performed. In the first study the hollow spheres were filled with 150 kBq/ml FDG and the background with 15 kBq/ml. In the second study we used 120 kBq/ml in the spheres and 50 kBq/ml in the background. Both a low and a high object-to-background ratio are studied this way. Multiple scans were acquired on a standalone PET and a PET-CT until 1% of the initial concentration remained. Activity concentration in the spheres and background was measured from the reconstructed images and compared to the actual concentration. For standalone PET, emission scans were reconstructed using hot Tx (emission contaminated) and cold Tx (not contaminated). Uniformity within the spheres was investigated by profile analysis. For PET-CT, the concentration in the big spheres (> 16 mm) was recovered. For the smaller spheres, recovery was insufficient due to partial volume effects. For standalone PET the recoveries of the spheres (> 16 mm) were 20% (first study) and 13% (second study) lower than the actual concentration. Using hot Tx, underestimation of activity concentration was up to > 50%. Nonuniformities within the biggest spheres were up to 35%, 12%, and 5% (first study), using standalone PET with hot Tx, cold Tx, and using PET-CT, respectively. Due to contamination of AC by emission photons, standalone PET results in a bias in the activity concentration and uniformity. Especially when patients get follow-up PET scans on both standalone PET and PET-CT, this may lead to misinterpretation.

[Comparison of noise equivalent count rate and image quality for two-dimensional and three-dimensional PET scans]

The aim of this study was to investigate the correlation between noise equivalent count (NEC) rates and the signal-to-noise ratio (S/N) in reconstructed images. The NEC rates were determined using uniform 20 cm and 70 cm tall, 20 cm diameter cylinders filled with 11C. The phantoms were scanned in both two-dimensional and three-dimensional modes. The reconstructed image noise was evaluated using FBP and OSEM algorithms (4 iterations and 8 subsets). The images were filtered to a final image resolution of 6.5 mm. From the reconstructed image sets, averages and standard deviations of images were generated, from which the average image S/N (=average/standard deviation) was calculated within an 18 cm central ROI. The S/N of a central slice and an end slice was compared with the NEC. The NEC was found to have a linear relationship to the image S/N of all slices, depending on differences in noise properties specific to the reconstruction algorithm. In two-dimensional mode, although the image S/N of the central slice and the edge slice showed a linear relationship with the NEC, in three-dimensional mode, the S/N of the central slice did not show a relationship with the NEC. The linear relationship was also found in both two- and three-dimensional acquisition modes, as well as for the different activity distributions. These results indicate that the NEC is not only a measure for comparing the count rate performance of imaging systems. However, an absolute evaluation is impossible to depend on reconstruction algorithm, slice number, and phantom type.

Mapping the relative contribution of gray matter activity vs. volume in brain PET: a new approach

Interpretation of brain positron emission tomography (PET) in terms of function vs. structure is ambiguous owing to the partial volume effect (PVE). Therefore, observed differences in tracer distribution could reflect differences in either activity or volume, a problem that applies principally to gray matter (GM) since white matter (WM) virtually always has uniform activity. To assess the contribution of GM volume vs. activity, we implemented a method to directly compare PET images with underlying structure, and applied it to resting-state (18)Fluoro-deoxy-glucose-PET (FDG) of healthy subjects. Methods. Average GM and WM PVE-corrected mean FDG uptake values were applied onto co-registered segmented magnetic resonance imaging data sets to generate a "virtual PET" in which activity is proportional to GM volume and resolution set to that of PET. The raw PET and virtual PET values were then compared across the sample of subjects, first voxel-wise to detect clusters with significant activity-volume mismatch, and second within regions-of-interest (ROI) to quantify mismatches between unsmoothed voxel values. Results. Relative to volume, there was significant hyperactivity of most GM structures of the dorsal brain-except the thalamus-and significant hypoactivity of the temporal lobe, hippocampal region, and cerebellum, consistent across the voxel- and ROI-based analyses. Conclusion. As applied to normals, our method documented the expected contribution of functional activity independently of local differences in GM volume in the normal pattern of FDG uptake, and disclosed marked heterogeneities in functional activity per unit GM volume among structures. This generic method should find applications in pathological states as well as for other PET and SPECT radiotracers.

Scatter modelling and compensation in emission tomography

In nuclear medicine, clinical assessment and diagnosis are generally based on qualitative assessment of the distribution pattern of radiotracers used. In addition, emission tomography (SPECT and PET) imaging methods offer the possibility of quantitative assessment of tracer concentration in vivo to quantify relevant parameters in clinical and research settings, provided accurate correction for the physical degrading factors (e.g. attenuation, scatter, partial volume effects) hampering their quantitative accuracy are applied. This review addresses the problem of Compton scattering as the dominant photon interaction phenomenon in emission tomography and discusses its impact on both the quality of reconstructed clinical images and the accuracy of quantitative analysis. After a general introduction, there is a section in which scatter modelling in uniform and non-uniform media is described in detail. This is followed by an overview of scatter compensation techniques and evaluation strategies used for the assessment of these correction methods. In the process, emphasis is placed on the clinical impact of image degradation due to Compton scattering. This, in turn, stresses the need for implementation of more accurate algorithms in software supplied by scanner manufacturers, although the choice of a general-purpose algorithm or algorithms may be difficult.