本程序来源于《MATLAB图像处理:能力提高与应用案例》这本书籍,这个书籍是适合初学matlab的新手加深对于matlab的学习与理解。
本段为主程序:
clear data
disp('input video');
% 读入视频图像wzw.avi并显示
avi = aviread('wzw.avi');
video = {avi.cdata};
for a = 1:length(video)
imagesc(video{a});
axis image off
drawnow;
end;
disp('output video');
% 调用tracking()函数对运动目标进行跟踪
tracking(video);
以下为子程序:
function d = tracking(video)
% 功能:跟踪视频中的运动目标并显示
% 输入:video-待跟踪的视频
% 输出:d-差值图像序列
% 读入视频图像
if ischar(video)
avi = aviread(video);
pixels = double(cat(4,avi(1:2:end).cdata))/255;
clear avi
else
pixels = double(cat(4,video{1:2:end}))/255;
clear video
end
% 将RGB图像转换成灰度图像
nFrames = size(pixels,4);
for f = 1:nFrames
pixel(:,:,f) = (rgb2gray(pixels(:,:,:,f)));
end
[rows,cols]=size(pixel(:,:,1));
nrames=f;
% 将相邻两帧图像进行作差,并将差值图像转换为二值图像
for l = 2:nrames
d(:,:,l)=(abs(pixel(:,:,l)-pixel(:,:,l-1)));
k=d(:,:,l);
bw(:,:,l) = im2bw(k, .2);
bw1=bwlabel(bw(:,:,l));
imshow(bw(:,:,l))
hold on
% 标记运动物体的位置并显示
cou=1;
for h=1:rows
for w=1:cols
if(bw(h,w,l)>0.5)
toplen = h;
if (cou == 1)
tpln=toplen;
end
cou=cou+1;
break
end
end
end
disp(toplen);
coun=1;
for w=1:cols
for h=1:rows
if(bw(h,w,l)>0.5)
leftsi = w;
if (coun == 1)
lftln=leftsi;
coun=coun+1;
end
break
end
end
end
disp(leftsi);
disp(lftln);
widh=leftsi-lftln;
heig=toplen-tpln;
widt=widh/2;
disp(widt);
heit=heig/2;
with=lftln+widt;
heth=tpln+heit;
wth(l)=with;
hth(l)=heth;
disp(heit);
disp(widh);
disp(heig);
rectangle('Position',[lftln tpln widh heig],'EdgeColor','r');
disp(with);
disp(heth);
plot(with,heth, 'r*');
drawnow;
hold off
end
本程序利用帧间差分法能够实现对单个目标的运动进行检测,并且用方框提示运动的目标。对于上述代码所提及视频,建议自己拍摄一个小视频,时间在5秒左右适宜,太大的话会引起内存溢出,视频的格式有严格的要求 ,须为 avi中的zjmedia uncompress rgb24,可以通过winavi来转换,如果视频出现报错disp(toplen),那就是视频不符合本程序的要求,程序代码经过本人多次的检验,能够正确的在matlab7.0中运行。