clear?all
videofile=‘people.avi‘;%讀取視頻文件
%videofile=‘carreco.avi‘;
%videofile=‘viptraffic.avi‘;
info=mmfileinfo（videofile）;%獲取視頻文件信息
cols=info.Video.Width;%列數
rows=info.Video.Height;%行數
%創建視頻對象，讀取視頻文件
hreader=vision.VideoFileReader（videofile‘ImageColorSpace‘‘RGB‘‘VideoOutputDataType‘‘single‘）;
%創建光流對象用于檢測運動方向和速度
hflow=vision.OpticalFlow（‘OutputValue‘‘Horizontal?and?vertical?components?in?complex?form‘‘ReferenceframeDelay‘3‘Method‘‘Horn-Schunck‘）;
%創建兩個均值對象，用于分析光流矢量
hmean1=vision.Mean;
hmean2=vision.Mean（‘RunningMean‘true）;
%創建均值濾波對象
hfilter=vision.MedianFilter;
%創建形態學關閉對象，填充分割以后的汽車空洞
hclose=vision.MorphologicalClose（‘Neighborhood‘strel（‘line‘545））;
%創建BLOB分析對象，用于從視頻中分割汽車
hblob=vision.BlobAnalysis;
hblob.CentroidOutputPort=false;?
hblob.AreaOutputPort=true;
hblob.BoundingBoxOutputPort=true;
hblob.OutputDataType=‘double‘;
hblob.MinimumBlobArea=250;
hblob.MaximumBlobArea=3600;
hblob.MaximumCount=80;
%創建形態學腐蝕對象，移除不需要的對象
herode=vision.MorphologicalErode（‘Neighborhood‘strel（‘square‘2））;
%創建形狀嵌入對象，在視頻中添加形狀，框出汽車邊界。
hshape1=vision.ShapeInserter（‘BorderColor‘‘Custom‘‘CustomBorderColor‘[0?1?0]）;
hshape2=vision.ShapeInserter（‘Shape‘‘Lines‘‘BorderColor‘‘Custom‘‘CustomBorderColor‘[255?255?0]）;
htext=vision.TextInserter（‘Text‘‘%4d‘‘Location‘[1?1]‘Color‘[1?1?1]‘FontSize‘12）;
sz=get（0‘ScreenSize‘）;
pos=[（sz（3）-4*（cols+75））/2?（sz（4）-rows）/2?cols+60?rows+80];
hvideo1=vision.VideoPlayer（‘Name‘‘Original?Video‘‘Position‘pos）;
%psos（1）=pos（1）+cols+5;
psos（1）=pos（1）+cols+15;%行數加上15個像素，實際上那個是在進行平移
hvideo2=vision.VideoPlayer（‘Name‘‘Motion?Vector‘‘Position‘pos）;
pos（1）=pos（1）+cols+15;
hvideo3=vision.VideoPlayer（‘Name‘‘Thresholded?Video‘‘Position‘pos）;
pos（1）=pos（1）+cols+15;
hvideo4=vision.VideoPlayer（‘Name‘‘Result?Video‘‘Position‘pos）;

%從視頻中檢測追蹤汽車
%顯示光流矢量像素點
[xposypos]=meshgrid（1:5:cols1:5:rows）;%grid?網格，坐標方格
xpos=xpos（:）;
ypos=ypos（:）;
locs=sub2ind（[rowscols]yposxpos）;

%循環處理視頻每一幀，直到文件結束
while?~isDone（hreader）
????pause（0.3）;
frame=step（hreader）;
gray=rgb2gray（frame）;

%計算光流場矢量，返回一個復數矩陣，分別代表每個像素點的U，V
%光流基本方程，dit（I）+It=0;
flow=step（hflowgray）;%就算階躍響應，實際上是在計算光流失量
%繪制光流圖系統默認每隔5行5列選擇一個像素點，繪制他的光流圖，20表示將光流幅值放大20倍，lines?每行對應一條曲線，
lines=[xposyposxpos+20*real（flow（locs））ypos+20*imag（flow（locs））];
%將光流矢量添加到視頻幀上
vector=step（hshape2framelines）;%實際上就是求取一個階躍響應
magnitude=flow.*conj（flow）;%conj求取復數的共軛
threshold=0.5*step（hmean2step（hmean2magnitude））;
carobj=step（hfiltermagnitude>=threshold）;
carobj=step（hclosestep（herodecarobj））;
[areabbox]=step（hblobcarobj）;
grow=30;
idx=bbox（:1）>grow;
ratio=zeros（length（idx）1）;
ratio（idx）=single（area（idx1））./single（bbox（idx3）.*bbox（idx4））;
flag=ratio>0.45;

count=int32（sum（flag））;
%count=count+1;
%bbox（~flag:）=int32（-1）;
result=step（hshape1framebbox）;
result（grow:grow+1::）=1;
result（1:151:30:）=0;
result=step（htextresultcount）;

%顯示處理結果
step（hvideo1frame）;
step（hvideo2vector）;
step（hvideo3car