Dataset Information

A scatter and randoms weighted (SRW) iterative PET reconstruction.

ABSTRACT: In positron emission tomography (PET) imaging, the main function of scatter and randoms corrections is to improve contrast and quantitative accuracy. Both corrections are essential and critically important. Several iterative reconstruction schemes incorporating scatter and randoms corrections have been developed over the years. In this work, the authors propose a new method to incorporate the scatter and randoms corrections into the iterative image reconstruction, which has shown promising results in regards to improving reconstruction performance and image quality as compared to the standard methods.The authors describe a scatter and randoms weighted (SRW) iterative PET reconstruction algorithm. The SRW method is based on the estimation of the trues fraction (TF) within the prompts. Once the TF is estimated, it is then incorporated into the weighting component of the system matrix, and the net result is a scatter and randoms weighting in the sensitivity image similar to the attenuation correction weighting. Although using the measured prompts in the TF estimation was demonstrated to achieve the fastest convergence at high statistics, it is not reliable in low counts situations due to the sparse and noisy nature of the measured prompts. Therefore, a mean estimation of the prompts derived from the forward projection of the reconstructed prompts image was introduced into the TF estimation. A contrast phantom was scanned, and the data were reconstructed using the standard and the SRW methods.The contrast vs noise, precision vs accuracy in contrast, absolute error vs number of iterations comparisons, and standard deviation image over different realizations of the same object were evaluated at low counts situations, and it was observed that the SRW method outperforms the standard approaches such as the scatter and randoms data precorrection and the ordinary Poisson methods. The image intensity (activity) outside the object can also be minimized using the SRW method. In addition, further improvement in accuracy, precision, convergence, and noise properties can be achieved by further improving the TF and the prompts estimate.The authors have developed a practical scatter and randoms weighting scheme in the sensitivity image for iterative PET reconstructions. Our proposed SRW method has a number of advantages over the conventional methods, and it has shown promising results with additional optimization for various applications to be further investigated.

SUBMITTER: Cheng JC

PROVIDER: S-EPMC3125084 | biostudies-other | 2011 Jun

REPOSITORIES: biostudies-other

ACCESS DATA

Similar Datasets

Project description:PurposeTo assess the applicability of the Fluorine-18 performance specifications defined by EANM Research Ltd (EARL), in Gallium-68 multi-centre PET-CT trials using conventional (ordered subset expectation maximisation, OSEM) and advanced iterative reconstructions which include the systems' point spread function (PSF) and a Bayesian penalised likelihood algorithm (BPL) commercially known as Q.CLEAR. The possibility of standardising the two advanced reconstruction methods was examined.MethodsThe NEMA image quality phantom was filled with Gallium-68 and scanned on a GE PET-CT system. PSF and BPL with varying post-reconstruction Gaussian filter width (2-6.4 mm) and penalisation factor (200-1200), respectively, were applied. The average peak-to-valley ratio from six profiles across each sphere was estimated to inspect any edge artefacts. Image noise was assessed using background variability and image roughness. Six GE and Siemens PET-CT scanners provided Gallium-68 images of the NEMA phantom using both conventional and advanced reconstructions from which the maximum, mean and peak recoveries were drawn. Fourteen patients underwent 68Ga-PSMA PET-CT imaging. BPL (200-1200) reconstructions of the data were compared against PSF smoothed with a 6.4-mm Gaussian filter.ResultsA Gaussian filter width of approximately 6 mm for PSF and a penalisation factor of 800 for BPL were needed to suppress the edge artefacts. In addition, those reconstructions provided the closest agreement between the two advanced iterative reconstructions and low noise levels with the background variability and the image roughness being lower than 7.5% and 11.5%, respectively. The recoveries for all methods generally performed at the lower limits of the EARL specifications, especially for the 13- and 10-mm spheres for which up to 27% (conventional) and 41% (advanced reconstructions) lower limits are suggested. The lesion standardised uptake values from the clinical data were significantly different between BPL and PSF smoothed with a Gaussian filter of 6.4 mm wide for all penalisation factors except for 800 and 1000.ConclusionIt is possible to standardise the advanced reconstruction methods with the reconstruction parameters being also sufficient for minimising the edge artefacts and noise in the images. For both conventional and advanced reconstructions, Gallium-68 specific recovery coefficient limits were required, especially for the smallest phantom spheres.

Project description:BackgroundPositron-emission tomography (PET) simulators are frequently used for development and performance evaluation of segmentation methods or quantitative uptake metrics. To date, most PET simulation tools are based on Monte Carlo simulations, which are computationally demanding. Other analytical simulation tools lack the implementation of time of flight (TOF) or resolution modelling (RM). In this study, a fast and easy-to-use PET simulation-reconstruction package, SiMulAtion and ReconsTruction (SMART)-PET, is developed and validated, which includes both TOF and RM. SMART-PET, its documentation and instructions to calibrate the tool to a specific PET/CT system are available on Zenodo. SMART-PET allows the fast generation of 3D PET images. As input, it requires one image representing the activity distribution and one representing the corresponding CT image/attenuation map. It allows the user to adjust different parameters, such as reconstruction settings (TOF/RM), noise level or scan duration. Furthermore, a random spatial shift can be included, representing patient repositioning. To evaluate the tool, simulated images were compared with real scan data of the NEMA NU 2 image quality phantom. The scan was acquired as a 60-min list-mode scan and reconstructed with and without TOF and/or RM. For every reconstruction setting, ten statistically equivalent images, representing 30, 60, 120 and 300 s scan duration, were generated. Simulated and real-scan data were compared regarding coefficient of variation in the phantom background and activity recovery coefficients (RCs) of the spheres. Furthermore, standard deviation images of each of the ten statistically equivalent images were compared.ResultsSMART-PET produces images comparable to actual phantom data. The image characteristics of simulated and real PET images varied in similar ways as function of reconstruction protocols and noise levels. The change in image noise with variation of simulated TOF settings followed the theoretically expected behaviour. RC as function of sphere size agreed within 0.3-11% between simulated and actual phantom data.ConclusionsSMART-PET allows for rapid and easy simulation of PET data. The user can change various acquisition and reconstruction settings (including RM and TOF) and noise levels. The images obtained show similar image characteristics as those seen in actual phantom data.

Project description:PURPOSE:In 90 Y microsphere radioembolization (RE), accurate post-therapy imaging-based dosimetry is important for establishing absorbed dose versus outcome relationships for developing future treatment planning strategies. Additionally, accurately assessing microsphere distributions is important because of concerns for unexpected activity deposition outside the liver. Quantitative 90 Y imaging by either SPECT or PET is challenging. In 90 Y SPECT model based methods are necessary for scatter correction because energy window-based methods are not feasible with the continuous bremsstrahlung energy spectrum. The objective of this work was to implement and evaluate a scatter estimation method for accurate 90 Y bremsstrahlung SPECT/CT imaging. METHODS:Since a fully Monte Carlo (MC) approach to 90 Y SPECT reconstruction is computationally very demanding, in the present study the scatter estimate generated by a MC simulator was combined with an analytical projector in the 3D OS-EM reconstruction model. A single window (105 to 195-keV) was used for both the acquisition and the projector modeling. A liver/lung torso phantom with intrahepatic lesions and low-uptake extrahepatic objects was imaged to evaluate SPECT/CT reconstruction without and with scatter correction. Clinical application was demonstrated by applying the reconstruction approach to five patients treated with RE to determine lesion and normal liver activity concentrations using a (liver) relative calibration. RESULTS:There was convergence of the scatter estimate after just two updates, greatly reducing computational requirements. In the phantom study, compared with reconstruction without scatter correction, with MC scatter modeling there was substantial improvement in activity recovery in intrahepatic lesions (from > 55% to > 86%), normal liver (from 113% to 104%), and lungs (from 227% to 104%) with only a small degradation in noise (13% vs. 17%). Similarly, with scatter modeling contrast improved substantially both visually and in terms of a detectability index, which was especially relevant for the low uptake extrahepatic objects. The trends observed for the phantom were also seen in the patient studies where lesion activity concentrations and lesion-to-liver concentration ratios were lower for SPECT without scatter correction compared with reconstruction with just two MC scatter updates: in eleven lesions the mean uptake was 4.9 vs. 7.1 MBq/mL (P = 0.0547), the mean normal liver uptake was 1.6 vs. 1.5 MBq/mL (P = 0.056) and the mean lesion-to-liver uptake ratio was 2.7 vs. 4.3 (P = 0.0402) for reconstruction without and with scatter correction respectively. CONCLUSIONS:Quantitative accuracy of 90 Y bremsstrahlung imaging can be substantially improved with MC scatter modeling without significant degradation in image noise or intensive computational requirements.

Dataset Information

A scatter and randoms weighted (SRW) iterative PET reconstruction.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets