clc; clear all; close all;
% 原始视频
targetavi = 'traffic.avi';
% 检测结果视频
resultavi = 'e:\\result.avi';
%% 读取视频
mov = VideoReader(targetavi);
fnum = mov.NumberOfFrames;
%% 建立结果视频
aviobj = avifile(resultavi);
aviobj.Quality = 100;
%%aviobj.Fps = 25;
%%aviobj.compression = 'Indeo5';
%% 帧间差分法
figure(1);
for i = 2 : fnum
x = read(mov, i-1);
y = read(mov, i);
subplot(1, 2, 1); imshow(x, []); title(sprintf('第%d帧视频', i-1), 'FontWeight', 'Bold', 'Color', 'r');
% 灰度化
if ndims(x) == 3
m = rgb2gray(x);
else
m = x;
end
if ndims(y) == 3
n = rgb2gray(y);
else
n = y;
end
% 中值滤波
m = medfilt2(m);
n = medfilt2(n);
% 数据类型转换
q = im2double(n);
w = im2double(m);
% 差分
c = q-w;
% 阈值，此值可以调节
t = 80/256;
% 阈值分割
c(abs(c)>=t)=255;
c(abs(c)<t) = 0;
c = logical(c);
x1 = x(:, :, 1); x2 = x(:, :, 2); x3 = x(:, :, 3);
x1(c) = 255; x2(c) = 0; x3(c) = 0;
xc = cat(3, x1, x2, x3);
subplot(1, 2, 2); imshow(xc, []); title(sprintf('第%d帧视频识别结果', i-1), 'FontWeight', 'Bold', 'Color', 'r');
f = getframe(gcf);
f = frame2im(f);
% 生成视频
aviobj=addframe(aviobj,f);
z=read(mov,2)
h=rgb2gray(z)
s=read(mov,fnum-1)
l=rgb2gray(s)
j=im2double(h)
k=im2double(l)
g=k-j
t=30
if g<=t
text(7,8,'无故障')
else
text(7,8,'故障')
end
end
%% 关闭视频句柄
%aviobj = close(aviobj);
