视音频信号处理实验3.docVIP

浙江理工大学 《视音频信号处理》 实验报告 2012～2013学年第一学期 学 院 信息学院 班 级 10数媒(1)班 姓 名 黄芳 学 号 任课教师 马文娟 （上课时间：周三1-2节） 数字媒体技术专业 2012年 12 月 《视音频信号处理》实 验 报 告 学号： 姓名：黄芳 班级：10数媒1班 一、实验名称 视频数据处理 二、实验内容 1、 请编写图像矩阵产生程序，具体要求为： a) 共产生两个图像矩阵A[256][256],和B[256][256] b) 图像矩阵中的点A(I,j) = i*256+j,B(I,j) = i*256+j,； c) 特别地，当25=i25+16, 20=j20+16, A(i,j) = i - 20; d) 特别地，当32=i32+16, 16=j16+16, B(i,j) = i - 27; e) 把两幅图像显示出来。 2、以题1中得到的图像数据基础上，以图像A为参考帧，请用三步搜索方法计算图像B中16*16图像块MBlock(32,16)的运动矢量。 三、实验步骤 1、编写图像矩阵产生程序 主要代码： hfView.cpp void CHfView::OnDraw(CDC* pDC) { CHfDoc* pDoc = GetDocument(); ASSERT_VALID(pDoc); char A[256][256],B[256][256]; int j; FILE *fa, *fb; CDC * dc; dc=GetDC(); for(int i=0;i256;i++)//生成2个图像矩阵A和B，并把两幅图像显示出来 { for(j=0; j256;j++) { if((i41i=25)(j36j=20)) {A[i][j]=i-20; dc-SetPixel(i,j,i-20);} else {A[i][j]=i*256+j; dc-SetPixel(i,j,i*256+j);} if((i48i=32)(j32j=16)) {B[i][j]=i-27; dc-SetPixel(i+300,j,i-27);} else {B[i][j]=i*256+j; dc-SetPixel(i+300,j,i*256+j);} } } //把生成的像素矩阵存入A.raw和B.raw文件中 if(!(fa=fopen(A.raw,wb))) printf(error to open A.raw\n); if(!(fb=fopen(B.raw,wb))) printf(error to open B.raw\n); fwrite(A,sizeof(A),1,fa); fwrite(B,sizeof(B),1,fb); fclose(fa); fclose(fb); } 结果： 2、以题1中得到的图像数据基础上，以图像A为参考帧，请用三步搜索方法计算图像B中16*16图像块MBlock(32,16)的运动矢量。 代码： #includestdio.h #includestdlib.h struct Block{ public: int x,y; char Array[16][16]; }; int SAD(Block Block1,Block Block2){ int i=0,j=0,sad=0; for(i=0;i16;i++) for(j=0;j16;j++) sad+=abs(Block2.Array[i][j]-Block1.Array[i][j]); return sad; } Block GetBlock(int x,int y,char tmp[256][256]) { int i,j; Block Blocktmp; Blocktmp.x=x; Blocktmp.y=y; for(i=0;i16;i++) for(j=0;j16;j++) Blocktmp.Array[i][j]=tmp[i+x][j+y]; return Blocktmp; } void ThreeStepSearch(int x,int y,Block Block1,char A[256][256],char B[256][256]){ int sad=9999999,temp,flag=0; int a,b; int c=x,d=y; int a0,a1,b0,b1; Block Block2; Block2=GetBl