运动估计算法MATLAB课程设计.docx
《运动估计算法MATLAB课程设计.docx》由会员分享，可在线阅读，更多相关《运动估计算法MATLAB课程设计.docx（26页珍藏版）》请在冰豆网上搜索。
运动估计算法MATLAB课程设计
课程设计任务书
学生姓名:
专业班级：
指导教师:
工作单位：
题目：
运动估计算法
初始条件：
ﻩ ＭＡＴLAB软件平台
要求完成的主要任务：
1.设计任务
设计视频压缩系统中的运动估计算法:
全搜索法（FＳ:
FullSearch）和三步法（TSS:
Ｔhreｅ Sｔep Sｅaｒch）,比较二种方法的搜索点和每帧的峰值信噪比（PＳNR:
peaksｉgnalｔonoｉｓeratio）
2.设计要求
编制算法代码;对视频进行运动估计;计算PSNR
时间安排:
答辩时间2013年1月24日。
指导教师签名:
年 月 日
系主任签名：
年 月 日
摘要1
1 运动估计算法概念ﻩ２
1.1运功估计算法基本思想2
1.２运动估计算法实验原理2
2设计原理和方法ﻩ3
２.2 新三步法ﻩ3
2．3全搜索法ﻩ４
2.4峰值信噪比ﻩ５
3 运动估计算法的MＡTLＡB编程6
3.2三步法9
3．３全搜索法指标ﻩ１１
３.4三步法指标１1
3.4仿真结果分析ﻩ１1
4小结与体会11
参考文献ﻩ11
附录ﻩ１2
摘要
运动估计的基本思想是尽可能准确地获得序列图像帧间的运动位移,即运动矢量。
因为运动估计越准确，预测补偿的图像质量越高，补偿的残差就越小，补偿编码所需位数越少,需要传输的比特率就越小。
利用得到的运动矢量在帧间进行运动补偿。
补偿残差经过变换、量化、编码后与运动矢量一起经过熵编码,然后以比特流形式发送出去。
在视频编码和处理系统中,运动估计和运动补偿技术对降低视频序列时间冗余度、提高编码效率起着非常关键的作用。
运动估计的准确程度将直接决定视频编码器的编码效率。
关键词：
运动估计、运动补偿技术、位移（运动）矢量
Abstrａcｔ
The bａｓicｉdｅaiｓthatthe mｏtionestimaｔionａｓａccuraｔｅasｐoｓsｉｂlｅｔheimaｇesｅｑuｅnce interfrａｍe moｔiondiｓplaceｍent,i．e. theｍotｉｏn vector.Motionｅstimatioｎmoreａccｕｒateｐrｅｄｉｃｔiｏnｃｏｍｐｅnsａtion， thehｉｇheｒtｈeｉmagequalitｙ iscompｅｎｓａｔedresiｄｕalsｉｓsmalｌer,leｓｓcoｍpensａtionｃodingｂitｓ ｒequired,thｅ ｓmaｌlerthetransmiｓsionｂitｒａte.Ｐerｆｏrｍing motｉoｎｃoｍｐensａｔｉｏnuｓingｔhemoｔiｏnvectoｒoｂtainedinｔheｉntｅrfｒａｍe．Cｏmpensatｉｏｎresiduａｌsthrougｈ tｒaｎsformation， qｕａntizaｔion,entｒopy-cｏｄedtｏgｅｔherwith tｈemoｔｉonveｃtｏrisencｏdｅｄ,ａnｄ tｈeｎsｅnt outinｔhｅform ｏfabiｔstreａｍ．
Iｎ vｉdeoｃｏdｉｎgand ｐroceｓsingsｙstem,tｈemｏtｉonｅstimａtionand motiｏncｏmｐensａtiontoreｄucｅｔｈｅ ｔempoｒａlreｄuｎdancyoｆvideｏ seｑuenｃｅtｏ improvethe codinｇefficiencyplayｓａcｒｕciaｌrole． Thｅdegreeofacｃurａcyｏftｈｅ ｍoｔioｎ esｔimaｔｉoｎwilldｉｒeｃtlydetｅrminｅｔheｅncｏｄingeｆｆicieｎcyｏｆthevｉｄeｏ encｏder.
Ｋeｙwords：
MotionｅstimatioｎMotiｏncompｅｎｓatioｎｔｅchniquesTｈevｅctorof diｓpｌａcｅment（movemｅｎt）
多媒体信息处理
1运动估计算法概念
视频原始图像中存在着大量的信息冗余,如时间冗余、空间冗余、信息熵冗余、谱间冗余、几何结构冗余、视觉冗余和知识冗余等等。
运动估计是视频压缩编码中的核心技术之一，采用运动估计和运动补偿技术可以消除视频信号的时间冗余以提高编码效率。
如何提高运动估计的效率，使运动估计算法的搜索过程更健壮、更快速、更高效成为目前研究的热点。
运动估计的基本思想是尽可能准确地获得序列图像帧间的运动位移，即运动矢量。
因为运动估计越准确,预测补偿的图像质量越高,补偿的残差就越小，补偿编码所需位数越少，需要传输的比特率就越小。
利用得到的运动矢量在帧间进行运动补偿。
补偿残差经过变换、量化、编码后与运动矢量一起经过熵编码,然后以比特流形式发送出去。
运动估计算法多种多样,大体上可以把它们分成四类：
块匹配法、递归估计法、贝叶斯估计法和光流法。
其中块匹配运动估计算法因其具有算法简单、便于ＶＬSI实现等优点得到广泛应用。
所以本文将重点介绍块匹配运动估计算法,并对各种块匹配算法在计算速度和估计精度上进行简单比较。
1.１运功估计算法基本思想
运动估计的基本思想是尽可能准确地获得序列图像帧间的运动位移，即运动矢量。
运动估计算法的目标是效率和准确性。
由于在成象的场景中一般有多个物体作不同的运动，如果直接按照不同类型的运动将图像分割成复杂的区域是比较困难的。
最直接和不受约束的方法是在每个像素都指定运动矢量，这就是所谓基于像素表示法。
这种表示法是对任何类型图像都是适用的,但是它需要估计大量的未知量,并且它的解时常在物理上是不正确,除非在估计过程中施加适当的物理约束。
这在具体实现时是不可能的,通常采用基于块的物体运动表示法。
因为运动估计越准确，预测补偿的图像质量越高，补偿的残差就越小，补偿编码所需位数越少，需要传输的比特率就越小。
1．２运动估计算法实验原理
在帧间预测编码中,由于活动图像邻近帧中的景物存在着一定的相关性。
因此，可将活动图像分成若干块或宏块,并设法搜索出每个块或宏块在邻近帧图像中的位置，并得出两者之间的空间位置的相对偏移量，得到的相对偏移量就是通常所指的运动矢量,得到运动矢量的过程被称为运动估计。
运动矢量和经过运动匹配后得到的预测误差共同发送到解码端,在解码端按照运动矢量指明的位置,从已经解码的邻近参考帧图像中找到相应的块或宏块，和预测误差相加后就得到了块或宏块在当前帧中的位置。
运动估计的准确程度往往用补偿图像与原图像比较的PSNＲ来衡量表示。
2设计原理和方法
2.1 三步法
三步法是应用得相当广泛的一种次优的运动估计搜索算法它的搜索区间一般为[-7，7]即在候选区中与编码块相同坐标位置处为原点,将参考块在其上下左右距离为7的范围内按照一定规律移动移到一个位置就做匹配计算它总共进行了三步搜索在下一次搜索时步长减半以前一步搜索得到的最优点为中心。
下图为三步法的搜索示意图。
算法的中心思想是，采用一种由粗到细的搜索模式，从原点开始，按一定步长取周围8个点构成每次搜索的点群，然后进行匹配计算,利用上一步搜索得到的最小块误差MBＤ点作为当前搜索的中心位置，每做一步，搜索的步长减1。
步搜索算法搜索窗选取（-7，+7），最多只需要做25个位置的匹配计算，相对于全搜索来比，大大减少了匹配运算的复杂度，而且数据读取比较规则。
2.2新三步法
TＳS假定运动矢量分布特点是在搜索窗口中均匀分布,但事实证明运动矢量是偏置中心的，RenｘiangLi等人在ＴSS的基础上提出了一种增强运动矢量中心偏置搜索和减小补偿误差的新三步法。
NTSＳ 是对TSS的一个改进,对运动量比较小的视频序列如可视电话序列有比较好的性能。
对于绝大多数的视频序列，运动矢量的分布都是在中心位置上的概率最大，随着与中心位置的距离的增大,概率会急剧地下降，这也就是前面所说的运动矢量的中心偏移特性。
运动量比较小的视频序列的这一特性会更加明显。
NTＳS算法在最好的情况下只需要做1７个点的匹配，在最坏的情况下需要做33个点的匹配,由于运动矢量中心偏置在现实视频序列中是普遍存在的,在通常情况下，NTSＳ算法需要做33点匹配的概率比较小,因此,在低速率视频应用中,如视频电话或视频会议中，NＴSS算法的优点可以得到较好的发挥。
２．３全搜索法
全搜索法（Fuｌl Seａｒｃh Method, FＳ）也称为穷尽搜索法，是对（Ｍ＋２dｘ）×（N ＋2dy）搜索范围内所有可能的候选位置计算MAD（i,j）值,从中找出最小MAD，其对应偏移量即为所求运动矢量。
此算法虽计算量大,但最简单、可靠,找到的必为全局最优点。
ＦS算法描述如下:
从原点出发，按顺时针螺旋方向由近及远,在逐个像素处计算MAD值，直到遍历搜索范围内听有的点,然后在计算的所有点的MAD中找到最小值,该点所在位置即对应最佳运动矢量。
但是正因为它是穷尽搜索因此会产生巨大的计算量如［7，7]的搜索区间每个宏块1６＊16需计算225个MAD值,这就直接制约了编码的实时实现。
快速算法本质上是一种穷尽搜索法其计算量仍是相当巨大的。
全搜索算法是简单也是效果最好的一种匹配算法,通过的全搜索匹配得到的结果是全局最优的，但由于计算量很大，我们在编解码中往往不采用这种方法,而只把他作为与其他算法的一种比较。
ﻩ
图2－1程序流程图
2.4峰值信噪比
峰值信噪比（经常缩写为PSＮR）是一个表示信号最大可能功率和影响它的表示精度的破坏性噪声功率的比值的工程术语。
由于许多信号都有非常宽的动态范围,峰值信噪比常用对数分贝单位来表示。
PSNＲ是“PｅakSｉgnaltｏ NoiseRatiｏ”的缩写。
peak的中文意思是顶点。
而raｔio的意思是比率或比列的。
整个意思就是到达噪音比率的顶点信号，psnr是一般是用于最大值信号和背景噪音之间的一个工程项目。
通常在经过影像压缩之后，输出的影像通常都会有某种程度与原始影像不一样。
为了衡量经过处理后的影像品质，我们通常会参考PSNR值来认定某个处理程序够不够令人满意。
PＳNＲ计算公式如下:
Ｐeak就是指８bits表示法的最大值２55。
ＭＳE 指 Ｍean ＳquaｒeError（均方误差，各值相差的n次方和的平均值的n次平方根（这几个字应该没有））,I（角标n）指原始影像第ｎ个pixel值,Ｐ（角标n）指经处理后的影像第n个pｉxel值。
PＳＮR的单位为dB。
所以PSNR值越大,就代表失真越少。
PSNＲ是最普遍,最广泛使用的评鉴画质的客观量测法,不过许多实验结果都显示，PSNR 的分数无法和人眼看到的视觉品质完全一致，有可能PＳNＲ较高者看起来反而比ＰSNR较低者差。
这是因为人眼的视觉对于误差的敏感度并不是绝对的,其感知结果会受到许多因素的影响而产生变化（例如:
人眼对空间频率较低的对比差异敏感度较高，人眼对亮度对比差异的敏感度较色度高,人眼对一个区域的感知结果会受到其周围邻近区域的影响）。
3 运动估计算法的ＭATLＡB编程
3．1全搜索法
图3-1-1ｔｈｅｆiｒstfｒame 图3－1-2ｔhe ｓeｃondframｅ
图３－1－3 帧间差值 图3-1-4 DFD
图3-1-５恢复后的第二帧图像
图3－１-6 第一帧图像的运动矢量图
图3－１-7第二帧图像的运动矢量图
3．2三步法
图３-２-1thｅｆｉrｓtfｒａme 图3-2-２ｔhesecond fｒame
图３-2－3帧间差值 图3-2－４ＤFD
图3-２-５恢复后的第二帧图像
图３－2－6第一帧图像的运动矢量图
图3-2-７ 第二帧图像的运动矢量图
3.3全搜索法指标
全搜索算法：
耗时长为Ｅlａpsedtｉmeis1.０53320ｓecondｓ.
3.4三步法指标
三步法：
耗时长为Eｌaｐsedtｉmｅis0.88１349seconｄｓ.
3．4仿真结果分析
根据峰值信噪比可知全搜索算法的效果比三步法的效果要好；但全搜索所消耗的时间要比三步法所消耗的时间长,其运行较慢，因为全搜索算法的搜索点数要比三步法的搜索点数多很多。
４小结与体会
通过本次对全搜索块匹配算法课程设计，让我更深刻的理解了多媒体信息的获取、搜索和保存的基本方法。
基于块的运动估计,是先将视频图像分成一个个规则的图像块,然后对每个图像块估计运动矢量。
基于块的运动估计和运动补偿已经广泛应用于各种视频压缩编码标准。
因此，本次课程设计是非常有作用的，让我们清醒的认识到自己还有好大的差距,它让我了解了多媒体处理的基本原理。
课程设计中会遇到很多很多的困难，但在同学的帮助和自己思考下最终还是没有解决问题。
所以,本次课程设计我个人的收获不大,但是让我知道自己还在图像编程方面有很多欠缺的地方,不仅对均方误差（ＭSE）、绝对值误差（MAD）更深的理解，并且采用全搜索块匹配方法得到所有宏块的运动矢量，把c语言的知识又重新的过了一遍,同时还理解了些多媒体信息处理的基本原理。
参考文献
［1]罗军辉. MATLAB7.０在图像处理中的应用.机械出版社.２005
［2]王家文．MATＬＡB7.０编程基础.机械出版社.２0０５
［３］钟麟．MＡTLAＢ仿真技术与应用教程.国防工业出版社．20０４
[４]冈萨雷斯． 数字图像处理（MATLＡＢ版）．电子工业出版社.20０5
附录
ＦＳ:
clｅarall；
Ｉ1=ｉｍreaｄ（'clａｉre1.bmp'）;%readｔｈefirsｔｆｒamｅ
I2=ｉmｒｅad（'clａｉre2.bmp＇）; %ｒeａdthｅ ｓecｏndframe
ｄｍ＝７;
%¸给图像扩边，每个边都扩大dm大小
%%%％%%%%%%%%%%％％%%%％％%%%%%%%％%%%%%％％%%％%％%％%%%％%%%%%％%％％％％%%%%%%％％%％%％%％%%%%%%%%%%%
％ ｜－----|-----------－----－---－--－－--－--－|--－---|
% ｜ B | C |D |
% |－---－|－--－--－－-----－-－------－－--－----|--－---|
% | | | |
% | | | |
％ | E| A |Ｆ ｜---－＞ 图像
％ | ｜ | |
％ || | ｜
％ | | | |
% |-－---|－----------－---------－-----－－-－|-－----|
% ｜ G｜ Ｈ ｜I |
% |-----|----------------－---－----－----－｜－-----|
％%%%%％％%%%%%％%%%%%%％%％％%%%%%％%%%％％％%%%％％%%%％%%%％%%％%%%%%%%%%％%%%%%％%%%%%%％％％%%%%%%%%％
I１＝double（I1）；
I2=doｕbｌe（Ｉ2）;
[rowｎｕmcｏlnｕm］= size（I1）；
IＩ=zｅrｏｓ（roｗnuｍ＋２*dm,cｏlnuｍ+２*ｄｍ）;
II（dｍ+1：
ｄm＋roｗnum,dm+１：
dm+ｃoｌnｕm）=I1;
foｒi=1:
dm
IＩ（ｉ,dm＋1：
dm＋coｌｎｕm）=II（ｄm+１,dm＋１:
ｄｍ+ｃolnuｍ）;
II（rownum+dm+i，dm+１:
dm+ｃolｎuｍ）=II（dm+rownuｍ,dm+1：
ｄm+cｏlnuｍ）;
ｅnd
foｒj＝１:
dm
ＩI（1:
rｏｗnum+2＊dm，j）=II（1:
rｏwnum+2＊ｄm,ｄm+1）;
ＩI（1:
rowｎum+２*dm，ｃoｌnｕｍ+dm+j）＝II（1:
roｗnum+２*dm，dm＋colnum）;
end
%下面进行全搜索算法
ｂｌocksize=1６;
rowbｌockｓ=ｒownuｍ/blｏcksizｅ;
coｌｂloｃks=cｏlｎum/ｂlocksize;
A＝9999９999999９99９９999９; %为了找到最小的均方误差,A用于设定一个很大的初值
Eｉj=0;
xrecord=ｏｎｅｓ（16，１6）； ％xｒecoｒd，ｙrecord用于存放匹配快的块号，即运动矢量
yrecｏｒd＝ｏnes（1６，16）;
dｉｆｆ=zeｒｏｓ（256，２56）; %这幅图的大小为2５６*256
tic
forx=０:
（ｒoｗｂloｃkｓ－１） %x表示行中第几个子块
row＝x*blockｓizｅ；
forｙ=0:
（colbｌocks-１） %y表示列中第几个子块
ｃol=y*bｌocksizｅ;
% temｐx=x＊blockｓiｚｅ+1：
（x+1）＊ｂlｏｃksize;
% tempy=y*blｏckｓｉze+1：
（ｙ+1）*ｂlocｋsiｚe;
ｆoｒ p=-dm:
dm
fｏｒq=-dm：
dｍ%（ｐ,q）表示x,ｙ对应子块在前一帧所的搜索位置
Eiｊ=0;
Eij=sum（ｓｕｍ（（I2（row+1：
roｗ+bｌｏckｓizｅ，ｃｏl+1:
col+bloｃksizｅ）-II（ｒｏｗ+dm+p＋１:
ｒｏw+dm+p+bloｃksｉzｅ,col＋ｄm＋q+1：
col+ｄm+ｑ+blocksｉzｅ））.^２））／（ｂlocksｉｚe^２）;
ｉfＥij Ａ=Ｅｉｊ;
xrｅｃord（ｘ+1，y＋１）=p；
yｒｅｃord（x＋1,y＋1）=q；
enｄ
enｄ
end
A＝９99９9999９9９99９9999;
forｍx=1：
blocksize
ｆorny=1:
ｂｌｏcksｉzｅ
diff（ｒoｗ+mｘ，cｏl＋ｎy）＝I2（ｒow+mx，ｃｏl+ny）-II（row+mx+dm+xrｅcord（x+１,y+1）,coｌ+ｎy+dｍ+yrecord（x+1,y＋1））；
eｎd
ｅnd
end
end
toc
ｆigurｅ,imｓhow（I１,[]）;
title（'tｈe ｆｉrｓt frame'）;
fiｇure,imshow（I2，［］）；
ｔitle（＇the sｅcｏndfｒａme＇）;
ＩIＩI=I2-I1；
fｉｇｕre,imshｏw（IIII,[]）;
title（‘帧间差值'）;
fｉｇurｅ，imshow（diff，[]）;
tｉｔle（＇DFD'）；
％tｉtle（'利用全搜索算法匹配后的帧间差'）;
forx=０:
（ｒowblｏcks-1） ％
rｏｗ=x*ｂloｃｋｓize；
fｏｒｙ=０：
（colｂloｃkｓ-1） %
col=y*ｂlｏcｋｓizｅ;
III（roｗ+1:
roｗ+blockｓize,coｌ+１:
col+bｌocksｉze）=ＩI（row+dｍ+xｒｅｃord（ｘ＋1,y＋１）+1:
row+ｄm+xrｅcord（x+1,ｙ+1）+bｌocｋsize,ｃｏｌ+ｄm+ｙrｅcoｒd（x+１,y+1）＋1：
cｏｌ＋dm+yrecord（x+1,y+1）+blocksｉｚｅ）+diｆf（ｒow+１:
rｏｗ+ｂｌｏｃksiｚe，coｌ+1:
col+blｏckｓize）；
ｅnd
ｅｎd
%ＩII=I1＋aｂs（diff）;
ｆigｕrｅ,imshｏｗ（III,[]）；
title（'»Ö¸´ºóµÄµÚ¶þ֡ͼÏñ'）;
ERR＝ｄiff；
% ｆigure,imshoｗ（ERR,[]）;
% tｉtle（'ＤFＤ＇）;
ｎｕmberaｒraｙ=0：
１:
255;
forｍ＝1：
255
nｕmｂerarｒａy（m+１）=0；
eｎｄ;
zeronumbｅｒ＝0;
fｏｒ n=１:
roｗnuｍ
ｆorm=1:
cｏlnｕm
dｉf＝abｓ（ERＲ（ｍ,n））;
if（dｉｆ==０）
ｔｅmp=zｅronumber;
zeｒonumｂer=temp＋1;
eｌｓe
numbeｒaｒrａｙ（dif）=nｕmｂerａrray（dｉf）＋1;
eｎd;
end;
ｅｎd;
ｆigure；plｏt（0,zeronumber，＇k*'）;hoｌdｏn;plot（nuｍｂeｒａrrａy，'r＊'）,titｌe（'DＦD distribuｔion'）;holdｏｆf;
EＲR１=zeｒos（１6,１６）;
fori=０:
１５
forｊ=0：
15
ERＲ1（ｉ+1,ｊ+1）=ｒound（sｕｍ（ｓuｍ（ERR（i*blｏｃksize+1:
i*blocksiｚe+blocｋs
升级会员 