Uncertainty Quantification and Experimental Design for large-scale linear Inverse Problems under Gaussian Process Priors

09/08/2021
by   Cédric Travelletti, et al.
21

We consider the use of Gaussian process (GP) priors for solving inverse problems in a Bayesian framework. As is well known, the computational complexity of GPs scales cubically in the number of datapoints. We here show that in the context of inverse problems involving integral operators, one faces additional difficulties that hinder inversion on large grids. Furthermore, in that context, covariance matrices can become too large to be stored. By leveraging results about sequential disintegrations of Gaussian measures, we are able to introduce an implicit representation of posterior covariance matrices that reduces the memory footprint by only storing low rank intermediate matrices, while allowing individual elements to be accessed on-the-fly without needing to build full posterior covariance matrices. Moreover, it allows for fast sequential inclusion of new observations. These features are crucial when considering sequential experimental design tasks. We demonstrate our approach by computing sequential data collection plans for excursion set recovery for a gravimetric inverse problem, where the goal is to provide fine resolution estimates of high density regions inside the Stromboli volcano, Italy. Sequential data collection plans are computed by extending the weighted integrated variance reduction (wIVR) criterion to inverse problems. Our results show that this criterion is able to significantly reduce the uncertainty on the excursion volume, reaching close to minimal levels of residual uncertainty. Overall, our techniques allow the advantages of probabilistic models to be brought to bear on large-scale inverse problems arising in the natural sciences.

READ FULL TEXT

page 17

page 18

page 19

page 24

research
12/22/2022

A brief note on Bayesian D-optimality criterion

We consider finite-dimensional Bayesian linear inverse problems with Gau...
research
04/17/2023

Goal-oriented Uncertainty Quantification for Inverse Problems via Variational Encoder-Decoder Networks

In this work, we describe a new approach that uses variational encoder-d...
research
07/31/2020

Statistical guarantees for Bayesian uncertainty quantification in non-linear inverse problems with Gaussian process priors

Bayesian inference and uncertainty quantification in a general class of ...
research
05/26/2023

CUQIpy – Part I: computational uncertainty quantification for inverse problems in Python

This paper introduces CUQIpy, a versatile open-source Python package for...
research
07/15/2020

Faster Uncertainty Quantification for Inverse Problems with Conditional Normalizing Flows

In inverse problems, we often have access to data consisting of paired s...
research
12/13/2019

Solving Optimal Experimental Design with Sequential Quadratic Programming and Chebyshev Interpolation

We propose an optimization algorithm to compute the optimal sensor locat...
research
04/21/2018

A Cell-Division Search Technique for Inversion with Application to Picture-Discovery and Magnetotellurics

Solving inverse problems in natural sciences often requires a search pro...

Please sign up or login with your details

Forgot password? Click here to reset