Share Email Print

Proceedings Paper

Discrete tomography by Bayesian labeling with its efficient algorithm
Author(s): Yuan Cheng; C. Forbes Dewey
Format Member Price Non-Member Price
PDF $14.40 $18.00
cover GOOD NEWS! Your organization subscribes to the SPIE Digital Library. You may be able to download this paper for free. Check Access

Paper Abstract

In inverse reconstruction, there are often cases where the volume or image to be reconstructed takes only a finite number of possible values. By explicitly modeling such information, discrete tomography aims to achieve better reconstruction quality and accuracy for these cases. The paper attempts to develop a framework for a general discrete tomography problem. The approach starts with an explicit model of the discreteness using a Bayesian formula. Class label variables are defined to denote the probabilities of each object point belonging to one particular class. The reconstruction then becomes the problem of assigning labels to each object point in the volume (3D) or image (2D) to be reconstructed. Unsurprisingly, this Bayesian labeling process resembles a segmentation process whose goal is also to estimate a discrete-valued field from continuous-valued observations. An expectation-maximization (EM) algorithm is developed to estimate the class label variables. By introducing another set of variables, the EM algorithm iteratively alternates the estimations of these two sets of variables. A linear equation is finally derived, composed of two terms. One accounts for the effect of the discreteness, and the other represents the integral property of the projection in tomography. This linear equation reveals a very interesting relationship between the discrete tomography and ordinary tomography, suggesting that the ordinary tomography may be treated as a special case of discrete tomography where the discreteness term is neglected. Solving the linear equation is usually very computational. This paper, however, derives an efficient algorithm by using concepts developed previously in the rho-filtered layergram method for conventional tomography. With the proposed high-pass filter, the solution for the linear equation can be computed very efficiently in the Fourier domain. In case the class values are unknown in advance, another level of EM algorithm is invoked to estimate these class values. The paper also discusses a Markov random field (MRF) model for non-stationary a priori probability, encouraging the local regularity and smoothness in the reconstruction. The experimental results demonstrate that discrete tomography using the proposed method improves the reconstruction quality greatly, especially when fewer projections are given.

Paper Details

Date Published: 4 October 2000
PDF: 9 pages
Proc. SPIE 4121, Mathematical Modeling, Estimation, and Imaging, (4 October 2000); doi: 10.1117/12.402448
Show Author Affiliations
Yuan Cheng, Massachusetts Institute of Technology (United States)
C. Forbes Dewey, Massachusetts Institute of Technology (United States)

Published in SPIE Proceedings Vol. 4121:
Mathematical Modeling, Estimation, and Imaging
David C. Wilson; Hemant D. Tagare; Fred L. Bookstein; Francoise J. Preteux; Edward R. Dougherty, Editor(s)

© SPIE. Terms of Use
Back to Top