运动目标检测光流法

在MATLAB2014上输入一个视频可以检测出其中的运动目标并对之进行标记

function [uv]=lk（img1img2wsize）;
[fxfyft]=deriv（img1img2）;
u=zeros（size（img1））;
v=zeros（size（img2））;
half_w=floor（wsize/2）;
for i=half_w+1:size（fx1）-half_w
    for j=half_w+1:size（fx2）-half_w
        curx=fx（i-half_w:i+half_wj-half_w:j+half_w）;
        cury=fy（i-half_w:i+half_wj-half_w:j+half_w）;
        curt=ft（i-half_w:i+half_wj-half_w:j+half_w）;
        curx=curx‘;
        cury=cury‘;
        curt=curt‘;
        curx=curx（:）;
        cury=cury（:）;
        curt=-curt（:）;
        a=[curxcury];
        UV=pinv（a‘*a）*a‘*curt;
        u（ij）=UV（1）;
        v（ij）=UV（2）;
    end;
end;
u（isnan（u））=0;
v（isnan（v））=0;

function [fxfyft]=deriv（img1img2）;
fx=conv2（img10.25*[-11;-11]）+conv2（img20.25*[-11;-11]）;
fy=conv2（img10.25*[-1-1;11]）+conv2（img20.25*[-1-1;11]）;
ft=conv2（img10.25*ones（2））+conv2（img2-0.25*ones（2））;

fx=fx（1:size（fx1）-11:size（fx2）-1）;
fy=fy（1:size（fy1）-11:size（fy2）-1）;
ft=ft（1:size（ft1）-11:size（ft2）-1）;

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

     文件       1922  2018-04-10 17:56  optial_flow.m

     文件       1516  2018-04-13 10:19  optical_flow_detect.m

     文件        995  2012-05-25 16:22  lk.m

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

                 4433                    3