Public reporting burden for the collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, gathering and maintaining the data needed, and completing and reviewing the collection of information.
AbstractIn many applications of tomography, the fundamental quantities of interest in an image are geometric ones. In these instances, pixel based signal processing and reconstruction is at best inefficient, and at worst, nonrobust in its use of the available tomographic data. Classical reconstruction techniques such as Filtered Back-Projection tend to produce spurious features when data is sparse and noisy; and these "ghosts" further complicate the process of extracting what is often a limited number of rather simple geometric features. In this paper we present a framework that, in its most general form, is a statistically optimal technique for the extraction of specific geometric features of objects directly from the noisy projection data. We focus on the tomographic reconstruction of binary polygonal objects from sparse and noisy data. In our setting, the tomographic reconstruction problem is essentially formulated as a (finite dimensional) parameter estimation problem. In particular, the vertices of binary polygons are used as their defining parameters. Under the assumption that the projection data are corrupted by Gaussian white noise, we use the Maximum Likelihood (ML) criterion, when the number of parameters is assumed known, and the Minimum Description Length (MDL) criterion for reconstruction when the number of parameters is not known. The resulting optimization problems are nonlinear and thus are plagued by numerous extraneous local extrema, making their solution far from trivial. In particular, proper initialization of any iterative technique is essential for good performance.To this end, we provide a novel method to construct a reliable yet simple initial guess for the solution. This procedure is based on the estimated moments of the object, which may be conveniently obtained directly from the noisy projection data.2