Cauchy difference priors for edge-preserving Bayesian inversion with an application to X-ray tomography

arXiv1603.06135MaRDI QIDQ6271663FDOQ6271663

Lassi Roininen, Sari Lasanen, Markku Markkanen, Janne M. J. Huttunen

Publication date: 19 March 2016

Abstract: We study Cauchy-distributed difference priors for edge-preserving Bayesian statistical inverse problems. On the contrary to the well-known total variation priors, one-dimensional Cauchy priors are non-Gaussian priors also in the discretization limit. Cauchy priors have independent and identically distributed increments. One-dimensional Cauchy and Gaussian random walks are special cases of L'evy $alpha$-stable random walks with $alpha=1$ and $alpha=2$, respectively. Both random walks can be written in closed-form, and as priors, they provide smoothing and edge-preserving properties. We briefly discuss also continuous and discrete L'evy $alpha$-stable random walks, and generalize the methodology to two-dimensional priors. We apply the developed algorithm to one-dimensional deconvolution and two-dimensional X-ray tomography problems. We compute conditional mean estimates with single-component Metropolis-Hastings and maximum a posteriori estimates with Gauss-Newton-type optimization method. We compare the proposed tomography reconstruction method to filtered back-projection estimate and conditional mean estimates with Gaussian and total variation priors.

This page was built for publication: Cauchy difference priors for edge-preserving Bayesian inversion with an application to X-ray tomography

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6271663)