A Geometric Projection-Space Reconstruction Algorithm
Abstract
We present a method to reconstruct images from finite sets of noisy projections that may be available only over limited or sparse angles. The algorithm calculates the maximum a posteriori (MAP) estimate of the full sinogram (which is an image of the 2-D Radon transform of the object) from the available data. It is implemented using a primal-dual constrained optimization procedure that solves a partial differential equation in the primal phase with an efficient local relaxation algorithm and uses a simple Lagrange multiplier update in the dual phase. The sinogram prior probability is given by a Markov random field (MRF) that includes information about the mass, center of mass, and convex hull of the object, and about the smoothness, fundamental constraints, and periodicity of the 2-D Radon transform. The object is reconstructed using convolution back projection applied to the estimated sinogram. We show several reconstructed objects which are obtained from simulated limited-angle and sparse angle data using the described algorithm, and compare these results to images obtained using convolution back projection directly.
Document Details
- Document Type
- Technical Report
- Publication Date
- Dec 13, 1988
- Accession Number
- ADA458813
Entities
People
- Alan S. Willsky
- Jerry L. Prince
Organizations
- Massachusetts Institute of Technology