With the capability of handling complicated velocity models, reverse-time migration (RTM) has become a powerful imaging method. Improving imaging accuracy and computational efficiency are two significant but challenging tasks in the applications of RTM. Despite being the most popular numerical technique applied in RTM, finite-difference (FD) methods often suffer from undesirable numerical dispersion, leading to a noticeable loss of imaging resolution. A new and effective FD operator, called the high-order stereo operator, has been developed to approximate the partial differential operators in the wave equation, from which a numerical scheme called the three-step stereo method (TSM) has been developed and has shown effectiveness in suppressing numerical dispersion. Numerical results show that compared with the conventional numerical methods such as the Lax-Wendroff correction (LWC) scheme and the staggered-grid (SG) FD method, this new method significantly reduces numerical dispersion and computational cost. Tests on the impulse response and the 2D prestack Hess acoustic VTI model demonstrated that the TSM achieves higher image quality than the LWC and SG methods do, especially when coarse computation grids were used, which indicated that the new method can be a promising algorithm for large-scale anisotropic RTM.