A hybrid Gibbs sampler for edge-preserving tomographic reconstruction with uncertain view angles

by   Felipe Uribe, et al.

In computed tomography, data consist of measurements of the attenuation of X-rays passing through an object. The goal is to reconstruct the linear attenuation coefficient of the object's interior. For each position of the X-ray source, characterized by its angle with respect to a fixed coordinate system, one measures a set of data referred to as a view. A common assumption is that these view angles are known, but in some applications they are known with imprecision. We propose a framework to solve a Bayesian inverse problem that jointly estimates the view angles and an image of the object's attenuation coefficient. We also include a few hyperparameters that characterize the likelihood and the priors. Our approach is based on a Gibbs sampler where the associated conditional densities are simulated using different sampling schemes - hence the term hybrid. In particular, the conditional distribution associated with the reconstruction is nonlinear in the image pixels, non-Gaussian and high-dimensional. We approach this distribution by constructing a Laplace approximation that represents the target conditional locally at each Gibbs iteration. This enables sampling of the attenuation coefficients in an efficient manner using iterative reconstruction algorithms. The numerical results show that our algorithm is able to jointly identify the image and the view angles, while also providing uncertainty estimates of both. We demonstrate our method with 2D X-ray computed tomography problems using fan beam configurations.



There are no comments yet.


page 13

page 15

page 16

page 18


Randomized Iterative Reconstruction for Sparse View X-ray Computed Tomography

With the availability of more powerful computers, iterative reconstructi...

Enhancing Industrial X-ray Tomography by Data-Centric Statistical Methods

X-ray tomography has applications in various industrial fields such as s...

Continuous Herded Gibbs Sampling

Herding is a technique to sequentially generate deterministic samples fr...

Advantage of Machine Learning over Maximum Likelihood in Limited-Angle Low-Photon X-Ray Tomography

Limited-angle X-ray tomography reconstruction is an ill-conditioned inve...

UVTomo-GAN: An adversarial learning based approach for unknown view X-ray tomographic reconstruction

Tomographic reconstruction recovers an unknown image given its projectio...

Data-driven Method for 3D Axis-symmetric Object Reconstruction from Single Cone-beam Projection Data

In this paper we consider 3D axis-symmetric (AS) object reconstruction f...

Perturbed Gibbs Samplers for Synthetic Data Release

We propose a categorical data synthesizer with a quantifiable disclosure...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.