16 December 1992 Nonlinear multigrid methods of optimization in Bayesian tomographic image reconstruction
Author Affiliations +
Abstract
Bayesian estimation of transmission tomographic images presents formidable optimization tasks. Numerical solutions of this problem are limited in speed of convergence by the number of iterations required for the propagation of information across the grid. Edge-preserving prior models for tomographic images inject a nonlinear element into the Bayesian cost function, which limits the effectiveness of algorithms such as conjugate gradient, intended for linear problems. In this paper, we apply nonlinear multigrid optimization to Bayesian reconstruction of a two-dimensional function from integral projections. At each resolution, we apply Gauss-Seidel type iterations, which optimize locally with respect to individual pixel values. If the cost function is differentiable, the algorithm speeds convergence; if it is nonconvex and/or nondifferentiable, multigrid can yield improved estimates.
© (1992) COPYRIGHT Society of Photo-Optical Instrumentation Engineers (SPIE). Downloading of the abstract is permitted for personal use only.
Charles A. Bouman, Charles A. Bouman, Ken D. Sauer, Ken D. Sauer, } "Nonlinear multigrid methods of optimization in Bayesian tomographic image reconstruction", Proc. SPIE 1766, Neural and Stochastic Methods in Image and Signal Processing, (16 December 1992); doi: 10.1117/12.130838; https://doi.org/10.1117/12.130838
PROCEEDINGS
11 PAGES


SHARE
Back to Top