clear?all;
?
%讀入視頻數(shù)據(jù)
mov?=??VideoReader（‘C:\Users\Administrator\Desktop\畢設(shè)\test.avi‘）;
nb?=mov.NumberOfframes;%視頻幀數(shù)

T=10;%運(yùn)動(dòng)目標(biāo)檢測(cè)閾值??可變
P=0.05;%運(yùn)動(dòng)目標(biāo)占比閾值??可變
for?i=2:nb
????x=read（mov?i）;%讀取視頻幀
????y=read（mov?i-1）;
????img=y;
????m=rgb2gray（x）;%灰度圖轉(zhuǎn)換
????n=rgb2gray（y）;

????
????m=medfilt2（m）;%中值濾波
????n=medfilt2（n）;%中值濾波
????



????
????c=abs（m-n）;%差分
%?????c=medfilt2（c）;%中值濾波
????
????c（abs（c）>=T）=255;%運(yùn)動(dòng)目標(biāo)檢測(cè)
????c（abs（c）????c=logical（c）;%轉(zhuǎn)換為0，1
????

????c=bwareaopen（c20）;%刪除小面積對(duì)象
????se1=strel（‘square‘3）;
????se2=strel（‘square‘13）;??
?
?????c=imdilate（cse2）;%形態(tài)學(xué)膨脹處理
????run=c;
??????
?
????%運(yùn)動(dòng)目標(biāo)提取
????c=uint8（c）;
????sh（::1）=c.*img（::1）;%點(diǎn)乘操作?
????sh（::2）=c.*img（::2）;
????sh（::3）=c.*img（::3）;
????imshow（sh）;
????title（sprintf（‘運(yùn)動(dòng)目標(biāo)提取?第%d幀‘i-1））;

????[Lbw4?numbw4]?=?bwlabel（run）;??????%??%Lbw4為貼標(biāo)簽之后的矩陣，numbw4為標(biāo)簽個(gè)數(shù)
????stats?=?regionprops（Lbw4）;%獲取區(qū)域的某個(gè)屬性（面積、最小包圍矩形的坐標(biāo)長(zhǎng)寬等）的值??
?????for?i?=?1?:?numbw4??
??????tempBound?=?stats（i）.BoundingBox?;?
??????%計(jì)算運(yùn)動(dòng)目標(biāo)區(qū)域占比
??????moving=find（c==i）;
??????moving=length（moving）;
??????mn=size（c）;
??????all_area=mn（1）*mn（2）;
??????pro=moving/all_area;
??????if?pro?>P
??????????rectangle（‘position‘tempBound‘edg