1. A method of processing a sampled seismic trace representing information relating to a subsurface region, comprising: estimating Fourier coefficients of the seismic trace using a discrete Fourier transform; selecting the estimated Fourier coefficients by applying a coherence criterion to the estimated Fourier coefficients to reduce an amount of aliased energy represented by the Fourier coefficients; re-orthogonalizing, in k-space, the Fourier coefficients; and using the re-orthogonalized Fourier coefficients to produce a seismic image of the subsurface region.