EARL compliance measurements on the biograph vision Quadra PET/CT system with a long axial field of view

Background

Our aim was to determine sets of reconstruction parameters for the Biograph Vision Quadra (Siemens Healthineers) PET/CT system that result in quantitative images compliant with the European Association of Nuclear Medicine Research Ltd. (EARL) criteria. Using the Biograph Vision 600 (Siemens Healthineers) PET/CT technology but extending the axial field of view to 106 cm, gives the Vision Quadra currently an around fivefold higher sensitivity over the Vision 600 with otherwise comparable spatial resolution. Therefore, we also investigated how the number of incident positron decays—i.e., exposure—affects EARL compliance. This will allow estimating a minimal acquisition time or a minimal applied dose in clinical scans while retaining data comparability.

Methods

We measured activity recovery curves on a NEMA IEC body phantom filled with an aqueous ¹⁸F solution and a sphere to background ratio of 10–1 according to the latest EARL guidelines. Reconstructing 3570 image sets with varying OSEM PSF iterations, post-reconstruction Gaussian filter full width at half maximum (FWHM), and varying exposure from 59 kDecays/ml (= 3 s frame duration) to 59.2 MDecays/ml (= 1 h), allowed us to determine sets of parameters to achieve compliance with the current EARL 1 and EARL 2 standards. Recovery coefficients (RCs) were calculated for the metrics RC_max, RC_mean, and RC_peak, and the respective recovery curves were analyzed for monotonicity. The background’s coefficient of variation (COV) was also calculated.

Results

Using 6 iterations, 5 subsets and 7.8 mm Gauss filtering resulted in optimal EARL1 compliance and recovery curve monotonicity in all analyzed frames, except in the 3 s frames. Most robust EARL2 compliance and monotonicity were achieved with 2 iterations, 5 subsets, and 3.6 mm Gauss FWHM in frames with durations between 30 s and 10 min. RC_peak only impeded EARL2 compliance in the 10 s and 3 s frames.

Conclusions

While EARL1 compliance was robust over most exposure ranges, EARL2 compliance required exposures between 1.2 MDecays/ml to 11.5 MDecays/ml. The Biograph Vision Quadra’s high sensitivity makes frames as short as 10 s feasible for comparable quantitative images. Lowering EARL2 RC_max limits closer to unity would possibly even permit shorter frames.

Hybrid model of photon propagation based on the analytical and Monte Carlo methods for a dual-head PET system

The construction of photon propagation has a close relationship with the quality of reconstructed images. The classical Monte Carlo (MC) based method can model the photon propagation precisely, but it is time-consuming. The analytical method can often quickly construct a model, but its precision is a problem. How to fully exploit the advantages of the MC simulation and analytical model is an open problem. Inspired by the characteristics of the depth of interaction (DOI) detectors, which can help confirm the deposited position of a photon with DOI-encoding technology, we virtually discretize each crystal into several subcrystals to obtain the statistical distribution by MC-based simulation. Then, the statistical distribution is combined with a spatially variant solid-angle model. This combination strategy provides a hybrid model to describe photon propagation with relatively high accuracy and low computational cost. Three discretization schemes are compared to optimize the constructed photon propagation model. Several experiments are carried out to evaluate the performance of the proposed hybrid method. The metrics of full width at half maximum (FWHM), contrast recovery (CR), and coefficient of variation (COV) are adopted to quantitate the imaging results. The classical MC-based method is compared as a gold-standard reference. When a crystal is divided into two discretized positions, the convergent tendencies of CRs and COVs are consistent with that based on MC simulation method, respectively. In terms of FWHMs, the resolutions of using the MC-based model and the proposed hybrid model are 0.71 mm and 0.68 mm in the direction parallel to the detector head, respectively. This indicates the potential of the proposed method in positron emission tomography imaging.

Towards guidelines to harmonize textural features in PET: Haralick textural features vary with image noise, but exposure-invariant domains enable comparable PET radiomics

Purpose

Image texture is increasingly used to discriminate tissues and lesions in PET/CT. For quantification or in computer-aided diagnosis, textural feature analysis must produce robust and comparable values. Because statistical feature values depend on image count statistics, we investigated in depth the stability of Haralick features values as functions of acquisition duration, and for common image resolutions and reconstructions.

Methods

A homogeneous cylindrical phantom containing 9.6 kBq/ml Ge-68 was repeatedly imaged on a Siemens Biograph mCT, with acquisition durations ranging from three seconds to three hours. Images with 1.5, 2, and 4 mm isometrically spaced voxels were reconstructed with filtered back-projection (FBP), ordered subset expectation maximization (OSEM), and the Siemens TrueX algorithm. We analysed Haralick features derived from differently quantized (3 to 8-bit) grey level co-occurrence matrices (GLCMs) as functions of exposure E, which we defined as the product of activity concentration in a volume of interest (VOI) and acquisition duration. The VOI was a 50 mm wide cube at the centre of the phantom. Feature stability was defined for df/dE → 0.

Results

The most stable feature values occurred in low resolution FBPs, whereas some feature values from 1.5 mm TrueX reconstructions ranged over two orders of magnitude. Within the same reconstructions, most feature value-exposure curves reached stable plateaus at similar exposures, regardless of GLCM quantization. With 8-bit GLCM, median time to stability was 16 s and 22 s for FBPs, 18 s and 125 s for OSEM, and 23 s, 45 s, and 76 s for PSF reconstructions, with longer durations for higher resolutions. Stable exposures coincided in OSEM and TrueX reconstructions with image noise distributions converging to a Gaussian. In FBP, the occurrence of stable values coincided the disappearance of negatives image values in the VOI.

Conclusions

Haralick feature values depend strongly on exposure, but invariance exists within defined domains of exposure. Here, we present an easily replicable procedure to identify said stable exposure domains, where image noise does not substantially add to textural feature values. Only by imaging at predetermined feature-invariant exposure levels and by adjusting exposure to expected activity concentrations, can textural features have a quantitative use in PET/CT. The necessary exposure levels are attainable by modern PET/CT systems in clinical routine.

Technical Note: Transconvolution based equalization of positron energy effects for the use of 68 Ge/68 Ga phantoms in determining 18 F PET recovery

PURPOSE

Avoiding measurement variability from ¹⁸ F phantom preparation by using ⁶⁸ Ge/⁶⁸ Ga phantoms for the determination of ¹⁸ F recovery curves (RC) in clinical quality assurance measurements and for PET/CT site qualification in multicentre clinical trials.

METHODS

RCs were obtained from PET/CT measurements of seven differently sized phantom spheres filled either with ¹⁸ F or with ⁶⁸ Ga. RCs for the respective other isotope were then determined by two different methods: In the first method, images were convolved with positron range transconvolution functions derived from positron annihilation distributions found in literature. This method generated recasted images matching images using the respective other isotope. In the second method, the PET/CT system's isotope independent (intrinsic) point spread function was determined from said phantom measurements and convolved with numerical representations simulating hot spheres filled with the respective other isotope. These simulations included the isotope specific positron annihilation distributions. Recovered activity concentrations were compared between recasted images, simulated images, and the originally acquired images.

RESULTS

¹⁸ F and ⁶⁸ Ga recovery was successfully determined from image acquisitions of the respective opposite isotope as well as from the simulations. ⁶⁸ Ga RCs derived from ¹⁸ F data had a normalized root-mean-square deviation (NRMSD) from real ⁶⁸ Ga measurements of 0.019% when using the first method and of 0.008% when using the second method. ¹⁸ F RCs derived from ⁶⁸ Ga data had a NRMSD from real ¹⁸ F measurements of 0.036% when using the first method and of 0.038% when using the second method.

CONCLUSIONS

Applying the principles of transconvolution, ¹⁸ F RCs can be recalculated from ⁶⁸ Ga phantom measurements with excellent accuracy. The maximal additionally introduced error was below 0.4% of the error currently accepted for RCs in the site qualification of multicentre clinical trials by the EARL program of the European Association of Nuclear Medicine (EANM). Therefore, our methods legitimately allow for the use of long-lived solid state ⁶⁸ Ge/⁶⁸ Ga phantoms instead of manually prepared ¹⁸ F phantoms to characterize comparability of ¹⁸ F measurements across different imaging sites or of longitudinal ¹⁸ F measurements at a single PET/CT system.