运动估计算法MATLAB课程设计.docx

1、运动估计算法MATLAB课程设计课程设计任务书学生姓名: 专业班级： 指导教师: 工作单位： 题 目： 运动估计算法 初始条件： LAB软件平台要求完成的主要任务：1.设计任务设计视频压缩系统中的运动估计算法:全搜索法(F: Full Search）和三步法(TSS: hreSepSach),比较二种方法的搜索点和每帧的峰值信噪比（PNR: peak sgnal o noe ratio）2.设计要求 编制算法代码;对视频进行运动估计;计算PSNR时间安排: 答辩时间2013年1月24日。指导教师签名: 年 月 日系主任签名： 年 月 日摘 要 11运动估计算法概念1.1 运功估计算法基本思想

2、21. 运动估计算法实验原理 22 设计原理和方法3.2新三步法323 全搜索法2.4 峰值信噪比3运动估计算法的MTLB编程 63.2三步法 93全搜索法指标.4三步法指标 13.4仿真结果分析14 小结与体会 11参考文献11附录2摘 要运动估计的基本思想是尽可能准确地获得序列图像帧间的运动位移,即运动矢量。因为运动估计越准确，预测补偿的图像质量越高，补偿的残差就越小，补偿编码所需位数越少,需要传输的比特率就越小。利用得到的运动矢量在帧间进行运动补偿。补偿残差经过变换、量化、编码后与运动矢量一起经过熵编码,然后以比特流形式发送出去。在视频编码和处理系统中,运动估计和运动补偿技术对降低视频序

3、列时间冗余度、提高编码效率起着非常关键的作用。运动估计的准确程度将直接决定视频编码器的编码效率。关键词：运动估计、运动补偿技术、位移(运动）矢量AbstrcThebic da i that themtion estimaion ccura as osl he imae sunceinterfremoion diplaceent, ie.the otnvector. Motion stimatio more ccate rin nstion，the hhe te mage qualitis comped resials smaler, le copenstion oding itequired,

4、 thmaler the transmision it te. erringmoto oensn uing he moin vecto otained in he ntrfe Cmpensat residus througtasformation，qntizaion, entopy-ced tgher withte moon vetr is encd, nte sntout in h formf a bi streIvdeo dg androcesing sstem, te mton stimtion andmotin cmenstion to reuc empol reudancy o vi

5、deseuen timprove thecodin efficiency play ccia roleTh degree of acurcy f toioesimao will detly detrmin he ncing eiciecy the veencder.ewords： Motion stimatio Motin compatio chniques Te vctor ofdipcment (movemt)多媒体信息处理1 运动估计算法概念视频原始图像中存在着大量的信息冗余,如时间冗余、空间冗余、信息熵冗余、谱间冗余、几何结构冗余、视觉冗余和知识冗余等等。运动估计是视频压缩编码中的核心

6、技术之一，采用运动估计和运动补偿技术可以消除视频信号的时间冗余以提高编码效率。如何提高运动估计的效率，使运动估计算法的搜索过程更健壮、更快速、更高效成为目前研究的热点。运动估计的基本思想是尽可能准确地获得序列图像帧间的运动位移，即运动矢量。因为运动估计越准确,预测补偿的图像质量越高,补偿的残差就越小，补偿编码所需位数越少，需要传输的比特率就越小。利用得到的运动矢量在帧间进行运动补偿。补偿残差经过变换、量化、编码后与运动矢量一起经过熵编码,然后以比特流形式发送出去。运动估计算法多种多样,大体上可以把它们分成四类：块匹配法、递归估计法、贝叶斯估计法和光流法。其中块匹配运动估计算法因其具有算法简单、

7、便于SI实现等优点得到广泛应用。所以本文将重点介绍块匹配运动估计算法,并对各种块匹配算法在计算速度和估计精度上进行简单比较。1. 运功估计算法基本思想运动估计的基本思想是尽可能准确地获得序列图像帧间的运动位移，即运动矢量。运动估计算法的目标是效率和准确性。由于在成象的场景中一般有多个物体作不同的运动，如果直接按照不同类型的运动将图像分割成复杂的区域是比较困难的。最直接和不受约束的方法是在每个像素都指定运动矢量，这就是所谓基于像素表示法。这种表示法是对任何类型图像都是适用的,但是它需要估计大量的未知量,并且它的解时常在物理上是不正确,除非在估计过程中施加适当的物理约束。这在具体实现时是不可能的,

8、通常采用基于块的物体运动表示法。因为运动估计越准确，预测补偿的图像质量越高，补偿的残差就越小，补偿编码所需位数越少，需要传输的比特率就越小。1 运动估计算法实验原理在帧间预测编码中,由于活动图像邻近帧中的景物存在着一定的相关性。因此，可将活动图像分成若干块或宏块,并设法搜索出每个块或宏块在邻近帧图像中的位置，并得出两者之间的空间位置的相对偏移量，得到的相对偏移量就是通常所指的运动矢量,得到运动矢量的过程被称为运动估计。 运动矢量和经过运动匹配后得到的预测误差共同发送到解码端,在解码端按照运动矢量指明的位置,从已经解码的邻近参考帧图像中找到相应的块或宏块，和预测误差相加后就得到了块或宏块在当前帧

9、中的位置。运动估计的准确程度往往用补偿图像与原图像比较的PSN来衡量表示。2 设计原理和方法2.1三步法三步法是应用得相当广泛的一种次优的运动估计搜索算法它的搜索区间一般为-7，7 即在候选区中与编码块相同坐标位置处为原点,将参考块在其上下左右距离为7的范围内按照一定规律移动移到一个位置就做匹配计算它总共进行了三步搜索在下一次搜索时步长减半以前一步搜索得到的最优点为中心。下图为三步法的搜索示意图。算法的中心思想是，采用一种由粗到细的搜索模式，从原点开始，按一定步长取周围8个点构成每次搜索的点群，然后进行匹配计算,利用上一步搜索得到的最小块误差MB点作为当前搜索的中心位置，每做一步，搜索的步长减

10、1。步搜索算法搜索窗选取(-7，+7)，最多只需要做25个位置的匹配计算，相对于全搜索来比，大大减少了匹配运算的复杂度，而且数据读取比较规则。2.2 新三步法TS 假定运动矢量分布特点是在搜索窗口中均匀分布, 但事实证明运动矢量是偏置中心的，Reniang Li等人在SS的基础上提出了一种增强运动矢量中心偏置搜索和减小补偿误差的新三步法。NTS是对TSS 的一个改进,对运动量比较小的视频序列如可视电话序列有比较好的性能。对于绝大多数的视频序列，运动矢量的分布都是在中心位置上的概率最大，随着与中心位置的距离的增大,概率会急剧地下降，这也就是前面所说的运动矢量的中心偏移特性。运动量比较小的视频序列

11、的这一特性会更加明显。NTS算法在最好的情况下只需要做1个点的匹配，在最坏的情况下需要做33个点的匹配,由于运动矢量中心偏置在现实视频序列中是普遍存在的,在通常情况下，NTS算法需要做33点匹配的概率比较小,因此,在低速率视频应用中,如视频电话或视频会议中，NSS算法的优点可以得到较好的发挥。 全搜索法全搜索法(FulSehMethod,F）也称为穷尽搜索法，是对（ d ）（N2dy ）搜索范围内所有可能的候选位置计算MAD (i,j)值,从中找出最小MAD，其对应偏移量即为所求运动矢量。此算法虽计算量大,但最简单、可靠,找到的必为全局最优点。S算法描述如下:从原点出发，按顺时针螺旋方向由近及

12、远,在逐个像素处计算MAD值，直到遍历搜索范围内听有的点,然后在计算的所有点的MAD中找到最小值,该点所在位置即对应最佳运动矢量。但是正因为它是穷尽搜索因此会产生巨大的计算量如 7， 7的搜索区间每个宏块116 需计算225个MAD值,这就直接制约了编码的实时实现。快速算法本质上是一种穷尽搜索法其计算量仍是相当巨大的。全搜索算法是简单也是效果最好的一种匹配算法,通过的全搜索匹配得到的结果是全局最优的，但由于计算量很大，我们在编解码中往往不采用这种方法,而只把他作为与其他算法的一种比较。图21程序流程图2.4 峰值信噪比峰值信噪比（经常缩写为PSR）是一个表示信号最大可能功率和影响它的表示精度的

13、破坏性噪声功率的比值的工程术语。由于许多信号都有非常宽的动态范围,峰值信噪比常用对数分贝单位来表示。PSN是“Pak Sgnal tNoise Rati”的缩写。peak的中文意思是顶点。而raio的意思是比率或比列的。整个意思就是到达噪音比率的顶点信号，psnr是一般是用于最大值信号和背景噪音之间的一个工程项目。通常在经过影像压缩之后，输出的影像通常都会有某种程度与原始影像不一样。为了衡量经过处理后的影像品质，我们通常会参考PSNR 值来认定某个处理程序够不够令人满意。PN计算公式如下:eak 就是指 bits 表示法的最大值55。E指eanquae Error(均方误差，各值相差的n次方和

14、的平均值的n次平方根(这几个字应该没有),I（角标n)指原始影像第个pixel 值,(角标n)指经处理后的影像第n个pxel 值。PR 的单位为dB。所以PSNR值越大,就代表失真越少。PSN 是最普遍,最广泛使用的评鉴画质的客观量测法,不过许多实验结果都显示，PSNR的分数无法和人眼看到的视觉品质完全一致，有可能 PN 较高者看起来反而比 SNR 较低者差。这是因为人眼的视觉对于误差的敏感度并不是绝对的,其感知结果会受到许多因素的影响而产生变化(例如:人眼对空间频率较低的对比差异敏感度较高，人眼对亮度对比差异的敏感度较色度高,人眼对一个区域的感知结果会受到其周围邻近区域的影响)。3运动估计算

15、法的ATLB编程31全搜索法 图3-1-1 ist fame 图31-2 heeond fram 图13 帧间差值 图3-1-4DFD图3-1- 恢复后的第二帧图像图3-6第一帧图像的运动矢量图图3-7 第二帧图像的运动矢量图32三步法 图-1 th rt fme 图3-2- he secondfame 图-23 帧间差值 图3-2 FD图3- 恢复后的第二帧图像图26 第一帧图像的运动矢量图图3-2-第二帧图像的运动矢量图3.3全搜索法指标全搜索算法：耗时长为lpsed tme is 1.53320 econd.3.4三步法指标三步法：耗时长为Eased tm is 0.88349 seco

16、n.34仿真结果分析根据峰值信噪比可知全搜索算法的效果比三步法的效果要好；但全搜索所消耗的时间要比三步法所消耗的时间长,其运行较慢，因为全搜索算法的搜索点数要比三步法的搜索点数多很多。 小结与体会 通过本次对全搜索块匹配算法课程设计，让我更深刻的理解了多媒体信息的获取、搜索和保存的基本方法。基于块的运动估计,是先将视频图像分成一个个规则的图像块,然后对每个图像块估计运动矢量。基于块的运动估计和运动补偿已经广泛应用于各种视频压缩编码标准。因此，本次课程设计是非常有作用的，让我们清醒的认识到自己还有好大的差距,它让我了解了多媒体处理的基本原理。课程设计中会遇到很多很多的困难，但在同学的帮助和自己思

17、考下最终还是没有解决问题。所以,本次课程设计我个人的收获不大,但是让我知道自己还在图像编程方面有很多欠缺的地方,不仅对均方误差(SE)、绝对值误差(MAD)更深的理解，并且采用全搜索块匹配方法得到所有宏块的运动矢量，把c语言的知识又重新的过了一遍,同时还理解了些多媒体信息处理的基本原理。参考文献1罗军辉.MATLAB7.在图像处理中的应用. 机械出版社.0052王家文 MATB7.编程基础.机械出版社.0钟麟 MTLA仿真技术与应用教程.国防工业出版社20冈萨雷斯数字图像处理(MATL版)电子工业出版社.205附录:clar all；1=rea(clre1.bmp); %read e firs

18、 amI2=mad(clre2.bmp); %ed thecnd frame;%给图像扩边，每个边都扩大dm大小% -|-|-|% B | C | D |% |-|-|-|% | | | |% | | | | | E | A | -图像 | | | | | | | | | |% |-|-|-|% G I |% |-|-|% Idouble（I1)；I2=dobe（2);rowm clnm =size(I1）； I=zr（ronu*dm,clnu+*);II(d+1：mronum,dm+：dm+onm)=I1;fo i=1:dm I(,dm1：dmcom）=II(m+,dm:+olnu); II

19、(rownum+dm+i，dm+:dm+olu)=II（dm+rownu,dm+1：m+clnu);ndfo j:dm I(1:rnum+2dm，j)=II(1:rwnum+2m,m+1); I(1:rowum+*dm，on+dm+j）II(1:ronum+*dm，dmcolnum);end %下面进行全搜索算法ocksize=1;rowbock =ownu/blcksiz; coloks =clum/locksize;A999999999999999; %为了找到最小的均方误差,A用于设定一个很大的初值 Ej=0; xrecord=(16，6)； xecod，record用于存放匹配快的块号

20、，即运动矢量 yrecdnes(1，16);d=ze（256，56); %这幅图的大小为2*256 tic for x=:(olok） %x表示行中第几个子块 rowx*blockiz； for =0:(colbocks-) %y表示列中第几个子块 ol=y*bocksiz;% temx=xblocki+1：（x+1）lksize;% tempy=y*blckze+1：(+1）*locsie; op=-dm:dm f q=-dm：d %(,q)表示x,对应子块在前一帧所的搜索位置 Ei=0; Eij=sum(I2（row+1：ro+bckiz，l+1:col+bloksiz)-II(+dm+p

21、:w+dm+p+bloksz,colmq+1：col+m+blocksz).)（lockse); f ijA =; xrord(+1，y）=p； yord(x1,y1）=q； en en end A9999999999999; for x=1：blocksize or ny=1:cksz diff(o+m，cly)I2(ow+mx，l+ny）-II(row+mx+dm+xrcord(x+,y+1),co+y+d+yrecord(x+1,y1)； ed nd end endtoc igur,imhow(I,); title(tertframe); fiure,imshow（I2，)； itle（

22、thescnd fme）; II=I2-I1； fre,imshw（IIII,); title（帧间差值); fur，imshow(diff，); tle（DFD)； ttle(利用全搜索算法匹配后的帧间差); for x=:(owblcks-1） r=x*loize； f =：(collok-1) % col=y*lciz; III（ro+1:ro+blockize,co+:col+bocksze)=I(row+d+xord(1,y)+1:row+m+xrcord(x+1,+1)+bocsize,+m+rcod(x+,y+1)1：cdm+yrecord（x+1,y+1)+blocks)+di

23、f（ow+:r+ksie，co+1:col+blckize)； nd d %II=I1as(diff）; igr,imsh(III,）； title(); ERRiff；% igure,imsho(ERR,);% ttle（F); mberara=0：:255; for 1：255 nmerary(m+)=0； e; zeronumb0; fn=:ronu or m=1:clnm dfab(ER（,n）); if(d=) mp=zronumber; zeonumer=temp1; ee numbear(dif)=nmerrray(df）1; ed; end; d; igure；plt(0,zeronumber，k*）;hod n;plot（nuerry，r),tite(DDdistribuion);hold f; ER=zeos(6,); for i=: for =0：15 ER1(+1,+1)=ound（s(u(ERR(i*blksize+1:i*blocksie+blocs