Image reconstruction for positron emission tomography using fuzzy nonlinear anisotropic diffusion penalty

Iterative algorithms such as maximum likelihood-expectation maximization (ML-EM) become the standard for the reconstruction in emission computed tomography. However, such algorithms are sensitive to noise artifacts so that the reconstruction begins to degrade when the number of iterations reaches a certain value. In this paper, we have investigated a new iterative algorithm for penalized-likelihood image reconstruction that uses the fuzzy nonlinear anisotropic diffusion (AD) as a penalty function. The proposed algorithm does not suffer from the same problem as that of ML-EM algorithm, and it converges to a low noisy solution even if the iteration number is high. The fuzzy reasoning instead of a nonnegative monotonically decreasing function was used to calculate the diffusion coefficients which control the whole diffusion. Thus, the diffusion strength is controlled by fuzzy rules expressed in a linguistic form. The proposed method makes use of the advantages of fuzzy set theory in dealing with uncertain problems and nonlinear AD techniques in removing the noise as well as preserving the edges. Quantitative analysis shows that the proposed reconstruction algorithm is suitable to produce better reconstructed images when compared with ML-EM, ordered subsets EM (OS-EM), Gaussian-MAP, MRP, TV-EM reconstructed images.

[1]  Omer Demirkaya,et al.  Post-reconstruction filtering of positron emission tomography whole-body emission images and attenuation maps using nonlinear diffusion filtering. , 2004, Academic radiology.

[2]  Alfred O. Hero,et al.  Space-alternating generalized expectation-maximization algorithm , 1994, IEEE Trans. Signal Process..

[3]  Jerry M. Mendel,et al.  Generating fuzzy rules by learning from examples , 1992, IEEE Trans. Syst. Man Cybern..

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

[5]  Juan Ruiz,et al.  Fuzzy anisotropic diffusion for speckle filtering , 2001, 2001 IEEE International Conference on Acoustics, Speech, and Signal Processing. Proceedings (Cat. No.01CH37221).

[6]  Jeffrey A. Fessler,et al.  Grouped-coordinate ascent algorithms for penalized-likelihood transmission image reconstruction , 1997, IEEE Transactions on Medical Imaging.

[7]  Jeffrey A. Fessler,et al.  Globally convergent image reconstruction for emission tomography using relaxed ordered subsets algorithms , 2003, IEEE Transactions on Medical Imaging.

[8]  J.M.M. Anderson,et al.  Accelerated penalized weighted least-squares and maximum likelihood algorithms for reconstructing transmission images from PET transmission data , 2005, IEEE Transactions on Medical Imaging.

[9]  K. Rajan,et al.  Image reconstruction by conditional entropy maximisation for PET system , 2004 .

[10]  Jeffrey A. Fessler,et al.  Grouped coordinate descent algorithms for robust edge-preserving image restoration , 1997, Optics & Photonics.

[11]  Anand Rangarajan,et al.  Provably convergent OSEM-like reconstruction algorithm for emission tomography , 2002, SPIE Medical Imaging.

[12]  Jesse S. Jin,et al.  An adaptive nonlinear diffusion algorithm for filtering medical images , 2000, IEEE Transactions on Information Technology in Biomedicine.

[13]  O. Demirkaya Anisotropic diffusion filtering of PET attenuation data to improve emission images. , 2002, Physics in medicine and biology.

[14]  Jeffrey A. Fessler Penalized weighted least-squares image reconstruction for positron emission tomography , 1994, IEEE Trans. Medical Imaging.

[15]  Dimitri Van De Ville,et al.  Noise reduction by fuzzy image filtering , 2003, IEEE Trans. Fuzzy Syst..

[16]  Pao-Ta Yu,et al.  Fuzzy selection filters for image restoration with neural learning , 1999, IEEE Trans. Signal Process..

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

[18]  F. Hermansen,et al.  Noise reduction and convergence of Bayesian algorithms with blobs based on the Huber function and median root prior , 2004, Physics in medicine and biology.

[19]  Siegfried Gottwald,et al.  Fuzzy Sets and Fuzzy Logic , 1993 .

[20]  George J. Klir,et al.  Fuzzy sets and fuzzy logic , 1995 .

[21]  Ulla Ruotsalainen,et al.  Accelerated median root prior reconstruction for pinhole single-photon emission tomography (SPET). , 2003, Physics in medicine and biology.

[22]  N.V. Denisova Bayesian reconstruction in SPECT with entropy prior and iterative statistical regularization , 2004, IEEE Transactions on Nuclear Science.

[23]  Grant T. Gullberg,et al.  Total variation regulated EM algorithm , 1998, 1998 IEEE Nuclear Science Symposium Conference Record. 1998 IEEE Nuclear Science Symposium and Medical Imaging Conference (Cat. No.98CH36255).

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

[25]  Guillermo Sapiro,et al.  Robust anisotropic diffusion , 1998, IEEE Trans. Image Process..

[26]  Donald Geman,et al.  Stochastic Relaxation, Gibbs Distributions, and the Bayesian Restoration of Images , 1984, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[27]  C. Riddell,et al.  Diffusion regularization for iterative reconstruction in emission tomography , 2003, 2003 IEEE Nuclear Science Symposium. Conference Record (IEEE Cat. No.03CH37515).

[28]  Dan J. Kadrmas Statistically regulated and adaptive EM reconstruction for emission computed tomography , 2000 .

[29]  Ulla Ruotsalainen,et al.  Generalization of median root prior reconstruction , 2002, IEEE Transactions on Medical Imaging.

[30]  Yutaka Hata,et al.  Interactive segmentation of the cerebral lobes with fuzzy inference in 3T MR images , 2006, IEEE Transactions on Systems, Man, and Cybernetics, Part B (Cybernetics).

[31]  Alvaro R. De Pierro,et al.  A row-action alternative to the EM algorithm for maximizing likelihood in emission tomography , 1996, IEEE Trans. Medical Imaging.

[32]  Fabrizio Russo An image enhancement technique combining sharpening and noise reduction , 2002, IEEE Trans. Instrum. Meas..

[33]  Jitendra Malik,et al.  Scale-Space and Edge Detection Using Anisotropic Diffusion , 1990, IEEE Trans. Pattern Anal. Mach. Intell..

[34]  A Boudraa,et al.  Fuzzy clustering-based segmented attenuation correction in whole-body PET imaging. , 2002, Physics in medicine and biology.

[35]  Partha P Mondal,et al.  Fuzzy-rule-based image reconstruction for positron emission tomography. , 2005, Journal of the Optical Society of America. A, Optics, image science, and vision.

[36]  Gengsheng L. Zeng,et al.  Total variation regulated EM algorithm [SPECT reconstruction] , 1999 .

[37]  P. Green Bayesian reconstructions from emission tomography data using a modified EM algorithm. , 1990, IEEE transactions on medical imaging.