An algorithm for optical flow computation based on a quasi-interpolant operator