Probabilistic approach to limited-data computed tomography reconstruction

Zenith Purisha*, Carl Jidling, Niklas Wahlstrom, Thomas B. Schön, Simo Särkkä

*Corresponding author for this work

Research output: Contribution to journalArticleScientificpeer-review

Abstract

In this work, we consider the inverse problem of reconstructing the internal structure of an object from limited x-ray projections. We use a Gaussian process (GP) prior to model the target function and estimate its (hyper)parameters from measured data. In contrast to other established methods, this comes with the advantage of not requiring any manual parameter tuning, which usually arises in classical regularization strategies. Our method uses a basis function expansion technique for the GP which significantly reduces the computational complexity and avoids the need for numerical integration. The approach also allows for reformulation of come classical regularization methods as Laplacian and Tikhonov regularization as GP regression, and hence provides an efficient algorithm and principled means for their parameter tuning. Results from simulated and real data indicate that this approach is less sensitive to streak artifacts as compared to the commonly used method of filtered backprojection.

Original languageEnglish
Article number105004
Number of pages20
JournalInverse Problems
Volume35
Issue number10
DOIs
Publication statusPublished - Oct 2019
MoE publication typeA1 Journal article-refereed

Keywords

  • X-ray tomography
  • Probabilistic model
  • Gaussian processes
  • Limited data
  • computed tomography
  • Radiographs
  • Statistical inversion
  • Image-reconstruction
  • Bayesian-estimation
  • Line integrals
  • Representation

Cite this