function?x?=?Whl
%初始化背景所需要的幀數
N=15;%背景更新參數
mov?=?aviread（‘1.AVI‘）;%讀取AVI視頻文件
%獲取視頻幀數
dimension?=?size（mov）;
frameQTY?=?dimension（2）;
M=mov;
tmFr?=?rgb2gray（frame2im（mov（1）））;%變成灰度圖
s?=?size（tmFr）;
scrsz?=?get（0‘ScreenSize‘）;%獲取屏幕大小
h1?=?figure（‘position‘?[11s（2）*2s（1）*2]）;
cdataSum?=?int16（zeros（s））;
sedisk?=?strel（‘rectangle‘[22]）;
sedisk2?=?strel（‘rectangle‘[32]）;
count?=?zeros（s）;
for?c1?=?1:frameQTY;
????colorFrC?=?M（c1）;%獲取第i幀
????frC?=?rgb2gray（frame2im（colorFrC））;
????frC?=?im2frame（frC?gray（256））;
????f?=?int16（frC.cdata）;
????if（?c1<=?N）
????????%前N幀用于背景建模
????????movie（h1?frC?1?15?[0s（1）00]）;%播放第i幀
????????cdataSum?=?cdataSum?+?f;????
????else
????????%開始運動檢測，首先計算出背景
????????if（c1==N+1）cdB?=?cdataSum/N;
????????else??colorFrC1?=?M（c1-N）;%獲取第c1-N幀
??????????????frC1?=?rgb2gray（frame2im（colorFrC1））;
??????????????frC1?=?im2frame（frC1?gray（256））;
??????????????f1?=?int16（frC1.cdata）;
??????????????cdataSum=cdataSum+f-f1;%背景更新
??????????????cdB=cdataSum/N;
??????????????cdB?=?int16（cdB）;
????????end
????????f2?=?abs（?f?-?cdB?）;
????????f2?=?uint8（f2）;
????????bwf?=?f2;
????????bw?=?find（?f2?????????bwf（bw）?=?0;
????????bw?=?find（?f2?>=?45?）;????
????????bwf（bw）?=?255;?????
????????cdB?=?uint8（cdB）;???????
??