数据压缩实验一实验报告

1.实验原理

YUV到RGB的转换算法为 R=Y+1.4075(V−128) G=Y+0.3455(U−128)−0.7169(V−128) B=Y+1.779(U−128) 分析 此时的U、V实际上是Cb和Cr，为把(R-Y)和(B-Y)的动态范围控制在±350mV内，引入压缩系数，可以得到公式 Cr=0.713（R-Y） Cb=0.564（B-Y） 又为了便于数字化处理，还要将Cr、Cb引入350mV的偏置，又由于经过8比特量化，色差信号经过归一化处理后，动态范围为-0.5－0.5，让色差零电平对应码电平128。所以公式变为 V=0.713（R-Y）+128 U=0.564（B-Y）+128 经过方程式的计算，可以得到最后的公式

RGB和YUV彩色空间基本原理

RGB彩色空间存储方式是每一个像素点都按照B、G、R的顺序进行存储，每一帧占用空间为width * height * 3个字节的空间。YUV彩色空间存储方式为每一帧的YUV分开存储，每一帧Y占用的空间为width * height，每一帧U和V所占用的空间分别为width * height /4。

取样格式

读取的图片为YUV空间，取样格式为4：2：0，即V和U的水平和垂直采样点数均为Y的一半。但要生成的RGB文件的取样格式为4：4：4，这就要求我们对给定的Y保持不变，对给定的U和V进行上采样，使其可与RGB的采样空间相同。 根据将RGB转成YUV的实验中下采样的方法，我采用的方法是用指针指向yuv文件中的u、v，并开辟两个128*128的空间存*储文件中的u、v。同时还开辟了两个256*256的空间，将一个指向给定uv指针的内容赋值给新的四个像素块，详见具体代码。

2.实验的流程分析

程序初始化（打开两个文件、定义变量和缓冲区 等） 读取YUV文件，抽取YUV数据写入缓冲区 调用YUV2RGB的函数实现YUV到RGB数据的转 换 写RGB文件 程序收尾工作（关闭文件，释放缓冲区）

命令行参数的设置 project->properties->Debugging-> Working Directory&Command Arguments

3.重要代码及分析

重要代码通过高亮标出

yuv转为rgb的函数代码：yuv2rgb.cpp

#include "stdlib.h" #include "yuv2rgb.h" static float RGBYUV14075[256], RGBYUV03455[256]; static float RGBYUV07169[256], RGBYUV17790[256]; int YUV2RGB(int x_dim, int y_dim, void *bmp, void *y_in, void *u_in, void *v_in, int flip) { static int init_done = 0; long i, j, size; unsigned char *r, *g, *b; unsigned char *y, *u, *v; unsigned char *pu1, *pu2, *pv1, *pv2, *psu, *psv; unsigned char *y_buffer, *u_buffer, *v_buffer; unsigned char *sub_u_buf, *sub_v_buf; if (init_done == 0) { InitLookupTable(); init_done = 1; } // check to see if x_dim and y_dim are divisible by 2 if ((x_dim % 2) || (y_dim % 2)) return 1; size = x_dim * y_dim; // allocate memory y_buffer = (unsigned char *)y_in; sub_u_buf = (unsigned char *)u_in; sub_v_buf = (unsigned char *)v_in; u_buffer = (unsigned char *)malloc(size * sizeof(unsigned char)); v_buffer = (unsigned char *)malloc(size * sizeof(unsigned char)); if (!(u_buffer && v_buffer)) { if (u_buffer) free(u_buffer); if (v_buffer) free(v_buffer); return 2; } b = (unsigned char *)bmp; g = b + 1; r = b + 2; y = y_buffer; u = u_buffer; v = v_buffer; for (j = 0; j < y_dim / 2; j++) { psu = sub_u_buf + j * x_dim / 2;//128*128 psv = sub_v_buf + j * x_dim / 2; pu1 = u_buffer + 2 * j * x_dim; pu2 = u_buffer + (2 * j + 1) * x_dim; pv1 = v_buffer + 2 * j * x_dim; pv2 = v_buffer + (2 * j + 1) * x_dim; for (i = 0; i < x_dim / 2; i++) { *pu1 = *psu; *(pu1 + 1) = *psu; *pu2 = *psu; *(pu2 + 1) = *psu; *pv1 = *psv; *(pv1 + 1) = *psv; *pv2 = *psv; *(pv2 + 1) = *psv; psu++; psv++; pu1 += 2; pu2 += 2; pv1 += 2; pv2 += 2; /*为了将yuv的取样空间由4：2：0变为4：4：4,就要将一个像素块的值赋给四个像素，psu指针指向的是输入的u的空间，将psu指针指向的一个值赋给新开辟空间中的四个像素，比如将psu指向的第一个值赋给第一行的第一、二个值，即pu1和pu1+1及第二行中的第一二个值pu2和pu2+1。然后psv指针指向的是输入的v的空间。*/ } } for (i = 0; i < size; i++) { g = b + 1; r = b + 2; float tmp; /*r的数据类型为unsigned char，范围为0~255，如果直接进行计算则会出现溢出现象，比如计算得出256，则会变为0，因此通过设置浮点型中间变量tmp防止数据的溢出*/ tmp= (*y + RGBYUV14075[*v]); if (tmp > 255) { tmp = 255; } if (tmp < 0) { tmp = 0; } *r = (unsigned char)tmp; tmp = (*y - RGBYUV03455[*u] - RGBYUV07169[*v]); if (tmp > 255) { tmp = 255; } if (tmp < 0) { tmp = 0; } *g = (unsigned char)tmp; tmp = (*y + RGBYUV17790[*u]); if (tmp > 255) { tmp = 255; } if (tmp < 0) { tmp = 0; } *b = (unsigned char) tmp; b += 3; y++; u++; v++; } return 0; } /*采用部分查找表法，提高运行效率*/ void InitLookupTable() { int i; for (i = 0; i < 256; i++) RGBYUV14075[i] = (float)1.0475 * (i-128); for (i = 0; i < 256; i++) RGBYUV03455[i] = (float)0.3455 * (i-128); for (i = 0; i < 256; i++) RGBYUV07169[i] = (float)0.7169 * (i-128); for (i = 0; i < 256; i++) RGBYUV17790[i] = (float)1.779 * (i-128); }

主函数

#include<stdio.h> #include<stdlib.h> #include<malloc.h> #include"yuv2rgb.h" #define u_int8_t unsigned __int8 #define u_int unsigned __int32 #define u_int32_t unsigned __int32 #define FALSE false #define TRUE true int main(int argc,char** argv) { u_int frameWidth = 352; u_int frameHeight = 240; bool flip =TRUE; //unsigned int i; /* internal variables */ char* rgbFileName = NULL; char* yuvFileName = NULL; FILE* rgbFile = NULL; FILE* yuvFile = NULL; u_int8_t* rgbBuf = NULL; u_int8_t* yBuf = NULL; u_int8_t* uBuf = NULL; u_int8_t* vBuf = NULL; u_int32_t videoFramesWritten = 0; /* begin process command line */ /* point to the specified file names */ yuvFileName = argv[1]; rgbFileName = argv[2]; frameWidth = atoi(argv[3]); frameHeight = atoi(argv[4]); /* 打开yuv文件*/ yuvFile = fopen(yuvFileName, "rb"); if (yuvFile == NULL) { printf("cannot find yuv file\n"); exit(1); } else { printf("The input yuv file is %s\n", yuvFileName); } /* open the RAW file */ rgbFile = fopen(rgbFileName, "wb"); if (rgbFile == NULL) { printf("cannot find rgb file\n"); exit(1); } else { printf("The output rgb file is %s\n", rgbFileName); } /* get an input buffer for a frame */ rgbBuf = (u_int8_t*)malloc(frameWidth * frameHeight * 3); /* get the output buffers for a frame */ yBuf = (u_int8_t*)malloc(frameWidth * frameHeight); uBuf = (u_int8_t*)malloc((frameWidth * frameHeight) / 4); vBuf = (u_int8_t*)malloc((frameWidth * frameHeight) / 4); while ((fread(yBuf, 1, frameWidth * frameHeight , yuvFile)) &&(fread(uBuf, 1, frameWidth * frameHeight/4, yuvFile)) &&(fread(vBuf, 1, frameWidth * frameHeight/4, yuvFile))) { if (YUV2RGB(frameWidth, frameHeight, rgbBuf, yBuf, uBuf, vBuf, flip)) { printf("error"); return 0; } fwrite(rgbBuf, 1, frameWidth * frameHeight * 3, rgbFile); printf("\r...%d", ++videoFramesWritten); } printf("\n%u %ux%u video frames written\n", videoFramesWritten, frameWidth, frameHeight); /*释放开辟的内存空间和关闭打开的文件*/ if(rgbBuf!=NULL){free(rgbBuf);} if(yBuf!=NULL){free(yBuf);} if(uBuf!=NULL){free(uBuf);} if(vBuf!=NULL){free(vBuf);} if(yuvFile!=NULL){fclose(yuvFile);} if(rgbFile!=NULL){fclose(rgbFile);} return (0);

}

4.实验结果及分析

错误及调试

1.程序build后没有报错，但在运行过程中出现了breakpoint，经过F11进入main函数，F10进行分步运行后，在main函数中发现错误，释放空间时多次释放rgb空间，而没有释放y、u、v的空间。

2.程序运行后，出现的测试yuv图像所有像素的yuv值均相同，通过检查，发现是yuv指针都始终指向第一个像素点，指针没有移动，后来通过y++，u++，v++的语句移动指针。

3.没有设置tmp为中间变量时出现少量溢出错误，在最终呈现的test。yuv图像上出现蓝色的小色块。以下是对比图。 生成的test图像：

最后结果

2.