RGB转LAB色彩空间
1.原理
RGB无法直接转换成LAB，需要先转换成XYZ再转换成LAB，即：RGB——XYZ——LAB
因此转换公式分两部分：
（1）RGB转XYZ
假设r,g,b为像素三个通道，取值范围均为[0,255]，转换公式如下：
（1）
（2）
（3）
M=
0.4124,0.3576,0.1805
0.2126,0.7152,0.0722
0.0193,0.1192,0.9505
等同于如下公式：
X = var_R * 0.4124 + var_G * 0.3576 + var_B * 0.1805
Y = var_R * 0.2126 + var_G * 0.7152 + var_B * 0.0722
Z = var_R * 0.0193 + var_G * 0.1192 + var_B * 0.9505
上面的gamma函数，是用来对图象进行非线性色调编辑的，目的是提高图像对比度。这个函数不是唯一的，但是我在网上查到的基本都使用上式。
（2）XYZ转LAB
（5）
（6）
上面两个公式中，L*，a*，b*是最终的LAB色彩空间三个通道的值。X，Y，Z是RGB转XYZ后计算出来的值，Xn，Yn，Zn一般默认是95.047，100.0，108.883。
2.代码实现
（1）完全按照算法无优化实现（经@竹子小蜻蜓指正，原来代码有误，现已改正--206.04.12）
inline float gamma(float x) {return x>0.04045?pow((x+0.055f)/1.055f,2.4f):x/12.92;}; void RGBToLab(unsigned char*rgbImg,float*labImg) { float B=gamma(rgbImg[0]/255.0f); float G=gamma(rgbImg[1]/255.0f); float R=gamma(rgbImg[2]/255.0f); float X=0.412453*R+0.357580*G+0.180423*B; float Y=0.212671*R+0.715160*G+0.072169*B; float Z=0.019334*R+0.119193*G+0.950227*B;
float X/=0.95047;
float Y/=1.0;
float Z/=1.08883;
float FX = X > 0.008856f ? pow(X,1.0f/3.0f) : (7.787f * X +0.137931f); float FY = Y > 0.008856f ? pow(Y,1.0f/3.0f) : (7.787f * Y +0.137931f); float FZ = Z > 0.008856f ? pow(Z,1.0f/3.0f) : (7.787f * Z +0.137931f); labImg[0] = Y > 0.008856f ? (116.0f * FY - 16.0f) : (903.3f * Y); labImg[1] = 500.f * (FX - FY); labImg[2] = 200.f * (FY - FZ); }
上面是完全按照转换算法做的无优化的实现，里面涉及了大量的浮点运算，在PC上可能没什么问题，但是如果是在Android操作系统的移动端上，即使利用JNI，把转换算法写成C++版本进行加速，速度也不理想，因为这个操作时逐像素的，每个像素做几十次浮点运算，耗时还是十分巨大的。
（2）牺牲一些精度的快速实现
首先可以把gamma函数去掉，因为经过我的测试，这个函数带来的影响很小。
其次，M矩阵中的值都为小数，可以通过放大变成整数，接下来的计算就可以变成整数计算（在手机上浮点计算一般耗时是整形计算的几倍）。在这里，可以把M矩阵的值都扩大2^20倍变为整数，在计算得到X,Y,Z分量时再缩小回来。但是，XYZ的小数部分会消失，会引入误差，而在公式(5)（6）中（尤其是（6），涉及到指数运算和阈值运算，之前XYZ的误差在这里会被放大），这种误差还会继续传递给最终的L,a,b值。那怎么才能保证一定的准确率呢？
在这里我的方法是引入了一个查表机制。过程如下：
A.在计算XYZ的时候，参数扩大2^20倍，但最后缩小时只缩小2^18，这样计算出来的XYZ值范围是[0,1023]
B.建立一个table，用于将取值为[0,1023]的XYZ通过f(t)函数映射到中间结果，记为：LabTable(m)
C.将LabTable(m)代入公式(5)计算最终的Lab分量值
代码如下：
const static int big_shift=18; const static int HalfShiftValue=512; const static int shift=10; const static int offset=128<<shift; const static int ScaleLC = (16 * (1 << shift)); const static int ScaleLT = 116; const static int ScaleY= 903; const static int para1=500; const static int para2=200; const static int ThPara=9; void RGBToLab(unsigned char*rgbImg,int*labImg) { long long X=(rgbImg[0] * 199049 + rgbImg[1] * 394494 + rgbImg[2] * 455033 + 524288)>> (big_shift); long long Y=(rgbImg[0] * 75675 + rgbImg[1] * 749900 + rgbImg[2] * 223002 + 524288) >> (big_shift); long long Z=(rgbImg[0] * 915161 + rgbImg[1] * 114795 + rgbImg[2] * 18621 + 524288) >> (big_shift); labImg[0] = Y>ThPara?((ScaleLT * LabTable[Y] - ScaleLC + HalfShiftValue)>>shift):((ScaleY* Y)>>shift); labImg[1] = (para1*(LabTable[X] - LabTable[Y])+HalfShiftValue+offset)>>shift; labImg[2] = (para2*(LabTable[Y] - LabTable[Z])+HalfShiftValue+offset)>>shift; }
建立LabTable的代码如下：
int LabTable[1024]; for (int I = 0; I < 1024; I++) { if (I > 9) LabTable[I] = (int)(pow((float)I / 1020, 1.0F / 3) * (1 << 10) + 0.5 ); else LabTable[I] = (int)((29 * 29.0 * I / (6 * 6 * 3 * 1020) + 4.0 / 29) * (1 << 10) + 0.5 ); printf("%d,",LabTable[I]); }
参考文章：
http://www.easyrgb.com/index.php?X=MATH&H=07#text7
http://vcsos.com/Article/pageSource/120314/20120314114422.shtml
http://www.cnblogs.com/Imageshop/archive/2013/01/31/2888097.html
http://www.cnblogs.com/Imageshop/archive/2013/02/02/2889897.html
