This work presents a new method for motion estimation of tagged cardiac magnetic resonance sequences based on variational techniques. The variational method has been improved by adding a new term in the optical flow equation that incorporates tracking points with high stability of phase. Results were obtained through simulated and real data, and were validated by manual tracking and with respect to a reference state-of-the-art method: harmonic phase imaging (HARP). The error, measured in pixels per frame, obtained with the proposed variational method is one order of magnitude smaller than the one achieved by the reference method, and it requires a lower computational cost.