Premium
An algorithm for interpolation in the pyramid domain ‡
Author(s) -
Guitton Antoine,
Claerbout Jon
Publication year - 2010
Publication title -
geophysical prospecting
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.735
H-Index - 79
eISSN - 1365-2478
pISSN - 0016-8025
DOI - 10.1111/j.1365-2478.2010.00874.x
Subject(s) - frequency domain , algorithm , aliasing , interpolation (computer graphics) , domain (mathematical analysis) , pyramid (geometry) , filter (signal processing) , mathematics , spatial frequency , missing data , physics , optics , computer science , mathematical analysis , statistics , artificial intelligence , image (mathematics) , computer vision
With the pyramid transform, 2D dip spectra can be characterized by 1D prediction‐error filters (pefs) and 3D dip spectra by 2D pefs. These filters, contrary to pefs estimated in the frequency‐space domain (ω, x ) , are frequency independent. Therefore, one pef can be used to interpolate all frequencies. Similarly, one pef can be computed from all frequencies, thus yielding robust estimation of the filter in the presence of noise. This transform takes data from the frequency‐space domain (ω, x ) to data in a frequency‐velocity domain (ω, u =ω· x ) using a simple mapping procedure that leaves locations in the pyramid domain empty. Missing data in (ω, x ) ‐space create even more empty bins in (ω, u ) ‐space. We propose a multi‐stage least‐squares approach where both unknown pefs and missing data are estimated. This approach is tested on synthetic and field data examples where aliasing and irregular spacing are present.