clc;
clear?all;

targetavi=‘走過1.avi‘;?????????????????????????????%%原始視頻
resultavi=‘result1.avi‘;??????????????????????????????????%%檢測結果視頻

%%%%%%%%%%%%%%%%%%%%%%讀取視頻
mov=aviread（targetavi）;
%mov=mmreader（targetavi）;
fnum=size（mov2）;

%%%%%%%%%%%%%%%%%%%%%%建立結果視頻
aviobj?=?avifile（resultavi）;
aviobj.Quality?=?100;
aviobj.Fps?=?25;
aviobj.compression=‘Indeo5‘;
%aviobj.compression=‘None‘;


%%%%%%%%%%%%%%%%%%%%%%幀間差分法
for?i=4:fnum-3
????x1=mov（i-3）.cdata（:::）;
????x2=mov（i-2）.cdata（:::）;
????x3=mov（i-1）.cdata（:::）;
????y=mov（i）.cdata（:::）;
????z1=mov（i+1）.cdata（:::）;
????z2=mov（i+2）.cdata（:::）;
????z3=mov（i+3）.cdata（:::）;
????
????????m1=rgb2gray（x1）;
????????m1=medfilt2（m1）;%%
?????????m2=rgb2gray（x2）;
????????m2=medfilt2（m2）;%%
?????????m3=rgb2gray（x3）;
????????m3=medfilt2（m3）;%%

????????n=rgb2gray（y）;
????????n=medfilt2（n）;
????????
???????
?????????o1=rgb2gray（z1）;
????????o1=medfilt2（o1）;
?????????o2=rgb2gray（z2）;
????????o2=medfilt2（o2）;
???????o3=rgb2gray（z3）;
????????o3=medfilt2（o3）;

????????q=im2double（n）;
????????w1=im2double（m1）;
????????w2=im2double（m2）;
????????w3=im2double（m3）;
????????p1=im2double（o1）;
????????p2=im2double（o2）;
????????p3=im2double（o3）;
????????
????????c1=q-w1;%%前兩幀差分
????????c2=p3-q;
????????g1=q-w2;
????????g2=p2-q;
????????h1=q-w3;
????????h2=p1-q;
???????
????????t=40;?%%閾值，此值可以調節
????????t=t/256;%%轉化為double型數據
????????k1=find（abs（c1）>=t）;%%find函數作用是找到圖c中的值大于t的點坐標
????????k2=find（abs（c2）>=t）;
????????k3=find（abs（g1）>=t）;
????????k4=find（abs（g2）>=t）;
????????k5=find（abs（h1）>=t）;
????????k6=find（abs（h2）>=t）;
??????
????????c1（k1）=255;%%二值化的一
????????c2（k2）=255;
????????g1（k3）=255;
????????g2（k4）=255;
????????h1（k5）=255;
????????h2（k6）=255;
????????k1=find（abs（c1）????????k2=find（abs（c2）????????k3=find（abs（g1）????????k4=find（abs（g2）????????k5=find（abs（h1）????????k6=find（abs（h2）????????
????????c1（k1）=0;
????????c2（k2）=0;
???????g1（k3）=0;
???????g2（k4）=0;
??????h1（k5）=0;
???????h2（k6）=0;
????????imshow（c1）;
????????imshow（c2）;
????????imshow（g1）;
????????imshow（g2）;
????????imshow（h1）;
????????imshow（h2）;
???????
????????f1=c1&c2;
????????f2=g1&g2;
????????f3=h1&h2;
????????%BW1?=?im2bw（c?t）;%%二值化，t為閾值
????????%BW2?=?im2bw（b?t）;
????????%f=c&b;%此時f已變為logical變量
????????f=f1|f2|f3;
????????x1=double（f）;%將f轉換為double變量
????????
????????imshow（x1）;
????????
????????adata=cat（3x1x1x1）;%%灰度圖像生成視頻必加此句
????????%adata=cat（3bbb）;
????????%adata=cat（3ccc）;
????????aviobj?=?addframe（aviobjadata）;%%生成視頻，圖片變為視頻
????????
????????hold?on;
????drawnow;
???hold?off;
end
aviobj=close（aviobj）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????2835??2012-11-28?21:00??duozhen.m