project3_xiaohuavivian_DCT变换.docVIP

多媒体技术实验报告 项目名称： DCT变换 提交文档学生姓名： xiaohuavivian 提交文档学生学号： 教师评阅成绩： 教师评阅意见： . . 提交报告时间： 2014 年 12 月 10 日 实验题目: 掌握DCT变换及DCT反变换的原理及步骤 了解DCT变换后图像的特点 实验要求: 1） 读入给定的BMP格式图像，如Fig7.bmp, Fig8.bmp，要求通过控制台输入图像文件名 2） 对输入图像进行 N*N 分块,然后对每块作DCT变换，将变换 结果以BMP格式文件输出.要求通过控制台输入 N的具体值. 3） 观察步骤 2）输出的结果图像 4） 将步骤2)得到的变换结果进行 N*N DCT反变换，将反变换结果以BMP格式文件输出。注意，此处的 N一定与步骤 2）的 N 取值一致。 5）观察对比原始图像与步骤 4）结果图像的异同。 实验环境: Windows7、VC6.0 算法描述: 程序流程图 测试程序说明 测试程序为编写在主函数main()里的测试程序。运行进入控制台选择要进行DCT变换和DCT反变换的BMP图像，再输入N值。程序跳过头文件信息，包括位图文件头，位图信息头，调色板等，申请颜色表的空间，读颜色表进入内存，申请图像数据空间，变换后写入文件，最后进行图像对比。 源程序清单: 添加必要的注释 #include windows.h #include stdio.h #include iostream.h #include fstream.h #include math.h #define WIDTHBYTES(i) ((i+31)/32*4) #define PI 3.141592 #define SQRTMN sqrt(N*N) BITMAPFILEHEADER bf; //位图文件头 BITMAPINFOHEADER bi; //位图信息头 RGBQUAD *pColorTable; //调色板信息 unsigned char **ImageData; //原始图像信息 unsigned char **Image; //DCT反变换之后的图像信息 unsigned char **ImageDCT; //DCT变换之后的图像信息 double **DCTData; int N; char filename[20]; int s = 0; long LineBytes; long height; int numberblockHeight, numberblockWidth; double C(int a) { double value; if(a == 0) { value = sqrt(2)/2; } else { value = 1; } return value; } //DCT变换函数 double DCT(int u, int v ,int i, int j, unsigned char **image) { double value = 0, a,b; int m ,n; double para = 2*C(u)*C(v)/SQRTMN; for(m = 0; m N; m++) { for(n = 0; n N; n++) { a = cos(PI*(2*m+1)*u/(2*N)); b = cos(PI*(2*n+1)*v/(2*N)); value += a*b*image[i*N+m][j*N+n]; } } value = value * para; return value; } //DCT反变换函数 unsigned char IDCT(in