Repository logo
 

Efficient nonparametric bayesian inference for X-ray transforms

Published version
Peer-reviewed

Change log

Authors

Monard, F 
Nickl, R 
Paternain, GP 

Abstract

We consider the statistical inverse problem of recovering a function f:MR, where M is a smooth compact Riemannian manifold with boundary, from measurements of general X-ray transforms Ia(f) of f, corrupted by additive Gaussian noise. For M equal to the unit disk with flat' geometry and $a=0$ this reduces to the standard Radon transform, but our general setting allows for anisotropic media $M$ and can further model local attenuation' effects -- both highly relevant in practical imaging problems such as SPECT tomography. We propose a nonparametric Bayesian inference approach based on standard Gaussian process priors for f. The posterior reconstruction of f corresponds to a Tikhonov regulariser with a reproducing kernel Hilbert space norm penalty that does not require the calculation of the singular value decomposition of the forward operator Ia. We prove Bernstein-von Mises theorems that entail that posterior-based inferences such as credible sets are valid and optimal from a frequentist point of view for a large family of semi-parametric aspects of f. In particular we derive the asymptotic distribution of smooth linear functionals of the Tikhonov regulariser, which is shown to attain the semi-parametric Cram'er-Rao information bound. The proofs rely on an invertibility result for the `Fisher information' operator IaIa between suitable function spaces, a result of independent interest that relies on techniques from microlocal analysis. We illustrate the performance of the proposed method via simulations in various settings.

Description

Keywords

Inverse problem, Bernstein-von Mises theorem, MAP estimate, Tikhonov regulariser, Gaussian prior, Radon transform, semiparametric efficiency

Journal Title

Annals of Statistics

Conference Name

Journal ISSN

0090-5364

Volume Title

47

Publisher

Institute of Mathematical Statistics
Sponsorship
Engineering and Physical Sciences Research Council (EP/M023842/1)
European Research Council (647812)