数字视频3.docVIP

《数字视频处理》 课内实验报告 2015-2016学年 第2学期 专业： 智能科学与技术 班级： 智能1301班 学号：19） 姓名： 姚衎 时间： 2016年5月26日 实验三 序列图像的运动估计实验 一、实验目的 进一步熟悉MATLAB软件开发环境。学习在MATLAB编程环境下实现对标准序列图像的前后两帧图像间进行运动估计，并统计得到的运动矢量。 二、实验要求 利用MATLAB图像处理工具箱中的函数，在MATLAB编程环境下，用穷尽块匹配算法实现序列图像的帧间运动估计。将标准序列图像的前后两帧图像文件存入计算机，运行MATLAB集成开发环境；编程实现序列图像的整象素精度穷尽块匹配算法。 三、实验内容 （1）在MATLAB上录入程序代码： （附）程序源代码如下： A=imread(1.bmp); B=imread(2.bmp); f2=rgb2gray(A); f1=rgb2gray(B); fp=0; figure,imshow(f2),title(目标帧); figure,imshow(f1),title(锚定帧); N=16; R=16; height=256; width=256; for i=1:N:height-N+1 for j=1:N:width-N+1 MAD_min=256*N*N; dy=0; dx=0; for k=-R:1:R for l=-R:1:R if i+k1 MAD=256*N*N; elseif i+kheight-N MAD=256*N*N; elseif j+l1 MAD=256*N*N; elseif j+lwidth-N MAD=256*N*N; else MAD=sum(sum(abs(double(f1(i:i+N-1,j:j+N-1))-double(f2(i+k:i+k+N-1,j+l:j+l+N-1))))); end if MADMAD_min MAD_min=MAD; dy=k; dx=l; end; end; end; fp(i:i+N-1,j:j+N-1)=f2(i+dy:i+dy+N-1,j+dx:j+dx+N-1); iblk=floor((i-1)/N+1); jblk=floor((j-1)/N+1); mvx(iblk,jblk)=dx; mvy(iblk,jblk)=dy; end; end; figure,imshow(uint8(fp)),title(预测帧); [X,Y]=meshgrid(N/2:N:256-N/2); Y=256-Y; figure,quiver(X,Y,mvx,mvy),title(整像素精度的运动场); diff=abs(double(f1)-fp); figure,imshow(uint8(diff)),title(diff); PSNR=20*log10(255)+20*log10(256)-10*log10(sum(sum(diff.*diff))); （2）运行结果： 四、实验总结 通过这次的序列图像的运动估计实验，我利用了MATLAB图像处理工具箱中的函数，在MATLAB编程环境下，用穷尽块匹配算法实现了序列图像的帧间运动估计。本次实验之后，我进一步熟悉了MATLAB软件开发环境。同时，我还学习了在MATLAB编程环境下实现对标准序列图像的前后两帧图像间进行运动估计，并经过统计得到了图像的运动矢量。我相信，这次实验在整个《数字视频处理》这门课当中有着十分重要的作用，并会为我今后的学习夯实好基础。