A method, article of manufacture, and data set for simulating seismic prestack depth migrated images on the basis of a model of a selected GF-node, without the use of real or synthetic recorded data, provides a very efficient and flexible way to calculate a simulated depth migrated image as a function of parameters such as survey, overburden model, pulse, elastic wavefield, and local reflectivity structure. Important information needed for the method is the scattering wavenumber, calculated for example, by ray methods and other equivalent methods. Complex model geometry can be done in 2-D and 3-D.