%f1?and?f2?are??two?consecutive?frames?of?a?video?sequence
clc;
clear;
close?all



f1=imread（‘CT0010.jpg‘）;
f40=imread（‘CT0011.jpg‘）;

figure;
imshow（f1）
figure
imshow（f40）


fr_f1=rgb2gray（f1）;
fr_f40=rgb2gray（f40）;?????
Xn=double（fr_f1）;
Xnp1=double（fr_f40）;
?
%get?image?size?and?adjust?for?border??獲取圖像對邊界進行調整
[hw]=size（fr_f1）;??
hm5=h-5;?wm5=w-5;
z=zeros（hw）;?v1=z;?v2=z;

%initialize????????初始化
dsx2=v1;?dsx1=v1;?dst=v1;
alpha2=625;
imax=20;

%Calculate?gradients??計算梯度
dst（5:hm55:wm5）?=?（?Xnp1（6:hm5+16:wm5+1）-Xn（6:hm5+16:wm5+1）?+?Xnp1（6:hm5+15:wm5）-Xn（6:hm5+15:wm5）+?Xnp1（5:hm56:wm5+1）-Xn（5:hm56:wm5+1）?+Xnp1（5:hm55:wm5）-Xn（5:hm55:wm5））/4;
dsx2（5:hm55:wm5）?=?（?Xnp1（6:hm5+16:wm5+1）-Xnp1（5:hm56:wm5+1）?+?Xnp1（6:hm5+15:wm5）-Xnp1（5:hm55:wm5）+?Xn（6:hm5+16:wm5+1）-Xn（5:hm56:wm5+1）?+Xn（6:hm5+15:wm5）-Xn（5:hm55:wm5））/4;
dsx1（5:hm55:wm5）?=?（?Xnp1（6:hm5+16:wm5+1）-Xnp1（6:hm5+15:wm5）?+?Xnp1（5:hm56:wm5+1）-Xnp1（5:hm55:wm5）+?Xn（6:hm5+16:wm5+1）-Xn（6:hm5+15:wm5）?+Xn（5:hm56:wm5+1）-Xn（5:hm55:wm5））/4;


for?i=1:imax
???delta=（dsx1.*v1+dsx2.*v2+dst）./（alpha2+dsx1.^2+dsx2.^2）;
???v1=v1-dsx1.*delta;
???v2=v2-dsx2.*delta;
end;
u=z;?u（5:hm55:wm5）=v1（5:hm55:wm5）;
v=z;?v（5:hm55:wm5）=v2（5:hm55:wm5）;

xskip=round（h/32）;
[hsws]=size（u（1:xskip:h1:xskip:w））
us=zeros（hsws）;?vs=us;

N=xskip^2;
for?i=1:hs-1
??for?j=1:ws-1
?????hk=i*xskip-xskip+1;
?????hl=i*xskip;
?????wk=j*xskip-xskip+1;
?????wl=j*xskip;
?????us（ij）=sum（sum（u（hk:hlwk:wl）））/N;
?????vs（ij）=sum（sum（v（hk:hlwk:wl）））/N;
???end;
end;

figure（3）;
quiver（usvs）;

axis?ij;
axis?tight;
axis?equal;

%　光流法檢測運動物體的基本原理是：給圖像中的每一個像素點賦予一個速度矢量，這就形成了一個圖像運動場，
%??在運動的一個特定時刻，圖像上的點與三維物體上的點一一對應，這種對應關系可由投影關系得到，根據各個像素點的速度矢量特征，
%??可以對圖像進行動態分析。如果圖像中沒有運動物體，則光流矢量在整個圖像區域是連續變化的。
%??當圖像中有運動物體時，目標和圖像背景存在相對運動，運動物體所形成的速度矢量必然和鄰域背景速度矢量不同，從而檢測出運動物體及位置。
%??采用光流法進行運動物體檢測的問題主要在于大多數光流法計算耗時，實時性和實用性都較差。
%??但是光流法的優點在于光流不僅攜帶了運動物體的運動信息，而且還攜帶了有關景物三維結構的豐富信息，它能夠在不知道場景的任何信息的情況下，檢測出運動對象。?

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2015-04-22?09:15??光流法檢測運動目標\
?????文件???????15776??2012-05-17?21:25??光流法檢測運動目標\CT0010.jpg
?????文件???????15843??2012-05-17?21:25??光流法檢測運動目標\CT0011.jpg
?????文件???????15913??2012-05-17?21:25??光流法檢測運動目標\CT0012.jpg
?????文件???????15919??2012-05-17?21:25??光流法檢測運動目標\CT0013.jpg
?????文件???????15584??2012-05-17?21:25??光流法檢測運動目標\CT0014.jpg
?????文件??????190100??2007-11-06?22:57??光流法檢測運動目標\car0.png
?????文件??????191767??2007-11-06?22:57??光流法檢測運動目標\car1.png
?????文件?????????138??2013-03-09?16:46??光流法檢測運動目標\www.hslogic.com.txt
?????文件????????2363??2012-05-24?19:34??光流法檢測運動目標\光流法.m