matlab光流lk,Matlab数字视频处理 光流LK算法

Matlab数字视频处理 光流LK算法

Matlab数字视频处理 光流LK算法

正在学习数字视频处理，老师给了一些资料 拿出来给大家分享！ function [u, v] = LucasKanade(im1, im2, windowSize);

%LucasKanade lucas kanade algorithm, without pyramids (only 1 level);

%REVISION: NaN vals are replaced by zeros

[fx, fy, ft] = ComputeDerivatives(im1, im2);

u = zeros(size(im1));

v = zeros(size(im2));

halfWindow = floor(windowSize/2);

for i = halfWindow+1:size(fx,1)-halfWindow

for j = halfWindow+1:size(fx,2)-halfWindow

curFx = fx(i-halfWindow:i+halfWindow, j-halfWindow:j+halfWindow); curFy = fy(i-halfWindow:i+halfWindow, j-halfWindow:j+halfWindow); curFt = ft(i-halfWindow:i+halfWindow, j-halfWindow:j+halfWindow);

curFx = curFx';

curFy = curFy';

curFt = curFt';

curFx = curFx(:);

curFy = curFy(:);

curFt = -curFt(:);

A = [curFx curFy];

U = pinv(A'*A)*A'*curFt;

u(i,j)=U(1);

v(i,j)=U(2);

end;

end;

u(isnan(u))=0;

v(isnan(v))=0;

%u=u(2:size(u,1), 2:size(u,2));

%v=v(2:size(v,1), 2:size(v,2));