clear?all
close?all
clc
%?-----------------------?使用第一幀獲取圖片大小?-----------------------

addpath（‘./data‘）;
img_dir?=?dir（‘./data/*.bmp‘）;
fr?=?imread（img_dir（220）.name）;???????????%?讀第一幅圖片
fr_bw?=?rgb2gray（fr）;????????????%?灰度變換
[height?width]?=?size（fr_bw）;?
fg?=?zeros（height?width）;
bg_bw?=?zeros（height?width）;

%?-----------------------?變量初始化?-----------------------------------

C?=?4;???????????????????????%?組成混合高斯模型的單高斯數(shù)目
M?=?4;???????????????????????%?組成背景的數(shù)目
D?=?2.5;?????????????????????%?閥值（一般為2.5個(gè)標(biāo)準(zhǔn)差）
alpha?=?0.005;???????????????%?學(xué)習(xí)速率（一般為0到1）
thresh?=?0.75;???????????????%?前景閥值（0.25到0.47）
sd_init?=?6;?????????????????%?初始化標(biāo)準(zhǔn)差
update_rate?=?alpha*C;???????%?參數(shù)學(xué)習(xí)速率
rank?=?zeros（1C）;???????????%?優(yōu)先權(quán)
w?=?zeros（heightwidthC）;???%?初始化權(quán)值數(shù)組，預(yù)先分配內(nèi)存
mean?=?zeros（heightwidthC）;?????%?像素均值，預(yù)先分配內(nèi)存
pixel_sd?=?zeros（heightwidthC）;?%?像素標(biāo)準(zhǔn)差，先分配內(nèi)存
u_diff?=?zeros（heightwidthC）;????%?當(dāng)前像素與均值的差
B?=?0;??????????????????????????????????%?權(quán)值累加值

%?-------------------------初始化均值和權(quán)值----------------------------

pixel_depth?=?8;????????????????????????%??像素深度為8位
pixel_range?=?2^pixel_depth?-?1;?????????%??灰度范圍
for?i?=?1:height
????for?j?=?1:width
????????for?k?=?1:C
????????????mean（ijk）?=?pixel_range?*?rand（）;?????%均值???
????????????w（ijk）?=?1/C;???????????????????????
????????????pixel_sd（i?j?k）?=?sd_init;??%標(biāo)準(zhǔn)差????????????????
????????????
????????end
????end
end

%---------------------?處理每一幀?-----------------------------------
??figure（）
%for?n?=?1:numel（img_dir）
for?n?=?220:265
????
????fr?=?imread（img_dir（n）.name）;???????
????fr_bw?=?rgb2gray（fr）;????????
????
????%?計(jì)算每個(gè)像素與均值差
????for?m?=?1:C
????????u_diff（::m）?=?abs（double（fr_bw）?-?double（mean（::m）））;
????end
?????
???%?對(duì)于相應(yīng)條件，更新參數(shù)
????for?i=1:height
????????for?j=1:width?
????????????
????????????weight?=?w（ij:）;
????????????match?=?0;
????????????match_index?=?0;???????%設(shè)置匹配模型的位置
????????????
????????????for?k?=?1:C???????????????????????
????????????????if?（abs（u_diff（ijk））?<=?D*pixel_sd（ijk）?）???????%滿足條件?
????????????????????match?=?1;??????????????????????????
????????????????????
????????????????????%?更新均值、標(biāo)準(zhǔn)差等
????????????????????w（ijk）?=?（1-alpha）*w（ijk）?+?alpha;
????????????????????p?=?alpha/w（ijk）;??????????????????
????????????????????mean（ijk）?=?（1-p）*mean（ijk）?+?p*double（fr_bw（ij））;
????????????????????pixel_sd（ijk）?=sqrt（（1-p）*（pixel_sd（ijk）^2）?+?p*（（double（fr_bw（ij））?-?mean（ijk）））^2）;
????????????????else?????????????????????????????????????
????????????????????w（ijk）?=?（1-alpha）*w（ijk）;??????
????????????????end
????????????end??
????????????
????????????bg_bw（ij）=0;???%背景看做各個(gè)分量的加權(quán)平均
????????????for?k=1:C
????????????????bg_bw（ij）?=?bg_bw（ij）+?mean（ijk）*w（ijk）;
????????????end
????????????
????????????%?沒有一個(gè)分量滿足當(dāng)前像素，創(chuàng)建一個(gè)新的分量
????????????if?（match?==?0）
????????????????[min_w?min_w_index]?=?min（w（ij:））;??
????????????????mean（ijmin_w_index）?

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2020-02-05?18:34??實(shí)現(xiàn)了基于混合高斯模型的背景減除，在matlab下運(yùn)行\(zhòng)
?????文件????????5766??2015-11-09?18:23??實(shí)現(xiàn)了基于混合高斯模型的背景減除，在matlab下運(yùn)行\(zhòng)run.m