1993/01/01 by K. Sauer, Charles A. Bouman, C. Bouman · 6 citations
Computer Science · Engineering · Medicine · #Image and Signal Denoising Methods #Medical Imaging Techniques and Applications #Sparse and Compressive Sensing Techniques
paper · doi:10.1109/78.193196
openalex publication_date 1993/01/01 · openalex created_date 2025/10/10 · openalex updated_date 2026/07/29
A method for Bayesian reconstruction which relies on updates of single pixel values, rather than the entire image, at each iteration is presented. The technique is similar to Gauss-Seidel (GS) iteration for the solution of differential equations on finite grids. The computational cost per iteration of the GS approach is found to be approximately equal to that of gradient methods. For continuously valued images, GS is found to have significantly better convergence at modes representing high spatial frequencies. In addition, GS is well suited to segmentation when the image is constrained to be discretely valued. It is shown that Bayesian segmentation using GS iteration produces useful estimates at much lower signal-to-noise ratios than required for continuously valued reconstruction. The convergence properties of gradient ascent and GS for reconstruction from integral projections are analyzed, and simulations of both maximum-likelihood and maximum a posteriori cases are included.>