Scatter Correction in PET Using the Transport Equation

In the presence of a strong scatterer, it is inappropriate to perform scatter correction by subtracting scattered from measured events and work on unscattered events only. Take, for example, a strong scatterer that completely surrounds the reconstruction volume. In this case, all events are scattered and removal simply would leave nothing. We propose an algorithm for PET reconstruction that is based on modeling the scanning process by the full transport equation, incorporating (multiple) scatter and attenuation. Then, a solution to the full forward model is sought. The concept sounds simple, but it is hard to implement since the forward model becomes arbitrarily complex. Also, care has to be taken to use expectation maximization techniques to retain image quality.

[1]  C. C. Watson,et al.  New, faster, image-based scatter correction for 3D PET , 1999, 1999 IEEE Nuclear Science Symposium. Conference Record. 1999 Nuclear Science Symposium and Medical Imaging Conference (Cat. No.99CH37019).

[2]  G. C. Pomraning,et al.  Linear Transport Theory , 1967 .

[3]  L. Shepp,et al.  Maximum Likelihood Reconstruction for Emission Tomography , 1983, IEEE Transactions on Medical Imaging.

[4]  M. Schellmann,et al.  Parallelization and Runtime Prediction of the ListMode OSEM Algorithm for 3D PET Reconstruction , 2006, 2006 IEEE Nuclear Science Symposium Conference Record.

[5]  C. Bohm,et al.  Correction for Scattered Radiation in a Ring Detector Positron Camera by Integral Transformation of the Projections , 1983, Journal of computer assisted tomography.