clear;
[filenamepathname]=uigetfile（‘*.avi‘‘請選擇文件‘）;
fpath=[pathname?filename];
read_obj=VideoReader（fpath）;
fr_num?=?read_obj.NumberOfframes;%?視頻幀數

num_1=0;
for?n?=50:fr_num-1?
fr?=?read（read_obj?n）;????
%----動態監測----得到currentframe和grayframe_1----------------------
currentframe?=?read（read_obj?n）;
????I=0.3*currentframe（::1）+0.59*currentframe（::2）+0.11*currentframe（::3）;?%來源于基于視頻圖像的火災檢測設計與實現_孫凱軍
????grayframe_1=[round（I）];
????%----------------
%1---fr=imread（‘dde1.jpg‘）;
%?讀取該圖像作為背景?
fr_bw1?=?rgb2gray（fr）;???
%?將背景轉換為灰度圖像?
fr_size?=?size（fr）;?????????????
%取幀大小?
width?=?fr_size（2）;?
height?=?fr_size（1）;?%獲取原圖像的尺寸
fg?=?zeros（height?width）;???%前景讀取的第二張圖片獲得
bg_bw?=?zeros（height?width）;%背景讀取的第一張圖片獲得
fr_bw1?=?double（fr_bw1）;%256?400
%?---------------------?mog?variables?-----------------------------------?
C?=?3;???????????%?組成混合高斯的單高斯數目?（一般3-5）?
D?=?2.5;?????????%?閾值（一般2.5個標準差）?
alpha?=?0.01;????%?learning?rate?學習率決定更新速度（between?0?and?1）?（from?paper?0.01）?
thresh?=?0.25;???%?foreground?threshold?前景閾值（0.25?or?0.75?in?paper）?
sd_init?=?36;????%?initial?standard?deviation?初始化標準差（for?new?components）?var?=?36?in?paper?
w?=?zeros（heightwidthC）;?%256?1200?????????%?initialize?weights?array?初始化權值數組?
w（::1）?=?1;
w（::2:C）?=?2^-10;??????????????????%?第一個高斯分布的初始權重為1，其余分布的權重為0
mean?=?zeros（heightwidthC）;%256?1200????????%?pixel?means?像素均值?
mean（::1）?=?fr_bw1;????????????????%?第一個高斯分布的初始均值為參考幀的值，其余分布的均值為0s
sd?=?sd_init*ones（heightwidthC）;???%?pixel?standard?deviations?像素標準差?
matchcnt?=?ones（height?widthC）;????%?匹配的次數，初始值都設為1
u_diff?=?zeros（heightwidthC）;??????%?difference?of?each?pixel?from?mean?圖片與高斯均值的差?
?%?從第二幀開始處理
????fr?=?read（read_obj?n+2）;?????%?讀取n1-n2幀并轉化為灰度圖、雙精度
??
fr_bw?=?rgb2gray（fr）;??	???%?convert?frame?to?grayscale?轉換為灰度圖像
fr_bw?=?double（fr_bw）;?????%?將灰度圖值設置為雙精度
?%求導入進來的圖片與各個高斯均值的差
for?m=1:C??
	u_diff（::m）?=?abs（fr_bw?-?double（mean（::m）））;?????
end???????
%?update?gaussian?components?for?each?pixel?更新每個像素的背景模型?
%rank_ind?=?zeros（C1）;
for?i=1:height
	for?j=1:width??????????????
????????match?=?0;?%像素與高斯模型匹配的標識
????????match_ind?=?0;%為該像素最匹配的高斯模型的標號
????????for?k=1:C?%與第k個高斯模型進行比對，然后更新參數?????????????????
????????????if?（abs（u_diff（ijk））?<=?D*sd（ijk）?&&?（match?==?0））???????%?pixel?matches?component像素匹配了高斯中的第k個模型?			
????????????????match?=?1;???????????????????????????????
????????????????%?variable?to?signal?component?match?設置匹配標記?				
????????????????match_ind?=?k;
????????????????%?update?weights?mean?sd?p??更新權值，均值，標準差和參數學習率?????????????????????
??				p?=?alpha/w（ijk）;???????????%理應使用p?=?alpha/gaussian才對，這里勉強??????
????????????????w（ijk）?=?（1-alpha）*w（ijk）?+?alpha;?????????????????????
				%p?=?alpha/w（ijk）;???????????%理應使用p?=?alpha/gaussian才對，這里勉強???????
????????????????mean（ijk）?=?（1-p）*mean（ijk）?+?p*double（fr_bw（ij））;?
????????????????sd（ijk）?=???sqrt（（1-p）*（sd（ijk）^2）?+?p*（（double（fr_bw（ij））?-?mean（ijk）））^2）;
??????