背景差分法MATLAB实现

程序1
背景差分法MA TLAB实现
function temp3
d=60;
b='d\6\capfile3.avi';
e='.bmp';
for i=0:13
u=d+1;
m=int2str(d);
n=int2str(u);
s=strcat(b,m,e); %连接字符
m=imread(s);
m=rgb2gry(m); %灰度值
% imshow(m);
if(d>=61)
s=strcat(b,n,e);
n=imread(s);
n=rgb2gry(n);
q=im2double(m);
figure(18);imshow(q);
w=in2double(n);
figure(19);imshow(w)
c=q-w; %进行差分处理
figure(20);imshow(c)
figure(99);imshow(c);%绘制直方图
t=100;%此値可以调节（阈值调整）
t=t/255;
k=find(abs(c)>=t);
c(k)=1;
k=find(abs(c)<t);
c(k)=0;
figure;imshow(c);
sc=strel('square',5);%用平坦的5*5结构元素进行了膨胀c=imdilate(c,sc);
end
d=d+1;
end
程序2
图像处理：
f0=imread('capfile3.avi61.bmp');
g=rgb2gray(f0);%图像灰度化
imshow(g);
g2=imnoise(g,'salt & pepper',0.05);%加淑盐噪声
figure(3);imshow(g2);
g3=medfilt2(g2,[3,3],'symmtric');%以3*3窗口中值滤波
figure(4);imshow(g3);
g1=im2bw(g);%图像二值化
figure(2);imshow(g1);
程序3
% “高速公路汽车阴影”视屏的阴影去除的方法
clc;
clear all;
close all;
% 打开视频文件对话框
h1=figure('toolbar','name','Shadow Removal');
[FileName.PathName]=uigetfile({'*.avi','avi(*.avi)','*.*',
'All Files(*.*)'},'Open video sequence');
if isequal([FileName,PathName],[0,0])
return
else
pic_path=fullfile(PathName,FileName);
mov=aviread(pic_path);
end
n=100; % 读取视频帧数
s=1; % 起始为1
e=n; % 结束为n
for i=s:e
mov_gray(i).image=mov(i).cdate; % mov_gray(i).image 为第i帧图像
end
[ih,iw,id]=size(mov_gray(1).image); % 取得视频帧的大小
% 计算1到第n帧的各个像素的均值
for i=s:e
u=mov_gray(i).image;
h=u(:)';
t(i,:)=h(:);
end
tmp=median(t); % 对t数组各列取均值
mean=reshape(tmp,ih,iw,id); % 将n帧均值图像恢复成h*w*d的图像矩阵，得到背景
figure,imshow(uint8(mean));title('背景图像');
fn=18;
ok=mean;
tmp(1:3)=0;
di=20;
th=ones(1,1,3)*di; % 预设阈值
diff=double(ok)-double(mov_gray(fn).image); % 计算第fn帧图像与背景的差值diff1=abs(diff);
for ai=1:ih
for bi=1:iw
if abs(diff(ai,bi,1:3))>=th(1,1,:)
ok(ai,bi,1:3)=mov_gray(fn).image(ai,bi,1:3);
else
ok(ai,bi,1:3)=0;
diff1(ai,bi,1:3)=0;
end
end
end
figure,imshow(uint8(ok));title('前景图像');
impixelinfo;
g=rgb2hsv(mean); % 将背景图像由RGB格式转为HSV格式
figure,imshow(g);title('背景图像的hsv格式');
impixelinfo;
figure,imtooi(g);
g1=rgb2hsv(ok); % 将前景图像由RGB格式转为HSV格式
figure,imshow(g1);title('前景图像的hsv格式');
impixelinfo;
figure,imtooi(g1);
% 阈值分割
% 这个地方的阈值选取是手动的，根据经验得到的！
for ci=1:ih
if(g1(ci,di,3)/g(ci,di,3)>=0.4) &&(g1(ci,di,3)/g(ci,di,3)<=0.5)
&& abs(g1(ci,di,2)-g(ci,di,2))<=0.5 && abs(g1(ci,di,1)-g(ci,di,1))<=0.1;
h1(ci,di,1:3)=1;
else
h1(ci,di,1:3)=0;
end
end
end
h=hsv2rgb(h1);
figure,imshow(h);title('阴影去除结果')；。