A general approach to regularizing inverse problems with regional data using Slepian wavelets

DOI10.1088/1361-6420/AA9909MaRDI QIDQ4601437zbMATHOpenAlexWikidataFDO

Authors V. Michel, Frederik J. Simons

Publication date 16 January 2018

Published in Inverse Problems (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1708.04462

zbMATH Keywords

waveletsregularizationHilbert spaceinverse problemsnumerical exampleill-posed problemsSlepian functionsingular-value decomposition

Mathematics Subject Classification ID

Numerical methods for wavelets (65T60) Nonlinear ill-posed problems (47J06) Numerical solutions of ill-posed problems in abstract spaces; regularization (65J20) Numerical solution to inverse problems in abstract spaces (65J22)

Abstract: Slepian functions are orthogonal function systems that live on subdomains (for example, geographical regions on the Earth's surface, or bandlimited portions of the entire spectrum). They have been firmly established as a useful tool for the synthesis and analysis of localized (concentrated or confined) signals, and for the modeling and inversion of noise-contaminated data that are only regionally available or only of regional interest. In this paper, we consider a general abstract setup for inverse problems represented by a linear and compact operator between Hilbert spaces with a known singular-value decomposition (svd). In practice, such an svd is often only given for the case of a global expansion of the data (e.g. on the whole sphere) but not for regional data distributions. We show that, in either case, Slepian functions (associated to an arbitrarily prescribed region and the given compact operator) can be determined and applied to construct a regularization for the ill-posed regional inverse problem. Moreover, we describe an algorithm for constructing the Slepian basis via an algebraic eigenvalue problem. The obtained Slepian functions can be used to derive an svd for the combination of the regionalizing projection and the compact operator. As a result, standard regularization techniques relying on a known svd become applicable also to those inverse problems where the data are regionally given only. In particular, wavelet-based multiscale techniques can be used. An example for the latter case is elaborated theoretically and tested on two synthetic numerical examples.

Recommendations

• scientific article; zbMATH DE number 484354
• The regularized orthogonal functional matching pursuit for ill-posed inverse problems
• On the regularization of projection methods for solving ill-posed problems
• Operator-theoretic and regularization approaches to ill-posed problems
• Regularization without preliminary knowledge of smoothness and error behaviour

Cited in

(4)

• A unified approach to scalar, vector, and tensor Slepian functions on the sphere and their construction by a commuting operator
• Concentration estimates for band-limited spherical harmonics expansions via the large sieve principle
• The regularized orthogonal functional matching pursuit for ill-posed inverse problems
• Multi-physics adjoint modeling of Earth structure: combining gravimetric, seismic, and geodynamic inversions

This page was built for publication: A general approach to regularizing inverse problems with regional data using Slepian wavelets

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