Improved quantitation for PET/CT image reconstruction with system modeling and anatomical priors

Enhancing Amyloid PET Quantification: MRI-Guided Super-Resolution Using Latent Diffusion Models

Amyloid PET imaging plays a crucial role in the diagnosis and research of Alzheimer's disease (AD), allowing non-invasive detection of amyloid-β plaques in the brain. However, the low spatial resolution of PET scans limits the accurate quantification of amyloid deposition due to partial volume effects (PVE). In this study, we propose a novel approach to addressing PVE using a latent diffusion model for resolution recovery (LDM-RR) of PET imaging. We leverage a synthetic data generation pipeline to create high-resolution PET digital phantoms for model training. The proposed LDM-RR model incorporates a weighted combination of L1, L2, and MS-SSIM losses at both noise and image scales to enhance MRI-guided reconstruction. We evaluated the model's performance in improving statistical power for detecting longitudinal changes and enhancing agreement between amyloid PET measurements from different tracers. The results demonstrate that the LDM-RR approach significantly improves PET quantification accuracy, reduces inter-tracer variability, and enhances the detection of subtle changes in amyloid deposition over time. We show that deep learning has the potential to improve PET quantification in AD, effectively contributing to the early detection and monitoring of disease progression.

Impact of irregular waveforms on data-driven respiratory gated PET/CT images processed using MotionFree algorithm

OBJECTIVES

MotionFree® (AMF) is a data-driven respiratory gating (DDG) algorithm for image processing that has recently been introduced into clinical practice. The present study aimed to verify the accuracy of respiratory waveform and the effects of normal and irregular respiratory motions using AMF with the DDG algorithm.

METHODS

We used a NEMA IEC body phantom comprising six spheres (37-, 28-, 22-, 17-, 13-, and 10 mm diameter) containing 18F. The sphere-to-background ratio was 4:1 (21.2 and 5.3 kBq/mL). We acquired PET/CT images from a stationary or moving phantom placed on a custom-designed motion platform. Respiratory motions were reproduced based on normal (sinusoidal or expiratory-paused waveforms) and irregular (changed amplitude or shifted baseline waveforms) movements. The "width" parameters in AMF were set at 10-60% and extracted data during the expiratory phases of each waveform. We verified the accuracy of the derived waveforms by comparing those input from the motion platform and output determined using AMF. Quantitative accuracy was evaluated as recovery coefficients (RCs), improvement rate, and %change that were calculated based on sphere diameter or width. We evaluated statistical differences in activity concentrations of each sphere between normal and irregular waveforms.

RESULTS

Respiratory waveforms derived from AMF were almost identical to the input waveforms on the motion platform. Although the RCs in each sphere for expiratory-paused and ideal stationary waveforms were almost identical, RCs except the expiratory-paused waveform were lower than those for the stationary waveform. The improvement rate decreased more for the irregular, than the normal waveforms with AMF in smaller spheres. The %change was improved by decreasing the width of waveforms with a shifted baseline. Activity concentrations significantly differed between normal waveforms and those with a shifted baseline in spheres < 28 mm.

CONCLUSIONS

The PET images using AMF with the DDG algorithm provided the precise waveform of respiratory motions and the improvement of quantitative accuracy in the four types of respiratory waveforms. The improvement rate was the most obvious in expiratory-paused waveforms, and the most subtle in those with a shifted baseline. Optimizing the width parameter in irregular waveform will benefit patients who breathe like the waveform with the shifted baseline.

SPECT reconstruction with a trained regularizer using CT-side information: Application to 177Lu SPECT imaging

Improving low-count SPECT can shorten scans and support pre-therapy theranostic imaging for dosimetry-based treatment planning, especially with radionuclides like 177Lu known for low photon yields. Conventional methods often underperform in low-count settings, highlighting the need for trained regularization in model-based image reconstruction. This paper introduces a trained regularizer for SPECT reconstruction that leverages segmentation based on CT imaging. The regularizer incorporates CT-side information via a segmentation mask from a pre-trained network (nnUNet). In this proof-of-concept study, we used patient studies with 177Lu DOTATATE to train and tested with phantom and patient datasets, simulating pre-therapy imaging conditions. Our results show that the proposed method outperforms both standard unregularized EM algorithms and conventional regularization with CT-side information. Specifically, our method achieved marked improvements in activity quantification, noise reduction, and root mean square error. The enhanced low-count SPECT approach has promising implications for theranostic imaging, post-therapy imaging, whole body SPECT, and reducing SPECT acquisition times.

Investigation and optimization of PET-guided SPECT reconstructions for improved radionuclide therapy dosimetry estimates

Introduction

To investigate and optimize the SPECTRE (Single Photon Emission Computed Theranostic REconstruction) reconstruction approach, using the hybrid kernelised expectation maximization (HKEM) algorithm implemented in the software for tomographic image reconstruction (STIR) software library, and to demonstrate the feasibility of performing algorithm exploration and optimization in 2D. Optimal SPECTRE parameters were investigated for the purpose of improving SPECT-based radionuclide therapy (RNT) dosimetry estimates.

Materials and Methods

Using the NEMA IEC body phantom as the test object, SPECT data were simulated to model an early and late imaging time point following a typical therapeutic dose of 8 GBq of 177Lu. A theranostic 68Ga PET-prior was simulated for the SPECTRE reconstructions. The HKEM algorithm parameter space was investigated for SPECT-unique and PET-SPECT mutual features to characterize optimal SPECTRE parameters for the simulated data. Mean and maximum bias, coefficient of variation (COV %), recovery, SNR and root-mean-square error (RMSE) were used to facilitate comparisons between SPECTRE reconstructions and OSEM reconstructions with resolution modelling (OSEM_RM). 2D reconstructions were compared to those performed in 3D in order to evaluate the utility of accelerated algorithm optimization in 2D. Segmentation accuracy was evaluated using a 42% fixed threshold (FT) on the 3D reconstructed data.

Results

SPECTRE parameters that demonstrated improved image quality and quantitative accuracy were determined through investigation of the HKEM algorithm parameter space. OSEM_RM and SPECTRE reconstructions performed in 2D and 3D were qualitatively and quantitatively similar, with SPECTRE showing an average reduction in background COV % by a factor of 2.7 and 3.3 for the 2D case and 3D case respectively. The 42% FT analysis produced an average % volume difference from ground truth of 158% and 26%, for the OSEM_RM and SPECTRE reconstructions, respectively.

Conclusions

The SPECTRE reconstruction approach demonstrates significant potential for improved SPECT image quality, leading to more accurate RNT dosimetry estimates when conventional segmentation methods are used. Exploration and optimization of SPECTRE benefited from both fast reconstruction times afforded by first considering the 2D case. This is the first in-depth exploration of the SPECTRE reconstruction approach, and as such, it reveals several insights for reconstructing SPECT data using PET side information.

DeepPVC: prediction of a partial volume-corrected map for brain positron emission tomography studies via a deep convolutional neural network

Background

Partial volume correction with anatomical magnetic resonance (MR) images (MR-PVC) is useful for accurately quantifying tracer uptake on brain positron emission tomography (PET) images. However, MR segmentation processes for MR-PVC are time-consuming and prevent the widespread clinical use of MR-PVC. Here, we aimed to develop a deep learning model to directly predict PV-corrected maps from PET and MR images, ultimately improving the MR-PVC throughput.

Methods

We used MR T1-weighted and [¹¹C]PiB PET images as input data from 192 participants from the Alzheimer’s Disease Neuroimaging Initiative database. We calculated PV-corrected maps as the training target using the region-based voxel-wise PVC method. Two-dimensional U-Net model was trained and validated by sixfold cross-validation with the dataset from the 156 participants, and then tested using MR T1-weighted and [¹¹C]PiB PET images from 36 participants acquired at sites other than the training dataset. We calculated the structural similarity index (SSIM) of the PV-corrected maps and intraclass correlation (ICC) of the PV-corrected standardized uptake value between the region-based voxel-wise (RBV) PVC and deepPVC as indicators for validation and testing.

Results

A high SSIM (0.884 ± 0.021) and ICC (0.921 ± 0.042) were observed in the validation and test data (SSIM, 0.876 ± 0.028; ICC, 0.894 ± 0.051). The computation time required to predict a PV-corrected map for a participant (48 s without a graphics processing unit) was much shorter than that for the RBV PVC and MR segmentation processes.

Conclusion

These results suggest that the deepPVC model directly predicts PV-corrected maps from MR and PET images and improves the throughput of MR-PVC by skipping the MR segmentation processes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40658-022-00478-8.

Theranostic SPECT reconstruction for improved resolution: application to radionuclide therapy dosimetry

Background

SPECT-derived dose estimates in tissues of diameter less than 3× system resolution are subject to significant losses due to the limited spatial resolution of the gamma camera. Incorporating resolution modelling (RM) into the SPECT reconstruction has been proposed as a possible solution; however, the images produced are prone to noise amplification and Gibbs artefacts. We propose a novel approach to SPECT reconstruction in a theranostic setting, which we term SPECTRE (single photon emission computed theranostic reconstruction); using a diagnostic PET image, with its superior resolution, to guide the SPECT reconstruction of the therapeutic equivalent. This report demonstrates a proof in principle of this approach.

Methods

We have employed the hybrid kernelised expectation maximisation (HKEM) algorithm implemented in STIR, with the aim of producing SPECT images with PET-equivalent resolution. We demonstrate its application in both a dual ⁶⁸Ga/¹⁷⁷Lu IEC phantom study and a clinical example using ⁶⁴Cu/⁶⁷Cu.

Results

SPECTRE is shown to produce images comparable in accuracy and recovery to PET with minimal introduction of artefacts and amplification of noise.

Conclusion

The SPECTRE approach to image reconstruction shows improved quantitative accuracy with a reduction in noise amplification. SPECTRE shows great promise as a method of improving SPECT radioactivity concentrations, directly leading to more accurate dosimetry estimates in small structures and target lesions. Further investigation and optimisation of the algorithm parameters is needed before this reconstruction method can be utilised in a clinical setting.

Iterative framework for image registration and partial volume correction in brain positron emission tomography

Imprecise registration between positron emission tomography (PET) and anatomical magnetic resonance (MR) images is a critical source of error in MR imaging-guided partial volume correction (MR-PVC). Here, we propose a novel framework for image registration and partial volume correction, which we term PVC-optimized registration (PoR), to address imprecise registration. The PoR framework iterates PVC and registration between uncorrected PET and smoothed PV-corrected images to obtain precise registration. We applied PoR to the [¹¹C]PiB PET data of 92 participants obtained from the Alzheimer’s Disease Neuroimaging Initiative database and compared the registration results, PV-corrected standardized uptake value (SUV) and its ratio to the cerebellum (SUVR), and intra-region coefficient of variation (CoV) between PoR and conventional registration. Significant differences in registration of as much as 2.74 mm and 3.02° were observed between the two methods (effect size <  − 0.8 or > 0.8), which resulted in considerable SUVR differences throughout the brain, reaching a maximal difference of 62.3% in the sensory motor cortex. Intra-region CoV was significantly reduced using the PoR throughout the brain. These results suggest that PoR reduces error as a result of imprecise registration in PVC and is a useful method for accurately quantifying the amyloid burden in PET.

Concurrent imaging of vascularization and metabolism in a mouse model of paraganglioma under anti-angiogenic treatment

Rationale: Deregulation of metabolism and induction of vascularization are major hallmarks of cancer. Using a new multimodal preclinical imaging instrument, we explored a sequence of events leading to sunitinib-induced resistance in a murine model of paraganglioma (PGL) invalidated for the expression of succinate dehydrogenase subunit B (Sdhb^-/-).

Methods: Two groups of Sdhb^-/- tumors bearing mice were treated with sunitinib (6 weeks) or vehicle (3 weeks). Concurrent Positron Emission Tomography (PET) with 2′ -deoxy-2′-[¹⁸F]fluoro-D-glucose (FDG), Computed Tomography (CT) and Ultrafast Ultrasound Imaging (UUI) imaging sessions were performed once a week and ex vivo samples were analyzed by western blots and histology.

Results: PET-CT-UUI enabled to detect a rapid growth of Sdhb^-/- tumors with increased glycolysis and vascular development. Sunitinib treatment prevented tumor growth, vessel development and reduced FDG uptake at week 1 and 2 (W1-2). Thereafter, imaging revealed tumor escape from sunitinib treatment: FDG uptake in tumors increased at W3, followed by tumor growth and vessel development at W4-5. Perfused vessels were preferentially distributed in the hypermetabolic regions of the tumors and the perfused volume increased during escape from sunitinib treatment. Finally, initial changes in total lesion glycolysis and maximum vessel length at W1 were predictive of resistance to sunitinib.

Conclusion: These results demonstrate an adaptive resistance of Sdhb^-/- tumors to six weeks of sunitinib treatment. Early metabolic changes and delayed vessel architecture changes were detectable and predictable in vivo early during anti-angiogenic treatment. Simultaneous metabolic, anatomical and functional imaging can monitor precisely the effects of anti-angiogenic treatment of tumors.

Analysis of partial volume correction on quantification and regional heterogeneity in cardiac PET

BACKGROUND

The partial volume correction (PVC) of cardiac PET datasets using anatomical side information during reconstruction is appealing but not straightforward. Other techniques, which do not make use of additional anatomical information, could be equally effective in improving the reconstructed myocardial activity.

METHODS

Resolution modeling in combination with different noise suppressing priors was evaluated as a means to perform PVC. Anatomical priors based on a high-resolution CT are compared to non-anatomical, edge-preserving priors (relative difference and total variation prior). The study is conducted on ex vivo datasets from ovine hearts. A simulation study additionally clarifies the relationship between prior effectiveness and myocardial wall thickness.

RESULTS

Simple resolution modeling during data reconstruction resulted in over- and underestimation of activity, which hampers the absolute left ventricular quantification when compared to the ground truth. Both the edge-preserving and the anatomy-based PVC techniques improve the absolute quantification, with comparable results (Student t-test, P = .17). The relative tracer distribution was preserved with any reconstruction technique (repeated ANOVA, P = .98).

CONCLUSIONS

The use of edge-preserving priors emerged as optimal choice for quantification of tracer uptake in the left ventricular wall of the available datasets. Anatomical priors visually outperformed edge-preserving priors when the thinnest structures were of interest.

Partial volume correction for PET quantification and its impact on brain network in Alzheimer's disease

Amyloid positron emission tomography (PET) imaging is a valuable tool for research and diagnosis in Alzheimer’s disease (AD). Partial volume effects caused by the limited spatial resolution of PET scanners degrades the quantitative accuracy of PET image. In this study, we have applied a method to evaluate the impact of a joint-entropy based partial volume correction (PVC) technique on brain networks learned from a clinical dataset of AV-45 PET image and compare network properties of both uncorrected and corrected image-based brain networks. We also analyzed the region-wise SUVRs of both uncorrected and corrected images. We further performed classification tests on different groups using the same set of algorithms with same parameter settings. PVC has sometimes been avoided due to increased noise sensitivity in image registration and segmentation, however, our results indicate that appropriate PVC may enhance the brain network structure analysis for AD progression and improve classification performance.

Analysis of partial volume correction on quantification and regional heterogeneity in cardiac PET

BACKGROUND

The partial volume correction (PVC) of cardiac PET datasets using anatomical side information during reconstruction is appealing but not straightforward. Other techniques, which do not make use of additional anatomical information, could be equally effective in improving the reconstructed myocardial activity.

METHODS

Resolution modeling in combination with different noise suppressing priors was evaluated as a means to perform PVC. Anatomical priors based on a high-resolution CT are compared to non-anatomical, edge-preserving priors (relative difference and total variation prior). The study is conducted on ex vivo datasets from ovine hearts. A simulation study additionally clarifies the relationship between prior effectiveness and myocardial wall thickness.

RESULTS

Simple resolution modeling during data reconstruction resulted in over- and underestimation of activity, which hampers the absolute left ventricular quantification when compared to the ground truth. Both the edge-preserving and the anatomy-based PVC techniques improve the absolute quantification, with comparable results (Student t-test, P = .17). The relative tracer distribution was preserved with any reconstruction technique (repeated ANOVA, P = .98).

CONCLUSIONS

The use of edge-preserving priors emerged as optimal choice for quantification of tracer uptake in the left ventricular wall of the available datasets. Anatomical priors visually outperformed edge-preserving priors when the thinnest structures were of interest.

PETPVE12: an SPM toolbox for Partial Volume Effects correction in brain PET - Application to amyloid imaging with AV45-PET

Positron emission tomography (PET) allows detecting molecular brain changes in vivo. However, the accuracy of PET is limited by partial volume effects (PVE) that affects quantitative analysis and visual interpretation of the images. Although PVE-correction methods have been shown to effectively increase the correspondence of the measured signal with the true regional tracer uptake, these procedures are still not commonly applied, neither in clinical nor in research settings. Here, we present an implementation of well validated PVE-correction procedures as a SPM toolbox, PETPVE12, for automated processing. We demonstrate its utility by a comprehensive analysis of the effects of PVE-correction on amyloid-sensitive AV45-PET data from 85 patients with Alzheimer's disease (AD) and 179 cognitively normal (CN) elderly. Effects of PVE-correction on global cortical standard uptake value ratios (SUVR) and the power of diagnostic group separation were assessed for the region-wise geometric transfer matrix method (PVEc-GTM), as well as for the 3-compartmental voxel-wise "Müller-Gärtner" method (PVEc-MG). Both PVE-correction methods resulted in decreased global cortical SUVRs in the low to middle range of SUVR values, and in increased global cortical SUVRs at the high values. As a consequence, average SUVR of the CN group was reduced, whereas average SUVR of the AD group was increased by PVE-correction. These effects were also reflected in increased accuracies of group discrimination after PVEc-GTM (AUC=0.86) and PVEc-MG (AUC=0.89) compared to standard non-corrected SUVR (AUC=0.84). Voxel-wise analyses of PVEc-MG corrected data also demonstrated improved detection of regionally increased AV45 SUVR values in AD patients. These findings complement the growing evidence for a beneficial effect of PVE-correction in quantitative analysis of amyloid-sensitive PET data. The novel PETPVE12 toolbox significantly facilitates the application of PVE-correction, particularly within SPM-based processing pipelines. This is expected to foster the use of PVE-correction in brain PET for more widespread use. The toolbox is freely available at http://www.fil.ion.ucl.ac.uk/spm/ext/#PETPVE12.

Anatomically-aided PET reconstruction using the kernel method

This paper extends the kernel method that was proposed previously for dynamic PET reconstruction, to incorporate anatomical side information into the PET reconstruction model. In contrast to existing methods that incorporate anatomical information using a penalized likelihood framework, the proposed method incorporates this information in the simpler maximum likelihood (ML) formulation and is amenable to ordered subsets. The new method also does not require any segmentation of the anatomical image to obtain edge information. We compare the kernel method with the Bowsher method for anatomically-aided PET image reconstruction through a simulated data set. Computer simulations demonstrate that the kernel method offers advantages over the Bowsher method in region of interest quantification. Additionally the kernel method is applied to a 3D patient data set. The kernel method results in reduced noise at a matched contrast level compared with the conventional ML expectation maximization algorithm.

Impact of spillover from white matter by partial volume effect on quantification of amyloid deposition with [11C]PiB PET

High non-specific uptake of [¹¹C]Pittsburgh compound B ([¹¹C]PiB) in white matter and signal spillover from white matter, due to partial volume effects, confound radioactivity measured in positron emission tomography (PET) with [¹¹C]PiB. We aimed to reveal the partial volume effect in absolute values of kinetic parameters for [¹¹C]PiB, in terms of spillover from white matter. Dynamic data acquired in [¹¹C]PiB PET scans with five healthy volunteers and eight patients with Alzheimer's disease were corrected with region-based and voxel-based partial volume corrections. Binding potential (BP_ND) was estimated using the two-tissue compartment model analysis with a plasma input function. Partial volume corrections significantly decreased cortical BP_ND values. The degree of decrease in healthy volunteers (-52.7±5.8%) was larger than that in Alzheimer's disease patients (-11.9±4.2%). The simulation demonstrated that white matter spillover signals due to the partial volume effect resulted in an overestimation of cortical BP_ND, with a greater degree of overestimation for lower BP_ND values. Thus, an overestimation due to partial volume effects is more severe in healthy volunteers than in Alzheimer's disease patients. Partial volume corrections may be useful for accurately quantifying Aβ deposition in cortical regions.

Enhancement of molecular sensitivity in positron emission tomography with quantum correlation of γ-ray photons

Enhancement of molecular sensitivity in positron emission tomography (PET) has long been discussed with respect to imaging instrumentation and algorithms for data treatment. Here, the molecular sensitivity in PET is discussed on the basis of 2-dimensional coincident measurements of 511 keV γ ray photons resultant from two-photon annihilation. Introduction of an additional selection window based on the energy sum and difference of the coincidently measured γ ray photons, without any significant instrumental and algorithmic changes, showed an improvement in the signal-to-noise ratio (SNR) by an order of magnitude. Improvement of performance characteristics in the PET imaging system was demonstrated by an increase in the noise equivalent count rate (NECR) which takes both the SNR and the detection efficiency into consideration. A further improvement of both the SNR and the NECR is expected for the present system in real clinical and in-vivo environments, where much stronger positron sources are employed.

Anatomy-guided brain PET imaging incorporating a joint prior model

We proposed a maximum a posterior (MAP) framework for incorporating information from co-registered anatomical images into PET image reconstruction through a novel anato-functional joint prior. The characteristic of the utilized hyperbolic potential function is determinate by the voxel intensity differences within the anatomical image, while the penalization is computed based on voxel intensity differences in reconstructed PET images. Using realistic simulated (18)FDG PET scan data, we optimized the performance of the proposed MAP reconstruction with the joint prior (JP-MAP) and compared its performance with conventional 3D MLEM and 3D MAP reconstructions. The proposed JP-MAP reconstruction algorithm resulted in quantitatively enhanced reconstructed images, as demonstrated in extensive FDG PET simulation study. The proposed method was also tested on a 20 min Florbetapir patient study performed on the high-resolution research tomograph. It was shown to outperform conventional methods in visual as well as quantitative accuracy assessment (in terms of regional noise versus activity value performance). The JP-MAP method was also compared with another MR-guided MAP reconstruction method, utilizing the Bowsher prior and was seen to result in some quantitative enhancements, especially in the case of MR-PET mis-registrations, and a definitive improvement in computational performance.

Fat-constrained 18F-FDG PET reconstruction in hybrid PET/MR imaging

Fusion of information from PET and MR imaging can increase the diagnostic value of both modalities. This work sought to improve (18)F FDG PET image quality by using MR Dixon fat-constrained images to constrain PET image reconstruction to low-fat regions, with the working hypothesis that fatty tissue metabolism is low in glucose consumption.

METHODS

A novel constrained PET reconstruction algorithm was implemented via a modification of the system matrix in list-mode time-of-flight ordered-subsets expectation maximization reconstruction, similar to the way time-of-flight weighting is incorporated. To demonstrate its use in PET/MR imaging, we modeled a constraint based on fat/water-separating Dixon MR images that shift activity away from regions of fat tissue during PET image reconstruction. PET and MR imaging scans of a modified National Electrical Manufacturers Association/International Electrotechnical Commission body phantom simulating body fat/water composition and in vivo experiments on 2 oncology patients were performed on a commercial time-of-flight PET/MR imaging system.

RESULTS

Fat-constrained PET reconstruction visibly and quantitatively increased resolution and contrast between high-uptake and fatty-tissue regions without significantly affecting the images in nonfat regions.

CONCLUSION

The incorporation of MR tissue information, such as fat, in image reconstruction can improve the quality of PET images. The combination of a variety of potential other MR tissue characteristics with PET represents a further justification for merging MR data with PET data in hybrid systems.

Management of respiratory motion in PET/computed tomography: the state of the art

Combined PET/computed tomography (CT) is of value in cancer diagnosis, follow-up, and treatment planning. For cancers located in the thorax or abdomen, the patient’s breathing causes artifacts and errors in PET and CT images. Many different approaches for artifact avoidance or correction have been developed; most are based on gated acquisition and synchronization between the respiratory signal and PET acquisition. The respiratory signal is usually produced by an external sensor that tracks a physiological characteristic related to the patient’s breathing. Respiratory gating is a compensation technique in which time or amplitude binning is used to exclude the motion in reconstructed PET images. Although this technique is performed in routine clinical practice, it fails to adequately correct for respiratory motion because each gate can mix several tissue positions. Researchers have suggested either selecting PET events from gated acquisitions or performing several PET acquisitions (corresponding to a breath-hold CT position). However, the PET acquisition time must be increased if adequate counting statistics are to be obtained in the different gates after binning. Hence, other researchers have assessed correction techniques that take account of all the counting statistics (without increasing the acquisition duration) and integrate motion information before, during, or after the reconstruction process. Here, we provide an overview of how motion is managed to overcome respiratory motion in PET/CT images.

Region-Based Partial Volume Correction Techniques for PET Imaging: Sinogram Implementation and Robustness

Background/Purpose. Limited spatial resolution of positron emission tomography (PET) requires partial volume correction (PVC). Region-based PVC methods are based on geometric transfer matrix implemented either in image-space (GTM) or sinogram-space (GTMo), both with similar performance. Although GTMo is slower, it more closely simulates the 3D PET image acquisition, accounts for local variations of point spread function, and can be implemented for iterative reconstructions. A recent image-based symmetric GTM (sGTM) has shown improvement in noise characteristics and robustness to misregistration over GTM. This study implements the sGTM method in sinogram space (sGTMo), validates it, and evaluates its performance. Methods. Two 3D sphere and brain digital phantoms and a physical sphere phantom were used. All four region-based PVC methods (GTM, GTMo, sGTM, and sGTMo) were implemented and their performance was evaluated. Results. All four PVC methods had similar accuracies. Both noise propagation and robustness of the sGTMo method were similar to those of sGTM method while they were better than those of GTMo method especially for smaller objects. Conclusion. The sGTMo was implemented and validated. The performance of the sGTMo in terms of noise characteristics and robustness to misregistration is similar to that of the sGTM method and improved compared to the GTMo method.

Magnetic resonance-guided positron emission tomography image reconstruction

The resolution of positron emission tomography (PET) images is limited by the physics of positron-electron annihilation and instrumentation for photon coincidence detection. Model-based methods that incorporate accurate physical and statistical models have produced significant improvements in reconstructed image quality when compared with filtered backprojection reconstruction methods. However, it has often been suggested that by incorporating anatomical information, the resolution and noise properties of PET images could be further improved, leading to better quantitation or lesion detection. With the recent development of combined MR-PET scanners, we can now collect intrinsically coregistered magnetic resonance images. It is therefore possible to routinely make use of anatomical information in PET reconstruction, provided appropriate methods are available. In this article, we review research efforts over the past 20 years to develop these methods. We discuss approaches based on the use of both Markov random field priors and joint information or entropy measures. The general framework for these methods is described, and their performance and longer-term potential and limitations are discussed.

Performance Evaluation of Small Animal PET Scanners With Different System Designs

This study evaluated the image quality metrics of small animal PET scanners based upon measured single detector module positioning performance. A semi-analytical approach was developed to study PET scanner performance in the scenario of multiple realizations. Positron range blurring, scanner system response function (SRF) and statistical noise were included in the modeling procedure. The scanner sensitivity map was included in the system matrix during maximum likelihood expectation maximization (MLEM) reconstruction. Several image quality metrics were evaluated for octagonal ring PET scanners consisting of continuous miniature crystal element (cMiCE) detector modules with varying designs. These designs included 8 mm and 15 mm thick crystal detectors using conventional readout with the photosensors on the exit surface of the crystal and a 15 mm thick crystal detector using our proposed sensor-on-the-entrance (SES) design. For the conventional readout design, the results showed that there was a tradeoff between bias and variance with crystal thickness. The 15 mm crystal detector had better detection task performance, while quantitation task performance was degraded. On the other hand, our SES detector had similar detection efficiency as the conventional design using a 15 mm thick crystal and had similar intrinsic spatial resolution as the conventional design using an 8 mm thick crystal. The end result was that by using the SES design, one could improve scanner quantitation task performance without sacrificing detection task performance.

Modeling and Simulation of 4D PET-CT and PET-MR Images

The driving force in the research and development of new hybrid PET-CT/MR imaging scanners is the production of images with optimum quality, accuracy, and resolution. However, the acquisition process is limited by several factors. Key issues are the respiratory and cardiac motion artifacts that occur during an imaging session. In this article the necessary tools for modeling and simulation of realistic high-resolution four-dimensional PET-CT and PET-MR imaging data are described. Beyond the need for four-dimensional simulations, accurate modeling of the acquisition process can be included within the reconstruction algorithms assisting in the improvement of image quality and accuracy of estimation of physiologic parameters from four-dimensional hybrid PET imaging.

Image guided biodistribution and pharmacokinetic studies of theranostics

Image guided technique is playing an increasingly important role in the investigation of the biodistribution and pharmacokinetics of drugs or drug delivery systems in various diseases, especially cancers. Besides anatomical imaging modalities such as computed tomography (CT), magnetic resonance imaging (MRI), molecular imaging strategy including optical imaging, positron emission tomography (PET) and single-photon emission computed tomography (SPECT) will facilitate the localization and quantization of radioisotope or optical probe labeled nanoparticle delivery systems in the category of theranostics. The quantitative measurement of the bio-distribution and pharmacokinetics of theranostics in the fields of new drug/probe development, diagnosis and treatment process monitoring as well as tracking the brain-blood-barrier (BBB) breaking through by high sensitive imaging method, and the applications of the representative imaging modalities are summarized in this review.

A method for model-free partial volume correction in oncological PET

Background

As is well known, limited spatial resolution leads to partial volume effects (PVE) and consequently to limited signal recovery. Determination of the mean activity concentration of a target structure is thus compromised even at target sizes much larger than the reconstructed spatial resolution. This leads to serious size-dependent underestimates of true signal intensity in hot spot imaging. For quantitative PET in general and in the context of therapy assessment in particular it is, therefore, mandatory to perform an adequate partial volume correction (PVC). The goal of our work was to develop and to validate a model-free PVC algorithm for hot spot imaging.

Methods

The algorithm proceeds in two automated steps. Step 1: estimation of the actual object boundary with a threshold based method and determination of the total activity A measured within the enclosed volume V. Step 2: determination of the activity fraction B, which is measured outside the object due to the partial volume effect (spill-out). The PVE corrected mean value is then given by C_mean = (A+B)/V. For validation simulated tumours were used which were derived from real patient data (liver metastases of a colorectal carcinoma and head and neck cancer, respectively). The simulated tumours have characteristics (regarding tumour shape, contrast, noise, etc.) which are very similar to those of the underlying patient data, but the boundaries and tracer accumulation are exactly known. The PVE corrected mean values of 37 simulated tumours were determined and compared with the true mean values.

Results

For the investigated simulated data the proposed approach yields PVE corrected mean values which agree very well with the true values (mean deviation (± s.d.): (−0.8±2.5)%).

Conclusions

The described method enables accurate quantitative partial volume correction in oncological hot spot imaging.

PET image reconstruction with anatomical edge guided level set prior

Acquiring both anatomical and functional images during one scan, PET/CT systems improve the ability to detect and localize abnormal uptakes. In addition, CT images provide anatomical boundary information that can be used to regularize positron emission tomography (PET) images. Here we propose a new approach to maximum a posteriori reconstruction of PET images with a level set prior guided by anatomical edges. The image prior models both the smoothness of PET images and the similarity between functional boundaries in PET and anatomical boundaries in CT. Level set functions (LSFs) are used to represent smooth and closed functional boundaries. The proposed method does not assume an exact match between PET and CT boundaries. Instead, it encourages similarity between the two boundaries, while allowing different region definition in PET images to accommodate possible signal and position mismatch between functional and anatomical images. While the functional boundaries are guaranteed to be closed by the LSFs, the proposed method does not require closed anatomical boundaries and can utilize incomplete edges obtained from an automatic edge detection algorithm. We conducted computer simulations to evaluate the performance of the proposed method. Two digital phantoms were constructed based on the Digimouse data and a human CT image, respectively. Anatomical edges were extracted automatically from the CT images. Tumors were simulated in the PET phantoms with different mismatched anatomical boundaries. Compared with existing methods, the new method achieved better bias-variance performance. The proposed method was also applied to real mouse data and achieved higher contrast than other methods.

Image reconstruction for PET/CT scanners: past achievements and future challenges

PET is a medical imaging modality with proven clinical value for disease diagnosis and treatment monitoring. The integration of PET and CT on modern scanners provides a synergy of the two imaging modalities. Through different mathematical algorithms, PET data can be reconstructed into the spatial distribution of the injected radiotracer. With dynamic imaging, kinetic parameters of specific biological processes can also be determined. Numerous efforts have been devoted to the development of PET image reconstruction methods over the last four decades, encompassing analytic and iterative reconstruction methods. This article provides an overview of the commonly used methods. Current challenges in PET image reconstruction include more accurate quantitation, TOF imaging, system modeling, motion correction and dynamic reconstruction. Advances in these aspects could enhance the use of PET/CT imaging in patient care and in clinical research studies of pathophysiology and therapeutic interventions.

Improved quantification in multiple-pinhole SPECT by anatomy-based reconstruction using microCT information

PURPOSE

The aim of this study was to evaluate the potential of anatomy-based reconstruction, using microCT information, to improve quantitative accuracy in multiple-pinhole SPECT.

METHODS

Multiple-pinhole SPECT and microCT images were acquired with the Micro Deluxe Phantom using both hot and cold rod inserts. The phantoms were filled with 3.7 MBq/ml of (99m)Tc. To improve microCT contrast, the phantoms were also filled with contrast agent. Emission images were reconstructed using a one-step-late (OSL) modification of the ordered subsets expectation maximization (OSEM) algorithm for incorporation of microCT information, to encourage smoothing within but not across boundaries. To allow quantification, the OSL OSEM algorithm takes into account imperfect camera motion, collimator response, angular variation of the sensitivity, intrinsic camera resolution, attenuation and scatter. For comparison, the emission images were also reconstructed by OSEM using post-reconstruction filtering and by OSL OSEM using a quadratic prior and an edge-preserving prior. In each rod of the phantoms the recovery coefficient (RC), defined as measured divided by the true activity concentration, was expressed as a function of the noise. Different noise levels were obtained by varying the amount of spatial filtering during or after reconstruction and by the use of binominal deviates.

RESULTS

Compared to conventional OSEM using post-reconstruction filtering and compared to OSL OSEM using a quadratic prior, our study demonstrated that the use of anatomical information during reconstruction significantly improved the quantitative accuracy in both cold and hot rods with a diameter larger than or equal to 2.4 mm. When compared to the edge-preserving prior, the anatomical prior performs significantly better for hot rods with a diameter ≥ 2.4 mm. For the 4.0-mm hot rods for example, the RC averaged over the different noise levels was 0.67 ± 0.02 when multiple-pinhole SPECT images were reconstructed using anatomical information, compared to 0.54 ± 0.08, 0.60 ± 0.04 and 0.64 ± 0.02 when OSEM in combination with a post-reconstruction filter, OSL OSEM using a quadratic prior and OSL OSEM using a median root prior was used, respectively. For the 4.0-mm cold rods, the RC averaged over the different noise levels was 0.61 ± 0.03 when the multiple-pinhole SPECT images were reconstructed using anatomical information, compared to 0.54 ± 0.07, 0.53 ± 0.08 and 0.60 ± 0.03 when OSEM in combination with a post-reconstruction filter, OSL OSEM using a quadratic prior and OSL OSEM using a median root prior was used, respectively.

CONCLUSION

Anatomy-based reconstruction using microCT information has the potential to improve quantitative accuracy in multiple-pinhole SPECT.

Attenuation-emission alignment in cardiac PET/CT based on consistency conditions

PURPOSE

In cardiac PET and PET/CT imaging, misaligned transmission and emission images are a common problem due to respiratory and cardiac motion. This misalignment leads to erroneous attenuation correction and can cause errors in perfusion mapping and quantification. This study develops and tests a method for automated alignment of attenuation and emission data.

METHODS

The CT-based attenuation map is iteratively transformed until the attenuation corrected emission data minimize an objective function based on the Radon consistency conditions. The alignment process is derived from previous work by Welch et al. ["Attenuation correction in PET using consistency information," IEEE Trans. Nucl. Sci. 45, 3134-3141 (1998)] for stand-alone PET imaging. The process was evaluated with the simulated data and measured patient data from multiple cardiac ammonia PET/CT exams. The alignment procedure was applied to simulations of five different noise levels with three different initial attenuation maps. For the measured patient data, the alignment procedure was applied to eight attenuation-emission combinations with initially acceptable alignment and eight combinations with unacceptable alignment. The initially acceptable alignment studies were forced out of alignment a known amount and quantitatively evaluated for alignment and perfusion accuracy. The initially unacceptable studies were compared to the proposed aligned images in a blinded side-by-side review.

RESULTS

The proposed automatic alignment procedure reduced errors in the simulated data and iteratively approaches global minimum solutions with the patient data. In simulations, the alignment procedure reduced the root mean square error to less than 5 mm and reduces the axial translation error to less than 1 mm. In patient studies, the procedure reduced the translation error by > 50% and resolved perfusion artifacts after a known misalignment for the eight initially acceptable patient combinations. The side-by-side review of the proposed aligned attenuation-emission maps and initially misaligned attenuation-emission maps revealed that reviewers preferred the proposed aligned maps in all cases, except one inconclusive case.

CONCLUSIONS

The proposed alignment procedure offers an automatic method to reduce attenuation correction artifacts in cardiac PET/CT and provides a viable supplement to subjective manual realignment tools.

A method for partial volume correction of PET-imaged tumor heterogeneity using expectation maximization with a spatially varying point spread function

Tumor heterogeneities observed in positron emission tomography (PET) imaging are frequently compromised by partial volume effects which may affect treatment prognosis, assessment or future implementations such as biologically optimized treatment planning (dose painting). This paper presents a method for partial volume correction of PET-imaged heterogeneous tumors. A point source was scanned on a GE Discovery LS at positions of increasing radii from the scanner's center to obtain the spatially varying point spread function (PSF). PSF images were fit in three dimensions to Gaussian distributions using least squares optimization. Continuous expressions were devised for each Gaussian width as a function of radial distance, allowing for generation of the system PSF at any position in space. A spatially varying partial volume correction (SV-PVC) technique was developed using expectation maximization (EM) and a stopping criterion based on the method's correction matrix generated for each iteration. The SV-PVC was validated using a standard tumor phantom and a tumor heterogeneity phantom and was applied to a heterogeneous patient tumor. SV-PVC results were compared to results obtained from spatially invariant partial volume correction (SINV-PVC), which used directionally uniform three-dimensional kernels. SV-PVC of the standard tumor phantom increased the maximum observed sphere activity by 55 and 40% for 10 and 13 mm diameter spheres, respectively. Tumor heterogeneity phantom results demonstrated that as net changes in the EM correction matrix decreased below 35%, further iterations improved overall quantitative accuracy by less than 1%. SV-PVC of clinically observed tumors frequently exhibited changes of +/-30% in regions of heterogeneity. The SV-PVC method implemented spatially varying kernel widths and automatically determined the number of iterations for optimal restoration, parameters which are arbitrarily chosen in SINV-PVC. Comparing SV-PVC to SINV-PVC demonstrated that similar results could be reached using both methods, but large differences result for the arbitrary selection of SINV-PVC parameters. The presented SV-PVC method was performed without user intervention, requiring only a tumor mask as input. Research involving PET-imaged tumor heterogeneity should include correcting for partial volume effects to improve the quantitative accuracy of results.

Impact on reader performance for lesion-detection/ localization tasks of anatomical priors in SPECT reconstruction

With increasing availability of multimodality imaging systems, high-resolution anatomical images can be used to guide the reconstruction of emission tomography studies. By measuring reader performance on a lesion detection task, this study investigates the improvement in image-quality due to use of prior anatomical knowledge, for example organ or lesion boundaries, during SPECT reconstruction. Simulated (67)Ga -citrate source and attenuation distributions were created from the mathematical cardiac-torso (MCAT) anthropomorphic digital phantom. The SIMIND Monte Carlo software was then used to generate SPECT projection data. The data were reconstructed using the De Pierro maximum a posteriori (MAP) algorithm and the rescaled-block-iterative (RBI) algorithm for comparison. We compared several degrees of prior knowledge about the anatomy: no knowledge about the anatomy; knowledge of organ boundaries; knowledge of organ and lesion boundaries; and knowledge of organ, lesion, and pseudo-lesion (non-emission uptake altering) boundaries. The MAP reconstructions used quadratic smoothing within anatomical regions, but not across any provided region boundaries. The reconstructed images were read by human observers searching for lesions in a localization receiver operating characteristic (LROC) study of the relative detection/localization accuracies of the reconstruction algorithms. Area under the LROC curve was computed for each algorithm as the comparison metric. We also had humans read images reconstructed using different prior strengths to determine the optimal trade-off between data consistency and the anatomical prior. Finally by mixing together images reconstructed with and without the prior, we tested to see if having an anatomical prior only some of the time changes the observer's detection/localization accuracy on lesions where no boundary prior is available. We found that anatomical priors including organ and lesion boundaries improve observer performance on the lesion detection/localization task. Use of just organ boundaries did not provide a statistically significant improvement in performance however. We also found that optimal prior strength depends on the level of anatomical knowledge, with a broad plateau in which observer performance is near optimal. We found no evidence that having anatomical priors use lesion boundaries only when available changes the observer's performance when they are not available. We conclude that use of anatomical priors with organ and lesion boundaries improves reader performance on a lesion-detection/localization task, and that pseudo-lesion boundaries do not hurt reader performance. However, we did not find evidence that a prior using only organ boundaries helps observer performance. Therefore we suggest prior strength should be tuned to the organ-only case, since a prior will likely not be available for all lesions.

Fast, accurate and shift-varying line projections for iterative reconstruction using the GPU

List-mode processing provides an efficient way to deal with sparse projections in iterative image reconstruction for emission tomography. An issue often reported is the tremendous amount of computation required by such algorithm. Each recorded event requires several back- and forward line projections. We investigated the use of the programmable graphics processing unit (GPU) to accelerate the line-projection operations and implement fully-3D list-mode ordered-subsets expectation-maximization for positron emission tomography (PET). We designed a reconstruction approach that incorporates resolution kernels, which model the spatially-varying physical processes associated with photon emission, transport and detection. Our development is particularly suitable for applications where the projection data is sparse, such as high-resolution, dynamic, and time-of-flight PET reconstruction. The GPU approach runs more than 50 times faster than an equivalent CPU implementation while image quality and accuracy are virtually identical. This paper describes in details how the GPU can be used to accelerate the line projection operations, even when the lines-of-response have arbitrary endpoint locations and shift-varying resolution kernels are used. A quantitative evaluation is included to validate the correctness of this new approach.

Improved attenuation correction via appropriate selection of respiratory-correlated PET data

PURPOSE

We propose a respiratory-correlated PET data processing method (called "BH-CT-based") based on breath-hold CT acquisition to reduce the smearing effect and improve the attenuation correction. The resulting images are compared with the ungated PET images acquired using a standard, free-breathing clinical protocol.

METHODS

The BH-CT-based method consisted of a list-mode acquisition with simultaneous respiratory signal recording. An additional breath-hold CT acquisition was also performed in order to define a tissue position from which PET events can be selected. A phantom study featured a 0.5-ml sphere (filled with 18F-fluorodeoxyglucose ((18)F-FDG) solution) pushed onto a rubber balloon (filled with (18)F-FDG solution and iodinated contrast agent). The feasibility of the BH-CT-based method was also assessed in two patients.

RESULTS

In the phantom study, the contrast-to-noise ratios (CNRs) were -1.6 for the Ungated volume and 5.1 for the BH-CT-based volume. For patients, CNRs were higher for BH-CT-based volumes than those for Ungated volumes (17.3 vs. 6.3 and 7.3 vs. 3.8, for patients 1 and 2, respectively). Bias-variance measurements were performed and yielded bias reduction of 40% with BH-CT-based.

CONCLUSION

The application of a BH-CT-based method decreases motion bias in PET images. This method resolves issues related to both PET-to-CT misregistration and erroneous attenuation correction and increases lesion detectability.

Joint model of motion and anatomy for PET image reconstruction

Anatomy-based positron emission tomography (PET) image enhancement techniques have been shown to have the potential for improving PET image quality. However, these techniques assume an accurate alignment between the anatomical and the functional images, which is not always valid when imaging the chest due to respiratory motion. In this article, we present a joint model of both motion and anatomical information by integrating a motion-incorporated PET imaging system model with an anatomy-based maximum a posteriori image reconstruction algorithm. The mismatched anatomical information due to motion can thus be effectively utilized through this joint model. A computer simulation and a phantom study were conducted to assess the efficacy of the joint model, whereby motion and anatomical information were either modeled separately or combined. The reconstructed images in each case were compared to corresponding reference images obtained using a quadratic image prior based maximum a posteriori reconstruction algorithm for quantitative accuracy. Results of these studies indicated that while modeling anatomical information or motion alone improved the PET image quantitation accuracy, a larger improvement in accuracy was achieved when using the joint model. In the computer simulation study and using similar image noise levels, the improvement in quantitation accuracy compared to the reference images was 5.3% and 19.8% when using anatomical or motion information alone, respectively, and 35.5% when using the joint model. In the phantom study, these results were 5.6%, 5.8%, and 19.8%, respectively. These results suggest that motion compensation is important in order to effectively utilize anatomical information in chest imaging using PET. The joint motion-anatomy model presented in this paper provides a promising solution to this problem.

Partial volume effect correction in PET using regularized iterative deconvolution with variance control based on local topology

Correcting positron emission tomography (PET) images for the partial volume effect (PVE) due to the limited resolution of PET has been a long-standing challenge. Various approaches including incorporation of the system response function in the reconstruction have been previously tested. We present a post-reconstruction PVE correction based on iterative deconvolution using a 3D maximum likelihood expectation-maximization (MLEM) algorithm. To achieve convergence we used a one step late (OSL) regularization procedure based on the assumption of local monotonic behavior of the PET signal following Alenius et al. This technique was further modified to selectively control variance depending on the local topology of the PET image. No prior 'anatomic' information is needed in this approach. An estimate of the noise properties of the image is used instead. The procedure was tested for symmetric and isotropic deconvolution functions with Gaussian shape and full width at half-maximum (FWHM) ranging from 6.31 mm to infinity. The method was applied to simulated and experimental scans of the NEMA NU 2 image quality phantom with the GE Discovery LS PET/CT scanner. The phantom contained uniform activity spheres with diameters ranging from 1 cm to 3.7 cm within uniform background. The optimal sphere activity to variance ratio was obtained when the deconvolution function was replaced by a step function few voxels wide. In this case, the deconvolution method converged in approximately 3-5 iterations for most points on both the simulated and experimental images. For the 1 cm diameter sphere, the contrast recovery improved from 12% to 36% in the simulated and from 21% to 55% in the experimental data. Recovery coefficients between 80% and 120% were obtained for all larger spheres, except for the 13 mm diameter sphere in the simulated scan (68%). No increase in variance was observed except for a few voxels neighboring strong activity gradients and inside the largest spheres. Testing the method for patient images increased the visibility of small lesions in non-uniform background and preserved the overall image quality. Regularized iterative deconvolution with variance control based on the local properties of the PET image and on estimated image noise is a promising approach for partial volume effect corrections in PET.