Scatter correction for positron emission mammography
In this paper we present a scatter correction method for a regularized list mode maximum likelihood reconstruction algorithm for the positron emission mammograph (PEM) that is being developed at our laboratory. The scatter events inside the object are modeled as additive Poisson random variables in the forward model of the reconstruction algorithm. The mean scatter sinogram is estimated using a Monte Carlo simulation program. With the assumption that the background activity is nearly uniform, the Monte Carlo scatter simulation only needs to run once for each PEM configuration. This saves computational time. The crystal scatters are modeled as a shift-invariant blurring in image domain because they are more localized. Thus, the useful information in the crystal scatters can be deconvolved in high-resolution reconstructions. The propagation of the noise from the estimated scatter sinogram into the reconstruction is analyzed theoretically. The results provide an easy way to calculate the required number of events in the Monte Carlo scatter simulation for a given noise level in the image. The analysis is also applicable to other scatter estimation methods, provided that the covariance of the estimated scatter sinogram is available.
- Research Organization:
- Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States)
- Sponsoring Organization:
- USDOE Director, Office of Science. Office of Biological and Environmental Research. Medical Sciences Division; National Institutes of Health (US)
- DOE Contract Number:
- AC03-76SF00098
- OSTI ID:
- 820250
- Report Number(s):
- LBNL-49398; PHMBA7; R&D Project: 860530; TRN: US0305727
- Journal Information:
- Physics in Medicine and Biology, Vol. 47, Issue 15; Other Information: Journal Publication Date: Aug. 7, 2002; PBD: 1 Apr 2002; ISSN 0031-9155
- Country of Publication:
- United States
- Language:
- English
Similar Records
Positron emission mammography with tomographic acquisition using dual planar detectors: initial evaluations
The influence of noise in full Monte Carlo ML-EM and dual matrix reconstructions in positron emission tomography