Consistent Inversion of Noisy Non-Abelian X-Ray Transforms

DOI10.1002/CPA.21942MaRDI QIDQ6318128arXivFDO

Authors François Monard, Richard Nickl, Gabriel P. Paternain

Publication date 2 May 2019

Published in arXiv

Full work available at URL https://doi.org/10.1002/CPA.21942

Mathematics Subject Classification ID

Statistics (62-XX) Numerical analysis (65-XX)

arXiv classification

math.APmath.STstat.TH

Abstract: For $M$ a simple surface, the non-linear statistical inverse problem of recovering a matrix field $Phi:Momathfrakso(n)$ from discrete, noisy measurements of the $SO(n)$-valued scattering data $C_{P}hi$ of a solution of a matrix ODE is considered ($ngeq2$). Injectivity of the map $Phimapsto{C}_{P}hi$ was established by [Paternain, Salo, Uhlmann; Geom.Funct.Anal. 2012]. A statistical algorithm for the solution of this inverse problem based on Gaussian process priors is proposed, and it is shown how it can be implemented by infinite-dimensional MCMC methods. It is further shown that as the number $N$ of measurements of point-evaluations of $C_{P}hi$ increases, the statistical error in the recovery of $Phi$ converges to zero in $L^2(M)$-distance at a rate that is algebraic in $1/N$, and approaches $1/sqrtN$ for smooth matrix fields $Phi$. The proof relies, among other things, on a new stability estimate for the inverse map $C_{P}hioPhi$. Key applications of our results are discussed in the case $n=3$ to polarimetric neutron tomography, see [Desai et al., Nature Sc.Rep. 2018] and [Hilger et al., Nature Comm. 2018]

Cited in

(2)

• Discrete to continuum limits in Bayesian inverse problems
• Uncertainty Quantification for Landweber Iteration with Randomized Normal-Operator Approximation

This page was built for publication: Consistent Inversion of Noisy Non-Abelian X-Ray Transforms

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