Convergent incremental optimization transfer algorithms: Application to tomography

Sangtae Ahn

doi:10.1109/TMI.2005.862740

Return

Convergent incremental optimization transfer algorithms: Application to tomography

2006-03-01

Sangtae Ahn

Jeffrey A. Fessler

Alfred O. Hero

DOI：10.1109/TMI.2005.862740

Original

Original request for help

Save

Abstract

En 中文

No convergent ordered subsets (OS) type image reconstruction algorithms for transmission tomography have been proposed to date. In contrast, in emission tomography, there are two known families of convergent OS algorithms: methods that use relaxation parameters [1], and methods based on the incremental expectation-maximization (EM) approach [2]. This paper generalizes the incremental EM approach [3] by introducing a general framework, incremental optimization transfer. The proposed algorithms accelerate convergence speeds and ensure global convergence without requiring relaxation parameters. The general optimization transfer framework allows the use of a very broad family of surrogate functions, enabling the development of new algorithms [4]. This paper provides the first convergent OS-type algorithm for (nonconcave) penalized-likelihood (PL) transmission image reconstruction by using separable paraboloidal surrogates (SPS) [5] which yield closed-form maximization steps. We found it is very effective to achieve fast convergence rates by starting with an OS algorithm with a large number of subsets and switching to the new transmission incremental optimization transfer (TRIOT) algorithm. Results show that TRIOT is faster in increasing the PL objective than nonincremental ordinary SPS and even OS-SPS yet is convergent.

Keywords:

incremental optimization transfer

maximum-likelihood estimation

penalized-likelihood estimation

statistical image reconstruction

transmission tomography