利用matlab的光流法实现显示光流的方向

基本实现光流法，可显示出两张图片中光流的方向。

%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

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2016-03-09 14:50  matlab\
     文件        2363  2016-03-09 15:55  matlab\aa.asv
     文件        2363  2016-03-09 14:40  matlab\aa.m
     文件      190100  2007-11-06 22:57  matlab\car0.png
     文件      191767  2007-11-06 22:57  matlab\car1.png
     文件       15776  2012-05-17 21:25  matlab\CT0010.jpg
     文件       15843  2012-05-17 21:25  matlab\CT0011.jpg
     文件       15913  2012-05-17 21:25  matlab\CT0012.jpg
     文件       15919  2012-05-17 21:25  matlab\CT0013.jpg
     文件       15584  2012-05-17 21:25  matlab\CT0014.jpg
     文件       24064  2016-03-04 10:33  matlab\Thumbs.db
     文件         138  2013-03-09 16:46  matlab\www.hslogic.com.txt