Bayesian Reconstruction for Emissiom Tomography via Deterministic Annealing

In emission tomography, a principled means of incorporating a piecewise smooth prior on the source f is via a mixed variable objective function E(f, l) defined on f and binary valued line processes l. MAP estimation on E(f, l) results in the difficult problem of minimizing an objective function that includes a nonsmooth Gibbs prior Φ* defined on the spatial derivatives of f. Previous approaches have used heuristic Gibbs potentials Φ that incorporate line processes, but only approximately. In this work, we present a continuation method in which the correct function Φ* is approached through a sequence of smooth Φ functions. Our continuation method is implemented using a GEM-ICM procedure. Simulation results show improvement using our continuation method relative to using Φ* alone, and to conventional EM reconstructions. Finally, we show a means of generalizing this formalism to the less restrictive case of piecewise linear instead of piecewise flat priors.