The determination of the P- and S-wave propagation directions is crucial for achieving wavefield decomposition, polarity reversal correction, and noise suppression in elastic reverse time migration (RTM). Compared with the conventional decoupled elastic wave equation, the first-order velocity-dilatation-rotation equations enable a more accurate computation of propagation directions for P- and S-waves. Moreover, compared with the Poynting vector, the optical flow vector signifies the wavefield propagation directions more accurately. To effectively enhance the accuracy of elastic wave imaging, we propose an elastic RTM based on first-order velocity-dilatation-rotation equations using the optical flow vector. Numerical tests illustrate that the proposed method, with or without noise, can better eliminate the migration artifacts and improve the imaging accuracy of the elastic RTM than conventional methods, achieving more accurate wavefield decomposition and superior S-wave polarity reversal correction.