Dual matrix ordered subsets reconstruction for accelerated 3D scatter compensation in single-photon emission tomography

Abstract.

Three-dimensional (3D) iterative maximum likelihood expectation maximization (ML-EM) algorithms for single-photon emission tomography (SPET) are capable of correcting image-degrading effects of non-uniform attenuation, distance-dependent camera response and patient shape-dependent scatter. However, the resulting improvements in quantitation, resolution and signal-to-noise ratio (SNR) are obtained at the cost of a huge computational burden. This paper presents a new acceleration method for ML-EM: dual matrix ordered subsets (DM-OS). DM-OS combines two acceleration methods: (a) different matrices for projection and back-projection and (b) ordered subsets of projections. DM-OS was compared with ML-EM on simulated data and on physical thorax phantom data, for both 180° and 360° orbits. Contrast, normalized standard deviation and mean squared error were calculated for the digital phantom experiment. DM-OS resulted in similar image quality to ML-EM, even for speed-up factors of 200 compared to ML-EM in the case of 120 projections. The thorax phantom data could be reconstructed 50 times faster (60 projections) using DM-OS with preservation of image quality. ML-EM and DM-OS with scatter compensation showed significant improvement of SNR compared to ML-EM without scatter compensation. Furthermore, inclusion of complex image formation models in the computer code is simplified in the case of DM-OS. It is thus shown that DM-OS is a fast and relatively simple algorithm for 3D iterative scatter compensation, with similar results to conventional ML-EM, for both 180° and 360° acquired data.