A method for processing seismic data includes receiving seismic data and a velocity model (c(x)) for a plurality of locations (x), scaling a dimension of the seismic data according to the velocity model (c(x)) to provide a velocity normalized seismic data, and generating a final image (S(x)) of the subsurface using the velocity normalized seismic data. The velocity normalized seismic data may be a reverse-time migration image (I(x,ξ)) corresponding to the plurality of locations (x) and a plurality of propagation distance offsets (ξ). The method may also include transforming the reverse-time migration image (I(x,ξ)) for the plurality of selected positions (x) to a wavenumber domain to provide velocity normalized wavenumber data (I(k,ψ)) and suppressing data components corresponding to non-physical or undefined reflection angles to provide enhanced wavenumber data (I′(k,ψ)) and using the enhanced wavenumber data (I′(k,ψ)) to generate the final image (S(x)). A corresponding apparatus is also disclosed herein.