I=imread（‘1bmpfile.bmp‘）;???????????????%?讀入第一幀作為背景幀
fr_bw?=?I;?????
[heightwidth]?=?size（fr_bw）;???????????%求每幀圖像大小
fg?=?zeros（height?width）;??????????????%定義前景和背景矩陣
bg_bw?=?zeros（height?width）;


C?=?3;??????????????????????????????????%?單高斯模型的個數（通常為3-5）
M?=?3;??????????????????????????????????%?代表背景的模型個數
D?=?2.5;????????????????????????????????%?偏差閾值
alpha?=?0.01;???????????????????????????%?學習率
thresh?=?0.25;??????????????????????????%?前景閾值
sd_init?=?15;????????????????????????????%?初始化標準差
w?=?zeros（heightwidthC）;??????????????%?初始化權重矩陣
mean?=?zeros（heightwidthC）;???????????%?像素均值
sd?=?zeros（heightwidthC）;?????????????%?像素標準差
u_diff?=?zeros（heightwidthC）;?????????%?像素與某個高斯模型均值的絕對距離
p?=?alpha/（1/C）;????????????????????????%?初始化p變量，用來更新均值和標準差
rank?=?zeros（1C）;??????????????????????%各個高斯分布的優先級（w/sd）




pixel_depth?=?8;????????????????????????%?每個像素8bit分辨率
pixel_range?=?2^pixel_depth?-1;?????????%?像素值范圍[0255]

for?i=1:height
????for?j=1:width
????????for?k=1:C
????????????
????????????mean（ijk）?=?rand*pixel_range;?????%初始化第k個高斯分布的均值
????????????w（ijk）?=?1/C;?????????????????????%?初始化第k個高斯分布的權重
????????????sd（ijk）?=?sd_init;????????????????%?初始化第k個高斯分布的標準差
????????????
????????end
????end
end


frame_num=23;%幀數
for?n?=?1:frame_num
????frame=strcat（num2str（n）‘bmpfile.bmp‘）;
????I1=imread（frame）;??%?依次讀入各幀圖像
????fr_bw?=I1;???????
????
????%?計算新像素與第m個高斯模型均值的絕對距離
????for?m=1:C
????????u_diff（::m）?=?abs（double（fr_bw）?-?double（mean（::m）））;
????end
?????
????%?更新高斯模型的參數
????for?i=1:height
????????for?j=1:width
????????????
????????????match?=?0;???????????????????????????????????????%匹配標記;
????????????for?k=1:C???????????????????????
????????????????if?（abs（u_diff（ijk））?<=?D*sd（ijk））???????%?像素與第k個高斯模型匹配
????????????????????
????????????????????match?=?1;???????????????????????????????%將匹配標記置為1
????????????????????
????????????????????%?更新權重、均值、標準差、p
????????????????????w（ijk）?=?（1-alpha）*w（ijk）?+?alpha;
????????????????????p?=?alpha/w（ijk）;??????????????????
????????????????????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）;
????????????????else?????????????????????????????????????????%?像素與第k個高斯模型不匹配
????????????????????w（ijk）?=?（1-alpha）*w（ijk）;???????????%略微減少權重
????????????????????
????????????????end
????????????end
????????????
??????????????????
????????????bg_bw（ij）=0;
????????????for?k=1:C
????????????????bg_bw（ij）?=?bg_bw（ij）+?mean（ijk）*w（ijk）;
????????????end
????????????
????????????%?像素值與任一高斯模型都不匹配，則創建新的模型
????????????if?（match?==?0）
????????????????[min_w?min_w_index]?=?min（w（ij:））;??????%尋找最小權重
????????????????mean（ijmin_w_index）?=?double（fr_bw（ij））;%初始化均值為當前觀測像素的均值
????????????????sd（ijmin_w_index）?=?sd_init;?????????????%初始化標準差為6
????????????????end

????????????rank?=?w（ij:）./sd（ij:）;???

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

?????文件??????87478??2004-06-24?11:53??混合高斯背景建模\10bmpfile.bmp

?????文件??????87478??2004-06-24?11:53??混合高斯背景建模\11bmpfile.bmp

?????文件??????87478??2004-06-24?11:53??混合高斯背景建模\12bmpfile.bmp

?????文件??????87478??2004-06-24?11:53??混合高斯背景建模\13bmpfile.bmp

?????文件??????87478??2004-06-24?11:53??混合高斯背景建模\14bmpfile.bmp

?????文件??????87478??2004-06-24?11:54??混合高斯背景建模\15bmpfile.bmp

?????文件??????87478??2004-06-24?11:54??混合高斯背景建模\16bmpfile.bmp

?????文件??????87478??2004-06-24?11:54??混合高斯背景建模\17bmpfile.bmp

?????文件??????87478??2004-06-24?11:54??混合高斯背景建模\18bmpfile.bmp

?????文件??????87478??2004-06-24?11:54??混合高斯背景建模\19bmpfile.bmp

?????文件??????87478??2004-06-24?11:53??混合高斯背景建模\1bmpfile.bmp

?????文件??????87478??2004-06-24?11:54??混合高斯背景建模\20bmpfile.bmp

?????文件??????87478??2004-06-24?11:54??混合高斯背景建模\21bmpfile.bmp

?????文件??????87478??2004-06-24?11:54??混合高斯背景建模\22bmpfile.bmp

?????文件??????87478??2004-06-24?11:54??混合高斯背景建模\23bmpfile.bmp

?????文件??????87478??2004-06-24?11:53??混合高斯背景建模\2bmpfile.bmp

?????文件??????87478??2004-06-24?11:53??混合高斯背景建模\3bmpfile.bmp

?????文件??????87478??2004-06-24?11:53??混合高斯背景建模\4bmpfile.bmp

?????文件??????87478??2004-06-24?11:53??混合高斯背景建模\5bmpfile.bmp

?????文件??????87478??2004-06-24?11:53??混合高斯背景建模\6bmpfile.bmp

?????文件??????87478??2004-06-24?11:53??混合高斯背景建模\7bmpfile.bmp

?????文件??????87478??2004-06-24?11:53??混合高斯背景建模\8bmpfile.bmp

?????文件??????87478??2004-06-24?11:53??混合高斯背景建模\9bmpfile.bmp

?????文件??????87478??2011-10-10?23:44??混合高斯背景建模\FG.bmp

?????文件???????4560??2011-10-10?23:53??混合高斯背景建模\gaussians.asv

?????文件???????4205??2011-10-12?12:33??混合高斯背景建模\gaussians.m

?????文件??????87294??2011-09-14?20:10??混合高斯背景建模\gau_pic23.jpg

?????目錄??????????0??2012-05-18?08:59??混合高斯背景建模

-----------?---------??----------?-----??----

??????????????2195531????????????????????28

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