Dataset Information

Harmonization of PET image reconstruction parameters in simultaneous PET/MRI

ABSTRACT:

Objective

Simultaneous PET/MRIs vary in their quantitative PET performance due to inherent differences in the physical systems and differences in the image reconstruction implementation. This variability in quantitative accuracy confounds the ability to meaningfully combine and compare data across scanners. In this work, we define image reconstruction parameters that lead to comparable contrast recovery curves across simultaneous PET/MRI systems.

Method

The NEMA NU-2 image quality phantom was imaged on one GE Signa and on one Siemens mMR PET/MRI scanner. The phantom was imaged at 9.7:1 contrast with standard spheres (diameter 10, 13, 17, 22, 28, 37 mm) and with custom spheres (diameter: 8.5, 11.5, 15, 25, 32.5, 44 mm) using a standardized methodology. Analysis was performed on a 30 min listmode data acquisition and on 6 realizations of 5 min from the listmode data. Images were reconstructed with the manufacturer provided iterative image reconstruction algorithms with and without point spread function (PSF) modeling. For both scanners, a post-reconstruction Gaussian filter of 3–7 mm in steps of 1 mm was applied. Attenuation correction was provided from a scaled computed tomography (CT) image of the phantom registered to the MR-based attenuation images and verified to align on the non-attenuation corrected PET images. For each of these image reconstruction parameter sets, contrast recovery coefficients (CRCs) were determined for the SUVmean, SUVmax and SUVpeak for each sphere. A hybrid metric combining the root-mean-squared discrepancy (RMSD) and the absolute CRC values was used to simultaneously optimize for best match in CRC between the two scanners while simultaneously weighting toward higher resolution reconstructions. The image reconstruction parameter set was identified as the best candidate reconstruction for each vendor for harmonized PET image reconstruction.

Results

The range of clinically relevant image reconstruction parameters demonstrated widely different quantitative performance across cameras. The best match of CRC curves was obtained at the lowest RMSD values with: for CRCmean, 2 iterations-7 mm filter on the GE Signa and 4 iterations-6 mm filter on the Siemens mMR, for CRCmax, 4 iterations-6 mm filter on the GE Signa, 4 iterations-5 mm filter on the Siemens mMR and for CRCpeak, 4 iterations-7 mm filter with PSF on the GE Signa and 4 iterations-7 mm filter on the Siemens mMR. Over all reconstructions, the RMSD between CRCs was 1.8%, 3.6% and 2.9% for CRC mean, max and peak, respectively. The solution of 2 iterations-3 mm on the GE Signa and 4 iterations-3 mm on Siemens mMR, both with PSF, led to simultaneous harmonization and with high CRC and low RMSD for CRC mean, max and peak with RMSD values of 2.8%, 5.8% and 3.2%, respectively.

Conclusions

For two commercially available PET/MRI scanners, user-selectable parameters that control iterative updates, image smoothing and PSF modeling provide a range of contrast recovery curves that allow harmonization in harmonization strategies of optimal match in CRC or high CRC values. This work demonstrates that nearly identical CRC curves can be obtained on different commercially available scanners by selecting appropriate image reconstruction parameters.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40658-021-00416-0.

SUBMITTER: Laforest R

PROVIDER: S-EPMC8571452 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Background During simultaneous PET/MRI, flexible MRI surface coils that lay on the patient are often omitted from PET attenuation correction processing, leading to quantification bias in PET images. Purpose To identify potential PET image quality improvement by using a recently developed lightweight MRI coil technology for the anterior array (AA) surface coil in both a phantom and in vivo study. Materials and Methods A phantom study and a prospective in vivo study were performed with a PET/CT scanner under three conditions: (a) no MRI surface coil (standard of reference), (b) traditional AA coil, and (c) lightweight AA coil. AA coils were not used in attenuation correction processing to emulate clinical PET/MRI. For the phantom study, PET images were reconstructed with and without time of flight (TOF) to assess quantification accuracy and uniformity. The in vivo study consisted of 10 participants (mean age, 66 years ± 10 [standard deviation]; six men) referred for a PET/CT oncologic examination who had undergone imaging between October 2019 and February 2020. Assessment of image quantification bias (defined as the standard error of the mean values) was conducted by comparing mean liver region of interest standardized uptake values with the no-coil standard of reference. A Wilcoxon signed-rank test was used to establish significance. Results For TOF and non-TOF, respectively, the phantom study revealed a mean PET quantification bias of -9.0% and -8.6% with the traditional AA coil and a mean PET quantification bias of -4.3% and -4.0% with the lightweight AA coil. The coefficients of variation reduced from 4.3% and 6.2% with the traditional AA coil to 2.1% and 2.7% with the lightweight AA coil, which demonstrated a homogeneity benefit from the lightweight coil that was greater with, versus without, TOF reconstruction. For the in vivo study, the mean liver standardized uptake value error was -5.9% with the traditional AA coil (P = .002 vs no coil) and -2.4% with the lightweight AA coil (P = .004 vs no coil). Conclusion The lightweight anterior array coil reduced PET image quantification bias by more than 50% compared with the traditional coil. Using the lightweight coil and performing time of flight-based reconstruction each reduced the variation of error. © RSNA, 2020 Online supplemental material is available for this article.

Project description:Respiratory and cardiac motion is the most serious limitation to whole-body PET, resulting in spatial resolution close to 1 cm. Furthermore, motion-induced inconsistencies in the attenuation measurements often lead to significant artifacts in the reconstructed images. Gating can remove motion artifacts at the cost of increased noise. This paper presents an approach to respiratory motion correction using simultaneous PET/MRI to demonstrate initial results in phantoms, rabbits, and nonhuman primates and discusses the prospects for clinical application.Studies with a deformable phantom, a free-breathing primate, and rabbits implanted with radioactive beads were performed with simultaneous PET/MRI. Motion fields were estimated from concurrently acquired tagged MR images using 2 B-spline nonrigid image registration methods and incorporated into a PET list-mode ordered-subsets expectation maximization algorithm. Using the measured motion fields to transform both the emission data and the attenuation data, we could use all the coincidence data to reconstruct any phase of the respiratory cycle. We compared the resulting SNR and the channelized Hotelling observer (CHO) detection signal-to-noise ratio (SNR) in the motion-corrected reconstruction with the results obtained from standard gating and uncorrected studies.Motion correction virtually eliminated motion blur without reducing SNR, yielding images with SNR comparable to those obtained by gating with 5-8 times longer acquisitions in all studies. The CHO study in dynamic phantoms demonstrated a significant improvement (166%-276%) in lesion detection SNR with MRI-based motion correction as compared with gating (P < 0.001). This improvement was 43%-92% for large motion compared with lesion detection without motion correction (P < 0.001). CHO SNR in the rabbit studies confirmed these results.Tagged MRI motion correction in simultaneous PET/MRI significantly improves lesion detection compared with respiratory gating and no motion correction while reducing radiation dose. In vivo primate and rabbit studies confirmed the improvement in PET image quality and provide the rationale for evaluation in simultaneous whole-body PET/MRI clinical studies.

Project description:Noninvasive methods are needed to explore the heterogeneous tumor microenvironment and its modulation by therapy. Hybrid PET/MRI systems are being developed for small-animal and clinical use. The advantage of these integrated systems depends on their ability to provide MR images that are spatially coincident with simultaneously acquired PET images, allowing combined functional MRI and PET studies of intratissue heterogeneity. Although much effort has been devoted to developing this new technology, the issue of quantitative and spatial fidelity of PET images from hybrid PET/MRI systems to the tissues imaged has received little attention. Here, we evaluated the ability of a first-generation, small-animal MRI-compatible PET scanner to accurately depict heterogeneous patterns of radiotracer uptake in tumors.Quantitative imaging characteristics of the MRI-compatible PET (PET/MRI) scanner were evaluated with phantoms using calibration coefficients derived from a mouse-sized linearity phantom. PET performance was compared with a commercial small-animal PET system and autoradiography in tumor-bearing mice. Pixel and structure-based similarity metrics were used to evaluate image concordance among modalities. Feasibility of simultaneous PET/MRI functional imaging of tumors was explored by following (64)Cu-labeled antibody uptake in relation to diffusion MRI using cooccurrence matrix analysis.The PET/MRI scanner showed stable and linear response. Activity concentration recovery values (measured and true activity concentration) calculated for 4-mm-diameter rods within linearity and uniform activity rod phantoms were near unity (0.97 ± 0.06 and 1.03 ± 0.03, respectively). Intratumoral uptake patterns for both (18)F-FDG and a (64)Cu-antibody acquired using the PET/MRI scanner and small-animal PET were highly correlated with autoradiography (r > 0.99) and with each other (r = 0.97 ± 0.01). On the basis of these data, we performed a preliminary study comparing diffusion MRI and radiolabeled antibody uptake patterns over time and visualized movement of antibodies from the vascular space into the tumor mass.The MRI-compatible PET scanner provided tumor images that were quantitatively accurate and spatially concordant with autoradiography and the small-animal PET examination. Cooccurrence matrix approaches enabled effective analysis of multimodal image sets. These observations confirm the ability of the current simultaneous PET/MRI system to provide accurate observations of intratumoral function and serve as a benchmark for future evaluations of hybrid instrumentation.

Project description:PURPOSE:Numerous image reconstruction methodologies for positron emission tomography (PET) have been developed that incorporate magnetic resonance (MR) imaging structural information, producing reconstructed images with improved suppression of noise and reduced partial volume effects. However, the influence of MR structural information also increases the possibility of suppression or bias of structures present only in the PET data (PET-unique regions). To address this, further developments for MR-informed methods have been proposed, for example, through inclusion of the current reconstructed PET image, alongside the MR image, in the iterative reconstruction process. In this present work, a number of kernel and maximum a posteriori (MAP) methodologies are compared, with the aim of identifying methods that enable a favorable trade-off between the suppression of noise and the retention of unique features present in the PET data. METHODS:The reconstruction methods investigated were: the MR-informed conventional and spatially compact kernel methods, referred to as KEM and KEM largest value sparsification (LVS) respectively; the MR-informed Bowsher and Gaussian MR-guided MAP methods; and the PET-MR-informed hybrid kernel and anato-functional MAP methods. The trade-off between improving the reconstruction of the whole brain region and the PET-unique regions was investigated for all methods in comparison with postsmoothed maximum likelihood expectation maximization (MLEM), evaluated in terms of structural similarity index (SSIM), normalized root mean square error (NRMSE), bias, and standard deviation. Both simulated BrainWeb (10 noise realizations) and real [18 F] fluorodeoxyglucose (FDG) three-dimensional datasets were used. The real [18 F]FDG dataset was augmented with simulated tumors to allow comparison of the reconstruction methodologies for the case of known regions of PET-MR discrepancy and evaluated at full counts (100%) and at a reduced (10%) count level. RESULTS:For the high-count simulated and real data studies, the anato-functional MAP method performed better than the other methods under investigation (MR-informed, PET-MR-informed and postsmoothed MLEM), in terms of achieving the best trade-off for the reconstruction of the whole brain and PET-unique regions, assessed in terms of the SSIM, NRMSE, and bias vs standard deviation. The inclusion of PET information in the anato-functional MAP method enables the reconstruction of PET-unique regions to attain similarly low levels of bias as unsmoothed MLEM, while moderately improving the whole brain image quality for low levels of regularization. However, for low count simulated datasets the anato-functional MAP method performs poorly, due to the inclusion of noisy PET information in the regularization term. For the low counts simulated dataset, KEM LVS and to a lesser extent, HKEM performed better than the other methods under investigation in terms of achieving the best trade-off for the reconstruction of the whole brain and PET-unique regions, assessed in terms of the SSIM, NRMSE, and bias vs standard deviation. CONCLUSION:For the reconstruction of noisy data, multiple MR-informed methods produce favorable whole brain vs PET-unique region trade-off in terms of the image quality metrics of SSIM and NRMSE, comfortably outperforming the whole image denoising of postsmoothed MLEM.