Ordered Subset Expectation Maximization Research Articles

Accuracy of holmium-166 SPECT/CT quantification over a large range of activities

BackgroundQuantitative imaging is a crucial step for dosimetry in radionuclide therapies. Traditionally, SPECT/CT imaging is quantified based on scanner-specific conversion factors or self-calibration, but recently absolute quantification methods have been introduced in commercial SPECT reconstruction software (Broad Quantification, Siemens Healthineers). In this phantom study we investigate the accuracy of three quantification methods for holmium-166 SPECT/CT imaging, and provide recommendations for clinical dosimetry.MethodsOne cylindrical phantom, filled with a homogeneous holmium-166-chloride activity concentration solution, was imaged at one time point to determine a scanner-specific conversion factor, and to characterize the spatial dependency of the activity concentration recovery. One Jaszczak phantom with six fillable spheres, 10:1 sphere-to-background ratio, was imaged over a large range of holmium-166 activities (61-3130 MBq). The images were reconstructed with either an ordered subset expectation maximization (OSEM, Flash3D-reconstruction; scanner-specific quantification or self-calibration quantification) or an ordered subset conjugate gradient (OSCG, xSPECT-reconstruction; Broad Quantification) algorithm. These three quantification methods were compared for the data of the Jaszczak phantom and evaluated based on whole phantom recovered activity, activity concentration recovery coefficients (ACRC), and recovery curves.ResultsThe activity recovery in the Jaszczak phantom was 28–115% for the scanner-specific, and 57–97% for the Broad Quantification quantification methods, respectively. The self-calibration-based activity recovery is inherently always 100%. The ACRC for the largest sphere (Ø60 mm, ~ 113 mL) ranged over (depending on the activity level) 0.22–0.89, 0.76–0.86, 0.39–0.72 for scanner-specific, self-calibration and Broad Quantification, respectively.ConclusionOf the three investigated quantification methods, the self-calibration technique produces quantitative SPECT images with the highest accuracy in the investigated holmium-166 activity range.

Image Quality of Cardiac Silicon Photomultiplier PET/CT Using an Infant Phantom of Extremely Low Birth Weight.

The lack of pediatrics-specific equipment for nuclear medicine imaging has resulted in insufficient diagnostic information for newborns, especially low-birth-weight infants. Although PET offers high spatial resolution and low radiation exposure, its use in newborns is limited. This study investigated the feasibility of cardiac PET imaging using the latest silicon photomultiplier (SiPM) PET technology in infants of extremely low birth weight (ELBW) using a phantom model. Methods: The study used a phantom model representing a 500-g ELBW infant with brain, cardiac, liver, and lung tissues. The cardiac tissue included a 3-mm-thick defect mimicking myocardial infarction. Organ tracer concentrations were calculated assuming 18F-FDG myocardial viability scans and 18F-flurpiridaz myocardial perfusion scans and were added to the phantom organs. Imaging was performed using an SiPM PET/CT scanner with a 5-min acquisition. The data acquired in list mode were reconstructed using 3-dimensional ordered-subsets expectation maximization with varying iterations. Image evaluation was based on the depiction of the myocardial defect compared with normal myocardial accumulation. Results: Increasing the number of iterations improved the contrast of the myocardial defect for both tracers, with 18F-flurpiridaz showing higher contrast than 18F-FDG. However, even at 50 iterations, both tracers overestimated the defect accumulation. A bull's-eye image can display the flow metabolism mismatch using images from both tracers. Conclusion: SiPM PET enabled cardiac PET imaging in a 500-g ELBW phantom with a 1-g heart. However, there were limitations in adequately depicting these defects. Considering the image quality and defect contrast,18F-flurpiridaz appears more desirable than 18F-FDG if only one of the two can be used.

Refining penalty parameter selection in whole-body PET image reconstruction for lung cancer patients using the cross-validation log-likelihood method

Objective.Penalty parameters in penalized likelihood positron emission tomography (PET) reconstruction are typically determined empirically. The cross-validation log-likelihood (CVLL) method has been introduced to optimize these parameters by maximizing a CVLL function, which assesses the likelihood of reconstructed images using one subset of a list-mode dataset based on another subset. This study aims to validate the efficacy of the CVLL method in whole-body imaging for cancer patients using a conventional clinical PET scanner.Approach.Fifteen lung cancer patients were injected with 243.7 ± 23.8 MBq of [18F]FDG and underwent a 22 min PET scan on a Biograph mCT PET/CT scanner, starting at 60 ± 5 min post-injection. The PET list-mode data were partitioned by subsampling without replacement, with 20 minutes of data for image reconstruction using an in-house ordered subset expectation maximization algorithm and the remaining 2 minutes of data for cross-validation. Two penalty parameters, penalty strengthβand Fair penalty function parameterδ, were subjected to optimization. Whole-body images were reconstructed, and CVLL values were computed across various penalty parameter combinations. The optimal image corresponding to the maximum CVLL value was selected by a grid search for each patient.Main results.Theδvalue required to maximize the CVLL value was notably small (⩽10-6in this study). The influences of voxel size and scan duration on image optimization were investigated. A correlation analysis revealed a significant inverse relationship between optimalβand scan count level, with a correlation coefficient of -0.68 (p-value = 3.5 × 10-5). The optimal images selected by the CVLL method were compared with those chosen by two radiologists based on their diagnostic preferences. Differences were observed in the selection of optimal images.Significance.This study demonstrates the feasibility of incorporating the CVLL method into routine imaging protocols, potentially allowing for a wide range of combinations of injected radioactivity amounts and scan durations in modern PET imaging.

Phantom test procedures for a new neuro-oncological amino acid PET tracer: [18F]fluciclovine.

Amino acid positron emission tomography (PET) examinations using anti-1-amino-3-[18F]-fluorocyclobutane-1-carboxylic acid ([18F]FACBC) were allowed for routine clinical use in July 2024. However, phantom test procedures for [18F]FACBC reconstruction parameters have not yet been established. The present study aimed to establish new phantom test procedures for [18F]FACBC brain PET imaging to determine optimal reconstruction parameters. Background (BG) activity as well as hot sphere and target-to-background ratios (TBRs) of [18F]FACBC were estimated based on brain activity and tumor-to-normal tissue ratios (TNR) in a Japanese clinical trial of [18F]FACBC. Phantom experiments proceeded under [18F]FACBC or L-[methyl-11C]-methionine ([11C]MET) conditions. The number of iterations and the Gaussian filter parameters were determined from the reconstruction parameters %contrastmean and coefficients of variation (CVs) in ordered subset expectation maximization (OSEM) and time-of-flight (TOF) with or without point-spread-function (PSF) correction. The amounts of activity in the hot spheres and BG were 1.1 and 5.5kBq/mL, respectively, and the TBR was 5.0 at the start of acquisition. The %contrastmean of all hot spheres was higher with [18F]FACBC than [11C]MET, and %contrastmean converged between 4 and 6 iterations in hot spheres with diameters < 10mm. We used four iterations for OSEM + TOF and five for OSEM + TOF + PSF correction for [18F]FACBC and [11C]MET images. The CV was higher for [18F]FACBC than [11C]MET. The optimal sizes of Gaussian filters for OSEM + TOF and OSEM + TOF + PSF correction of image reconstruction were 5mm for [18F]FACBC, and 4 and 3mm, respectively, for [11C]MET images. We estimated phantom activity and TBR based on brain activity in a Japanese clinical trial and established new phantom test procedures for [18F]FACBC. We recommend that the optimal reconstruction parameters for [18F]FACBC should be set to the same number of iterations as [11C]MET and that the FWHM of Gaussian filter should have a few mm higher than [11C]MET to reduce image noise from brain normal tissue.

The effect of reconstruction algorithms on semi-quantitative measurements in 18F-FDG PET/CT imaging.

This study was carried out to understand whether Q.Clear and ordered subset expectation maximization (OSEM), reconstruction algorithms used in fluorine-18-fluorodeoxyglucose (18F-FDG) positron emission tomography/computed tomography (PET/CT) applications, and parameters such as time of flight (TOF) and point spread function (PSF) cause different results in semi-quantitative measurements. Raw PET data of 264 patients who were referred to 18F-FDG PET/CT imaging with the purpose of evaluation of known or suspicious malignant disease were reconstructed separately with Q.Clear (GE Healthcare), a BPL, an OSEM algorithm, PSF (SharpIR®) and TOF (VUE Point FX®) methods. Each patient's liver, mediastinal blood pool, metabolic tumor volume (MTV), total lesion glycolysis (TLG), and standardized uptake values (SUV) (SUVmax, SUVmean, and SUVpeak) of a total of 264 lesions selected from the patients were performed. β350+ToF yielded higher measurement results than all other variables for all of the lesion SUVmax, lesion SUVmean, L/AP SUVmax, and L/AP SUVmean parameters. OSEM+ToF and OSEM+TOF+PSF algorithms yielded higher mean and median SUVmax values for the reference structures (liver and mediastinum) and for lesions SUVmax and SUVmean values were statistically significantly lower than the β350+ToF method. The method with the lowest mean value for the L/Liver SUVmax variable was OSEM+ToF 4iter16ss (mean=1.76), while the method with the highest mean value was β350+ToF (mean=2.26). β350+ToF was the reconstruction method with the highest ratios for L/AP SUVmax and SUVmean for both lesions below and above 1 cm. β350+ToF algorithm had also statistically significantly higher results for these variables compared to all other parameters in malignant lesions. When comparing 18F-FDG PET/CT images, the use of different reconstruction algorithms may lead to misleading results, especially in the evaluation of response to treatment of malignancies.

Comparison of Measurement and Prognostic Power of SUV Between High-Definition and Standard PET Imaging in Non-Small Cell Lung Cancer Patients.

This study aimed to evaluate the measurement and prognostic ability of the SUVmax of whole-body tumors (SUVmaxwb) in non-small cell lung cancer (NSCLC) patients, comparing high-definition (HD) PET imaging with standard-definition (SD) PET imaging. Methods: The study included 242 consecutive NSCLC patients who underwent baseline 18F-FDG PET/CT from April 2018 to January 2021. Two imaging techniques were used: HD PET (using ordered-subsets expectation maximization with point-spread function modeling and time-of-flight techniques and smaller voxels) and SD PET (with ordered-subsets expectation maximization and time-of-flight techniques). SUVmaxwb was determined by measuring all the tumor lesions in the whole body, and tumor-to-background ratio (TBR) was calculated using the background SUVmean of various body parts. Results: The patient cohort had an average age of 68.3 y, with 59.1% being female. During a median follow-up of 29.6 mo, 83 deaths occurred. SUVmaxwb was significantly higher in HD PET than SD PET, with respective medians of 17.4 and 11.8. The TBR of 1,125 tumoral lesions was also higher in HD PET. Univariate Cox regression analysis showed that SUVmaxwb from both HD and SD PET were significantly associated with overall survival. However, after adjusting for TNM (tumor, node, metastasis) stage, only SUVmaxwb from SD PET remained significantly associated with survival. Conclusion: HD PET imaging in NSCLC patients yields higher SUVmaxwb and TBR, enhancing tumor visibility. Despite this, its prognostic value is less significant than SD PET after adjusting clinical TNM stage. Thus, consideration should be given to using HD PET reconstruction to increase tumor visibility, and SD PET is recommended for NSCLC patient prognostication and therapeutic evaluation, as well as for the classification of lung nodules.

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

Positron Emission Tomography (PET) imaging after 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 for 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.

Optimizing scan time and bayesian penalized likelihood reconstruction algorithm in copper-64 PET/CT imaging: a phantom study

Objectives: The aim of this study was to evaluate Cu-64 PET phantom image quality using Bayesian Penalized Likelihood (BPL) and Ordered Subset Expectation Maximum with point-spread function modeling (OSEM-PSF) reconstruction algorithms. In the BPL, the regularization parameter β was varied to identify the optimum value for image quality. In the OSEM-PSF, the effect of acquisition time was evaluated to assess the feasibility of shortened scan duration. Methods: A NEMA IEC PET body phantom was filled with known activities of water soluble Cu-64. The phantom was imaged on a PET/CT scanner and was reconstructed using BPL and OSEM-PSF algorithms. For the BPL reconstruction, various β values (150, 250, 350, 450, and 550) were evaluated. For the OSEM-PSF algorithm, reconstructions were performed using list-mode data intervals ranging from 7.5 to 240 s. Image quality was assessed by evaluating the signal to noise ratio (SNR), contrast to noise ratio (CNR), and background variability (BV). Results: The SNR and CNR were higher in images reconstructed with BPL compared to OSEM-PSF. Both the SNR and CNR increased with increasing β, peaking at β = 550. The CNR for all β, sphere sizes and tumor-to-background ratios (TBRs) satisfied the Rose criterion for image detectability (CNR > 5). BPL reconstructed images with β = 550 demonstrated the highest improvement in image quality. For OSEM-PSF reconstructed images with list-mode data duration ≥ 120 s, the noise level and CNR were not significantly different from the baseline 240 s list-mode data duration. Conclusions: BPL reconstruction improved Cu-64 PET phantom image quality by increasing SNR and CNR relative to OSEM-PSF reconstruction. Additionally, this study demonstrated scan time can be reduced from 240 to 120 s when using OSEM-PSF reconstruction while maintaining similar image quality. This study provides baseline data that may guide future studies aimed to improve clinical Cu-64 imaging.

Comparison of quantitative whole body PET parameters on [68Ga]Ga-PSMA-11 PET/CT using ordered Subset Expectation Maximization (OSEM) vs. bayesian penalized likelihood (BPL) reconstruction algorithms in men with metastatic castration-resistant prostate cancer

BackgroundPSMA PET/CT is a predictive and prognostic biomarker for determining response to [177Lu]Lu-PSMA-617 in patients with metastatic castration resistant prostate cancer (mCRPC). Thresholds defined to date may not be generalizable to newer image reconstruction algorithms. Bayesian penalized likelihood (BPL) reconstruction algorithm is a novel reconstruction algorithm that may improve contrast whilst preventing introduction of image noise. The aim of this study is to compare the quantitative parameters obtained using BPL and the Ordered Subset Expectation Maximization (OSEM) reconstruction algorithms.MethodsFifty consecutive patients with mCRPC who underwent [68Ga]Ga-PSMA-11 PET/CT using OSEM reconstruction to assess suitability for [177Lu]Lu-PSMA-617 therapy were selected. BPL algorithm was then used retrospectively to reconstruct the same PET raw data. Quantitative and volumetric measurements such as tumour standardised uptake value (SUV)max, SUVmean and Molecular Tumour Volume (MTV-PSMA) were calculated on both reconstruction methods. Results were compared (Bland-Altman, Pearson correlation coefficient) including subgroups with low and high-volume disease burdens (MTV-PSMA cut-off 40 mL).ResultsThe SUVmax and SUVmean were higher, and MTV-PSMA was lower in the BPL reconstructed images compared to the OSEM group, with a mean difference of 8.4 (17.5%), 0.7 (8.2%) and − 21.5 mL (-3.4%), respectively. There was a strong correlation between the calculated SUVmax, SUVmean, and MTV-PSMA values in the OSEM and BPL reconstructed images (Pearson r values of 0.98, 0.99, and 1.0, respectively). No patients were reclassified from low to high volume disease or vice versa when switching from OSEM to BPL reconstruction.Conclusions[68Ga]Ga-PSMA-11 PET/CT quantitative and volumetric parameters produced by BPL and OSEM reconstruction methods are strongly correlated. Differences are proportional and small for SUVmean, which is used as a predictive biomarker. Our study suggests that both reconstruction methods are acceptable without clinical impact on quantitative or volumetric findings. For longitudinal comparison, committing to the same reconstruction method would be preferred to ensure consistency.

Novel Detector Configurations in Cone-Beam CT Systems: A Simulation Study

Cone-beam computed tomography (CBCT) has emerged in recent years as an adequate alternative to mammography and tomosynthesis due to the several advantages over traditional mammography, including its ability to provide 3D images, its reduced radiation dose, and its ability to image dense breasts more effectively and conduct more effective breast compressions, etc. Furthermore, CBCT is capable of providing images with high sensitivity and specificity, allowing a more accurate evaluation, even of dense breasts, where mammography and tomosynthesis may lead to a false diagnosis. Clinical and experimental CBCT systems rely on cesium iodine (CsI:Tl) scintillators for X-ray energy conversion. This study comprises an investigation among different novel CBCT detector technologies, consisting either of scintillators (BGO, LSO:Ce, LYSO:Ce, LuAG:Ce, CaF2:Eu, LaBr3:Ce) or semiconductors (Silicon, CZT) in order to define the optimum detector design for a future experimental setup, dedicated to breast imaging. For this purpose, a micro-CBCT system was adapted, using GATE v9.2.1, consisting of the aforementioned various detection schemes. Two phantom configurations were selected: (a) an aluminum capillary positioned at the center of the field of view in order to calculate the system’s spatial resolution and (b) a breast phantom consisting of spheres of different materials, such that their characteristics are close to the breast composition. Breast phantom contrast-to-noise ratios (CNRs) were extracted from the phantom’s tomographic images. The images were reconstructed with filtered back projection (FBP) and ordered subsets expectation-maximization (OSEM) algorithms. The semiconductors acted satisfactorily in low-density matter, while LYSO:Ce, LaBr3:Ce, and LuAG:Ce presented adequate CNRs for all the different spheres’ densities. The energy converters that are presented in this study were evaluated for their performance against the standard CsI:Tl crystal.

Development of Solutions Based on Computer Simulation for Quantitative SPECT/CT

Purpose: In modern nuclear medicine, an ordered subset expectation maximization (OSEM) iterative algorithm is used to obtain images by single photon emission computed tomography (SPECT) method. At the same time, the problem of choosing the stopping iteration number at which the solution will be as close as possible to the exact value has not been solved. The purpose of this work is to study the problem of stopping the OSEM algorithm. Material and methods: The studies were carried out using the method of mathematical modeling based on the software package “Virtual Platform for Simulation Testing of the SPECT/CT Method”. Computer experiments simulated SPECT/CT with 99mTc phosphotech of a patient with metastatic lesion of the posterior segment of the third rib. The stopping criterion of the iterative OSEM algorithm was studied to obtain a solution that is quantitatively closest to the exact value of the accumulated activity at the tumor. Results: In computer modeling experiments, the tumor dependent Pearson chi-square test was shown to provide a solution that should give correct values for the SUVmean score. However, further research is needed on this important problem to determine how strongly the number of iterations depends on the size of the tumor and the amount of activity accumulated in it. Conclusion: Based on the results obtained, we can conclude that in the SPECT method, in order to move to accurate quantitative estimates of the accumulated activity in the tumor, it is necessary to develop a tumor-oriented stopping criterion for iterative algorithm. If there are several tumors, the stopping iteration numbers for them may differ. Further research on this important issue is planned.

Phantom and clinical evaluation of the Bayesian penalised likelihood reconstruction algorithm Q.Clear without PSF correction in amyloid PET images

PurposeBayesian penalised likelihood (BPL) reconstruction, which incorporates point-spread-function (PSF) correction, provides higher signal-to-noise ratios and more accurate quantitation than conventional ordered subset expectation maximization (OSEM) reconstruction. However, applying PSF correction to brain PET imaging is controversial due to Gibbs artefacts that manifest as unpredicted cortical uptake enhancement. The present study aimed to validate whether BPL without PSF would be useful for amyloid PET imaging.MethodsImages were acquired from Hoffman 3D brain and cylindrical phantoms for phantom study and 71 patients administered with [18F]flutemetamol in clinical study using a Discovery MI. All images were reconstructed using OSEM, BPL with PSF correction, and BPL without PSF correction. Count profile, %contrast, recovery coefficients (RCs), and image noise were calculated from the images acquired from the phantoms. Amyloid β deposition in patients was visually assessed by two physicians and quantified based on the standardised uptake value ratio (SUVR).ResultsThe overestimated radioactivity in profile curves was eliminated using BPL without PSF correction. The %contrast and image noise decreased with increasing β values in phantom images. Image quality and RCs were better using BPL with, than without PSF correction or OSEM. An optimal β value of 600 was determined for BPL without PSF correction. Visual evaluation almost agreed perfectly (κ = 0.91–0.97), without depending on reconstruction methods. Composite SUVRs did not significantly differ between reconstruction methods.ConclusionGibbs artefacts disappeared from phantom images using the BPL without PSF correction. Visual and quantitative evaluation of [18F]flutemetamol imaging was independent of the reconstruction method. The BPL without PSF correction could be the standard reconstruction method for amyloid PET imaging, despite being qualitatively inferior to BPL with PSF correction for [18F]flutemetamol amyloid PET imaging.

Stability of radiomic features from positron emission tomography images: a phantom study comparing advanced reconstruction algorithms and ordered subset expectation maximization.

In this study, we compared the repeatability and reproducibility of radiomic features obtained from positron emission tomography (PET) images according to the reconstruction algorithm used-advanced reconstruction algorithms, such as HYPER iterative (IT), HYPER deep learning reconstruction (DLR), and HYPER deep progressive reconstruction (DPR), or traditional Ordered Subset Expectation Maximization (OSEM)-to understand the potential variations and implications of using advanced reconstruction techniques in PET-based radiomics. We used a heterogeneous phantom with acrylic spherical beads (4- or 8-mm diameter) filled with 18F. PET images were acquired and reconstructed using OSEM, IT, DLR, and DPR. Original and wavelet radiomic features were calculated using SlicerRadiomics. Radiomic feature repeatability was assessed using the Coefficient of Variance (COV) and intraclass correlation coefficient (ICC), and inter-acquisition time reproducibility was assessed using the concordance correlation coefficient (CCC). For the 4- and 8-mm diameter beads phantom, the proportion of radiomic features with a COV < 10% was equivocal or higher for the advanced reconstruction algorithm than for OSEM. ICC indicated that advanced methods generally outperformed OSEM in repeatability, except for the original features of the 8-mm beads phantom. In the inter-acquisition time reproducibility analysis, the combinations of 3 and 5min exhibited the highest reproducibility in both phantoms, with IT and DPR showing the highest proportion of radiomic features with CCC > 0.8. Advanced reconstruction methods provided enhanced stability of radiomic features compared with OSEM, suggesting their potential for optimal image reconstruction in PET-based radiomics, offering potential benefits in clinical diagnostics and prognostics.

Iterative Reconstruction of Micro Computed Tomography Scans Using Multiple Heterogeneous GPUs.

Graphics processing units (GPUs) facilitate massive parallelism and high-capacity storage, and thus are suitable for the iterative reconstruction of ultrahigh-resolution micro computed tomography (CT) scans by on-the-fly system matrix (OTFSM) calculation using ordered subsets expectation maximization (OSEM). We propose a finite state automaton (FSA) method that facilitates iterative reconstruction using a heterogeneous multi-GPU platform through parallelizing the matrix calculations derived from a ray tracing system of ordered subsets. The FSAs perform flow control for parallel threading of the heterogeneous GPUs, which minimizes the latency of launching ordered-subsets tasks, reduces the data transfer between the main system memory and local GPU memory, and solves the memory-bound of a single GPU. In the experiments, we compared the operation efficiency of OS-MLTR for three reconstruction environments. The heterogeneous multiple GPUs with job queues for high throughput calculation speed is up to five times faster than the single GPU environment, and that speed up is nine times faster than the heterogeneous multiple GPUs with the FIFO queues of the device scheduling control. Eventually, we proposed an event-triggered FSA method for iterative reconstruction using multiple heterogeneous GPUs that solves the memory-bound issue of a single GPU at ultrahigh resolutions, and the routines of the proposed method were successfully executed on each GPU simultaneously.

Impact of [11C]methionine PET with Bayesian penalized likelihood reconstruction on glioma grades based on new WHO 2021 classification.

The uptake of [11C]methionine in positron emission tomography (PET) imaging overlapped in earlier images of tumors. Bayesian penalized likelihood (BPL) reconstruction increases the quantitative values of tumors compared with conventional ordered subset-expectation maximization (OSEM). The present study aimed to grade glioma malignancy based on the new WHO 2021 classification using [11C]methionine PET images reconstructed using BPL. We categorized 32 gliomas in 28 patients as grades 2/3 (n = 15) and 4 (n = 17) based on the WHO 2021 classification. All [11C]methionine images were reconstructed using OSEM + time-of-flight (TOF) and BPL + TOF (β = 200). Maximum standardized uptake value (SUVmax) and tumor-to-normal tissue ratio (T/Nmax) were measured at each lesion. The mean SUVmax was 4.65 and 4.93 in grade 2/3 and 6.38 and 7.11 in grade 4, and the mean T/Nmax was 7.08 and 7.22 in grade 2/3 and 9.30 and 10.19 in grade 4 for OSEM and BPL, respectively. The BPL significantly increased these values in grade 4 gliomas. The area under the receiver operator characteristic (ROC) curve (AUC) for SUVmax was the highest (0.792) using BPL. The BPL increased mean SUVmax and mean T/Nmax in lesions with higher contrast such as grade 4 glioma. The discrimination power between grades 2/3 and 4 in SUVmax was also increased using [11C]methionine PET images reconstructed with BPL.

Feasibility of noise-reduction reconstruction technology based on non-local-mean principle in SiPM-PET/CT

Quantitative values of positron emission tomography (PET) images using non-local-mean in a silicon photomultiplier (SiPM)-PET/computed tomography (CT) system with phantom and clinical images.The evaluation was conducted on a National Electrical Manufacturers Association body phantom with micro-spheres (4, 5, 6, 8, 10, 13 mm) and clinical images using the SiPM-PET/CT system.The signal-to-background ratio of the phantom was set to 4, and all PET image data was obtained and reconstructed using three-dimensional ordered subset expectation maximization, time-of-flight, point-spread function, and a 4-mm Gaussian filter (GF) and clear adaptive low-noise method (CaLM) in mild, standard, and strong intensities.The evaluation included the standardized uptake value (SUV), percent contrast (QH), coefficient of variation of the background area (CVbackground) clinical imaging for SUV of lung nodules, liver signal-to-noise ratio (SNR), and visual evaluation.SUVmax for 8-mm sphere in phantom images at 2 min for GF and CaLM (mild, standard, strong) were 2.11, 2.32, 2.02, and 1.72; the QH, 8 mm was 27.33 %, 27.47 %, 21.81 %, and 16.09 %; and CVbackground was 12.78, 11.35, 7.86, and 4.71, respectively.CaLM demonstrated higher SUVmax in clinical images than GF for all lung nodule sizes. The average SUVmax for nodules with a diameter of ≤ 1 cm were 5.9 ± 2.4, 9.9 ± 4.9, 9.9 ± 5.0, and 9.9 ± 5.0 for GF and CaLM-mild, standard, and strong intensities, respectively.Liver SNRs were higher for CaLM (mild, standard, strong) compared to GF, with increasing CaLM intensity causing higher liver SNR. CaLM-mild and standard demonstrated suitability for diagnosis in visual evaluation.

Investigation of PET image quality with acquisition time/bed and enhancement of lesion quantification accuracy through deep progressive learning

ObjectiveTo improve the PET image quality by a deep progressive learning (DPL) reconstruction algorithm and evaluate the DPL performance in lesion quantification.MethodsWe reconstructed PET images from 48 oncological patients using ordered subset expectation maximization (OSEM) and deep progressive learning (DPL) methods. The patients were enrolled into three overlapped studies: 11 patients for image quality assessment (study 1), 34 patients for sub-centimeter lesion quantification (study 2), and 28 patients for imaging of overweight or obese individuals (study 3). In study 1, we evaluated the image quality visually based on four criteria: overall score, image sharpness, image noise, and diagnostic confidence. We also measured the image quality quantitatively using the signal-to-background ratio (SBR), signal-to-noise ratio (SNR), contrast-to-background ratio (CBR), and contrast-to-noise ratio (CNR). To evaluate the performance of the DPL algorithm in quantifying lesions, we compared the maximum standardized uptake values (SUVmax), SBR, CBR, SNR and CNR of 63 sub-centimeter lesions in study 2 and 44 lesions in study 3.ResultsDPL produced better PET image quality than OSEM did based on the visual evaluation methods when the acquisition time was 0.5, 1.0 and 1.5 min/bed. However, no discernible differences were found between the two methods when the acquisition time was 2.0, 2.5 and 3.0 min/bed. Quantitative results showed that DPL had significantly higher values of SBR, CBR, SNR, and CNR than OSEM did for each acquisition time. For sub-centimeter lesion quantification, the SUVmax, SBR, CBR, SNR, and CNR of DPL were significantly enhanced, compared with OSEM. Similarly, for lesion quantification in overweight and obese patients, DPL significantly increased these parameters compared with OSEM.ConclusionThe DPL algorithm dramatically enhanced the quality of PET images and enabled more accurate quantification of sub-centimeters lesions in patients and lesions in overweight or obese patients. This is particularly beneficial for overweight or obese patients who usually have lower image quality due to the increased attenuation.

The determination of the optimal threshold on measurement of thyroid volume using quantitative SPECT/CT for Graves' hyperthyroidism

PurposeTo investigate the optimal threshold for measuring thyroid volume in patients with Grave's hyperthyroidism (GH) by SPECT/CT.Materials and methodsA 53 mL butterfly-shaped hollow container made of two 45-degree transparent elbows was put into a NEMA IEC phantom tank. The butterfly-shaped container and the tank were then filled with Na99mTcO4 of different radioactive concentrations, respectively, which could simulate thyroid gland with GH by different target-to-background ratios (T/B) (200:1, 600:1, 1000:1). The different T/B of planar imaging and SPECT/CT were acquired by a Discovery NM/CT 670 Pro SPECT/CT. With Thyroid software (Version 4.0) of GE-Xeleris workstation, the region of the thyroid gland in planar imaging was delineated. The thyroid area and average long diameter of both lobes were substituted into the Allen formula to calculate the thyroid volume. The calculation error was compared with the actual volume. Q-Metrix software was used to perform CT-based attenuation correction, scatter correction, resolution recovery. Ordered-subsets expectation maximization was used to reconstruct SPECT data. 20%, 25%, 30%, 40%, 50%, 60% thresholds were selected to automatically delineate the volume of interest and compared with the real volume, which determinated the optimal threshold. We measured the thyroid volume of 40 GH patients using the threshold and compared the volumes obtained by planar imaging and ultrasound three-dimensional. The differences of the volumes with different T/B and thresholds were compared by the ANOVA and least significant difference t test. The volumes delineated by SPECT/CT were evaluated using ANOVA, least significant difference t test, correlation analysis and, linear regression and Bland–Altman concordance test plot. The differences and consistency of thyroid volume were compared among the above three methods.ResultsThere was no significant difference in the results between different T/B models (P > 0.05). The thyroid volume calculated by the planar imaging formula method was higher than the real volume, with an average overestimation of 22.81%. The volumes delineated by SPECT/CT threshold automatically decreased while the threshold increased. There were significant differences between groups with different thresholds (P < 0.001). With an average error of 3.73%, the thyroid volume analyzed by the threshold of 25% was close to the results of ultrasound measurement (P > 0.05). Thyroid volume measured by planar imaging method was significantly higher than ultrasound and SPECT/CT threshold automatic delineation method (P < 0.05). The agreement between the SPECT/CT 25% threshold and ultrasound (r = 0.956, b = 0.961) was better than that between the planar imaging and ultrasound (r = 0.590, b = 0.574). The Bland–Altman plot also showed that the thyroid volume measured by the 25% threshold automatic delineation method was in good agreement with the ultrasound measurement.ConclusionsThe T/B has no effect on the measurement of thyroid volume in GH patients; planar imaging method can significantly overestimate thyroid volume in GH patients, and 25% threshold automatic delineation method can obtain more accurate thyroid volume in GH patients.

Modular GAN: positron emission tomography image reconstruction using two generative adversarial networks.

The reconstruction of PET images involves converting sinograms, which represent the measured counts of radioactive emissions using detector rings encircling the patient, into meaningful images. However, the quality of PET data acquisition is impacted by physical factors, photon count statistics and detector characteristics, which affect the signal-to-noise ratio, resolution and quantitative accuracy of the resulting images. To address these influences, correction methods have been developed to mitigate each of these issues separately. Recently, generative adversarial networks (GANs) based on machine learning have shown promise in learning the complex mapping between acquired PET data and reconstructed tomographic images. This study aims to investigate the properties of training images that contribute to GAN performance when non-clinical images are used for training. Additionally, we describe a method to correct common PET imaging artefacts without relying on patient-specific anatomical images. The modular GAN framework includes two GANs. Module 1, resembling Pix2pix architecture, is trained on non-clinical sinogram-image pairs. Training data are optimised by considering image properties defined by metrics. The second module utilises adaptive instance normalisation and style embedding to enhance the quality of images from Module 1. Additional perceptual and patch-based loss functions are employed in training both modules. The performance of the new framework was compared with that of existing methods, (filtered backprojection (FBP) and ordered subset expectation maximisation (OSEM) without and with point spread function (OSEM-PSF)) with respect to correction for attenuation, patient motion and noise in simulated, NEMA phantom and human imaging data. Evaluation metrics included structural similarity (SSIM), peak-signal-to-noise ratio (PSNR), relative root mean squared error (rRMSE) for simulated data, and contrast-to-noise ratio (CNR) for NEMA phantom and human data. For simulated test data, the performance of the proposed framework was both qualitatively and quantitatively superior to that of FBP and OSEM. In the presence of noise, Module 1 generated images with a SSIM of 0.48 and higher. These images exhibited coarse structures that were subsequently refined by Module 2, yielding images with an SSIM higher than 0.71 (at least 22% higher than OSEM). The proposed method was robust against noise and motion. For NEMA phantoms, it achieved higher CNR values than OSEM. For human images, the CNR in brain regions was significantly higher than that of FBP and OSEM (p < 0.05, paired t-test). The CNR of images reconstructed with OSEM-PSF was similar to those reconstructed using the proposed method. The proposed image reconstruction method can produce PET images with artefact correction.

PET/CT Reconstruction and Its Impact on [Measures of] Metabolic Tumor Volume

In oncological imaging, the use of metabolic tumor volume (MTV) for further prognostic differentiation and the development of risk adapted strategies appears promising. The aim of this analysis was to evaluate ultra-high definition (UHD) and ordered subset expectation maximization (OSEM) PET/CT reconstructions for their potential impact on different methods of MTV measurement. We analyzed positron emission tomography combined with computed tomography (PET/CT) scans of 40 Hodgkin lymphoma patients before first-line treatment who had undergone fluorodeoxyglucose (FDG) PET/CT. The MTVs were determined taking an SUV of 4.0 (MTV4.0) as a fixed threshold or 41% of the single hottest voxel (MTV41%) as an adaptive threshold for automated lymphoma delineation in both UHD and OSEM reconstructions. We then compared the absolute and relative differences between MTV4.0 and MTV41% in UHD and OSEM reconstructions. The relative distribution of MTV4.0 and MTV41% in relation to the reconstruction method applied was recorded and respective differences were tested for statistical significance using the paired sample t-test. A comparison of MTV4.0 and MTV41% showed smaller relative and absolute differences in MTV between different reconstruction settings for the MTV4.0 method. Conversely, the absolute as well as the relative differences between MTVs obtained from different reconstructions settings were significantly greater when the MTV41% method was applied (p<0001). MTV4.0 brings higher robustness between different reconstruction settings, while with MTV41% the deviation between volumes obtained with different reconstruction settings is greater. For clinical routine and for multicenter settings, the MTV4.0 therefore appears most promising.