Fast estimation of scatter components using the ordered subset expectation maximization algorithm for scatter compensation

In this work we propose a method for scatter compensation in single photon emission computed tomography imaging, by which we can estimate the scatter components in projections in high speed with good accuracy. The method is that we first estimate the scatter components in projections based on scatter response kernels by one time of ordered subsets expectation maximization iterative reconstruction, and then subtract the estimated scatter components from the projections and complete reconstruction by filtered back-projection method. The principle is that the image corresponding to the scatter components in projections consists largely of low-frequency components of an activity distribution; these low-frequency components will converge faster than the high ones in iterative reconstruction. Therefore, we can estimate the low-frequency component image before the image converges with the high-frequency ones, and obtain the scatter components by re-projecting the low-frequency component image with scatter response kernels. The effects of the proposed method were compared with the dual- and triple-energy window methods using experimental measurements. The results show that good accuracy in estimated scatter components, good uniformity of scatter compensation at the center and the side of an object, and good noise property can be acquired by this method.