clc;
clear?all;
I=imread（‘0000.jpg‘）;???????????????????????????????????????????????????????????????????????%?讀入第一幀作為背景幀
fr_bw?=?rgb2gray（I）;?????
[heightwidth]?=?size（fr_bw）;???????????????????????????????????????????????????????????????%求每幀圖像大小
C=3;????????????????????????????????????????????????????????????????????????????????????????%?單高斯模型的個數（通常為3-5）
u0?=?zeros（heightwidthC）;
for?i=?1:C-1
????u0（::i）?=?double（I（::i））;%255*double（I（::i））./double（sum（I3）+（1e-10））;
end
u0（::C）=double（fr_bw）;?
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%設置初始參數%%%%%%%%%%%%%%%%%%%%%%
M?=?3;??????????????????????????????????????????????????????????????????????????????????????%?代表背景的模型個數
T?=?5;
D?=?2.5;????????????????????????????????????????????????????????????????????????????????????%?偏差閾值
alpha?=?0.001;??????????????????????????????????????????????????????????????????????????????%?學習率
thresh?=?0.5;???????????????????????????????????????????????????????????????????????????????%?前景閾值
sd0?=?15;???????????????????????????????????????????????????????????????????????????????????%?初始化標準差
w?=?ones（heightwidthC）/C;?????????????????????????????????????????????????????????????????%?初始化權重矩陣
sd?=?sd0*ones（heightwidthC）;??????????????????????????????????????????????????????????????%?像素標準差
p?=?alpha./w;???????????????????????????????????????????????????????????????????????????????%?初始化p變量，用來更新均值和標準差
rank?=?zeros（1C）;??????????????????????????????????????????????????????????????????????????%各個高斯分布的優先級（w/sd）
frame_num=200;??????????????????????????????????????????????????????????????????????????????%幀數
for?n?=?1:frame_num
????if?n<10
????????a=‘000‘;
????end
????if?n<100&&n>=10
????????a=‘00‘;
????end
????if?100<=n
????????a=‘0‘;
????end

????frame=strcat（anum2str（n）‘.jpg‘）;
????I1=imread（frame）;???????????????????????????????????????????????????????????????????????%?依次讀入各幀圖像
????fr_bw0?=rgb2gray（I1）;
%?????fr_bw0=medfilt2（fr_bw0[3?3]）;
????h=fspecial（‘average‘）;
????fr_bw0=imfilter（fr_bw0h）;
????for?i=?1:C-1
????????fr_bw（::i）?=?double（I1（::i））;%255*double（）./double（sum（I13）+（1e-10））;
????end
????fr_bw（::C）=double（fr_bw0）;
????for?i=?C+1:M
????????fr_bw（::i）?=?double（fr_bw0）;
????end
????%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%?計算新像素與高斯模型均值的絕對距離%%%%%%%%%%%%%
????u_diff?=?abs（double（fr_bw）?-?double（u0））;
?????
????%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%?更新高斯模型的參數%%%%%%%%%%%%%%%%%%%%%%%%%%%%
????match0=abs（u_diff）<=max（D*sdT）;????????????????????????????????????????????????????????%像素與高斯模型匹配則標注為1，否則為0
????match=max（match0[]3）;?????????????????????????????????????????????????????????????????%判斷是否有高斯分布與圖像匹配
????p=alpha./w;?????????????????????????????????????????????????????????????????????????????%更新p
????wmatch1=（1-alpha）*w.*match0?+?alpha.*match0;????????????????????????????????????????????%若更新匹配部分的w
????wmatch0?=?（1-alpha）*w.*（1-match0）;??????????????????????????????????????????????????????%若更新不匹配部分的w
?

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????5072??2012-01-08?05:23??GMM.m