BGR到CIELab的颜色空间转换

1、是8bit的BGR图片转换到8bit的CIELab图片,非归一化

2、转换公式用整形转换公式,而非浮点公式,这样存在误差


定义转换公式系数表:

static int _coeff[3][3];
static int _labTab[1024];

int bn = 1<<20;
//X,Y,Z的值从(0,1)放到(0,1020)的初始化表
static int _labTab[1024];
_coeff[0][0] = (int)(0.189828f * bn+0.5);
_coeff[0][1] = (int)(0.376219f * bn+0.5);
_coeff[0][2] = (int)(bn - _coeff[0][0] - _coeff[0][1]);//0.433953f * bn+0.5;
_coeff[1][0] = (int)(0.072169f * bn+0.5);
_coeff[1][1] = (int)(0.715160f * bn+0.5);
_coeff[1][2] = (int)(bn - _coeff[1][0] - _coeff[1][1]);//0.212671f * bn+0.5;
_coeff[2][0] = (int)(0.872766f * bn+0.5);
_coeff[2][1] = (int)(0.109477f * bn+0.5);
_coeff[2][2] = (int)(bn - _coeff[2][0] - _coeff[2][1]);

float c = 16.f/116.f;
for (int t=0; t<1024; ++t)
{
    if(t>9)
	_labTab[t] = (int)(powf(t/1020.f,1.f/3)*1020.f+0.5);
    else
	_labTab[t] = (int)((7.787*t/1020 + c)*1020 + 0.5);
}


转换:

	int rows = src.rows;
	int cols = src.cols;
	uchar* bgrData = src.ptr(0); //src为原始图片,cv::Mat类型
	for(int y=0; y<rows; ++y)
	{
		for(int x=0; x<cols; ++x)
		{
			uchar* BGR = bgrData + y*cols*3 + x*3;
			uchar* Lab = labData + y*cols*3 + x*3; //labData为转换的LAB数组 

			int R = BGR[2];
			int G = BGR[1];
			int B = BGR[0];

			//X,Y,Z==>(0,255)==>(0,1020)
			int X = (B * _coeff[0][0] + G * _coeff[0][1] + R * _coeff[0][2] + 524288)>>(20-2);
			int Y = (B * _coeff[1][0] + G * _coeff[1][1] + R * _coeff[1][2] + 524288)>>(20-2);
			int Z = (B * _coeff[2][0] + G * _coeff[2][1] + R * _coeff[2][2] + 524288)>>(20-2);
			int FX = _labTab[X];
			int FY = _labTab[Y];
			int FZ = _labTab[Z];

			//放大1024(1020)倍到(0,100)再放大到(0,255)
			uchar L = (296*FY - 41780)>>10;
			uchar a = ((500*(FX - FY) + 512)>>10) + 128;
			uchar b = ((200*(FY - FZ) + 512)>>10) + 128;

			Lab[0] = L;
			Lab[1] = a;
			Lab[2] = b;
		}
	}




  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值