Computing device, computer instructions and method process input seismic data d recorded in a first domain by seismic receivers that travel in water, the input seismic data d including pressure and particle motion measurements, including up-going and down-going wave-fields. A model p is generated in a second domain by solving an inverse problem for the input seismic data d, wherein applying an L transform to the model p describes the input data d. An L′ transform, which is different from the L transform, is then applied to the model p to obtain an output seismic data in the first domain, the output seismic data having a characteristic imparted by the transform L′. The characteristic is related to pressure wave-fields and/or particle motion wave-fields interpolated at positions in-between the input seismic receivers. An image of the surveyed subsurface is generated based on the output seismic dataset.