Iterative maximum a posteriori (MAP) restoration from partially-known blur for tomographic reconstruction

An iterative maximum a posteriori (MAP) algorithm is proposed for simultaneous signal-covariance estimation and restoration when only partial knowledge of the system response matrix (SRM) and the noisy-blurred sinogram of an image to be reconstructed are available. Convergence analysis is performed to ascertain that the proposed covariance estimator converges to the optimal one in the MAP sense. The superiority of the proposed algorithm, in comparison with the iterative linear minimum mean-squared-error (LMMSE) filter for incorrect SRM information, is experimentally verified.

[1]  Jeffrey A. Fessler.  Tomographic Reconstruction Using Information-Weighted Spline Smoothing , 1993, IPMI.

[2]  C T Chen,et al.  Superresolved tomography by convex projections and detector motion. , 1992, Journal of the Optical Society of America. A, Optics and image science.

[3]  Nikolas P. Galatsanos,et al.  MAP and regularized constrained total least-squares image restoration , 1994, Proceedings of 1st International Conference on Image Processing.

[4]  B. R. Hunt,et al.  Digital Image Restoration , 1977 .

[5]  Rabab K. Ward,et al.  Restoration of images distorted by systems of random impulse response , 1985 .

[6]  Roland T. Chin,et al.  Iterative Wiener filters for image restoration , 1991, IEEE Trans. Signal Process..

[7]  Nikolas P. Galatsanos,et al.  Regularized constrained total least-squares image restoration , 1994, Other Conferences.

[8]  T. J. Hebert,et al.  A union of deterministic and stochastic methods for image reconstruction , 1992, IEEE Conference on Nuclear Science Symposium and Medical Imaging.