A method for imaging a subsurface formation using at least two seismic sources and receivers includes: a) creating a wave field in the formation, b) obtaining a set of responses for each receiver and performing wave field separation on them, c) obtaining a virtual source signal for at least a first receiver k by cross-correlating at least part of the response of a second receiver m with at least part of the response of the first receiver k and summing the cross-correlated responses from the sources, making receiver k a virtual source, d) for each virtual source, estimating an amplitude radiation pattern by performing a transform of the response data, e) deconvolving at least one virtual source signal with at least part of the estimated radiation pattern for that receiver to produce an amplitude-corrected virtual source signal, and f) using that signal to create an image of the subsurface formation.