Horn-Schunck 法实现的光流法

Horn-Schunck 法实现的光流法

function [fx fy ft] = computeDerivatives（im1 im2）

if size（im21）==0
    im2=zeros（size（im1））;
end

% Horn-Schunck original method
fx = conv2（im10.25* [-1 1; -1 1]‘same‘） + conv2（im2 0.25*[-1 1; -1 1]‘same‘）;
fy = conv2（im1 0.25*[-1 -1; 1 1] ‘same‘） + conv2（im2 0.25*[-1 -1; 1 1] ‘same‘）;
ft = conv2（im1 0.25*ones（2）‘same‘） + conv2（im2 -0.25*ones（2）‘same‘）;

% derivatives as in Barron
% fx= conv2（im1（1/12）*[-1 8 0 -8 1]‘same‘）;
% fy= conv2（im1（1/12）*[-1 8 0 -8 1]‘‘same‘）;
% ft = conv2（im1 0.25*ones（2）‘same‘） + conv2（im2 -0.25*ones（2）‘same‘）;
% fx=-fx;fy=-fy;

% An alternative way to compute the spatiotemporal derivatives is to use simple finite difference masks.
% fx = conv2（im1[1 -1]）;
% fy = conv2（im1[1; -1]）;
% ft= im2-im1;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        759  2009-01-05 18:05  hs\computeDerivatives.m

     文件       1081  2009-01-22 13:04  hs\gaussFilter.m

     文件       2756  2010-09-23 15:50  hs\HS.asv

     文件       2756  2010-09-23 15:50  hs\HS.m

     文件     236079  2009-01-22 13:14  hs\HSsnapshot.JPG

     文件       1347  2010-09-23 16:10  hs\plotFlow.asv

     文件       1347  2010-09-23 16:20  hs\plotFlow.m

     文件        944  2009-01-22 13:04  hs\smoothImg.m

     文件       6252  2010-09-24 15:35  hs\Untitled.asv

     文件       6252  2010-09-24 15:36  hs\Untitled.m

     文件      78886  2008-03-03 11:56  hs\yos10.tif

     文件      78908  2008-03-03 11:56  hs\yos9.tif

     目录          0  2010-09-24 15:35  hs

----------- ---------  ---------- -----  ----

               417367                    13