数据压缩实验1|彩色空间转换

最新推荐文章于 2021-11-21 20:43:03 发布

mengjiuye

最新推荐文章于 2021-11-21 20:43:03 发布

阅读量129

点赞数

本文链接：https://blog.csdn.net/mengjiuye/article/details/115190729

版权

数据压缩2|彩色空间转换

实验目的
实验内容
实验步骤
- rgb转yuv
- yuv转rgb
实验总结

实验目的

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的一半。

实验步骤

rgb转yuv

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

RGB文件存储方式：按每个像素BGR分量依次存放，BGR BGR ……BGR
YUV文件存储方式：按照全部像素的Y数据块、U数据块和V数据块依次存放。先Y，后V，中间是U。

利用转化公式将R，G，B的每个点的值运算后存到对应的Y，U，V点中，转化公式如下：
Y＝0.2990R+0.5870G+0.1140B
U＝-0.1684R-0.3316G+0.5B
V＝0.5R-0.4187G-0.0813B

在工程-设置-调试中设置程序文件所在的工作目录和主函数所需的命令参数，若之后需要修改图片的长宽，修改会更为便捷。
在这里插入图片描述

主要程序：
查表法，靠存储量换取速度，存储事先计算好的结果。

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

分配内存空间

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;

将rgb转换成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 ++;
	}
}

下采样

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

main函数

#include <stdio.h>
#include <stdlib.h>
#include <malloc.h>
#include "rgb2yuv.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;

	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;

	rgbFileName = argv[1];
	yuvFileName = argv[2];

	frameWidth = atoi(argv[3]);
	frameHeight = atoi(argv[4]);
	
	rgbFile = fopen(rgbFileName, "rb");
	if (rgbFile == NULL)
	{
		printf("cannot find rgb file\n");
		exit(1);
	}
	else
	{
		printf("The input rgb file is %s\n", rgbFileName);
	}

	yuvFile = fopen(yuvFileName, "wb");
	if (yuvFile == NULL)
	{
		printf("cannot find yuv file\n");
		exit(1);
	}
	else
	{
		printf("The output yuv file is %s\n", yuvFileName);
	}

	rgbBuf = (u_int8_t*)malloc(frameWidth * frameHeight * 3);

	yBuf = (u_int8_t*)malloc(frameWidth * frameHeight);
	uBuf = (u_int8_t*)malloc((frameWidth * frameHeight) / 4);
	vBuf = (u_int8_t*)malloc((frameWidth * frameHeight) / 4);

	if (rgbBuf == NULL || yBuf == NULL || uBuf == NULL || vBuf == NULL)
	{
		printf("no enought memory\n");
		exit(1);
	}

	while (fread(rgbBuf, 1, frameWidth * frameHeight * 3, rgbFile)) 
	{
		if(RGB2YUV(frameWidth, frameHeight, rgbBuf, yBuf, uBuf, vBuf, flip))
		{
			printf("error");
			return 0;
					}

		for (i = 0; i < frameWidth*frameHeight; i++)
		{
			if (yBuf[i] < 16) yBuf[i] = 16;
			if (yBuf[i] > 235) yBuf[i] = 235;
		}

		for (i = 0; i < frameWidth*frameHeight/4; i++)
		{
			if (uBuf[i] < 16) uBuf[i] = 16;
			if (uBuf[i] > 240) uBuf[i] = 240;

			if (vBuf[i] < 16) vBuf[i] = 16;
			if (vBuf[i] > 240) vBuf[i] = 240;
		}

		fwrite(yBuf, 1, frameWidth * frameHeight, yuvFile);
		fwrite(uBuf, 1, (frameWidth * frameHeight) / 4, yuvFile);
		fwrite(vBuf, 1, (frameWidth * frameHeight) / 4, yuvFile);

		printf("\r...%d", ++videoFramesWritten);
	}

	printf("\n%u %ux%u video frames written\n", 
		videoFramesWritten, frameWidth, frameHeight);

	/* cleanup */

	fclose(rgbFile);
	fclose(yuvFile);

	return(0);
}

实验结果：
在这里插入图片描述
原图与转换后的图像基本相同

yuv转rgb

将U,V进行扩展，变为原来的4倍。利用转化公式将Y，U，V的每个点的值运算后存到对应的R，G，B点中，转化公式如下：
R= Y +1.4075V ;
G= Y-0.3455U- 0.7169V;
B= Y+1.7790U ;
主要代码：
将U,V进行扩展：

for (j = 0; j < yy/2; j ++)
{   psu=u_buffer +j*x/2;
    psv=v_buffer +j*x/2;
	pu1 = exp_u_buf + 2 * j * x;
	pu2 = exp_u_buf + (2 * j + 1) * x;
	pv1 = exp_v_buf + 2 * j * x;
	pv2 = exp_v_buf + (2 * j + 1) * x;
	for (i = 0; i < x/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转换成rgb：

if (flip) {
	for (j = 0; j < yy; j ++)
	{   y = y_buffer + (yy- j - 1) * x;
		u = exp_u_buf + (yy - j - 1) * x;
		v = exp_v_buf + (yy - j - 1) * x;

		for (i = 0; i < x; i ++) {
			g = b + 1;
			r = b + 2;

			
			float r0, g0, b0;
			r0 = ( (*y)+RGBYUV14075[*v]); 
			g0 = ((*y)-RGBYUV03455[*u]- RGBYUV07169[*v]);
			b0 = (  (*y)+RGBYUV1779[*u] );
			*r = (r0 > 0 ? (r0 > 255 ? 255 : (unsigned char)r0) : 0);
			*g = (g0 > 0 ? (g0 > 255 ? 255 : (unsigned char)g0) : 0);
			*b = (b0 > 0 ? (b0 > 255 ? 255 : (unsigned char)b0) : 0);
			b += 3;
			y ++;
			u ++;
			v ++;
		}
	}
} else {
	for (i = 0; i < size; i++)
	{
		u= v + 1;
		y = v + 2;
		*r = (unsigned char)( (*y)+RGBYUV14075[(*v)-128]); 
		*g = (unsigned char)((*y)-RGBYUV03455[(*u)-128]- RGBYUV07169[(*v)-128]);
		*b = (unsigned char)(  (*y)+RGBYUV1779[(*u)-128] );
		v += 3;
		 r++;
		g ++;
		b ++;
	}
}
free(exp_u_buf );
free(exp_v_buf);


return 0;
}

实验结果：
在这里插入图片描述

实验总结

注意数据格式的转换。
可以采用查表法靠存储量换取速度。

mengjiuye

关注

0
点赞
踩
0

收藏

觉得还不错? 一键收藏
0
评论
数据压缩实验1|彩色空间转换

数据压缩2|彩色空间转换实验目的实验内容实验步骤实验目的1．学会从计算和程序的角度分析问题通过完成本实验，理解计算思维，即从问题出发，通过逐步分析和分解，把原问题转化为可用程序方式解决的问题。在此过程中设计出一个解决方案。2．进一步理解彩色空间的概念并掌握不同彩色空间转换的基本方程。3．通过逐步设计程序，掌握编程细节：如查找表的设计，内存分配，对 U 和 V 信号进行下采样，文件读写过程等。掌握程序调试的基本方法。实验内容1．掌握彩色空间转换的基本思想及转换公式（1）YUV与RGB空间的相互
复制链接

扫一扫