Seismic data are processed to reduce or eliminate aliasing due, for example to sparse or irregular sampling. An iterative method includes an inhibiting function used in conjunction with a function evaluating a magnitude of Fourier coefficients that together act to reduce the effects of aliased energies and preferentially select true energies. Computational steps are conducted primarily in k-space, without returning to x-space, thereby reducing computational costs.