On some optimization techniques in image reconstruction from projections
Applied Numerical Mathematics - Applications of optimization
SNARK09 - A software package for reconstruction of 2D images from 1D projections
Computer Methods and Programs in Biomedicine
Hi-index | 0.00 |
In previous publications we proposed, in the area of Positron Emission Tomography (PET), the use of a Maximum A posteriori Probability (MAP) optimization criterion with a particular normal prior which enforces smoothness on the resulting reconstructions. Two iterative algorithms were previously proposed to optimize this functional: one was known to converge to the desired reconstruction but was slow, while for the other it was found experimentally that, although it always appeared to converge an order of magnitude faster than the first one on examples realistic for PET, it could be made to diverge on artificial examples. In this paper we present an algorithm which is as fast as the second of these previously proposed algorithms, but it shares with the first the desirable property that it is guaranteed to converge to the reconstruction which is optimal according to our MAP criterion. We demonstrate the behavior of the algorithm on a variety of examples.