25 June 1999 Fast iterative methods applied to tomography models with general Gibbs priors
Author Affiliations +
Abstract
Recently, we present a very fast method (BSREM) for solving regularized problems in emission tomography that is convergent to Maximum a Posteriori (MAO) solutions for convex priors. The method generalizes in a natural way the Expectation Maximization (EM) algorithm and some of its extensions to the MAP case. It consists of decomposing the likelihood function in blocks containing sets of projections plus a block that corresponds to the prior function, and iterating for each block using scaled gradient directions, resembling a smoothed iteration algorithm. In the general nonconvex case, it can be proven that the algorithm converges to a critical point, instead of the sought global maximum. In spite of this, BSREM is so fast and flexible, that it can be used to search for global optima when the model is not convex. In this article, we present an implementation of BSREM, that works as the local optimization method for each step of a Graduated Convexity approach as compared with BSREM itself without any convexifaction. We illustrate the behavior of the method with applications to emission tomography.
© (1999) COPYRIGHT Society of Photo-Optical Instrumentation Engineers (SPIE). Downloading of the abstract is permitted for personal use only.
Alvaro R. De Pierro, Alvaro R. De Pierro, Michel E. Beleza Yamagishi, Michel E. Beleza Yamagishi, } "Fast iterative methods applied to tomography models with general Gibbs priors", Proc. SPIE 3816, Mathematical Modeling, Bayesian Estimation, and Inverse Problems, (25 June 1999); doi: 10.1117/12.351307; https://doi.org/10.1117/12.351307
PROCEEDINGS
5 PAGES


SHARE
RELATED CONTENT

New imaging algorithm in diffusion tomography
Proceedings of SPIE (August 17 1997)
Semi-infinite positron emission tomography
Proceedings of SPIE (October 08 1995)
Fast Bayesian estimation methods in emission tomography
Proceedings of SPIE (September 21 1998)
Analytic inversion of a general model in diffuse tomography
Proceedings of SPIE (December 28 1992)

Back to Top