Seismic data interpolation with f-p domain spectra weighting function

Seismic trace interpolation can be formulated as an underdetermined least squares inverse problem. In order to force interpolation to follow the directions of seismic events, the weighting function is designed to build the interpolated energy along the desired direction. Usually, such a weighting function is built in the Fourier domain. However, when Fourier spectra are aliased, especially in the case of up-sampling seismic traces (creating more output than input), the design of the weighting function faces problems. In this paper, we discuss these problems and a proposed solution using the f-p domain.