1. A method for estimating a velocity field using an image sequence having more than two frames, the method comprising: receiving an input image sequence comprising at least an integer number M image frames individually including multidimensional image data corresponding to a plurality of pixel locations at M different times, where M is greater than 2 ; and providing a set of equations including: M−1 displaced frame difference equations, where one of the displaced frame difference equations includes a motion compensated prediction (MCP) equation, and where M−2 of the displaced frame difference equations include motion compensated interpolation (MCI) equations, and 2x(M−2) conservative velocity constraint (CVC) equations for solving a velocity field at times t 2 , t M−1 ; using at least one processor, solving the equation set using a set of iteration equations and the image data of the M image frames to determine a velocity field describing velocity vectors at pixel locations at a first time t 1 .