A greed algorithm for seismic data interpolation using the approximate irregular discrete Fourier transform

In this paper, an alternative greed algorithm for interpolation/regularization of irregularly-sampled seismic data in the Fourier domain is described in connection with the approximate irregular discrete Fourier transform (AIDFT). The greed algorithm is used to fill in empty bins, generated by defective sampling, under statistical and/or physical constraints, so as to achieve an acceptable Fourier spetrum. Much like in other implementations, a least square norm Fourier spectrum is the input for the process. Here, this least square initial solution is provided by the AIDFT. The greed algorithm proposed is an iterative procedure that consists in, step by step, correcting for survey’s footprints of main Fourier components.