Joint-MAP Tomographic Reconstruction with Patch Similarity Based Mixture Prior Model

Tomographic reconstruction from noisy projections do not yield adequate results. Mathematically, this tomographic reconstruction represents an ill-posed problem due to information missing caused by the presence of noise. Maximum a posteriori (MAP) or Bayesian reconstruction methods offer possibilities to improve the image quality as compared with analytical methods in particular by introducing a prior to guide the reconstruction and regularize the noise. With an aim to achieve robust utilization of continuity/connectivity information and overcome the heuristic weight update for other nonlocal prior methods, this paper proposes a novel patch similarity based mixture (PSM) prior model for tomographic reconstruction. This prior is defined by a weighted Gaussian distance between neighborhood intensities. The weight quantifies the similarity between local neighborhoods and is computed using a maximization entropy constraint. This prior is then introduced within a joint image/weight MAP computed tomography reco...

[1]  Jeffrey A. Fessler,et al.  Edge-preserving tomographic reconstruction with nonlocal regularization , 1998, IEEE Transactions on Medical Imaging.

[2]  Song-Chun Zhu,et al.  FRAME: filters, random fields, and minimax entropy towards a unified theory for texture modeling , 1996, Proceedings CVPR IEEE Computer Society Conference on Computer Vision and Pattern Recognition.

[3]  Cong Nie,et al.  Bayesian statistical reconstruction for low-dose X-ray computed tomography using an adaptive-weighting nonlocal prior , 2009, Comput. Medical Imaging Graph..

[4]  Michael Knaup,et al.  GPU-based parallel-beam and cone-beam forward- and backprojection using CUDA , 2008, 2008 IEEE Nuclear Science Symposium Conference Record.

[5]  Wojciech Zbijewski,et al.  Suppression of intensity transition artifacts in statistical x-ray computer tomography reconstruction through radon inversion initialization. , 2003, Medical physics.

[6]  J. Hornegger,et al.  Fast GPU-Based CT Reconstruction using the Common Unified Device Architecture (CUDA) , 2007, 2007 IEEE Nuclear Science Symposium Conference Record.

[7]  Yann Gousseau,et al.  The TVL1 Model: A Geometric Point of View , 2009, Multiscale Model. Simul..

[8]  Guy Gilboa,et al.  Nonlocal Operators with Applications to Image Processing , 2008, Multiscale Model. Simul..

[9]  K. Lange.  Convergence of EM image reconstruction algorithms with Gibbs smoothing. , 1990, IEEE transactions on medical imaging.

[10]  Jeffrey A. Fessler,et al.  Ieee Transactions on Image Processing: to Appear Globally Convergent Algorithms for Maximum a Posteriori Transmission Tomography , 2022 .

[11]  Siwei Lyu.  An implicit Markov random field model for the multi-scale oriented representations of natural images , 2009, 2009 IEEE Conference on Computer Vision and Pattern Recognition.

[12]  Song-Chun Zhu,et al.  Filters, Random Fields and Maximum Entropy (FRAME): Towards a Unified Theory for Texture Modeling , 1998, International Journal of Computer Vision.

[13]  M. Glas,et al.  Principles of Computerized Tomographic Imaging , 2000 .

[14]  Jeffrey A. Fessler,et al.  Efficient and accurate likelihood for iterative image reconstruction in x-ray computed tomography , 2003, SPIE Medical Imaging.

[15]  Laurent D. Cohen,et al.  Non-local Regularization of Inverse Problems , 2008, ECCV.

[16]  Florence Tupin,et al.  Iterative Weighted Maximum Likelihood Denoising With Probabilistic Patch-Based Weights , 2009, IEEE Transactions on Image Processing.

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

[18]  F. Beekma,et al.  Ordered subset reconstruction for x-ray CT. , 2001, Physics in medicine and biology.

[19]  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.

[20]  Chih-Jen Lin,et al.  Iterative Scaling and Coordinate Descent Methods for Maximum Entropy , 2009, ACL.

[21]  Anand Rangarajan,et al.  Bayesian image reconstruction for transmission tomography using deterministic annealing , 2003, J. Electronic Imaging.

[22]  Stan Z. Li,et al.  Markov Random Field Modeling in Image Analysis , 2001, Computer Science Workbench.

[23]  Ken D. Sauer,et al.  A unified approach to statistical tomography using coordinate descent optimization , 1996, IEEE Trans. Image Process..

[24]  M. Bertero,et al.  Ill-posed problems in early vision , 1988, Proc. IEEE.

[25]  Hakan Erdogan,et al.  Monotonic algorithms for transmission tomography , 1999, IEEE Transactions on Medical Imaging.

[26]  Jianhua Ma,et al.  Nonlocal Prior Bayesian Tomographic Reconstruction , 2008, Journal of Mathematical Imaging and Vision.

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

[28]  D. Hunter,et al.  Optimization Transfer Using Surrogate Objective Functions , 2000 .

[29]  Chi Zhang,et al.  Fast and Robust Deconvolution-Based Image Reconstruction for Photoacoustic Tomography in Circular Geometry: Experimental Validation , 2010, IEEE Photonics Journal.

[30]  Jean-Michel Morel,et al.  A Review of Image Denoising Algorithms, with a New One , 2005, Multiscale Model. Simul..

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

[32]  Gerhard Winkler,et al.  Image Analysis, Random Fields and Markov Chain Monte Carlo Methods: A Mathematical Introduction , 2002 .