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