%% 帧间差分法
close all;% 清除所有工作空间的变量
clear all;% 关闭所有窗口
clc; % 清除命令窗口的内容
n=1; %变量
mov=VideoReader('boy-running.mp4'); %读取视频
nNUM = mov.NumberOfFrame; %得到视频总帧数
background=double(rgb2gray(read(mov,1)));% read读取帧图像，rgb2gray()把图像转换为灰度图像
for k=n+1:nNUM % n+1到总帧数循环域
CurrentImage=double(rgb2gray(read(mov,k))); % 当前帧的灰度图像
Formerimage=double(rgb2gray(read(mov,k-1))); % 前一帧的灰度图像
ID=imadd(CurrentImage,Formerimage); % 帧间差分=运动物体前后所在区域
th=graythresh(ID); % 最大类间方差法找到图片的一个合适的阈值
BW=im2bw(ID,th); %通过上条阈值二值化图像
% figure;
% subplot(121);imshow(ID);
% subplot(122);imshow(BW);
% 背景更新
alpha=0.1;
CurrentBack=background.*BW + (alpha.*CurrentImage + (1-alpha).*background).*(1-BW);
% 前后两帧有变化区域不更新，无变化区域 更新到背景中
background=CurrentBack; % 背景更新完毕
%---------弄半天，前面只是在给背景更新。。。。-----------------------------
cut=abs(CurrentImage-background);
cut=uint8(cut);
cth =graythresh(cut);
cBW=im2bw(cut,cth);
% figure;
% subplot(121);imshow(BW);
% subplot(122);imshow(cBW);
se=[1 1 1;1 1 1;1 1 1]; % 结构元素
BWCuterode=imerode(cBW,se);
BWCuterodil=bwmorph(BWCuterode,'dilate',3); % 膨胀
% BWCuterodil=bwareaopen(BWCuterodil,2000,8); % 去掉连通区域面积小于2000的区域
% figure;
% subplot(131);imshow(read(mov,k));
% subplot(132);imshow(BWCuterodil);
% 找出图形中的八连通区域，视为目标所在区域
[L,nm]=bwlabel(BWCuterodil,8);
% 找出最大区域
stats=regionprops(L); % 测量标记矩阵L中每一个标注区域的一系列属性 ，返回stats是一个结构数组
% 整合
Ar=cat(1,stats.Area); % 1-按 列 连接所有输入数组，列数相同；2-按 行
% 最大区域
[~,ind_max_ar]=max(Ar);
% 提取目标
rect=stats(ind_max_ar).BoundingBox; % 最小外接矩形,形式 [x y w h]
% frame=double(rgb2gray(read(mov,k)));
frame=read(mov,k);
imshow(frame);
rectangle('position',rect,'EdgeColor','r');
pause(0.001);
end
%% 需要解决的问题是---最后的几帧图像识别错误，已经没有目标，但是出现 虚景目标
% 其次，目标定位 的边框目标太大，或者太小
评论0