We retrieve reflected plane waves by applying seismic interferometry to the recorded ground motion from a cluster of earthquakes. We employ upgoing/downgoing P/S wavefield decomposition, time windows, time reversal, and multi-dimensional deconvolution (MDD) to improve the quality of the extraction of reflected waves with seismic interferometry. Because MDD interferometry requires the separation of wavefields depending on the direction of wave propagation, almost no studies apply this technique to earthquake data observed at the surface to extract body waves. The wavefield separation and seismic interferometry based on MDD allow us to reconstruct PP, PS, SP, and SS reflected waves without unwanted crosstalk between P and S waves. From earthquake data, we obtain PP, PS, and SS reflected plane waves that reflect off the same reflector, and estimate P-and S-wave velocities.