回想当初学习时做的一个帧差法检测运动物体
核心程序基本流程：
1.VideoReader读取视频文件
2.获取视频长度
3.构建循环
4.读取连续两帧
5.计算运动的四个边界点
6.画出当前帧
7.在当前帧中绘制框
demo:
mov=VideoReader('traffic.avi'); %读取视频文件
N=mov.NumFrames; %读取视频长度
%%
figure
for i=1:N-1
frame=read(mov,i+1);
Pframe=read(mov,i); % 读取第一帧视频帧
subplot(1,2,1)
imshow(Pframe);
[Xmin,Ymin,Xmax,Ymax]=FrameDiff(frame,Pframe); %输出运动目标的四个边界点
subplot(1,2,2)
imshow(Pframe); %画出当前帧
rectangle('position',[Xmin,Ymin,Xmax-Xmin,Ymax-Ymin],'edgecolor','g'); %画出跟踪框
FrameDiff
FrameDiff的操作流程：
1.前后两帧图像灰度化
2.中值滤波去噪
3.图像均衡化
4.将图像转化为图像双精度
5.检测帧差
6.中值滤波
7.设置帧差阈值
8.转化为逻辑阵
9.找出边缘像素点
function [Xmin,Ymin,Xmax,Ymax]=FrameDiff(frame,Pframe)
%%
x=rgb2gray(frame);
y=rgb2gray(Pframe); %前后两帧图像灰度化
%%
x=medfilt2(x);
y=medfilt2(y);%中值滤波去噪
%%
x = histeq(x);
y = histeq(y);% 图像均衡化
%%
n=im2double(x);
p=im2double(y); %将图像转化为图像双精度
%%
c=n-p; %检测帧差
c=medfilt2(c); %中值滤波
t=80/256; %设置帧差阈值
c(abs(c)>=t)=255;
c(abs(c)<t)=0;
c=logical(c); %将c转为逻辑阵
[a,b]=find(c>0);
Xmax=max(b);
Xmin=min(b);
Ymax=max(a);
Ymin=min(a); %找出边缘像素点
