%%?用再次改進的Guass迭代法求背景
%?導入視頻得到用平均值求得的背景
run（‘pingjunzhi‘）

%?初始化背景
%?這里用平均值求得的背景作為初始化的背景
%?首先，將所有的點值轉換為灰度，并計算其均值、方差
for?i?=?1:100
????mov_final（i）.cdata?=?rgb2gray（mov_final（i）.cdata）;
end

back=rgb2gray（back）;

for?i?=?1:100
????mov_final（i）.cdata?=?im2double（mov_final（i）.cdata）;
end

back=im2double（back）;

for?i=1:vidHeight
????for?j=1:vidWidth
????????sum_temp=0;
????????for?k=1:100
????????????sum_temp=sum_temp+（mov_final（k）.cdata（ij）-back（ij））.^2;
????????end
????????squaresigma（ij）=sum_temp./100;
????end
end

%?參數初始化
alfa=0.05;
yuzhi1=0.5;
yuzhi2=2.5;

%?開始對前兩幅圖片做處理
mov（1:nframes）?=?...
????struct（‘cdata‘?zeros（vidHeight?vidWidth?3?‘uint8‘）...
???????????‘colormap‘?[]）;

%對第一幅圖進行處理
mov（1）.cdata=read（mov_original??1）;
mov（1）.cdata=rgb2gray（mov（1）.cdata）;
mov（1）.cdata=im2double（mov（1）.cdata）;
P=1./sqrt（2.*pi.*squaresigma）.*（-（mov（1）.cdata-back）.^2./（2.*squaresigma））;
flag=ones（vidHeightvidWidth）;
%干擾點處理
for?i=2:vidHeight-1
????for?j=2:vidWidth-1
????????P_n（1）=mov（1）.cdata（i-1j-1）;
????????P_n（2）=mov（1）.cdata（i-1j）;
????????P_n（3）=mov（1）.cdata（i-1j+1）;
????????P_n（4）=mov（1）.cdata（ij-1）;
????????P_n（5）=mov（1）.cdata（ij+1）;
????????P_n（6）=mov（1）.cdata（i+1j-1）;
????????P_n（7）=mov（1）.cdata（i+1j）;
????????P_n（8）=mov（1）.cdata（i+1j+1）;
????????if?P_n????????????flag（ij）=1;
????????else
????????????flag（ij）=0;
????????end
????end
end
%計算新的背景
for?i=?vidHeight
????for?j=?vidWidth
????????if?P（ij）<=yuzhi1&&flag（ij）==1
????????????m（ij1）=1;
????????????back（ij）=alfa*mov（1）.cdata（ij）+（1-alfa）*back（ij）;
????????elseif?P（ij）>yuzhi1
????????????m（ij1）=0;
????????end
????end
end
%計算新的方差
for?i=1:vidHeight
????for?j=1:vidWidth
????????if?P（ij）<=yuzhi1&&flag（ij）==1
????????????squaresigma（ij）=alfa*（mov（1）.cdata（ij）-back（ij））^2+（1-alfa）*squaresigma（ij）;
????????end
????end
end
????
%對第二幅圖進行處理
mov（2）.cdata=read（mov_original??2）;
mov（2）.cdata=rgb2gray（mov（2）.cdata）;
mov（2）.cdata=im2double（mov（2）.cdata）;
P=1./sqrt（2.*pi.*squaresigma）.*（-（mov（2）.cdata-back）.^2./（2.*squaresigma））;
flag=ones（vidHeightvidWidth）;
%干擾點處理
for?i=2:vidHeight-1
????for?j=2:vidWidth-1
????????P_n（1）=mov（2）.cdata（i-1j-1）;
????????P_n（2）=mov（2）.cdata（i-1j）;
????????P_n（3）=mov（2）.cdata（i-1j+1）;
????????P_n（4）=mov（2）.cdata（ij-1）;
????????P_n（5）=mov（2）.cdata（ij+1）;
????????P_n（6）=mov（2）.cdata（i+1j-1）;
????????P_n（7）=mov（2）.cdata（i+1j）;
????????P_n（8）=mov（2）.cdata（i+1j+1）;
????????if?P_n????????????flag（ij）=1;
????????else
????????????flag（ij）=0;
????????end
????end
end
%計算新的背景
for?i=?vidHeight
????for?j=?vidWidth
????????if?P（ij）<=yuzhi1&&flag（ij）==1
????????????m（ij2）=1;
????????????back（ij）=alfa*mov（2）.cdata（ij）+（1-alfa）*back（ij）;
????????elseif?P（ij）>yuzhi1
????????????m（ij2）=0;
????????end
????end
end
%計算新的方差
for?i=1:vidHeight
????for?j=1:vidWidth
????????if?P（ij）<=yuzhi1&&flag（ij）==1
????????????squaresigma（ij）=alfa*（mov（2）.cdata（ij）-back（ij））^2+（1

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????2862??2011-05-03?01:41??trafficflow\test1.m

?????文件???????4874??2011-05-01?20:46??trafficflow\adadvanced_Guass.m

?????文件???????5263??2011-05-01?13:41??trafficflow\Advanced_Guass.m

?????文件???????1225??2011-05-01?20:47??trafficflow\pingjunzhi.m

?????文件???????3132??2011-05-03?00:39??trafficflow\test1.asv

?????文件??????33935??2011-05-01?23:41??trafficflow\agback1.png

?????文件??????33935??2011-05-02?00:01??trafficflow\back_gray.png

?????文件??????33935??2011-05-02?00:11??trafficflow\aagback.png

?????文件??????33935??2011-05-02?00:11??trafficflow\agback.png

?????文件???????1578??2011-05-02?21:31??trafficflow\test2.m

?????文件????????259??2011-05-02?00:45??trafficflow\test2.asv

?????文件???22984316??2011-05-01?18:50??trafficflow\B0.avi

?????文件??????14169??2011-05-02?04:56??trafficflow\a.jpg

?????文件??????12816??2011-05-02?04:58??trafficflow\b.jpg

?????文件???????1683??2011-05-02?09:10??trafficflow\tracking.asv

?????文件???????2576??2011-05-03?02:46??trafficflow\test5.m

?????文件???????2846??2011-05-02?20:46??trafficflow\goodrestest1.m

?????文件????????841??2011-05-02?22:48??trafficflow\test_BG.m

?????文件????????780??2011-05-02?22:34??trafficflow\test_BG.asv

?????文件?????304182??2011-05-02?23:19??trafficflow\back_mhf.bmp

?????文件???????1359??2011-05-02?23:06??trafficflow\goodrestest1.rar

?????文件???????2933??2011-05-03?00:15??trafficflow\test3.m

?????文件??????12734??2011-05-02?23:23??trafficflow\增強對比度二值化結果.bmp

?????文件?????102454??2011-05-02?23:24??trafficflow\原始圖像.bmp

?????文件???????2917??2011-05-03?00:44??trafficflow\test3.asv

?????文件????????692??2011-05-03?00:10??trafficflow\和差.png

?????文件????????224??2011-05-03?00:10??trafficflow\前后幀差.png

?????文件????????540??2011-05-03?00:10??trafficflow\背景差.png

?????文件????????360??2011-05-03?00:16??trafficflow\公共部分.png

?????文件????????407??2011-05-03?00:35??trafficflow\檢測.png

............此處省略17個文件信息