MATLAB视频运动目标检测参考源码

MATLAB视频运动目标检测参考源码

将视频截取成帧

fileName = 'q.avi'; %此处输入视频名称，要带后缀，注意好像只有wmv和avi格式

obj = VideoReader(fileName);

numFrames = obj.NumberOfFrames;% 帧的总数

for k = 1 :numFrames% 读取数据

frame = read(obj,k);

imshow(frame);

imwrite(frame,strcat(num2str(k),'.jpg'),'jpg');% 保存帧

End

灰度二值化图像

for i=1:numberofframe%numberofframe为图像张数

imageName=strcat(num2str(i),'.jpg');

I = imread(imageName);

imggray = rgb2gray(I);%灰度化图像

imwrite(imggray,strcat(num2str(i),'huidu.jpg'),'jpg');% 保存帧

i2=im2bw(I);% i2是二值图像，不需要求阈值

imwrite(i2,strcat(num2str(i),'erzhi.jpg'),'jpg');% 保存帧

end

背景差分法

filename = '103.jpg';%读原图

I = imread(filename);

X = rgb2gray(I);%读，灰度原图

filename = 'YCBJ.jpg';%读取背景图

M = imread(filename);

Y = rgb2gray(M);%读，灰度背景图

Z = imabsdiff(X,Y);%两图相减

level = graythresh(Z);

Z = im2bw(Z,level);%二值化图

% 膨胀腐蚀算子，相当于闭操作

se = strel('disk',2);

Q=imerode(Z,se);%腐蚀

Q=imdilate(Q,se);%膨胀

lab1='原图';

lab4='腐蚀膨胀后';

subplot(221);imshow(I);title(lab1);

subplot(222);imshow(M);title(lab2);

subplot(223);imshow(Z);title(lab3);

subplot(224);imshow(Q);title(lab4);

三帧差分法

filename1='110.jpg';%读取前一帧

I1 = imread(filename1);

O = rgb2gray(I1);%读，灰度前一帧

filename2='115.jpg';%读取中间帧图像

I2 = imread(filename2);

P = rgb2gray(I2);%读，灰度当前帧

filename3='120.jpg';%读取后一帧图像

I3 = imread(filename3);

Q = rgb2gray(I3);%读，灰度后一帧

A = imabsdiff(P,O);

level1 = graythresh(A);

A = im2bw(A,level1);%二值化前两帧差图

subplot(131);imshow(A);

B = imabsdiff(Q,P);

level = graythresh(B);

B = im2bw(B,level);%二值化后两帧差图

subplot(132);imshow(B);

for i=1:480

for j=1:720

C(i,j)=A(i,j)*B(i,j);

end

end

subplot(133);imshow(C);%二AB两图相与

帧差分循环读写

for i=1:numberofframe%numberofframe为帧数

filename=strcat(num2str(i),'.jpg');

I = imread(filename);

X = rgb2gray(I);%读，灰度原图

filename = 'YCBJ.jpg';

M = imread(filename);

Y = rgb2gray(M);%读，灰度背景图

Z = imsubtract(Y,X);%两图相减

level = graythresh(Z);

Z = im2bw(Z,level);%二值化图

% 膨胀腐蚀算子，相当于闭操作

imwrite(Z,strcat(num2str(i),'Chafenghou.jpg'),'jpg');% 保存帧

end

掩模处理

D = imread('zhan.jpg');

%掩模处理

for j=1:720

for i=1:480

if(j<-2.39*i+765)

D(i,j)=0;

end

if(j>0.8*i+460)

D(i,j)=0;

end

if(i<140)

D(i,j)=0;

end

end

end

图像对称翻转

AA=imread('b.jpg');%文件名是你自己的哦

for k=1:3

B(:,:,k)=fliplr(AA(:,:,k));%左右翻转

end

imwrite(B,'obj.jpg','jpg');% 保存帧

简单的判断摔倒

BJ=imread('YCBJ.jpg');% 读取有物体彩色图片

fileName = '6.avi';

obj = VideoReader(fileName);

numFrames = obj.NumberOfFrames;% 帧的总数

for k = 1 :numFrames% 读取数据

T = read(obj,k);

A=imsubtract(BJ,T);

level = graythresh(A);

B = im2bw(A,level);%二值

for j=1:720

for i=1:480

if(j<-2.39*i+765)

B(i,j)=0;

end

if(j>0.8*i+460)

B(i,j)=0;

end

if(i<140)

B(i,j)=0;

end

end

end

% imwrite(B,strcat(num2str(k),'.jpg'),'jpg');% 保存帧

figure(1);

imshow(B);

[bwl,ln]=bwlabel(B,8);

STATS=regionprops(bwl,'Area');

a=0;

for n=1:ln

if(STATS(n).Area>500)

abc=(bwl==n);

[c,r]=find(abc==1);

a=a+1;

xbegin=min(r);

ybegin=min(c);

xlength=max(r)-xbegin;

ylength=max(c)-ybegin;

p=ylength/xlength;

if(p<=1)

hold on

rectangle('Position',[xbegin,ybegin,xlength,ylength],'EdgeColor','R','linewidth',3);

hold off

end

end

end

end

改善后的判断摔倒

BJ=imread('YCBJ.jpg');% 读取有物体彩色图片

fileName = '13.wmv';

obj = VideoReader(fileName);

numFrames = obj.NumberOfFrames;% 帧的总数

H=zeros(1,1);

S=zeros(1,1);

for k = 1:numFrames% 读取数据

T = read(obj,k);

A=imsubtract(BJ,T);

level = graythresh(A);

B = im2bw(A,level);%二值

for j=1:720

for i=1:480

if(i<140)

B(i,j)=0;

end

end

end

% imwrite(B,strcat(num2str(k),'.jpg'),'jpg');% 保存帧

figure(1);

imshow(T);

[bwl,ln]=bwlabel(B,8);

STATS=regionprops(bwl,'Area');

a=0;

for n=1:ln

if(STATS(n).Area>50)

abc=(bwl==n);

[c,r]=find(abc==1);

xbegin=min(r);

ybegin=min(c);

xlength=max(r)-xbegin;

ylength=max(c)-ybegin;

p=ylength/xlength;

if(p<=1)

x=max(c);

y = 39.916*x-10881;

y1=y*0.8566;

y2=y*1.1434;

if(STATS(n).Area<y2&&STATS(n).Area>y1)

a=a+1;

H(k,a)=max(c);

S(k,a)=STATS(n).Area;

hold on

rectangle('Position',[xbegin,ybegin,xlength,ylength],'EdgeColor','R','linewidth',3);

hold off

end

end

end

end

End

将目标区域值存入矩阵数据分析

BJ=imread('YCBJ.jpg');% 读取有物体彩色图片

fileName = '1.avi';

obj = VideoReader(fileName);

numFrames = obj.NumberOfFrames;% 帧的总数

H=zeros(1,1);

S=zeros(1,1);

for k = 1:numFrames% 读取数据

T = read(obj,k);

A=imsubtract(BJ,T);

level = graythresh(A);

B = im2bw(A,level);%二值

imshow(B);

% imwrite(B,strcat(num2str(k),'.jpg'),'jpg');% 保存帧

[bwl,ln]=bwlabel(B,8);

STATS=regionprops(bwl,'Area');

a=0;

for n=1:ln

if(STATS(n).Area>400)

abc=(bwl==n);

[c,r]=find(abc==1);

a=a+1;

H(k,a)=max(c);%将联通区域底部坐标放到矩阵H中

S(k,a)=STATS(n).Area;%将联通区域面积坐标放到矩阵S中

end

end

end