github:https://github.com/MichaelBeechan
CSDN:https://blog.csdn.net/u011344545
===================================================================
function [fx, fy, ft] = computeDerivatives(im1, im2)
if size(im2,1)==0
im2=zeros(size(im1));
end
% Horn-Schunck original method
fx = conv2(im1,0.25* [-1 1; -1 1],'same