Computer instructions, computing device and method for processing seismic data under-sampled in an angle domain, the seismic data corresponding to a reverse time migration, three-dimensional, angle domain common image gather (ADCIG). The method includes receiving the seismic data; calculating, based on the seismic data, shot and receiver wave-fields with an RTM wave propagation engine; applying a wave-fields decomposition algorithm to obtain a propagation direction for the shot and receiver wave-fields; forming the ADCIG by applying an image condition to the shot and receiver wave-fields; determining that specular energies of the ADCIG are under-sampled around a reflection angle; during the step of forming the ADCIG, extrapolating the specular energies to a neighborhood of the reflection angle; and generating an image of a subsurface that is being surveyed based on the extrapolated specular energies.