Data-driven, energy-based method for estimation of scattered events in positron emission tomography.

Objective. Scattered events add bias in the reconstructed positron emission tomography (PET) images. Our objective is the accurate estimation of the scatter distribution, required for an effective scatter correction. Approach. In this paper, we propose a practical energy-based (EB) scatter estimation method that uses the marked difference between the energy distribution of the non-scattered and scattered events in the presence of randoms. In contrast to previous EB methods, we model the unscattered events using data obtained from measured point sources. Main results. We demonstrate feasibility using Monte Carlo simulated as well as experimental data acquired on the long axial field-of-view (FOV) PennPET EXPLORER scanner. Simulations show that the EB scatter estimated sinograms, for all phantoms, are in excellent agreement with the ground truth scatter distribution, known from the simulated data. Using the standard NEMA image quality (IQ) phantom we find that both the EB and single scatter simulation (SSS) provide good contrast recovery values. However, the EB correction gives better lung residuals. Significance. Application of the EB method on measured data showed, that the proposed method can be successfully translated to real-world PET scanners. When applied to a 20 cm diameter ×20 cm long cylindrical phantom the EB and SSS algorithms demonstrated very similar performance. However, on a larger 35 cm × 30 cm long cylinder the EB can better account for increased multiple scattering and out-of-FOV activity, providing more uniform images with 12%-36% reduced background variability. In typical PET ring sizes, the EB estimation can be performed in a matter of a few seconds compared to the several minutes needed for SSS, leading to efficiency advantages over the SSS implementation. as well.