A process comprising: obtaining seismic data by deploying a seismic energy source and seismic receivers, actuating the source, and detecting seismic energy reflected from a plurality of subsurface positions resulting therefrom at the receivers; forming a plurality of traces by digitally sampling the seismic energy detected at the receivers and indexing the seismic energy with respect to time; sorting the traces to form a plurality of shot gathers; applying a depth migration technique to the shot gathers to obtain the source wave field and the receiver wave field and to formulate an imaging condition; and characterized in that the imaging condition is represented by: I x → = ∫ R ω 2 c 2 x → p ̅ s x → ω p r x → ω − ∇ p ̅ s x → ω ∇ p r x → ω dω and that the method further comprises: generating two images according to the cross-correlation imaging condition (I 0) and the gradient-based imaging condition (I 1) for each shot gather as represented by: I x → = I 0 x → + I 1 x → I 0 x → = ∫ R ω 2 c 2 x → p ̅ s x → ω p r x → ω dω I 1 x → = ∫ R − ∇ p ̅ s x → ω ∇ p r x → ω) dω wherein p s ( x, ω) is the complex conjugate of the source wave field; p r ( x, ω) is the receiver wave field, ω is the circular frequency, c is the phase velocity of wave propagation, and x denotes an image point; computing a general domain parameter α at each subsurface position wherein computing the general domain parameter α is expressed as α = c 0 2 c s c r cos θ = I 1 / I 0, where θ is the local reflection angle, c 0 is the phase velocity along the axis of symmetry of seismic wave propagation, c s and c r are the phase velocities of the source and receiver wave fields; mapping the images according to the general domain parameter to form common image gathers; applying amplitude correction to common image gathers by decomposing by mapping to the reflection angle domain, or by normalizing by source-side illumination; and storing the images and common image gathers on non-transitory computer-readable media. obtaining seismic data by deploying a seismic energy source and seismic receivers, actuating the source, and detecting seismic energy reflected from a plurality of subsurface positions resulting therefrom at the receivers; forming a plurality of traces by digitally sampling the seismic energy detected at the receivers and indexing the seismic energy with respect to time; sorting the traces to form a plurality of shot gathers; applying a depth migration technique to the shot gathers to obtain the source wave field and the receiver wave field and to formulate an imaging condition; and characterized in that the imaging condition is represented by: I x → = ∫ R ω 2 c 2 x → p ̅ s x → ω p r x → ω − ∇ p ̅ s x → ω ∇ p r x → ω dω and that the method further comprises: generating two images according to the cross-correlation imaging condition (I 0) and the gradient-based imaging condition (I 1) for each shot gather as represented by: I x → = I 0 x → + I 1 x → I 0 x → = ∫ R ω 2 c 2 x → p ̅ s x → ω p r x → ω dω I 1 x → = ∫ R − ∇ p ̅ s x → ω ∇ p r x → ω) dω wherein p s ( x, ω) is the complex conjugate of the source wave field; p r ( x, ω) is the receiver wave field, ω is the circular frequency, c is the phase velocity of wave propagation, and x denotes an image point; computing a general domain parameter α at each subsurface position wherein computing the general domain parameter α is expressed as α = c 0 2 c s c r cos θ = I 1 / I 0, where θ is the local reflection angle, c 0 is the phase velocity along the axis of symmetry of seismic wave propagation, c s and c r are the phase velocities of the source and receiver wave fields; mapping the images according to the general domain parameter to form common image gathers; applying amplitude correction to common image gathers by decomposing by mapping to the reflection angle domain, or by normalizing by source-side illumination; and storing the images and common image gathers on non-transitory computer-readable media. generating two images according to the cross-correlation imaging condition (I 0) and the gradient-based imaging condition (I 1) for each shot gather as represented by: I x → = I 0 x → + I 1 x → I 0 x → = ∫ R ω 2 c 2 x → p ̅ s x → ω p r x → ω dω I 1 x → = ∫ R − ∇ p ̅ s x → ω ∇ p r x → ω) dω wherein p s ( x, ω) is the complex conjugate of the source wave field; p r ( x, ω) is the receiver wave field, ω is the circular frequency, c is the phase velocity of wave propagation, and x denotes an image point; computing a general domain parameter α at each subsurface position wherein computing the general domain parameter α is expressed as α = c 0 2 c s c r cos θ = I 1 / I 0, where θ is the local reflection angle, c 0 is the phase velocity along the axis of symmetry of seismic wave propagation, c s and c r are the phase velocities of the source and receiver wave fields; mapping the images according to the general domain parameter to form common image gathers; applying amplitude correction to common image gathers by decomposing by mapping to the reflection angle domain, or by normalizing by source-side illumination; and storing the images and common image gathers on non-transitory computer-readable media.