Improved iterative reconstruction method for Compton imaging using median filter

A Compton camera is a device for imaging a radio-source distribution without using a mechanical collimator. Ordered-subset expectation-maximization (OS-EM) is widely used to reconstruct Compton images. However, the OS-EM algorithm tends to over-concentrate and amplify noise in the reconstructed image. It is, thus, necessary to optimize the number of iterations to develop high-quality images, but this has not yet been achieved. In this paper, we apply a median filter to an OS-EM algorithm and introduce a median root prior expectation-maximization (MRP-EM) algorithm to overcome this problem. In MRP-EM, the median filter is used to update the image in each iteration. We evaluated the quality of images reconstructed by our proposed method and compared them with those reconstructed by conventional algorithms using mathematical phantoms. The spatial resolution was estimated using the images of two point sources. Reproducibility was evaluated on an ellipsoidal phantom by calculating the residual sum of squares, zero-mean normalized cross-correlation, and mutual information. In addition, we evaluated the semi-quantitative performance and uniformity on the ellipsoidal phantom. MRP-EM reduces the generated noise and is robust with respect to the number of iterations. An evaluation of the reconstructed image quality using some statistical indices shows that our proposed method delivers better results than conventional techniques.