Elastic wave imaging has been a significant challenge in the exploration industry due to the complexities in wave physics and numerical implementation. We have separated the governing equations for P- and S-wave propagation without the assumptions of homogeneous Lamé parameters to capture the mode conversion between the two body waves in an isotropic, constant-density medium. The resulting set of two coupled second-order equations for P- and S-potentials clearly demonstrates that mode conversion only occurs at the discontinuities of the shear modulus. Applying the Born approximation to the new equations, we derive the PP, PS, SP, and SS imaging conditions from the first gradients of waveform matching objective functions. The resulting images are consistent with the physical perturbations of the elastic parameters, and, hence, they are automatically free of the polarity reversal artifacts in the converted images. When implementing elastic reverse time migration (RTM), we find that scalar wave equations can be used to back propagate the recorded P-potential, as well as individual components in the vector field of the S-potential. Compared with conventional elastic RTM, the proposed elastic RTM implementation using acoustic propagators not only simplifies the imaging condition, it but also reduces the computational cost and the artifacts in the images. We have determined the accuracy of our method using 2D and 3D numerical examples.