Relaxation algorithms for matrix completion, with applications to seismic travel-time data interpolation
Autor: | C. W. Ulberg, Robert Baraldi, Aleksandr Y. Aravkin, Kenneth C. Creager, Rajiv Kumar |
---|---|
Rok vydání: | 2018 |
Předmět: |
Rank (linear algebra)
FOS: Physical sciences 010103 numerical & computational mathematics 01 natural sciences Theoretical Computer Science Physics - Geophysics Wavelet FOS: Mathematics 0101 mathematics Mathematics - Optimization and Control Mathematical Physics Mathematics Smoothness (probability theory) Matrix completion Applied Mathematics Constrained optimization Computer Science Applications Geophysics (physics.geo-ph) 010101 applied mathematics Optimization and Control (math.OC) Signal Processing Curve fitting Relaxation (approximation) Algorithm Interpolation 62F35 65K10 |
DOI: | 10.48550/arxiv.1808.04515 |
Popis: | Travel time tomography is used to infer the underlying three-dimensional wavespeed structure of the Earth by fitting seismic travel time data collected at surface stations. Data interpolation and denoising techniques are important pre-processing steps that use prior knowledge about the data, including parsimony in the frequency and wavelet domains, low-rank structure of matricizations, and local smoothness. We show how local smoothness structure can be combined with low rank constraints using level-set optimization formulations, and develop a new relaxation algorithm that can efficiently solve these joint problems. In the seismology setting, we use the approach to interpolate missing stations and de-noise observed stations. The new approach is competitive with alternative algorithms, and offers new functionality to interpolate observed data using both smoothness and low rank structure in the presence of data fitting constraints. Comment: 22 pages, 10 figures |
Databáze: | OpenAIRE |
Externí odkaz: |