Energy dependence of nonstationary scatter subtraction-restoration in high resolution PET (original) (raw)
1999, IEEE Transactions on Medical Imaging
In previous works, the determination of object and detector scatter kernels from line source measurements was described and their application in scatter correction was investigated. It was also shown that low energy data contains a large fraction of useful events (true and detector scatter events). In the present work, data acquired in multispectral mode was summed in single broad energy windows of lower energy thresholds varying from 129 to 516 keV in steps of 43 keV and a constant upper energy threshold of 645 keV. Line-source projections were fitted by extracting the object and detector scatter kernels as a function of energy threshold. These kernels were then used to process scatter by the nonstationary convolution subtraction-restoration method in phantom images. After scatter correction, the sensitivity is found to increase by up to 64% at the lower threshold of 129 keV, relative to the conventional photopeak energy window (344-645 keV). Whereas contrast and spatial resolution are degraded as the energy discriminator is lowered, such degradation is fully recovered by the scatter correction. As a result of scatter correction, the noise increases insignificantly in hot regions but substantially in cold regions, in proportion of the amount of scatter.
Sign up for access to the world's latest research.
checkGet notified about relevant papers
checkSave papers to use in your research
checkJoin the discussion with peers
checkTrack your impact
Related papers
Nonstationary scatter subtraction-restoration in high-resolution PET
PubMed, 1996
Although removal of object scatter has been shown to improve both contrast and quantitation accuracy, subtraction of detector scatter leads to marginal contrast enhancement and negligible resolution recovery at the expense of reduced sensitivity and increased statistical noise. Since detector scatter has correct information about radioactivity but slightly erroneous information about source location, we suggest that this component should be restored to preserve sensitivity and improve resolution. Methods: A scatter correction model that consecutively removes object scatter and restores detector scatter is proposed. The scatter components are processed in the spatial domain using nonstationary scatter kernels. The detector scatter restoration kernel is obtained by piecewise inversion in the Fourier space. The model was tested using line source and hot spot phantom measurements. Results: Object scatter subtraction increased contrast substantively with no effect on resolution. Detector scatter restoration recovered resolution almost completely with modest contrast enhancement in small lesions. Spillover effects were reduced to less than 5% for hot spots > or = 3 x FWHM, at the expense of moderate noise amplification. Conclusion: While subtraction of object scatter is necessary for contrast enhancement and quantitation accuracy, restoration of detector scatter preserves sensitivity and improves quantitation accuracy by reducing spillover effects in high-resolution PET.
Object and detector scatter-function dependence on energy and position in high resolution PET
IEEE Transactions on Nuclear Science, 1995
We have shown in previous works that distinct non-stationary analytical scatter kernels can be extracted from line source measurements and used to independently subtract object scatter and subtract or restore detector scatter in high resolution PET. In this work, the dependence of the scatter components on energy threshold and source position I was investigated. Line source measurements were acquired in multispectral mode using the Sherbrooke PET simulator. Scatter parameters were extracted from data summed in energy windows with a lower threshold varying from 129 keV to 516 keV in steps of 42 keV, and a fixed upper threshold of 644 keV. Decreasing the lower threshold from 344 keV to 129 keV increases the trues by only 25%. but increases object scatter by 136% and almost triples detector scatter. A gain in efficiency by a factor of 2 or more would result from recovering the latter by restoration in the broad window. The intensity and shape of the scatter functions for both object and detector are shown to have a significant dependence on energy and position. This dependence needs to be taken into accwnt in the design of kernels for accurate scatter correction over a broad energy range.
Energy dependence of scatter components in multispectral PET imaging
IEEE Transactions on Medical Imaging, 1995
High resolution images in PET based on small individual detectors are obtained at the cost of low sensitivity and increased detector scatter. These limitations can be partially overcome by enlarging discrimination windows to include more low-energy events and by developing more efficient energydependent methods to correct for scatter radiation from all sources. The feasibility of multispectral scatter correction was assessed by decomposing response functions acquired in multiple energy windows into four basic components: object, collimator and detector scatter, and trues. The shape and intensity of these components are different and energy-dependent. They are shown to contribute to image formation in three ways: useful (true), potentially useful (detector scatter), and undesirable (object and collimator scatter) information to the image over the entire energy range. With the Sherbrooke animal PET system, restoration of detector scatter in every energy window would allow nearly 90% of all detected events to participate in image formation. These observations suggest that multispectral acquisition is a promising solution for increasing sensitivity in high resolution PET. This can be achieved without loss of image quality if energydependent methods are made available to preserve useful events as potentially useful events are restored and undesirable events removed.
Novel Scatter Compensation of List-Mode PET Data Using Spatial and Energy Dependent Corrections
IEEE Transactions on Medical Imaging, 2000
With the widespread use of positron emission tomography (PET) crystals with greatly improved energy resolution (e.g., 11.5% with LYSO as compared to 20% with BGO) and of list-mode acquisitions, the use of the energy of individual events in scatter correction schemes becomes feasible. We propose a novel scatter approach that incorporates the energy of individual photons in the scatter correction and reconstruction of list-mode PET data in addition to the spatial information presently used in clinical scanners. First, we rewrite the Poisson likelihood function of list-mode PET data including the energy distributions of primary and scatter coincidences and show that this expression yields an MLEM reconstruction algorithm containing both energy and spatial dependent corrections. To estimate the spatial distribution of scatter coincidences we use the single scatter simulation (SSS). Next, we derive two new formulae which allow estimation of the 2-D (coincidences) energy probability density functions (E-PDF) of primary and scatter coincidences from the 1-D (photons) E-PDFs associated with each photon. We also describe an accurate and robust object-specific method for estimating these 1-D E-PDFs based on a decomposition of the total energy spectra detected across the scanner into primary and scattered components. Finally, we show that the energy information can be used to accurately normalize the scatter sinogram to the data. We compared the performance of this novel scatter correction incorporating both the position and energy of detected coincidences to that of the traditional approach modeling only the spatial distribution of scatter coincidences in 3-D Monte Carlo simulations of a medium cylindrical phantom and a large, nonuniform NCAT phantom. Incorporating the energy information in the scatter correction decreased bias in the activity distribution estimation by 20% and 40% in the cold regions of the large NCAT phantom at energy resolutions 11.5% and 20% at 511 keV, respectively, compared to when using the spatial information alone.
Simultaneous Attenuation and Scatter Corrections in Small Animal PET Imaging
2006
The aim of this work is to simultaneously correct for attenuation and scatter in PET by analytically assessing the distribution of the scattered photons using the emission and the transmission images, the probability of scatter interactions and the detection efficiency. Above the usual lower energy threshold of 350 keV, the attenuated photons are dominantly those which have undergone a Compton scattering. By considering that each pixel in the image is the measurement of the transmitted photons through the subject, added to the contribution from the other sources by means of their scatter at this position, a simple equation is established accounting for the primaries by simultaneously correcting the data for attenuation and scatter for all emitting sources. This new method was applied in Monte Carlo simulated and measured data with the Sherbrooke small animal PET scanner in line sources, hot spot phantoms, and in rat hearts and tumors, and was compared to the approach developed by Watson et al.
Scatter degrades the contrast and quantitative accuracy of positron emission tomography (PET) images, and most methods for estimating and correcting scattered coincidences in PET subtract scattered events from the measured data. Compton scattering kinematics can be used to map out the locus of possible scattering locations. These curved lines (2D) or surfaces (3D), which connect the coincidence detectors, encompass the surface (2D) or volume (3D) where the decay occurs. In the limiting case where the scattering angle approaches zero, the scattered coincidence approaches the true coincidence. Therefore, both true and scattered coincidences can be considered similarly in a generalized scatter maximum-likelihood expectation-maximization reconstruction algorithm. The proposed method was tested using list-mode data obtained from a GATE simulation of a Jaszczak-type phantom. For scatter fractions from 10% to 60%, this approach reduces noise and improves the contrast recovery coefficients by 0.5-3.0% compared with reconstructions using true coincidences and by 3.0-24.5% with conventional reconstruction methods. The results demonstrate that this algorithm is capable of producing images entirely from scattered photons, eliminates the need for scatter corrections, increases image contrast, and reduces noise. This could be used to improve diagnostic quality and/or to reduce patient dose and radiopharmaceutical cost.
Adaptive correction of scatter and random events for 3-D backprojected PET data
IEEE Transactions on Nuclear Science, 2001
Backprojected images offer compact and accurate storage of three-dimensional (3-D) positron emission tomography (PET) data and so are potentially very suited to high-resolution 3-D PET scanners. However, there are very few methods available for correcting backprojected images for the effects of scatter and random events. Scatter correction methods usually involve experimental measurement of scatter response functions or accurate simulations, and are thus limited by the accuracy of the model used. In contrast, the approach proposed here uniquely adapts a general form of a combined scatter and random response function. This method inherently does not rely upon accurate determination of the scatter and random response functions but does, however, rely upon knowledge of the attenuating object's outline. The technique is able to adapt to different imaging situations and account for detector scatter, and in particular is readily implemented in backprojection space. The method has been implemented for a small-volume HIDAC (II) PET scanner.
Model-based scatter correction for fully 3D PET
Physics in Medicine and Biology, 1996
A method is presented that directly calculates the mean number of scattered coincidences in data acquired with fully 3D positron emission tomography (PET). This method uses a transmission scan, an emission scan, the physics of Compton scatter, and a mathematical model of the scanner in a forward calculation of the number of events for which one photon has undergone a single Compton interaction. The distribution of events for which multiple Compton interactions have occurred is modelled as a linear transformation of the single-scatter distribution. Computational efficiency is achieved by sampling at rates no higher than those required by the scatter distribution and by implementing the algorithm using look-up tables. Evaluation studies in phantoms with large scatter fractions show that the method yields images with quantitative accuracy equivalent to that of slice-collimated PET in clinically useful times.
Loading Preview
Sorry, preview is currently unavailable. You can download the paper by clicking the button above.