用c++实现颜色空间rgb，grey，luv和lab的互转_luv与lab色度公式之间的转换-CSDN博客

1 rgb转grey，rgb转luv，rgb转lab

1. 1 rgb转grey

void RgbToGrey(unsigned char *rgb, double *grey)
{
double R = ((double) rgb[0]) / (double)255.0;
double G = ((double) rgb[1]) / (double)255.0;
double B = ((double) rgb[2]) / (double)255.0;
    double y = (double)0.212671 * R + (double)0.715160 * G + (double)0.072169 * B;


if (y > (double)0.008856) grey[0] = (double)116.0 * pow(y, (double)1.0 / (double)3.0) - (double)16.0;
    else grey[0] = (double)903.3 * y;
}

1. 2 rgb转luv

void RgbToLuv(unsigned char *rgb, double *luv)
{
double L, u, v;


double R = ((double) rgb[0]) / (double)255.0;
double G = ((double) rgb[1]) / (double)255.0;
double B = ((double) rgb[2]) / (double)255.0;


    double x = (double)0.412453 * R + (double)0.357580 * G + (double)0.180423 * B;
    double y = (double)0.212671 * R + (double)0.715160 * G + (double)0.072169 * B;
    double z = (double)0.019334 * R + (double)0.119193 * G + (double)0.950227 * B;


if (y > 0.008856) L = (double)116.0 * pow(y, (double)1.0 / (double)3.0) - (double)16.0;
    else L = (double)903.3 * y;


double sum = x + 15 * y + 3 * z;


if(sum != 0) u = 4 * x / sum, v = 9 * y / sum;
else u = 4.0, v = (double)9.0 / (double)15.0;


luv[0] = L;
    luv[1] = 13 * L * (u - (double)0.19784977571475);
    luv[2] = 13 * L * (v - (double)0.46834507665248);
}

1. 3 rgb转grey

void RgbToLab(unsigned char *rgb, double *lab)
{
double R = ((double) rgb[0]) / (double)255.0;
double G = ((double) rgb[1]) / (double)255.0;
double B = ((double) rgb[2]) / (double)255.0;


    double X =  0.412453 * R + 0.357580 * G + 0.180423 * B;
    double Y =  0.212671 * R + 0.715160 * G + 0.072169 * B;
    double Z =  0.019334 * R + 0.119193 * G + 0.950227 * B;


    double xr = X / 0.950456, yr = Y / 1.000, zr = Z / 1.088854;
    if(yr > 0.008856) lab[0] = 116.0 * pow(yr, 1.0 / 3.0) - 16.0;
else lab[0] = 903.3 * yr;


double fxr, fyr, fzr;


    if(xr > 0.008856) fxr = pow(xr, 1.0 / 3.0);
else fxr = 7.787 * xr + 16.0 / 116.0;


if(yr > 0.008856) fyr = pow(yr, 1.0 / 3.0);
else fyr = 7.787 * yr + 16.0 / 116.0;


if(zr > 0.008856) fzr = pow(zr, 1.0 / 3.0);
else fzr = 7.787 * zr + 16.0 / 116.0;


lab[1] = 500.0 * (fxr - fyr);
lab[2] = 200.0 * (fyr - fzr);
}

2 grey转rgb，grey转luv，grey转lab

2.1 grey转rgb

void GreyToRgb(double *grey, unsigned char *rgb)
{
if(grey[0] < (double)0.1) rgb[0] = rgb[1] = rgb[2] = 0;
else
{
double y;
if(grey[0] <= (double)7.9996) y = 255 * grey[0] / (double)903.3;
else y = (grey[0] + (double)16.0) / (double)116.0, y  = 255 * y * y * y;


rgb[0] = (unsigned char)((y < 0) ? 0 : ((y > 255) ? 255 : y));
rgb[2] = rgb[1] = rgb[0];
}
}

2.2 grey转luv

void GreyToLuv(double *grey, double *luv)
{
luv[0] = grey[0], luv[1] = luv[2] = 0;
}

2.3 grey转lab

void GreyToLab(double *grey, double *lab)
{
lab[0] = grey[0], lab[1] = lab[2] = 0;
}

3 luv转rgb，luv转grey

3.1 luv转rgb

void LuvToRgb(double *luv, unsigned char *rgb)
{
if(luv[0] < 0.1) rgb[0] = rgb[1] = rgb[2] = 0;
else
{
double x, y, z, u, v;


if(luv[0] <= (double)7.9996) y = luv[0] / (double)903.3;
else y = (luv[0] + (double)16.0) / (double)116.0, y  = y * y * y;


u = luv[1] / (13 * luv[0]) + (double)0.19784977571475;
v = luv[2] / (13 * luv[0]) + (double)0.46834507665248;
x = 9 * u * y / (4 * v);
z = (12 - 3 * u - 20 * v) * y / (4 * v);


double R = (double)3.240479 * x - (double)1.537150 * y - (double)0.498535 * z;
double G = (double)-0.969256 * x + (double)1.875992 * y + (double)0.041556 * z;
double B = (double)0.055648 * x - (double)0.204043 * y + (double)1.057311 * z;
R *= 255, G *= 255, B *= 255;


rgb[0] = (unsigned char)((R < 0) ? 0 : ((R > 255) ? 255 : R));
rgb[1] = (unsigned char)((G < 0) ? 0 : ((G > 255) ? 255 : G));
rgb[2] = (unsigned char)((B < 0) ? 0 : ((B > 255) ? 255 : B));
}
}

3.2 luv转grey

void LuvToGrey(double *luv, double *grey)
{
grey[0] = luv[0];
}

4 lab转rgb，lab转grey

4.1 lab转rgb

void LabToRgb(double *lab, unsigned char *rgb)
{
    double X, Y, Z;
    double P = (lab[0] + 16.0) / 116.0;


if(lab[0] > 7.9996) Y = 1.000 * P * P * P;
else Y = 1.000 * lab[0] / 903.3;


    double yr = Y / 1.000, fy;
if(yr > 0.008856) fy = pow(yr, 1.0 / 3.0);
else fy = 7.787 * yr + 16.0 / 116.0;


double fx = lab[1] / 500.0 + fy, fz = fy - lab[2] / 200.0;


if(fx > 0.2069) X = 0.950456 * fx * fx * fx;
else X = 0.950456 / 7.787 * (fx - 16.0 / 116.0);


if(fz > 0.2069) Z = 1.088854 * fz * fz * fz;
else Z = 1.088854 / 7.787 * (fz - 16.0 / 116.0);


    double R = 3.240479 * X - 1.537150 * Y - 0.498535 * Z;
    double G = -0.969256 * X + 1.875992 * Y + 0.041556 * Z;
    double B = 0.055648 * X - 0.204043 * Y + 1.057311 * Z;
R *= 255, G *= 255, B *= 255;


rgb[0] = (unsigned char)((R < 0) ? 0 : ((R > 255) ? 255 : R));
rgb[1] = (unsigned char)((G < 0) ? 0 : ((G > 255) ? 255 : G));
rgb[2] = (unsigned char)((B < 0) ? 0 : ((B > 255) ? 255 : B));
}

4.2 lab转grey

void LabToGrey(double *lab, double *grey)
{
 grey[0] = lab[0];
}

5 参考文献

[1]http://blog.sina.com.cn/s/blog_5309cefc0101cdeg.html 作者：RescuEe 题目：用c++实现颜色空间rgb，grey，luv和lab的互转

[2]http://blog.sina.com.cn/s/blog_46e3af5b0100u7wt.html 作者：孙志海题目：RGB与Lab颜色空间转换的程序实现及效果(VC++6.0)