实验：彩色空间转换

实验目的

1. 学会从计算和程序的角度分析问题
   通过完成本实验，理解计算思维，即从问题出发，通过逐步分析和分解，把原问题转化为可用程序方式解决的问题。在此过程中设计出一个解决方案。
2. 进一步理解彩色空间的概念并掌握不同彩色空间转换的基本方程。
3. 通过逐步设计程序，掌握编程细节：如查找表的设计，内存分配，对 U 和 V 信号进行下采样，文件读写过程等。掌握程序调试的基本方法。

实验内容

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的一半。
2. 掌握由问题到程序的实现过程
   （1）掌握所用的程序语言
   只有熟悉所使用的程序语言，精通语言中各种结构（包括其形式和意义），才可能实现将所要解决的问题转换为程序。。
   （2）学会写程序
   首先确定适用的程序结构，保证写出的程序结构良好清晰、易于阅读和理解。还应注意当有些条件或要求改变时，程序是否容易修改去满足新的要求。
   （3）检查程序错误的能力
   程序错误通常包括语言错误、逻辑错误和算法错误。要熟悉所用工具和编程环境帮助迅速找到错误并改正。

实验步骤

基本要求（必做）：

编写RGB转化为YUV程序，重点掌握函数定义，部分查找表的初始化和调用，缓冲区分配。将得到的RGB文件转换为YUV文件，用YUV Viewer播放器观看，验证是否正确。

核心代码

#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;
	}

	// 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_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;

	// convert RGB to YUV
	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 ++;
			}
		}
	} else {
		for (i = 0; i < size; 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 ++;
		}
	}

	// subsample UV
	for (j = 0; j < y_dim/2; j ++)
	{
		psu = sub_u_buf + j * x_dim / 2;
		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 ++)
		{
			*psu = (*pu1 + *(pu1+1) + *pu2 + *(pu2+1)) / 4;
			*psv = (*pv1 + *(pv1+1) + *pv2 + *(pv2+1)) / 4;
			psu ++;
			psv ++;
			pu1 += 2;
			pu2 += 2;
			pv1 += 2;
			pv2 += 2;
		}
	}

	free(u_buffer);
	free(v_buffer);

	return 0;
}


void InitLookupTable()
{
	int i;

	for (i = 0; i < 256; i++) RGBYUV02990[i] = (float)0.2990 * i;
	for (i = 0; i < 256; i++) RGBYUV05870[i] = (float)0.5870 * i;
	for (i = 0; i < 256; i++) RGBYUV01140[i] = (float)0.1140 * i;
	for (i = 0; i < 256; i++) RGBYUV01684[i] = (float)0.1684 * i;
	for (i = 0; i < 256; i++) RGBYUV03316[i] = (float)0.3316 * i;
	for (i = 0; i < 256; i++) RGBYUV04187[i] = (float)0.4187 * i;
	for (i = 0; i < 256; i++) RGBYUV00813[i] = (float)0.0813 * i;
}

验证

用YUVviewerPlus验证原RGB图和生成的YUV图

提高要求（可选）：

编写将YUV转换为RGB的程序。将给定的实验数据用该程序转换为RGB文件。并与原RGB文件进行比较，如果有误差，分析误差来自何处。

核心代码

#include <iostream>
#include "stdlib.h"
#include "yuv2rgb.h"

static float YUVRGB14075[256];
static float YUVRGB07169[256], YUVRGB03455[256];
static float YUVRGB17790[256];

int YUV2RGB(void* y_in,void* u_in,void* v_in,void* rgb_out, int W, int H)
{
	static int init_done = 0;
	long i, j;
	float rt, gt, bt; 
	unsigned char *r, *g, *b;
	unsigned char *y, *u, *v;
	b = (unsigned char *)rgb_out;
	y = (unsigned char *)y_in;
	u = (unsigned char *)u_in;
	v = (unsigned char *)v_in;

	if (init_done == 0)
	{
		InitLookupTable();
		init_done = 1;
	}

	for(j = 0;j < H; j++)
	{
		for (i = 0;i < W; i++)
		{
			g = b+1;
			r = b+2;

			bt = (*y) + YUVRGB17790[*u] - YUVRGB17790[128];
			gt = (*y) - YUVRGB07169[*v] - YUVRGB03455[*u] + YUVRGB03455[128] + YUVRGB07169[128];
			rt = (*y) + YUVRGB14075[*v] - YUVRGB14075[128];

			//bt = (*y) + YUVRGB17790[*u];
			//gt = (*y) - YUVRGB07169[*v] - YUVRGB03455[*u];
			//rt = (*y) + YUVRGB14075[*v];

			if (rt < 0)    rt=0;
			if (rt > 255)  rt=255;
			if (gt < 0)    gt=0;
			if (gt > 255)  gt=255;
			if (bt < 0)    bt=0;
			if (bt > 255)  bt=255;

			*b = (unsigned char)bt;
			*g = (unsigned char)gt;
			*r = (unsigned char)rt;

			b += 3;
			y++;
			if(i % 2 == 0 && j % 2 == 0)
			{
				u++;
				v++;
			}		
		}
	}
	return 0;
}

void InitLookupTable()
{
	int i;
	for (i = 0 ; i < 256 ; i++)		YUVRGB14075[i] = (float)1.4075 * i;
	for (i = 0 ; i < 256 ; i++)   	YUVRGB07169[i] = (float)0.7169 * i;
	for (i = 0 ; i < 256 ; i++)  	YUVRGB03455[i] = (float)0.3455 * i;
	for (i = 0 ; i < 256 ; i++)  	YUVRGB17790[i] = (float)1.7790 * i;
}

验证

用YUVviewerPlus查看三个图，分别为原YUV、实验二RGB2YUV的结果、实验一YUV2RGB的结果

放大后发现细节部分有像素色块，才像是因为U、V分量的上取样产生了误差。

实验总结

• 实验中用部分查找表提前存储了YUV与RGB转化公式中需要用到的大量计算数据，在函数中可以直接调用查找表，避免实施运算。
• 学会用YUVviewerPlus查看RGB图像。
• 深入理解RGB\YUV文件的存储结构。