Data interpolation is an important step for seismic data analysis
because many processing tasks, such as multiple attenuation and
migration, are based on regularly sampled seismic data. Failed
interpolations may introduce artifacts and eventually lead to
inaccurate final processing results. In this paper, we generalize
seismic data interpolation as a basis pursuit problem and propose an
iteration framework for recovering missing data. The method is based
on nonlinear iteration and sparse transform. A modified Bregman
iteration is used for solving the constrained minimization problem
based on compressed sensing. The new iterative strategy guarantees
fast convergence by using a fixed threshold value. We also propose a
generalized velocity-dependent (VD) formulation of the seislet
transform as an effective sparse transform, in which the nonhyperbolic
normal moveout equation serves as a bridge between local slope
patterns and moveout parameters in the common-midpoint domain. It can
also be reduced to the traditional VD-seislet if special heterogeneity
parameter is selected. The generalized VD-seislet transform predicts
prestack reflection data in offset coordinates, which provides a high
compression of reflection events. The method was applied to synthetic
and field data examples and the results show that the generalized
VD-seislet transform can reconstruct missing data with the help of the
modified Bregman iteration even for nonhyperbolic reflections under
complex conditions, such as VTI media or aliasing.