数据压缩实验一：彩色空间转换

1．彩色空间转换的基本思想及转换公式

（1）YUV与RGB空间的相互转换
由电视原理可知，亮度和色差信号的构成如下：
Y＝0.2990R+0.5870G+0.1140B
R-Y＝0.7010R-0.5870G-0.1140B
B-Y＝-0.2990R-0.5870G+0.8860B
为了使色差信号的动态范围控制在0.5之间，需要进行归一化，对色差信号引入压缩系数。归一化后的色差信号为：
U＝-0.1684R-0.3316G+0.5B
V＝0.5R-0.4187G-0.0813B
（2） 码电平分配及数字表达式
⚫ 亮电平信号量化后码电平分配
在对分量信号进行8比特均匀量化时，共分为256个等间隔的量化级。为了防止信号变动造成过载，在256级上端留20级，下端留16级作为信号超越动态范围的保护带。
⚫ 色差信号量化后码电平分配
色差信号经过归一化处理后，动态范围为-0.5－0.5，让色差零电平对应码电平128，
色差信号总共占225个量化级。在256级上端留15级，下端留16级作为信号超越动态范围的保护带。
（3）色度格式
4:2:0格式是指色差信号U，V的取样频率为亮度信号取样频率的四分之一，在水平方向和垂直方向上的取样点数均为Y的一半。
（4）yuv和rgb存储方式
YUV
在内存中yuv（4:2:0）格式的文件，以Y（0帧）U（0帧）V（0帧）、Y（1帧）U（1帧）V（1帧）……的方式存储在连续的内存中，其中U、V的数据量是Y的1/4
一帧的数据量 = 水平像素数 * 垂直像素数 * 量化比特数 * 1.5
RGB
和yuv以帧为单位存储不同，内存中rgb是以像素为单位存储的，每个像素都对应了1字节的R、1字节的G、1字节的B，一个像素用3个字节表示
一帧的数据量 = 水平像素数 * 垂直像素数 * 量化比特数 * 3

二、实验流程

rgb2yuv
1.程序初始化（打开两个文件、定义变量和缓冲区等）
2.读取RGB文件，抽取RGB数据写入缓冲区
3.调用RGB2YUV的函数实现RGB到YUV数据的转换
具体过程：
1）从主函数中获得输入的rgb文件指针
2）定义相关中间变量及指针（开辟缓冲区）
3）色差信号下采样
4）数据转换算法
5）释放缓冲区
4.写RGB文件
4.写YUV文件
5.程序收尾工作（关闭文件，释放缓冲区）
yuv2rgb
1.程序初始化（打开两个文件、定义变量和缓冲区等）
2.读取YUV文件，抽取YUV数据写入缓冲区
3.调用YUV2RGB的函数实现YUV到RGB数据的转换
具体过程：
1）从主函数中获得输入的yuv文件指针
2）定义相关中间变量及指针（开辟缓冲区）
3）u、v上采样
4）数据转换算法
5）释放缓冲区
4.写RGB文件
5.程序收尾工作（关闭文件，释放缓冲区）

三、实验

rgb2yuv

rgb2yuv.h

#ifndef YUV2RGB_H_ 
#define YUV2RGB_H_ 
int RGB2YUV (int x_dim, int y_dim, void *bmp, void *y_out, void *u_out, void *v_out, int flip);
void InitLookupTable();
#endif

rgb2yuv.cpp

#include "stdlib.h"
#include "rgb2yuv.h"

static float RGBYUV02990[256], RGBYUV05870[256], RGBYUV01140[256];
static float RGBYUV01684[256], RGBYUV03316[256];
static float RGBYUV04187[256], RGBYUV00813[256];
int RGB2YUV (int x_dim, int y_dim, void *bmp, void *y_out, void *u_out, void *v_out, 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;
	}

	if ((x_dim % 2) || (y_dim % 2)) return 1;
	size = x_dim * y_dim;

	y_buffer = (unsigned char *)y_out;
	sub_u_buf = (unsigned char *)u_out;
	sub_v_buf = (unsigned char *)v_out;
	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;
	y = y_buffer;
	u = u_buffer;
	v = v_buffer;

	if (!flip) {
   
		for (j = 0; j < y_dim; j ++)
		{
   
			y = y_buffer + (y_dim - j - 1) * x_dim;
			u = u_buffer + (y_dim - j - 1) * x_dim;
			v = v_buffer + (y_dim - j - 1) * x_dim;

			for (i = 0; i < x_dim; i ++) 
			{
   
				g = b + 1;
				r = b + 2;
				*y = (unsigned char)(  RGBYUV02990[*r] + RGBYUV05870[*g] + RGBYUV01140[*b]);
				*u = (unsigned char)(- RGBYUV01684[*r] - RGBYUV03316[*g] + (*b)/2          + 128);
				*v = (unsigned char)(  (*r)/2          - RGBYUV04187[*g] - RGBYUV00813[*b] + 128);
				b += 3;
				y ++;
				u ++;
				v ++;
			}
		}