Ordered Subset Expectation Maximization Research Articles

Assessing the impact of different penalty factors of the Bayesian reconstruction algorithm Q.Clear on in vivo low count kinetic analysis of [11C]PHNO brain PET-MR studies

IntroductionQ.Clear is a Bayesian penalised likelihood (BPL) reconstruction algorithm available on General Electric (GE) Positron Emission Tomography (PET)-Computed Tomography (CT) and PET-Magnetic Resonance (MR) scanners. This algorithm is regulated by a β value which acts as a noise penalisation factor and yields improvements in signal to noise ratio (SNR) in clinical scans, and in contrast recovery and spatial resolution in phantom studies. However, its performance in human brain imaging studies remains to be evaluated in depth. This pilot study aims to investigate the impact of Q.Clear reconstruction methods using different β value versus ordered subset expectation maximization (OSEM) on brain kinetic modelling analysis of low count brain images acquired in the PET-MR.MethodsSix [11C]PHNO PET-MR brain datasets were reconstructed with Q.Clear with β100–1000 (in increments of 100) and OSEM. The binding potential relative to non-displaceable volume (BPND) were obtained for the Substantia Nigra (SN), Striatum (St), Globus Pallidus (GP), Thalamus (Th), Caudate (Cd) and Putamen (Pt), using the MIAKAT™ software. Intraclass correlation coefficients (ICC), repeatability coefficients (RC), coefficients of variation (CV) and bias from Bland–Altman plots were reported. Statistical analysis was conducted using a 2-way ANOVA model with correction for multiple comparisons.ResultsWhen comparing a standard OSEM reconstruction of 6 iterations/16 subsets and 5 mm filter with Q.Clear with different β values under low counts, the bias and RC were lower for Q.Clear with β100 for the SN (RC = 2.17), Th (RC = 0.08) and GP (RC = 0.22) and with β200 for the St (RC = 0.14), Cd (RC = 0.18)and Pt (RC = 0.10). The p-values in the 2-way ANOVA model corroborate these findings. ICC values obtained for Th, St, GP, Pt and Cd demonstrate good reliability (0.87, 0.99, 0.96, 0.99 and 0.96, respectively). For the SN, ICC values demonstrate poor reliability (0.43).ConclusionBPND results obtained from quantitative low count brain PET studies using [11C]PHNO and reconstructed with Q.Clear with β < 400, which is the value used for clinical [18F]FDG whole-body studies, demonstrate the lowest bias versus the typical iterative reconstruction method OSEM.

Initial clinical experience with dedicated multi-pinhole (MPH) collimator for 99mTc-HMPAO brain perfusion SPECT.

Dedicated multi-pinhole (MPH) collimators have been successfully tested in selected clinical investigations. The aim of our work was to report initial experiences with an MPH collimator set designed for brain perfusion single photon emission tomography (SPECT). Ten patients underwent sequential technetium-99m-hexamethylpropyleneamineoxime (99mTc-HMPAO) SPECT with a dual-head SPECT camera equipped with conventional low-energy parallel hole collimators (LEHR), and with a triple-head system equipped with MPH collimators. Low-energy parallel hole collimators data were reconstructed by filtered back projection (FBP), ordered subset expectation maximization (OSEM), software for tomographic image reconstruction (STIR). In addition, both the parallel hole data and MPH data were reconstructed by Tera-TomoTM 3D iterative reconstruction denoted LEHR_TT3D and MPH_TT3D, respectively. Five medical experts visually compared the reconstructed images of the five data sets and defined a ranking sequence from the lowest (1) to the highest (5) image quality. Results were compared using the Friedman test. P values below 0.05 were considered significant. Low-energy parallel hole collimators acquisition resulted in 5 million, while MPH acquisition in 13 million total counts with 30 and 34 minutes of acquisition time, respectively. Mean rank coefficients of the reconstruction methods were 1.96±0.52, 2.66±0.46, 2.86±0.60, 3.62±0.55, 3.9±0.68 for FBP, STIR, LEHR_TT3D, LEHR_OSEM, MPH_TT3D respectively. The differences between MPH_TT3D-FBP (P<0.01); MPH_TT3D-STIR (P<0.05); LEHR_OSEM-FBP (P<0.01) were significant. Image quality provided by MPH collimator is comparable to that provided by conventional LEHR imaging. Higher sensitivity has the potential to shorten acquisition time or to reduce the amount of administered activity.

Evaluation of the spline reconstruction technique for preclinical PET imaging

Background and ObjectiveThe Spline Reconstruction Technique (SRT) is a fast algorithm based on a novel numerical implementation of an analytic representation of the inverse Radon transform. The purpose of this study is to provide a comparison between SRT, Filtered Back-Projection (FBP), Ordered Subset Expectation Maximization 2D (2D-OSEM), and the Tera-Tomo 3D algorithm, using phantom data at various acquisition durations as well as small-animal data obtained from the Mediso nanoScan® PET/CT scanner. MethodsFor this purpose, the “NEMA NU 4–2008 standards” protocol was employed at five different realizations and acquisition durations. In addition to the image quality metrics described by the NEMA protocol, Cold Region Contrast was also considered as a figure-of-merit. Furthermore, Cold Region Contrast was measured in the myocardial infarction region of six male Wistar rats. The volumetric defect quantification was assessed with dedicated computer software. Lastly, plots of Recovery Coefficient and Spill-Over Ratio as a function of the Percentage Standard Deviation were generated, after smoothing the phantom reconstructions with four different Gaussian filters. Statistical significance was determined by employing the Kruskal-Wallis test or One-way Analysis of Variance depending on the normality of the variable's distribution. ResultsThe present study revealed that, at the expense of slightly increased noise in the reconstructed images, SRT resulted in higher Recovery Coefficient values for small hot regions of interest, when compared with FBP and 2D-OSEM at all acquisition durations. Furthermore, SRT reconstructed images exhibit higher Recovery Coefficient values, for all hot regions of interest, when compared to the other 2D algorithms at short acquisition durations. In both phantom and animal studies, SRT achieved a significant improvement over 2D-OSEM for the Spill-Over Ratio and the Cold Region Contrast. These advantages were maintained even after comparing the algorithms at equal noise levels. The Tera-Tomo 3D algorithm (4 subsets, iterations≥ 13) performed significantly better compared to the other algorithms for all figures-of-merit. No statistically significant differences regarding the myocardial defect size were observed between the algorithms investigated. ConclusionsOverall, SRT appears that could be useful for the quantification of small hot regions of interest, cold regions of interest, as well as in low-count imaging applications.

Assessment of different quantification metrics of [18F]-NaF PET/CT images of patients with abdominal aortic aneurysm

BackgroundWe aim to assess the spill-in effect and the benefit in quantitative accuracy for [18F]-NaF PET/CT imaging of abdominal aortic aneurysms (AAA) using the background correction (BC) technique. MethodsSeventy-two datasets of patients diagnosed with AAA were reconstructed with ordered subset expectation maximization algorithm incorporating point spread function (PSF). Spill-in effect was investigated for the entire aneurysm (AAA), and part of the aneurysm excluding the region close to the bone (AAAexc). Quantifications of PSF and PSF+BC images using different thresholds (% of max. SUV in target regions-of-interest) to derive target-to-background (TBR) values (TBRmax, TBR90, TBR70 and TBR50) were compared at 3 and 10 iterations. ResultsTBR differences were observed between AAA and AAAexc due to spill-in effect from the bone into the aneurysm. TBRmax showed the highest sensitivity to the spill-in effect while TBR50 showed the least. The spill-in effect was reduced at 10 iterations compared to 3 iterations, but at the expense of reduced contrast-to-noise ratio (CNR). TBR50 yielded the best trade-off between increased CNR and reduced spill-in effect. PSF+BC method reduced TBR sensitivity to spill-in effect, especially at 3 iterations, compared to PSF (P-value ≤ 0.05). ConclusionTBR50 is robust metric for reduced spill-in and increased CNR.

Optimization of image reconstruction method of cerebral blood flow perfusion imaging with digital CZT SPECT.

The purpose of this study is to evaluate the effects of filtered back projection (FBP), ordered subset expectation maximisation (OSEM), and different filters on cadmium zinc telluride single-photon emission computed tomography [CZT single-photon emission computed tomography (SPECT)] cerebral blood perfusion image quality to optimise the image reconstruction method. Under routine clinical conditions, tomographic imaging was performed on the phantom and patients. Image processing included image reconstruction using FBP and OSEM, and the filtering method used Butterworth (Bw) and Gaussian (Gs) filters. Visual and semi-quantitative parameters [integral uniformity, root mean square (RMS) noise and contrast and contrast-to-noise ratio (CNR)] were used to evaluate image quality to optimise image reconstruction parameters. One-way and two-way analysis of variance were used to process phantom and clinical data. In the tomographic images of the phantom, the semi-quantitative analysis showed that the integral uniformity of FBP+Bw was better than that of OSEM+Bw and OSEM+Gs (P < 0.05), and that the RMS noise of FBP+Bw was lower than that of OSEM+Bw and OSEM+Gs (P < 0.001). The contrast of FBP+Bw and OSEM+Bw in the cold area diameter ≥2 cm group was higher than that of OSEM+Gs (P < 0.001), whereas the CNR of FBP+Bw was higher than that of OSEM+Bw and OSEM+Gs (P < 0.001); the contrast of OSEM+Bw cold area diameter <2 cm was higher than that of FBP+Bw (P < 0.01). The semi-quantitative analysis results of the clinical images were consistent with the phantom's. In CZT SPECT cerebral blood flow perfusion imaging, it is suggested that the image postprocessing method of FBP+Bw (fc = 0.40; n = 10) should be used routinely in clinical application, and if there are uncertain small lesions in the processed image, it is suggested to use the reconstruction method of OSEM+Bw (EM-equivalent iterations = 60; fc = 0.45; n = 10) instead.

Optimization of 99m Tc whole-body SPECT/CT image quality: A phantom study.

PurposeInvestigate the impact of acquisition time and reconstruction parameters on single‐photon emission computed tomography/computed tomography (SPECT/CT) image quality with the ultimate aim of finding the shortest possible acquisition time for clinical whole‐body SPECT/CT (WB‐SPECT/CT) while maintaining image qualityMethodsThe National Electrical Manufacturers Association (NEMA) image quality measurements were performed on a SPECT/CT imaging system using a NEMA International Electrotechnical Commission (IEC) phantom with spherical inserts of varying diameter (10–37 mm), filled with 99mTc in activity sphere‐to‐background concentration ratio of 8.5:1. A gated acquisition was acquired and binned data were summed to simulate acquisitions of 15, 8, and 3 s per projection angle. Images were reconstructed on a Hermes (HERMES Medical Solutions AB, Stockholm, Sweden) workstation using eight subsets and between 4 and 24 iterations of the three‐dimensional (3D) ordered subset expectation maximization (OSEM) algorithm. Reconstructed images were post‐smoothed with 3D Gaussian filter ranging from 0 to 12 mm full‐width at half maximum (FWHM). Contrast recovery, background variability, and contrast‐to‐noise ratio were evaluatedResultsAs expected, the spheres were more clearly defined as acquisition time and count statistics improved. The optimal iteration number and Gaussian filter were determined from the contrast recovery convergence and level of noise. Convergence of contrast recovery was observed at eight iterations while 12 iterations yielded stabilized values at all acquisition times. In addition, it was observed that applying 3D Gaussian filter of 8–12 mm FWHM suppressed the noise and mitigated Gibbs artifacts. Background variability was larger for small spheres than larger spheres and the noise decreased when acquisition time became longer. A contrast‐to‐noise ratio >5 was reached for the two smallest spheres of 10 and 13 mm at acquisition times of 8 sConclusionOptimized reconstruction parameters preserved image quality with reduce acquisition time in present study. This study suggests an optimal protocol for clinical 99mTc SPECT/CT can be reached at 8 s per projection angle, with data reconstructed using 12 iterations and eight subset of the 3D OSEM algorithm and 8 mm Gaussian post‐filter.

Quantitative evaluation of PSMA PET imaging using a realistic anthropomorphic phantom and shell-less radioactive epoxy lesions

BackgroundPositron emission tomography (PET) with prostate specific membrane antigen (PSMA) have shown superior performance in detecting metastatic prostate cancers. Relative to [18F]fluorodeoxyglucose ([18F]FDG) PET images, PSMA PET images tend to visualize significantly higher-contrast focal lesions. We aim to evaluate segmentation and reconstruction algorithms in this emerging context. Specifically, Bayesian or maximum a posteriori (MAP) image reconstruction, compared to standard ordered subsets expectation maximization (OSEM) reconstruction, has received significant interest for its potential to reach convergence with minimal noise amplifications. However, few phantom studies have evaluated the quantitative accuracy of such reconstructions for high contrast, small lesions (sub-10 mm) that are typically observed in PSMA images. In this study, we cast 3 mm–16-mm spheres using epoxy resin infused with a long half-life positron emitter (sodium-22; 22Na) to simulate prostate cancer metastasis. The anthropomorphic Probe-IQ phantom, which features a liver, bladder, lungs, and ureters, was used to model relevant anatomy. Dynamic PET acquisitions were acquired and images were reconstructed with OSEM (varying subsets and iterations) and BSREM (varying β parameters), and the effects on lesion quantitation were evaluated.ResultsThe 22Na lesions were scanned against an aqueous solution containing fluorine-18 (18F) as the background. Regions-of-interest were drawn with MIM Software using 40% fixed threshold (40% FT) and a gradient segmentation algorithm (MIM’s PET Edge+). Recovery coefficients (RCs) (max, mean, peak, and newly defined “apex”), metabolic tumour volume (MTV), and total tumour uptake (TTU) were calculated for each sphere. SUVpeak and SUVapex had the most consistent RCs for different lesion-to-background ratios and reconstruction parameters. The gradient-based segmentation algorithm was more accurate than 40% FT for determining MTV and TTU, particularly for lesions le 6 mm in diameter (R2 = 0.979–0.996 vs. R2 = 0.115–0.527, respectively).ConclusionAn anthropomorphic phantom was used to evaluate quantitation for PSMA PET imaging of metastatic prostate cancer lesions. BSREM with β = 200–400 and OSEM with 2–5 iterations resulted in the most accurate and robust measurements of SUVmean, MTV, and TTU for imaging conditions in 18F-PSMA PET/CT images. SUVapex, a hybrid metric of SUVmax and SUVpeak, was proposed for robust, accurate, and segmentation-free quantitation of lesions for PSMA PET.

A Novel Method to Suppress the Effect of Subdiaphragmatic Activity in 99mTc Myocardial Perfusion SPECT and Evaluation of Its Usefulness Using a Myocardial Phantom.

Background: Smoothing in 99mTc myocardial perfusion single-photon emission computed tomography (SPECT) often increases myocardial artifacts due to subdiaphragmatic activity near the heart. To reduce these artifacts, we developed a new process flow, masking on unsmoothed images (MUS), that includes the extraction of the myocardium by masking before smoothing. Methods: This study evaluated the relationships between matrix sizes and distances to the subdiaphragmatic activity using the MUS method compared to conventional methods using a combination of image reconstruction methods (filtered back-projection [FBP] and ordered subset expectation maximization [OSEM]) with or without corrections (attenuation [AC], scatter [SC], and resolution recovery [RR]) using a myocardial phantom. The results were compared for two matrix sizes (pixel sizes) (128×128 [3.3 mm] and 64×64 [6.6 mm]); four subdiaphragmatic activity distances (5, 10, 15, and 20 mm); and three reconstruction methods (FBP without correction; OSEM with RR; and OSEM with AC, SC, and RR). Results: In the conventional method, increasing distance resulted in interference with myocardial perfusion SPECT evaluation however, the artifacts were less apparent when the MUS method was applied. The images converted to 64×64 did not show the same effect as the 128×128 images, even when RR was used. The MUS method was useful for acquisition at 128×128, along with the use of RR in the reconstruction process. Conclusion: MUS mitigated the effects of subdiaphragmatic activity on myocardial perfusion SPECT, particularly combined with 128×128 acquisitions and iterative reconstruction with RR.

Fast and Accurate Thoracic SPECT Image Reconstruction

In Single-Photon Emission Computed Tomography (SPECT), the reconstructed image has insufficient contrast, poor resolution and inaccurate volume of the tumor size due to physical degradation factors. Generally, nonstationary filtering of the projection or the slice is one of the strategies for correcting the resolution and therefore improving the quality of the reconstructed SPECT images. This paper presents a new 3D algorithm that enhances the quality of reconstructed thoracic SPECT images and reduces the noise level with the best degree of accuracy. The suggested algorithm is composed of three steps. The first one consists of denoising the acquired projections using the benefits of the complementary properties of both the Curvelet transform and the Wavelet transforms to provide the best noise reduction. The second step is a simultaneous reconstruction of the axial slices using the 3D Ordered Subset Expectation Maximization (OSEM) algorithm. The last step is post-processing of the reconstructed axial slices using one of the newest anisotropic diffusion models named Partial Differential Equation (PDE). The method is tested on two digital phantoms and clinical bone SPECT images. A comparative study with four algorithms reviewed on state of the art proves the significance of the proposed method. In simulated data, experimental results show that the plot profile of the proposed model keeps close to the original one compared to the other algorithms. Furthermore, it presents a notable gain in terms of contrast to noise ratio (CNR) and execution time. The proposed model shows better results in the computation of contrast metric with a value of 0.68 ± 7.2 and the highest signal to noise ratio (SNR) with a value of 78.56 ± 6.4 in real data. The experimental results prove that the proposed algorithm is more accurate and robust in reconstructing SPECT images than the other algorithms. It could be considered a valuable candidate to correct the resolution of bone in the SPECT images.

Impact of Reduced Acquisition Time on Bone Single-photon Emission Computed Tomography Images in Oncology Patients.

Background:The use of resolution recovery (RR) in bone and myocardial perfusion imaging is becoming increasingly popular in nuclear medicine departments. RR produces reconstructed images that show improved spatial resolution and signal-to-noise ratio compared with conventional single-photon emission computed tomography (SPECT) images.Objective:To evaluate the impact of the ordered subset expectation maximization (OSEM) RR modality on preserving noise, signal-to-noise ratio (SNR), and contrast-to-noise ratio (CNR) for short SPECT acquisition.Methods:This prospective study was conducted on 80 patients. Full SPECT acquisition was performed as a standardized protocol, while reduced acquisition was achieved with the Poisson resampling method. Noise, SNR, and CNR were measured for different reconstruction parameters for the same image levels. The impact of surface area and body mass index was also measured for the same reconstruction parameters.Results:The results show significantly higher SNR and CNR for the Evolution for Bone protocol compared to the other two reconstruction protocols for full and reduced SPECT acquisition. With the shortening of the SPECT acquisition, an increase in the value of noise was recorded. SNR and CNR decreased with the reduction in SPECT acquisition.Conclusion:The Evolution for Bone protocol for all three analyzed acquisition protocols had the lowest noise values. The highest SNR and CNR were recorded in the Evolution for Bone protocol for the three acquisition protocols and SPECT acquisition time can be reduced from 20 to 10 min for bone SPECT.

Effect of Bayesian penalty likelihood algorithm on 18F-FDG PET/CT image of lymphoma.

ObjectiveRecently, a new Bayesian penalty likelihood (BPL) reconstruction algorithm has been applied in PET, which is expected to provide better image resolution than the widely used ordered subset expectation maximization (OSEM). The purpose of this study is to compare the differences between these two algorithms in terms of image quality and effects on clinical diagnostics and quantification of lymphoma.MethodsA total of 246 FDG-positive lesions in 70 patients with lymphoma were retrospectively analyzed by using BPL and OSEM + time-of-flight + point spread function algorithms. Visual analysis was used to evaluate the effects of different reconstruction algorithms on clinical image quality and diagnostic certainty. Quantitative analysis was used to compare the differences between pathology and lesion size.ResultsThere were significant differences in lesion-related SUVmax, total-lesion-glycolysis (TLG), and signal-to-background ratio (SBR) (P < 0.01). The variation Δ SUVmax% and Δ SBR% caused by the two reconstruction algorithms were negatively correlated with tumor diameter, while Δ MTV% and Δ TLG% were positively correlated with tumor diameter. In the grouped analysis based on pathology, there were significant differences in lesion SUVmax, lesion SUVmean, and SBR. In non-Hodgkin’s lymphoma (diffuse large B cells and follicular lymphoma), diversities were significantly found in SUVmax, SUVmean, SBR, and TLG of the lesions (P < 0.05). According to the grouped analysis based on lesion size, for lesions smaller than 1 cm and 2 cm, there was a significant difference in SUVmean, SUVmax, SBR, and MTV, but not in lesions larger than or equal to 2 cm (P > 0.05), and the liver background SUVmean (P > 0.05) remained unchanged.ConclusionBPL reconstruction algorithm could effectively improve clinical image quality and diagnostic certainty. In quantitative analysis, there were no significant differences among different pathological groups, but there were significant diversities in lesion sizes. Especially for small lesions, lesion SUVmax increased and SBR was significantly improved, which may better assist in the diagnosis of small lesions of lymphoma.

Data-driven respiratory signal estimation from temporally finely sampled projection data in conventional cardiac perfusion SPECT imaging.

The aim of this work was to revisit the data-driven approach of axial center-of-mass (COM) measurements to recover a surrogate respiratory signal from finely sampled (100ms) single photon emission computed tomography (SPECT) projection data derived from list-mode acquisitions. For our initial evaluation, we acquired list-mode projection data from an anthropomorphic cardiac phantom mounted on a Quasar respiratory motion platform simulating 15mm amplitude respiratory motion. We also selected 302 consecutive patients (138 males, 164 females) with list-mode acquisitions, external respiratory motion tracking, and written consent to evaluate the clinical efficacy of our data-driven approach. Linear regression, Pearson's correlation coefficient (r), and standard error of the estimates (SEE) between the respiratory signals obtained with a visual tracking system (VTS) and COM measurements were calculated for individual projection data sets and for the patient group as a whole. Both the VTS- and COM-derived respiratory signals were used to estimate and correct respiratory motion. The reconstruction for six-degree of freedom rigid-body motion estimation was done in two ways: (1) using three iterations of ordered-subsets expectation-maximization (OSEM) with four subsets (16 projection angles per subset), or 12 iterations of maximum-likelihood expectation-maximization (MLEM). Respiratory motion compensation was done employing either OSEM with 16 subsets (four projection angles per subset) and five iterations or MLEM and 80 iterations, using the two respiratory estimates, respectively. Polar map quantification was also performed, calculating the percentage count difference (%Diff) between polar maps without and with respiratory motion included. Average % Diff was calculated in 17 segments (defined according to ASNC Guidelines). Paired t-tests were used to determine significance (p-values). The r-value calculated when comparing the VTS and COM respiratory signals varied widely between -0.01 and 0.96 with an average of 0.70, while the SEE varied between 0.80 and 6.48mm with an average of 2.05mm for our patient set, while the same values for the one anthropomorphic phantom acquisition are 0.91 and 1.11mm, respectively. A comparison between the respiratory motion estimates for VTS and COM in the S-I direction yielded an r=0.90 (0.94), and an SEE of 1.56mm (1.20mm) for OSEM (MLEM), respectively. Bland-Altman plots and calculated intraclass correlation coefficients also showed excellent agreement between the VTS and COM respiratory motion estimates. Average S-I respiratory estimates for the VTS (COM) were 9.04 (9.2mm) and 9.01mm (9.14mm) for the OSEM and MLEM, respectively. The paired t-test approached significance when comparing VTS and COM estimated respiratory signals with p-values of 0.069 and 0.051 for OSEM and MLEM. The respiratory estimates from the anthropomorphic cardiac phantom experiment using the VTS (COM) were 12.62 (14.10mm) and 12.55mm (14.29mm) for OSEM and MLEM, respectively. Polar map quantification yielded average % Diff consistently better when employing VTS-derived respiratory estimates to correct for respiration compared to the COM-derived estimates. The results indicate that our COM method has the potential to provide an automated data-driven correction of cardiac respiratory motion without the drawbacks of our VTS methodology. However, it is not generally equivalent to the VTS method in extent of correction.

68Ga-PSMA11 PET/CT for biochemically recurrent prostate cancer: Influence of dual-time and PMT- vs SiPM-based detectors

Objectives68Ga-PSMA11 PET/CT is excellent for evaluating biochemically recurrent prostate cancer (BCR PC). Here, we compared the positivity rates of dual-time point imaging using a PET/CT scanner (DMI) with silicon photomultiplier (SiPM) detectors and a PET/CT scanner (D690) with photomultiplier tubes (PMT), in patients with BCR PC. MethodsFifty-eight patients were prospectively recruited and randomized to receive scans on DMI followed by D690 or vice-versa. Images from DMI were reconstructed using the block sequential regularized expectation maximization (BSREM) algorithm and images from D690 were reconstructed using ordered subset expectation maximization (OSEM), according to the vendor's recommendations. Two readers independently reviewed all images in randomized order, recorded the number and location of lesions, as well as standardized uptake value (SUV) measurements. ResultsTwenty-eight patients (group A) had DMI as first scanner followed by D690, while 30 patients (group B) underwent scans in reversed order. Mean PSA was 30±112.9 (range 0.3–600.66) ng/mL for group A and 41.5 ± 213.2 (range 0.21–1170) ng/mL for group B (P = 0.796). The positivity rate in group A was 78.6% (22/28 patients) vs. 73.3% (22/30 patients) in group B. Although the performance of the two scanners was equivalent on a per-patient basis, DMI identified 5 additional sites of suspected recurrent disease when used as first scanner. The second scan time point did not reveal additional abnormal uptake. ConclusionsThe delayed time point in 68Ga-PSMA11 PET/CT did not show a higher positivity rate. SiPM-based PET/CT identified additional lesions. Further studies with larger cohorts are needed to confirm these results.

A NIM PET/CT phantom for evaluating the PET image quality of micro-lesions and the performance parameters of CT

BackgroundThe commonly used NEMA IEC Body phantom has a number of defects, hindering its application for detecting micro-lesions and measuring the performance parameters of computed tomography (CT). This study aimed to propose a PET/CT phantom designed by National Institute of Metrology (NIM), China, which is capable of simultaneously testing the performance of PET and CT systems, and to evaluate the quality of imaging.MethodsThe phantom developed in the present study, the NIM PET/CT phantom, is composed of a PET imaging module and a CT imaging module, and these modules are connected together through bolts, which can simultaneously measure the imaging performance of PET and CT systems. Hot spheres were filled with 4:1 sphere-to-background activity concentration using 18F-fluorodeoxyglucose (18F-FDG), and cold spheres were filled with non-radioactive water. We compared the results of imaging obtained from the NIM PET/CT phantom and the NEMA IEC Body phantom to assess their diagnostic efficacy. In order to evaluate the generalization ability of the NIM PET/CT phantom, three different PET/CT systems were used to scan on the same scanning protocol. To evaluate the effects of image reconstruction algorithms on image quality assessment, ordered subset expectation maximization (OSEM), OSEM-point-spread function (PSF), OSEM-TOF, and OSEM-PSF-TOF algorithms were employed.ResultsThe imaging quality of the NIM PET/CT phantom and the NEMA IEC Body phantom was relatively consistent. The NIM PET/CT phantom could detect 7 mm spheres without influencing the imaging quality. It was found that PSF reconstruction exhibited to reduce the speed of convergence, the contrast and background variability of spheres (13–28 mm) were significantly improved after two iterations. In addition to improve the image contrast and background variability, TOF could markedly improve the overall image quality and instrument detection limit. TOF-PSF could noticeably reduce noise level, enhance imaging details, and improve quality of imaging.ConclusionsThe results showed that in comparison with the NEMA IEC Body phantom, the NIM PET/CT phantom outperformed in evaluating the PET image quality of micro-lesions and the performance parameters of CT.

Determining the minimal ultra-low dose CT for reliable attenuation correction of 18F-FDG PET-CT: a phantom study.

To investigate minimal required sub milli-Sievert (mSv) ultra-low dose CT and corresponding tube current and voltage for reliable attenuation correction and semi- quantitation in 18F-FDG PET-CT in an effort for radiation dose reduction. Methods: We performed a PET-CT investigational study using a NEMA torso phantom containing six spheres (diameter: 10, 13, 17, 22, 28, 37 mm) filled with a fixed concentration of 60 kBq/ml and a background of 15 kBq/ml of 18F-FDG. Two sets of PET images, separated by 2 hours, were acquired for 3 minutes in a single bed position using 3-D mode with and without time-of-flight in a GE D-690 scanner. Several sets of CT images were acquired for attenuation correction with different combinations of tube voltage (80, 100, 120 kVp) and effective mAs (tube current-time product divided by pitch), using the maximum beam collimation (64 x 0.625 mm). The lowest CT acquisition technique available on this scanner is 10 mA, 0.4 s and 1.375 for the tube current, tube rotation time and pitch, respectively. The CT radiation dose was estimated based on the computed tomography dose index volume (CTDIvol) measurements performed following the standard methodology and the Imaging Performance Assessment of CT Scanners (ImPACT) calculator. Each of the CT techniques was used for attenuation correction to the same PET acquisition, using ordered-subset expectation maximum (OSEM) algorithm with 24 subsets and 2 iterations. The maximal and average radioactivity (kBq/ml) and standardized uptake values (SUV) of the spheres were measured. The minimal ultra-low dose CT for attenuation correction was determined by reproducible SUV measurements (±10%) compared to our reference CT protocol of 100 kVp and 80 mA for 0.5 s rotation. Results: The minimal ultra-low dose of CT for reproducible quantification in all spheres (<10% relative difference) was determined to be 0.3 mSv for a combination of 100 kVp and 10 mA at 0.5 s rotation, 0.984 helical pitch (0.26 mGy measured CTDIvol) . Based on these results we could confidently determine the CT parameters for reliable attenuation correction of PET images while significantly reducing the associated radiation dose. Conclusion: Our phantom study provided guidance in using ultra-low dose CT for precise attenuation correction and semi-quantification of 18F-FDG PET imaging, which can further reduce CT dose and radiation exposure to patients in clinical PET-CT studies. Clinical application: Based on the data, we can further reduce the radiation dose to sub-mSv using an ultra-low dose CT protocol for reliable attenuation correction in clinical 18F-FDG PET-CT studies.

Initial results of a mouse brain PET insert with a staggered 3-layer DOI detector

Objective. Small animal positron emission tomography (PET) requires a submillimeter resolution for better quantification of radiopharmaceuticals. On the other hand, depth-of-interaction (DOI) information is essential to preserve the spatial resolution while maintaining the sensitivity. Recently, we developed a staggered 3-layer DOI detector with 1 mm crystal pitch and 15 mm total crystal thickness, but we did not demonstrate the imaging performance of the DOI detector with full ring geometry. In this study we present initial imaging results obtained for a mouse brain PET prototype developed with the staggered 3-layer DOI detector. Approach. The prototype had 53 mm inner diameter and 11 mm axial field-of-view. The PET scanner consisted of 16 DOI detectors each of which had a staggered 3-layer LYSO crystal array (4/4/7 mm) coupled to a 4 × 4 silicon photomultiplier array. The physical performance was evaluated in terms of the NEMA NU 4 2008 protocol. Main Results. The measured spatial resolutions at the center and 15 mm radial offset were 0.67 mm and 1.56 mm for filtered-back-projection, respectively. The peak absolute sensitivity of 0.74% was obtained with an energy window of 400–600 keV. The resolution phantom imaging results show the clear identification of a submillimetric rod pattern with the ordered-subset expectation maximization algorithm. The inter-crystal scatter rejection using a narrow energy window could enhance the resolvability of a 0.75 mm rod significantly. Significance. In an animal imaging experiment, the detailed mouse brain structures such as cortex and thalamus were clearly identified with high contrast. In conclusion, we successfully developed the mouse brain PET insert prototype with a staggered 3-layer DOI detector.

Feasibility of ultra-high-speed acquisition in xSPECT bone algorithm: a phantom study with advanced bone SPECT-specific phantom.

Although xSPECT Bone (xB) provides quantitative single-photon emission computed tomography (SPECT) high-resolution images, patients' burden remains high due to long acquisition time; therefore, this study aimed to investigate the feasibility of shortening the xB acquisition time using a custom-designed phantom. A custom-designed xSPECT bone-specific (xSB) phantom with simulated cortical and spongious bones was developed based on the thoracic bone phantom. Both standard- and ultra-high-speed (UHS) xB acquisitions were performed in a male patient with lung cancer. In this phantom study, SPECT was acquired for 3, 6, 9, 12, and 30min. The clinical SPECT acquisition time per rotation was 9 and 3min for standard and UHS, respectively. SPECT images were reconstructed using ordered subset expectation maximization with three-dimensional resolution recovery (Flash3D; F3D) and xB algorithms. Quantitative SPECT value (QSV) and coefficient of variation (CV) were measured using the volume of interests (VOIs) placed at the center of the vertebral body and hot sphere. A linear profile was plotted on the spinous process at the center of the xSB phantom; then, the full width at half maximum (FWHM) was measured. The standardized uptake value (SUV) and standard deviation from the first thoracic to the fifth lumbar vertebrae in clinical standard- and UHS-xB images were measured using a 1-cm3 VOI. The QSV of F3D images was underestimated even in large regions, whereas those of xB images were close to actual radioactivity concentration. The CV was similar or lower for xB images than that for F3D images but was not decreased with increasing acquisition time for both reconstruction images. The FWHM of xB images was lower than those of F3D images at all acquisition times. The mean SUV values from the first thoracic to fifth lumbar vertebrae for standard- and UHS-xB images were 6.73 ± 0.64 and 6.19 ± 0.87, respectively, showing a strong positive correlation. Results of this phantom study suggest that xB imaging can be obtained in only one-third of the acquisition time without compromising the image quality. The SUV of UHS-xB images can be similar to that of standard-xB images in terms of clinical interpretation.

Stochastic EM methods with variance reduction for penalised PET reconstructions

Expectation-maximisation (EM) is a popular and well-established method for image reconstruction in positron emission tomography (PET) but it often suffers from slow convergence. Ordered subset EM (OSEM) is an effective reconstruction algorithm that provides significant acceleration during initial iterations, but it has been observed to enter a limit cycle. In this work, we investigate two classes of algorithms for accelerating OSEM based on variance reduction for penalised PET reconstructions. The first is a stochastic variance reduced EM algorithm, termed as SVREM, an extension of the classical EM to the stochastic context that combines classical OSEM with variance reduction techniques for gradient descent. The second views OSEM as a preconditioned stochastic gradient ascent, and applies variance reduction techniques, i.e., SAGA and SVRG, to estimate the update direction. We present several numerical experiments to illustrate the efficiency and accuracy of the approaches. The numerical results show that these approaches significantly outperform existing OSEM type methods for penalised PET reconstructions, and hold great potential.

Iterative reconstruction incorporating background correction improves quantification of [18F]-NaF PET/CT images of patients with abdominal aortic aneurysm

BackgroundA confounding issue in [18F]-NaF PET/CT imaging of abdominal aortic aneurysms (AAA) is the spill in contamination from the bone into the aneurysm. This study investigates and corrects for this spill in contamination using the background correction (BC) technique without the need to manually exclude the part of the AAA region close to the bone. MethodsSeventy-two (72) datasets of patients with AAA were reconstructed with the standard ordered subset expectation maximization (OSEM) algorithm incorporating point spread function (PSF) modelling. The spill in effect in the aneurysm was investigated using two target regions of interest (ROIs): one covering the entire aneurysm (AAA), and the other covering the aneurysm but excluding the part close to the bone (AAAexc). ROI analysis was performed by comparing the maximum SUV in the target ROI (SUVmax(T)), the corrected cSUVmax (SUVmax(T) − SUVmean(B)) and the target-to-blood ratio (TBR = SUVmax(T)/SUVmean(B)) with respect to the mean SUV in the right atrium region. ResultsThere is a statistically significant higher [18F]-NaF uptake in the aneurysm than normal aorta and this is not correlated with the aneurysm size. There is also a significant difference in aneurysm uptake for OSEM and OSEM + PSF (but not OSEM + PSF + BC) when quantifying with AAA and AAAexc due to the spill in from the bone. This spill in effect depends on proximity of the aneurysms to the bone as close aneurysms suffer more from spill in than farther ones. ConclusionThe background correction (OSEM + PSF + BC) technique provided more robust AAA quantitative assessments regardless of the AAA ROI delineation method, and thus it can be considered as an effective spill in correction method for [18F]-NaF AAA studies.

Impact of reconstruction parameters on spatial resolution and its comparison between cadmium-zinc-telluride SPECT/CT and conventional SPECT/CT.

To evaluate the impact of reconstruction parameters on the spatial resolution of the tomographic image in single photon emission computed tomography (SPECT)/computed tomography (CT), and compare spatial resolution between a new polyvalent whole-body Cadmium-Zinc-Telluride camera (CZT-SPECT/CT) and a conventional dual-head Anger camera (conventional SPECT/CT). Spatial resolution was evaluated with four-line sources filled with 99mTc in tomographic images reconstructed by varying reconstruction parameters. Ordered-subset expectation maximization (OSEM) algorithm was performed with varying iterations (1-20), the number of subsets was fixed at 10. Butterworth filter, Gauss filter and no-filter were selected, respectively. Computed tomography-based attenuation correction (CTAC), scatter correction, resolution recovery and no correction (NC) were adopted for image correction. Filtered back projection (FBP) with Butterworth filter and CTAC was performed in image reconstruction. Spatial resolution was expressed by the full width at half-maximum (FWHM) value. The impact of reconstruction parameters on the spatial resolution was identical in both cameras: FWHM values decreased with the increase of iterations and converged uniformly when the number of iterations was over 4. FWHM values decreased with the increase of cutoff frequency of the Butterworth filter and increased with the increase of the Gauss filter. scatter correction and resolution recovery improved spatial resolution, whereas CTAC had a negligible effect on spatial resolution when reconstructed by OSEM. FWHM was generally lower with OSEM reconstruction than FBP reconstruction. On the whole, under the same reconstruction conditions, CZT-SPECT/CT had a lower FWHM value than conventional SPECT/CT. The spatial resolution was improved with the increase of iterations. Increasing the cutoff frequency of the Butterworth filter and decreasing the Gauss filter enhanced spatial resolution. The spatial resolution was better reconstructed by OSEM associated with attenuation correction, scatter correction and resolution recovery than FBP. CZT-SPECT/CT had better spatial resolution than conventional SPECT/CT.