基于块的全搜索运动估计算法实现实验报告

《基于块的全搜索运动估计算法实现实验报告》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、数字视频处理实验报告学院：通信与信息工程学院系班：电信科0901班姓名：学号：时间：2012年11月23号一、实验名称：基于块的全搜索运动估计算法实现二、实验目的：1、掌握运动估计算法的实现原理。2、掌握运动估计算法的研究现状及多种计算方法。3、学习基于块的全搜索运动估计算法，研究分析其Matlab实现程序过程，并补充完成程序，对实验结果进行分析比较。三、实验要求三、实验要求1、对实验程序motionEstAnalysis.m进行分析，完成主程序流程图。函数流程图：2、编写补充完成部分不全程序代码，调试程序使其能正确运行(1)motionEstES()%Computesmoti

2、onvectorsusingexhaustivesearchmethod（全搜索法计算运动矢量）%%Input%imgP:Theimageforwhichwewanttofindmotionvectors（当前图像）%imgI:Thereferenceimage（参考图像）%mbSize:Sizeofthemacroblock（宏块尺寸）%p:Searchparameter(readliteraturetofindwhatthismeans)（搜索参数）%%Ouput%motionVect:themotionvectorsforeachintegralmacroblockini

3、mgP（当前图像中每一个积分宏块的运动矢量）%EScomputations:Theaveragenumberofpointssearchedforamacroblock（每个宏块搜索的平均点数）%%WrittenbyArohBarjatyafunction[BlockCenter,motionVect,EScomputations]=motionEstES(imgP,imgI,mbSize,p)%定义函数文件motionEstES.m，imgP、imgI、mbSize、p为传入参数，BlockCenter、motionVect、EScomputations为返回参数[rowco

4、l]=size(imgI);%将参考图像的行数赋值给row，列数赋值给colblockcenter=zeros(2,row*col/mbSize^2);vectors=zeros(2,row*col/mbSize^2);%定义全0的矢量矩阵的大小costs=ones(2*p+1,2*p+1)*65537;%定义最小绝对差矩阵的大小computations=0;%搜索点数赋初值为0%westartofffromthetopleftoftheimage（从图像左上角开始）%wewillwalkinstepsofmbSize（以宏块尺寸为步长）%foreverymarcoblockt

5、hatwelookatwewilllookfor%aclosematchppixelsontheleft,right,topandbottomofit（对于每一个宏块，在它的上下左右找到与搜索参数p最匹配的像素）mbCount=1;%搜索的宏块数赋初值为1%1为循环起始值，mbSize为步长值，row-mbSize+1为循环终止值fori=1:mbSize:row-mbSize+1forj=1:mbSize:col-mbSize+1%theexhaustivesearchstartshere（全搜索开始）%wewillevaluatecostfor(2p+1)blocksver

6、tically%and(2p+1)blockshorizontaly（我们将计算水平方向上(2p+1)个块的最小绝对差和垂直方向上(2p+1)个块的最小绝对差）%misrow(vertical)index（m为行指数）%niscol(horizontal)index（n为列指数）%thismeanswearescanninginrasterorderform=-p:p%水平方向上位移矢量范围forn=-p:p%垂直方向上位移矢量范围%补充下面程序%row/Vertco-ordinateforrefblock（参考块的行(垂直方向)的范围）refBlkVer=i+m;%col/H

7、orizontalco-ordinate（参考块的列(水平方向)的范围）refBlkHor=j+n;%如果参考块的行列范围的任意一个在已经搜索过的宏块之外，则继续下一步的搜索if(refBlkVer<1

8、

9、refBlkVer+mbSize-1>row...

10、

11、refBlkHor<1

12、

13、refBlkHor+mbSize-1>col)continue;endcosts(m+p+1,n+p+1)=costFuncMAD(imgP(i:i+mbSize-1,j:j+mbSize-1),...imgI