【跌倒检测】基于matlab帧差法跌倒检测【含Matlab源码 320期】

一、简介

摄像机采集的视频序列具有连续性的特点。如果场景内没有运动目标，则连续帧的变化很微弱，如果存在运动目标，则连续的帧和帧之间会有明显地变化。
帧间差分法(Temporal Difference)就是借鉴了上述思想。由于场景中的目标在运动，目标的影像在不同图像帧中的位置不同。该类算法对时间上连续的两帧或三帧图像进行差分运算，不同帧对应的像素点相减，判断灰度差的绝对值，当绝对值超过一定阈值时，即可判断为运动目标，从而实现目标的检测功能。

1.2 三帧差分法


综上所述，帧间差分法的原理简单，计算量小，能够快速检测出场景中的运动目标。但由实验结果可以看出，帧间差分法检测的目标不完整，内部含有“空洞”，这是因为运动目标在相邻帧之间的位置变化缓慢，目标内部在不同帧图像中相重叠的部分很难检测出来。帧间差分法通常不单独用在目标检测中，往往与其它的检测算法结合使用。

二、源代码

clear all clc
global Bili;
A=[];
B=imread('46.jpg');
B=rgb2gray(B);
B=medfilt2(B);        %中值滤波
fileName = '倾斜跌倒1.avi'; 
obj = VideoReader(fileName);           %读取录制的视频
vidFrames=read(obj);                 %读取所有帧图像
numFrames=get(obj, 'NumberOfFrames');%总帧数    
for k = 63:250                   %帧数循环，人为选取33帧作为背景帧
mov(k).cdata=vidFrames(:,:,:,k);
mov(k).cdata=rgb2gray(mov(k).cdata);
mov(k).cdata=medfilt2(mov(k).cdata);
R=imabsdiff(mov(k).cdata,B);
z=im2bw(R,graythresh(R));
z=bwmorph(z,'erode',3);                       % 3次腐蚀处理掉噪点
P = bwmorph(z,'dilate',3);                    %3次膨胀膨胀
P= bwareaopen(P,50);     
figure(1),subplot(121);imshow(vidFrames(:,:,:,k));
figure(1);subplot(122);imshow(P);
x=sum(P(:));
s=x/(480*640);
if s>=0.03
   [m1,n1] = find(P == 1); 
%第一级检测,角度
   top1= min(m1);%纵向
   bottom1= max(m1);
   left1= min(n1);%横向
   right1= max(n1);
   height1=bottom1-top1+1;%纵向高；
   width1 = right1-left1+1;%横向宽
   rectangle('Position',[left1,top1,width1,height1],'EdgeColor','r');
   line([left1,right1],[top1,bottom1],'color','r','LineWidth',1);%标记对角线
   line([left1,right1],[bottom1,top1],'color','r','LineWidth',1);
   tan=(bottom1-top1)/(right1-left1);
   tan=atan(tan);
   J=tan*180/pi;
%第二级检测，质心高度比Rz
   [rectx,recty,area,perimeter] = minboundrect(n1,m1,'p');
   line(rectx(:),recty(:),'color','w','LineWidth',1);
   line([rectx(1),rectx(3)],[recty(1),recty(3)],'color','b','LineWidth',1);
   line([rectx(2),rectx(4)],[recty(2),recty(4)],'color','b','LineWidth',1);
   line([0,640],[max(bottom1),max(bottom1)],'color','r','LineWidth',2);%水平线
   line([(rectx(1)+rectx(3))/2,(rectx(1)+rectx(3))/2],[(recty(1)+recty(3))/2,max(bottom1)],'color','b','LineWidth',2);%画高
   h=max(bottom1)-(recty(1)+recty(3))/2;%高度
   Rz=h/360;
%第三级检测，判断宽高比
   width=sqrt((rectx(2)-rectx(1))^2+(recty(2)-recty(1))^2);
   length=sqrt((rectx(1)-rectx(4))^2+(recty(1)-recty(4))^2);
   Bili=width/length;
%判断
    if J<=45
        if(Rz<0.5||Bili>2)
%           imwrite(vidFrames(:,:,:,k),strcat(num2str(k),'.jpg'),'jpg');% 保存帧
            load chirp
            sound(y,Fs)
        end
    end
end
end

三、运行结果

四、备注

2014a