Spatiotemporal Besov Priors for Bayesian Inverse Problems

Abstract

Fast development in science and technology has driven the need for proper statistical tools to capture special data features such as abrupt changes or sharp contrast. Many inverse problems in data science require spatiotemporal solutions derived from a sequence of time-dependent objects with these spatial features, e.g., the dynamic reconstruction of computerized tomography (CT) images with edges. Conventional methods based on Gaussian processes (GP) often fall short in providing satisfactory solutions since they tend to offer oversmooth priors. Recently, the Besov process (BP), defined by wavelet expansions with random coefficients, has emerged as a more suitable prior for Bayesian inverse problems of this nature. While BP excels in handling spatial inhomogeneity, it does not automatically incorporate temporal correlation inherited in the dynamically changing objects. In this paper, we generalize BP to a novel spatiotemporal Besov process (STBP) by replacing the random coefficients in the series expansion with stochastic time functions as Q-exponential process (Q-EP) which governs the temporal correlation structure. We thoroughly investigate the mathematical and statistical properties of STBP. Simulations, two limited-angle CT reconstruction examples, a highly non-linear inverse problem involving Navier-Stokes equation, and a spatiotemporal temperature imputation problem are used to demonstrate the advantage of the proposed STBP compared with the classic STGP and a time-uncorrelated approach.

64 pages, 17 figures

References in corpus (10)

• Rates of contraction of posterior distributions based on Gaussian process priors
• Convergence rates of posterior distributions for noniid observations
• Bayesian inverse problems with Gaussian priors
• Adaptive Bayesian estimation using a Gaussian random field with inverse Gamma bandwidth
• Geometric MCMC for Infinite-Dimensional Inverse Problems
• Rates of contraction for posterior distributions in $\bolds{L^r}$-metrics, $\bolds{1\le r\le\infty}$
• Posterior Consistency for Bayesian Inverse Problems through Stability and Regression Results
• Multivariate type G Matérn stochastic partial differential equation random fields
• A TV-Gaussian prior for infinite-dimensional Bayesian inverse problems and its numerical implementations
• MAP Estimators for Piecewise Continuous Inversion