Methods and systems for separating P-S wave field data are described. Slowness values for the PP mode and the PS mode are estimated and are, typically, unequal based on aliased and/or irregularly sampled data. A calculation, in the space-time domain, based on a matrix of equations, generates separated P-wave and S-wave data. The separated P-wave data and S-wave data are output for further imaging.