视音频实验3.doc

任课教师：张海翔 《视音频信号处理》（2013-2014学年第1学期） 实 验 报 告 学号： 姓名： 班级： 实验三 视频数据处理 一、实验内容1、请编写图像矩阵产生程序，具体要求为： 共产生两个图像矩阵A[256][256],和B[256][256] 图像矩阵中的点A(I,j) = i*256+j,B(I,j) = i*256+j,； 特别地，当25=i25+16, 20=j20+16, A(i,j) = i - 20; 特别地，当32=i32+16, 16=j16+16, B(i,j) = i - 27; （256*256中有16*16的块） 以题1中得到的图像数据基础上，以图像A为参考帧，请用全搜索方法计算图像B中16*16图像块MBlock(32,16)的运动矢量。（-7，+4） （在256*256上循环检查完全相等的） 以题1中得到的图像数据基础上，以图像A为参考帧，请用三步搜索方法计算图像B中16*16图像块MBlock(32,16)的运动矢量。 给定一个8*8的DCT系数矩阵C[8][8], 12 -8 -1 0 0 0 0 0 0 0 -1 0 0 0 0 0 -1 -1 0 0 0 0 0 0 0 -1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 请编写程序将DCT系数按Zigzag扫描顺序转换为系数序列，并用游程编码对系数序列进行编码，以连续‘0’的个数，非零系数值的格式输出的中间符号。 中间符号用printf打出 二、实验代码 1、图形矩阵产生程序 #includestdio.h #includestdlib.h int main(){ char A[256][256],B[256][256]; FILE *filea,*fileb; if(!(filea=fopen(a.raw,wb))){ printf(can not open file.\n); } if(!(fileb=fopen(b.raw,wb))){ printf(can not open file.\n); } for(int i=0;i256;i++){ for(int j=0;j256;j++){//当=i25+16, 20=j20+16, A(i,j) = i - 20; if((i=25i41)(j=20j36)) {A[i][j]=i-20;} else A[i][j] = i*256+j; } } for(int i=0;i256;i++){ for(int j=0;j256;j++){//当=i25+16, 20=j20+16, A(i,j) = i - 20; if((i=32i48)(j=16j32)) {B[i][j]=i-27;} else B[i][j] = i*256+j; } } fwrite(A,sizeof(A),1,filea); fwrite(B,sizeof(B),1,fileb); fclose(filea); fclose(fileb); return 0; } 全搜索法 #includestdio.h #includestdlib.h int main(){ int A[256][256],B[256][256]; FILE *filea,*fileb; if(!(filea=fopen(a.raw,wb))){ printf(can not open file.\n); } if(!(fileb=fopen(b.raw,wb))){ printf(can not open file.\n); } for(int i=0;i256;i++){ for(int j=0;j256;j++){//当=i25+16, 20=j20+16, A(i,j) = i - 20; if(25=ii4120=jj36) {A[i][j]=i-20;} else A[i][j] = i*256+j; }} for(int i=0;i256;i++){ for(int j=0;j256;j++){//当=i25+16, 20=j20+16, A(i,j) = i - 20; if(32=ii4816=jj32) {B[i][j]=i-27;} else B[i][j] = i*256+j; }} for(int m=0;m256;m++){ for(int n=0;n2