Transmission Tomography Reconstruction Using Compound Gauss-Markov Random Fields and Ordered Subsets

Emission tomography images are degraded due to the presence of noise and several physical factors, like attenuation and scattering. To remove the attenuation effect from the emission tomography reconstruction, attenuation correction factors (ACFs) are used. These ACFs are obtained from a transmission scan and it is well known that they are homogeneous within each tissue and present abrupt variations in the transition between tissues. In this paper we propose the use of compound Gauss Markov random fields (CGMRF) as prior distributions to model homogeneity within tissues and high variations between regions. In order to find the maximum a posteriori (MAP) estimate of the reconstructed image we propose a new iterative method, which is stochastic for the line process and deterministic for the reconstruction. We apply the ordered subsets (OS) principle to accelerate the image reconstruction. The proposed method is tested and compared with other reconstruction methods.

[1]  H. Malcolm Hudson,et al.  Accelerated image reconstruction using ordered subsets of projection data , 1994, IEEE Trans. Medical Imaging.

[2]  Ken D. Sauer,et al.  A generalized Gaussian image model for edge-preserving MAP estimation , 1993, IEEE Trans. Image Process..

[3]  Anand Rangarajan,et al.  An accelerated convergent ordered subsets algorithm for emission tomography , 2004, Physics in medicine and biology.

[4]  Aggelos K. Katsaggelos,et al.  SPECT Image Reconstruction Using Compound Prior Models , 2002, Int. J. Pattern Recognit. Artif. Intell..

[5]  Paul Kinahan,et al.  Postinjection single photon transmission tomography with ordered-subset algorithms for whole-body PET imaging , 2002 .

[6]  Charles L. Byrne,et al.  Accelerating the EMML algorithm and related iterative algorithms by rescaled block-iterative methods , 1998, IEEE Trans. Image Process..

[7]  Donald Geman,et al.  Stochastic relaxation, Gibbs distributions, and the Bayesian restoration of images , 1984 .

[8]  Improvement of reconstructed images in ordered subset-Bayesian reconstruction method , 1998, 1998 IEEE Nuclear Science Symposium Conference Record. 1998 IEEE Nuclear Science Symposium and Medical Imaging Conference (Cat. No.98CH36255).

[9]  Aggelos K. Katsaggelos,et al.  Restoration of Severely Blurred High Range Images Using Stochastic and Deterministic Relaxation Algorithms in Compound Gauss Markov Random Fields , 1997, EMMCVPR.

[10]  R. Munro.  The Scientific Basis of Medicine , 1885, Glasgow medical journal.

[11]  H. Stark,et al.  Tomographic image reconstruction using the theory of convex projections. , 1988, IEEE transactions on medical imaging.

[12]  Aggelos K. Katsaggelos,et al.  Bayesian multichannel image restoration using compound Gauss-Markov random fields , 2003, IEEE Trans. Image Process..

[13]  John W. Woods,et al.  Simulated annealing in compound Gaussian random fields , 1990, IEEE Trans. Inf. Theory.

[14]  Avinash C. Kak,et al.  Principles of computerized tomographic imaging , 2001, Classics in applied mathematics.

[15]  Brian D. Ripley,et al.  Spatial Statistics: Ripley/Spatial Statistics , 2005 .

[16]  Soojin Lee.  Accelerated deterministic annealing algorithms for transmission CT reconstruction using ordered subsets , 2002 .

[17]  Hakan Erdogan,et al.  Monotonic algorithms for transmission tomography , 2002, 5th IEEE EMBS International Summer School on Biomedical Imaging, 2002..

[18]  Aggelos K. Katsaggelos,et al.  Bayesian Reconstruction for Transmission Tomography with Scale Hyperparameter Estimation , 2005, IbPRIA.

[19]  Selection of projection set and the order of calculation in ordered subsets expectation maximization method , 1997, 1997 IEEE Nuclear Science Symposium Conference Record.