Interpolation of geophysical data using continuous global surfaces
- 1 November 2002
- journal article
- Published by Society of Exploration Geophysicists in Geophysics
- Vol. 67 (6), 1810-1822
- https://doi.org/10.1190/1.1527081
Abstract
A wide class of interpolation methods, including thin‐plate and tension splines, kriging, sinc functions, equivalent‐source, and radial basis functions, can be encompassed in a common mathematical framework involving continuous global surfaces (CGSs). The difficulty in applying these techniques to geophysical data sets has been the computational and memory requirements involved in solving the large, dense matrix equations that arise. We outline a three‐step process for reducing the computational requirements: (1) replace the direct inversion techniques with iterative methods such as conjugate gradients; (2) use preconditioning to cluster the eigenvalues of the interpolation matrix and hence speed convergence; and (3) compute the matrix–vector product required at each iteration with a fast multipole or fast moment method. We apply the new methodology to a regional gravity compilation with a highly heterogeneous sampling density. The industry standard minimum‐curvature algorithms and several scale‐dependent CGS methods are unable to adapt to the varying data density without introducing spurious artifacts. In contrast, the thin‐plate spline is scale independent and produces an excellent fit. When applied to an aeromagnetic data set with relatively uniform sampling, the thin‐plate spline does not significantly improve results over a standard minimum‐curvature algorithm.Keywords
This publication has 19 references indexed in Scilit:
- Smooth fitting of geophysical data using continuous global surfacesGeophysics, 2002
- Fast Solution of the Radial Basis Function Interpolation Equations: Domain Decomposition MethodsSIAM Journal on Scientific Computing, 2001
- Fast Evaluation of Radial Basis Functions: Moment-Based MethodsSIAM Journal on Scientific Computing, 1998
- Splines — more than just a smooth interpolatorGeoderma, 1994
- Computation of Thin-Plate SplinesSIAM Journal on Scientific and Statistical Computing, 1991
- A fast algorithm for particle simulationsJournal of Computational Physics, 1987
- Interpolation of scattered data: Distance matrices and conditionally positive definite functionsConstructive Approximation, 1986
- GMRES: A Generalized Minimal Residual Algorithm for Solving Nonsymmetric Linear SystemsSIAM Journal on Scientific and Statistical Computing, 1986
- Numerical Procedures for Surface Fitting of Scattered Data by Radial FunctionsSIAM Journal on Scientific and Statistical Computing, 1986
- Comparing splines and krigingComputers & Geosciences, 1984