Quantitative and Qualitative Improvement of Low-Count [68Ga]Citrate and [90Y]Microspheres PET Image Reconstructions Using Block Sequential Regularized Expectation Maximization Algorithm

Y-90 PET/MR imaging optimization with a Bayesian penalized likelihood reconstruction algorithm

Positron Emission Tomography (PET) imaging after90 Y liver radioembolization is used for both lesion identification and dosimetry. Bayesian penalized likelihood (BPL) reconstruction algorithms are an alternative to ordered subset expectation maximization (OSEM) with improved image quality and lesion detectability. The investigation of optimal parameters for90 Y image reconstruction of Q.Clear, a commercial BPL algorithm developed by General Electric (GE), in PET/MR is a field of interest and the subject of this study. The NEMA phantom was filled at an 8:1 sphere-to-background ratio. Acquisitions were performed on a PET/MR scanner for clinically relevant activities between 0.7 and 3.3 MBq/ml. Reconstructions with Q.Clear were performed varying the β penalty parameter between 20 and 6000, the acquisition time between 5 and 20 min and pixel size between 1.56 and 4.69 mm. OSEM reconstructions of 28 subsets with 2 and 4 iterations with and without Time-of-Flight (TOF) were compared to Q.Clear with β = 4000. Recovery coefficients (RC), their coefficient of variation (COV), background variability (BV), contrast-to-noise ratio (CNR) and residual activity in the cold insert were evaluated. Increasing β parameter lowered RC, COV and BV, while CNR was maximized at β = 4000; further increase resulted in oversmoothing. For quantification purposes, β = 1000-2000 could be more appropriate. Longer acquisition times resulted in larger CNR due to reduced image noise. Q.Clear reconstructions led to higher CNR than OSEM. A β of 4000 was obtained for optimal image quality, although lower values could be considered for quantification purposes. An optimal acquisition time of 15 min was proposed considering its clinical use.

Reduction of [68Ga]Ga-DOTA-TATE injected activity for digital PET/MR in comparison with analogue PET/CT

BACKGROUND

New digital detectors and block-sequential regularized expectation maximization (BSREM) reconstruction algorithm improve positron emission tomography (PET)/magnetic resonance (MR) image quality. The impact on image quality may differ from analogue PET/computed tomography (CT) protocol. The aim of this study is to determine the potential reduction of injected [68Ga]Ga-DOTA-TATE activity for digital PET/MR with BSREM reconstruction while maintaining at least equal image quality compared to the current analogue PET/CT protocol.

METHODS

NEMA IQ phantom data and 25 patients scheduled for a diagnostic PET/MR were included. According to our current protocol, 1.5 MBq [68Ga]Ga-DOTA-TATE per kilogram (kg) was injected. After 60 min, scans were acquired with 3 (≤ 70 kg) or 4 (> 70 kg) minutes per bedposition. PET/MR scans were reconstructed using BSREM and factors β 150, 300, 450 and 600. List mode data with reduced counts were reconstructed to simulate scans with 17%, 33%, 50% and 67% activity reduction. Image quality was measured quantitatively for PET/CT and PET/MR phantom and patient data. Experienced nuclear medicine physicians performed visual image quality scoring and lesion counting in the PET/MR patient data.

RESULTS

Phantom analysis resulted in a possible injected activity reduction of 50% with factor β = 600. Quantitative analysis of patient images revealed a possible injected activity reduction of 67% with factor β = 600. Both with equal or improved image quality as compared to PET/CT. However, based on visual scoring a maximum activity reduction of 33% with factor β = 450 was acceptable, which was further limited by lesion detectability analysis to an injected activity reduction of 17% with factor β = 450.

CONCLUSION

A digital [68Ga]Ga-DOTA-TATE PET/MR together with BSREM using factor β = 450 result in 17% injected activity reduction with quantitative values at least similar to analogue PET/CT, without compromising on PET/MR visual image quality and lesion detectability.

Quantitative Assessment of Myocardial Ischemia With Positron Emission Tomography

Recent advances in positron emission tomography (PET) technology and reconstruction techniques have now made quantitative assessment using cardiac PET readily available in most cardiac PET imaging centers. Multiple PET myocardial perfusion imaging (MPI) radiopharmaceuticals are available for quantitative examination of myocardial ischemia, with each having distinct convenience and accuracy profile. Important properties of these radiopharmaceuticals ( 15 O-water, 13 N-ammonia, 82 Rb, 11 C-acetate, and 18 F-flurpiridaz) including radionuclide half-life, mean positron range in tissue, and the relationship between kinetic parameters and myocardial blood flow (MBF) are presented. Absolute quantification of MBF requires PET MPI to be performed with protocols that allow the generation of dynamic multiframes of reconstructed data. Using a tissue compartment model, the rate constant that governs the rate of PET MPI radiopharmaceutical extraction from the blood plasma to myocardial tissue is calculated. Then, this rate constant ( K1 ) is converted to MBF using an established extraction formula for each radiopharmaceutical. As most of the modern PET scanners acquire the data only in list mode, techniques of processing the list-mode data into dynamic multiframes are also reviewed. Finally, the impact of modern PET technologies such as PET/CT, PET/MR, total-body PET, machine learning/deep learning on comprehensive and quantitative assessment of myocardial ischemia is briefly described in this review.

A novel figure of merit to investigate 68Ga PET/CT image quality based on patient weight and lesion size using Q.Clear reconstruction algorithm: A phantom study

INTRODUCTION

Q.Clear is a Bayesian penalised-likelihood algorithm that uses a β-value for positron emission tomography(PET)/computed tomography(CT) image reconstruction(IR). Our study proposes a novel figure of merit, named CRBV, to compare the Q.Clear performances using 68Ga PET/CT image with the ordered-subset-expectation-maximization(OSEM) algorithm and to identify the optimal β-values for these images using two phantoms mimicking normal and overweight patients.

METHODS

NEMA IQ phantom with or without a ring of water-filled plastic bags (NEMAstd and NEMAow, respectively) was acquired and reconstructed with OSEM and Q.Clear at various β-values and minutes/bed position(min/bp). Contrast recovery(CR), background variability(BV) and CRBV were calculated. Highest CRBV values were used to identify optimal β-value ranges.

RESULTS

Q.Clear with 250 ≤ β ≤ 800 improved CRBV compared to OSEM for all the investigated spheres and acquisition setups. Outside of this range, Q.Clear still outperformed OSEM with few exceptions depending on spheres diameters and phantoms(e.g.,β-value = 1600 for diameters ≤ 17 mm using the NEMAow phantom). Regarding the CRBV performance for IR optimization, for the 4 min/bp NEMAstd IR, β-values = 300 ÷ 350 allowed to simultaneously optimize all diameters(except for the 10 mm); for the NEMAow IR, β-values = 350 ÷ 500 were needed for diameters > 20 mm, while β-values = 200 ÷ 250 were selected for the remaining diameters. For the 2 min/bp, β-value = 500 was suitable for diameters > 17 mm in both NEMAstd and NEMAow IR, while for smaller diameters β-value = 200 and β-values = 250 ÷ 350 were obtained for NEMAstd and NEMAow, respectively.

CONCLUSION

Almost all tested β-values of Q.Clear improved the CRBV compared to OSEM. In both phantoms, simulating normal and over-weight patients, optimal β-values were found according to lesion sizes and investigated acquisition times.

Phantom-based evaluation of yttrium-90 datasets using biograph vision quadra

PURPOSE

The image quality characteristics of two NEMA phantoms with yttrium-90 (⁹⁰Y) were evaluated on a long axial field-of-view (AFOV) PET/CT. The purpose was to identify the optimized reconstruction setup for the imaging of patients with hepatocellular carcinoma after ⁹⁰Y radioembolization.

METHODS

Two NEMA phantoms were used, where one had a 1:10 sphere to background activity concentration ratio and the second had cold background. Reconstruction parameters used are as follows: iterations 2 to 8, Gaussian filter 2- to 6-mm full-width-at-half-maximum, reconstruction matrices 440 × 440 and 220 × 220, high sensitivity (HS), and ultra-high sensitivity (UHS) modes. 50-, 40-, 30-, 20-, 10-, and 5-min acquisitions were reconstructed. The measurements included recovery coefficients (RC), signal-to-noise ratio (SNR), background variability, and lung error which measures the residual error in the corrections. Patient data were reconstructed with 20-, 10-, 5-, and 1-min time frames and evaluated in terms of SNR.

RESULTS

The RC for the hot phantom was 0.36, 0.45, 0.53, 0.63, 0.68, and 0.84 for the spheres with diameters of 10, 13, 17, 22, 28, and 37 mm, respectively, for UHS 2 iterations, a 220 × 220 matrix, and 50-min acquisition. The RC values did not differ with acquisition times down to 20 min. The SNR was the highest for 2 iterations, measured 11.7, 16.6, 17.6, 19.4, 21.9, and 27.7 while the background variability was the lowest (27.59, 27.08, 27.36, 26.44, 30.11, and 33.51%). The lung error was 18%. For the patient dataset, the SNR was 19%, 20%, 24%, and 31% higher for 2 iterations compared to 4 iterations for 20-, 10-, 5-, and 1-min time frames, respectively.

CONCLUSIONS

This study evaluates the NEMA image quality of a long AFOV PET/CT scanner with ⁹⁰Y. It provides high RC for the smallest sphere compared to other standard AFOV scanners at shorter scan times. The maximum patient SNR was for 2 iterations, 20 min, while 5 min delivers images with acceptable SNR.

Positron Range Corrections and Denoising Techniques for Gallium-68 PET Imaging: A Literature Review

Gallium-68 (68Ga) is characterized by relatively high positron energy compared to Fluorine-18 (18F), causing substantial image quality degradation. Furthermore, the presence of statistical noise can further degrade image quality. The aim of this literature review is to identify the recently developed positron range correction techniques for 68Ga, as well as noise reduction methods to enhance the image quality of low count 68Ga PET imaging. The search engines PubMed and Scopus were employed, and we limited our research to published results from January 2010 until 1 August 2022. Positron range correction was achieved by using either deblurring or deep learning approaches. The proposed techniques improved the image quality and, in some cases, achieved an image quality comparable to 18F PET. However, none of these techniques was validated in clinical studies. PET denoising for 68Ga-labeled radiotracers was reported using either reconstruction-based techniques or deep learning approaches. It was demonstrated that both approaches can substantially enhance the image quality by reducing the noise levels of low count 68Ga PET imaging. The combination of 68Ga-specific positron range correction techniques and image denoising approaches may enable the application of low-count, high-quality 68Ga PET imaging in a clinical setting.

Efficient Strike Artifact Reduction Based on 3D-Morphological Structure Operators from Filtered Back-Projection PET Images

Positron emission tomography (PET) can provide functional images and identify abnormal metabolic regions of the whole-body to effectively detect tumor presence and distribution. The filtered back-projection (FBP) algorithm is one of the most common images reconstruction methods. However, it will generate strike artifacts on the reconstructed image and affect the clinical diagnosis of lesions. Past studies have shown reduction in strike artifacts and improvement in quality of images by two-dimensional morphological structure operators (2D-MSO). The morphological structure method merely processes the noise distribution of 2D space and never considers the noise distribution of 3D space. This study was designed to develop three-dimensional-morphological structure operators (3D MSO) for nuclear medicine imaging and effectively eliminating strike artifacts without reducing image quality. A parallel operation was also used to calculate the minimum background standard deviation of the images for three-dimensional morphological structure operators with the optimal response curve (3D-MSO/ORC). As a result of Jaszczak phantom and rat verification, 3D-MSO/ORC showed better denoising performance and image quality than the 2D-MSO method. Thus, 3D MSO/ORC with a 3 × 3 × 3 mask can reduce noise efficiently and provide stability in FBP images.

Effects of New Bayesian Penalized Likelihood Reconstruction Algorithm on Visualization and Quantification of Upper Abdominal Malignant Tumors in Clinical FDG PET/CT Examinations

Purpose

This study evaluated the effects of new Bayesian penalized likelihood (BPL) reconstruction algorithm on visualization and quantification of upper abdominal malignant tumors in clinical FDG PET/CT examinations, comparing the results to those obtained by an ordered subset expectation maximization (OSEM) reconstruction algorithm. Metabolic tumor volume (MTV) and texture features (TFs), as well as SUV-related metrics, were evaluated to clarify the BPL effects on quantification.

Materials and Methods

A total of 153 upper abdominal lesions (82 liver metastatic and 71 pancreatic cancers) were included in this study. FDG PET/CT images were acquired with a GE Discovery 710 scanner equipped with a time-of-flight system. Images were reconstructed using OSEM and BPL (beta 700) algorithms. In 58 lesions <1.5 cm in greatest diameter (small-lesion group), visual image quality of each lesion was evaluated using a four-point scale. SUVmax was obtained for quantitative metrics. Visual scores and SUVmax were compared between OSEM and BPL images. In 95 lesions >2.0 cm in greatest diameter (larger-lesion group), SUVmax, SUVpeak, MTV, and six TFs were compared between OSEM and BPL images. In addition to the size-based analyses, an increase of SUVmax with BPL was evaluated according to the original SUVmax in OSEM images.

Results

In the small-lesion group, both visual score and SUVmax were significantly higher in the BPL than OSEM images. The increase in visual score was observed in 20 (34%) of all 58 lesions. In the larger-lesion group, no statistical difference was observed in SUVmax, SUVpeak, or MTV between OSEM and BPL images. BPL increased high gray-level zone emphasis and decreased low gray-level zone emphasis among six TFs compared to OSEM with statistical significance. No statistical differences were observed in other TFs. SUVmax-based analysis demonstrated that BPL increased and decreased SUVmax in lesions with low (<5) and high (>10) SUVmax in original OSEM images, respectively.

Conclusion

This study demonstrated that BPL improved conspicuity of small or low-count upper abdominal malignant lesions in clinical FDG PET/CT examinations. Only two TFs represented significant differences between OSEM and BPL images of all quantitative metrics in larger lesions.

Yttrium-90 quantitative phantom study using digital photon counting PET

BACKGROUND

PET imaging of 90Y-microsphere distribution following radioembolisation is challenging due to the count-starved statistics from the low branching ratio of e+/e- pair production during 90Y decay. PET systems using silicon photo-multipliers have shown better 90Y image quality compared to conventional photo-multiplier tubes. The main goal of the present study was to evaluate reconstruction parameters for different phantom configurations and varying listmode acquisition lengths to improve quantitative accuracy in 90Y dosimetry, using digital photon counting PET/CT.

METHODS

Quantitative PET and dosimetry accuracy were evaluated using two uniform cylindrical phantoms specific for PET calibration validation. A third body phantom with a 9:1 hot sphere-to-background ratio was scanned at different activity concentrations of 90Y. Reconstructions were performed using OSEM algorithm with varying parameters. Time-of-flight and point-spread function modellings were included in all reconstructions. Absorbed dose calculations were carried out using voxel S-values convolution and were compared to reference Monte Carlo simulations. Dose-volume histograms and root-mean-square deviations were used to evaluate reconstruction parameter sets. Using listmode data, phantom and patient datasets were rebinned into various lengths of time to assess the influence of count statistics on the calculation of absorbed dose. Comparisons between the local energy deposition method and the absorbed dose calculations were performed.

RESULTS

Using a 2-mm full width at half maximum post-reconstruction Gaussian filter, the dosimetric accuracy was found to be similar to that found with no filter applied but also reduced noise. Larger filter sizes should not be used. An acquisition length of more than 10 min/bed reduces image noise but has no significant impact in the quantification of phantom or patient data for the digital photon counting PET. 3 iterations with 10 subsets were found suitable for large spheres whereas 1 iteration with 30 subsets could improve dosimetry for smaller spheres.

CONCLUSION

The best choice of the combination of iterations and subsets depends on the size of the spheres. However, one should be careful on this choice, depending on the imaging conditions and setup. This study can be useful in this choice for future studies for more accurate 90Y post-dosimetry using a digital photon counting PET/CT.

Impact of image reconstruction method on dose distributions derived from 90Y PET images: phantom and liver radioembolization patient studies

PET images acquired after liver ⁹⁰Y radioembolization therapies are typically very noisy, which significantly challenges both visualization and quantification of activity distributions. To improve their noise characteristics, regularized iterative reconstruction algorithms such as block sequential regularized expectation maximization (Q.Clear for GE Healthcare, USA) have been proposed. In this study, we aimed to investigate the effects which different reconstruction algorithms may have on patient images, with reconstruction parameters initially narrowed down using phantom studies. Moreover, we evaluated the impact of these reconstruction methods on voxel-based dose distribution in phantom and patient studies (lesions and healthy livers). The International Electrotechnical Commission (IEC)/NEMA phantom, containing six spheres, was filled with ⁹⁰Y and imaged using a GE Discovery 690 PET/CT scanner with time-of-flight enabled. The images were reconstructed using Q.Clear (with β parameter ranging from 0 to 8000) and ordered subsets expectation maximization. The image quality and quantification accuracy were evaluated by computing the hot ([Formula: see text]) and cold ([Formula: see text]) contrast recovery coefficients, background variability (BV) and activity bias. Next, dose distributions and dose volume histograms were generated using MIM® software's SurePlan LiverY90 toolbox. Subsequently, parameters optimized in these phantom studies were applied to five patient datasets. Dose parameters, such as D_max, D_mean, D₇₀, and V_100Gy, were estimated, and their variability for different reconstruction methods was investigated. Based on phantom studies, the β parameter values optimized for image quality and quantification accuracy were 2500 and 300, respectively. When all investigated reconstructions were applied to patient studies, D_mean, D₅₀, D₇₀, and V_100Gy showed coefficients of variation below 8%; whereas the variability of D_max was up to 30% for both phantom and patient images. Although β = 300-1000 would provide accurate activity quantification for a region of interest, when considering activity/dose voxelized distribution, higher β value (e.g. 4000-5000) would provide the greatest accuracy for dose distributions. In this ⁹⁰Y radioembolization PET/CT study, the β parameter in regularized iterative (Q.Clear) reconstruction was investigated for image quality, accurate quantification and dose distributions based on phantom experiments and then applied to patient studies. Our results indicate that more accurate dose distribution can be achieved from smoother PET images, reconstructed with larger β values than those yielding the best activity quantifications but noisy images. Most importantly, these results suggest that quantitative measures, which are commonly used in clinics, such as SUV_max or SUV_peak( equivalent of D_max), should not be employed for ⁹⁰Y PET images, since their values would highly depend on the image reconstruction.