A method for motion estimation of sequences of images is proposed, wherein for consecutive frames (f1, f2) respective corresponding consecutive Fourier transformed frames (F1, F2) are determined, and wherein motion parameters for translation, rotation and/or for scaling are derived based on a phase relationship between said respective corresponding consecutive Fourier transformed frames (F1, F2), and in particular based on translational, rotational and/or scaling properties of the Fourier transforming process.