色差公式发展的三个重要的阶段：1976年以前（CIELAB和CIELUV的采用）、1976年到2001年（CIEDE2000色差公式的推荐）、2001年以后。
国际照明委员会1998年成立了技术委员会(TC)1-47，主要任务是基于色相和明度的工业色差修正。
经过技术委员会成员的紧密协作，在2001年被推荐了名为CIEDE2000的色差公式。
它包含了5个对CIELAB的校正：明度权重函数、彩度权重函数、色相权重函数、彩度差和色相差的交互项（用来改进蓝色的性能）、调整CIELAB a*的因子（用来改善灰色的性能）。
结果表明，CIEDE2000公式对高级的色差公式如CIE94或CMC都有了很大地改进。
国际照明委员会(CIE，Commission Internationale de L'Eclairage / International Commissionon Illumination)
负责国际颜色标准的研究开发，其中一个非常重要的任务就是建立色差公式，以向用色工业提供物体色判断的标准，比如，用一个色差公式的数值来判断是否合格。自从第一个CIE色度系统研究出来以后，有40多个色差公式已经被研究出来.
主要的色差公式有：FCM色差公式、LABHNU色差公式、JPC79色差公式、ATDN色差公式、住友方法、CIELAB色差公式的改良式、BFD色差公式、LCD色差公式、CMC(1：c)色差公式、CIE94色差公式等。CIEDE2000色差公式。
为了进一步改善工业色差评价的视觉一致性，CIE专门成立了工业色差评价的色相和明度相关修正技术委员会TC1-47（Hue and Lightness Dependent Correction to Industrial Colour Difference Evaluation）,经过该技术委员会对现有色差公式和视觉评价数据的分析与测试，在2000年提出了一个新的色彩评价公式，并于2001年得到了国际照明委员会的推荐，称为CIE2000色差公式，简称CIEDE2000，色差符合为△E00 。CIEDE2000是到目前为止最新的色差公式，该公式与CIE94相比要复杂的多，同时也大大提高了精度。
由于人眼与色差计在色度图中不同位置的色彩评估方法和形状上的差异，引起了许多色彩评定时测量数据与目测结论不匹配的问题，CIE 2000色差计算公式的提出，就是为了解决这个问题和难题。
在CIE Lab（L*a*b*色空间）中评估色差的常用参数ΔE*ab和Δa*b*分别定义为圆圈和矩形。CIE 2000的色差参数ΔE00却是一个主轴在饱和度方向上的与人眼辨别临界区相接近的椭圆。在低饱和度的区域，重量系数SL，SC和SH都接近与1，使得整个椭圆更象个圆形。在高饱和度的区域，重量系数SC将比其他两个系数SL和SH变得更大，因此在饱和度方向上椭圆被拉得更狭长（饱和度的敏感度更低）。
CIEDE2000色差公式
色差计算工具，https://www.colortell.com/colorde
CIEDE2000色差公式如下：标准色和样品色。
引入不同的参数来调节三个色差△L‘、△C‘ab、△H‘ab在总色差中的权重，以提高色差计算结果与目视评判的一致性。
三个常量参数KL，KC和KH，用户可以根据不同的测量对象或者色彩品质控制要求来自定义它们的数值，以获得更灵活的色差计算方式。
计算步骤：
1、计算CIELAB公式中的L*、a*、b*、C*ab
2、计算L‘、a‘、C‘、h‘
G 表示 CIELAB 颜色空间的轴的调整因子。
3、计算△L‘、△C‘ab、△H‘ab
最新色差公式CIEDE2000 - 道客巴巴 https://www.doc88.com/p-1436669686686.html
4、计算SL、SC、SH、T和RT(由RC计算得来)
彩度C,色相角h。
T是色相角的函数，T和h之间的关系如下图：
CIEDE2000 色差公式宽容量椭圆如下图所示：
在色相角270°左右，颜色的描述范围【宽容量椭圆】有些失真
matlab程序：
function [DE,RT,hm,Dt,Rot]=CIEDE2000(LABREF,LAB,K) % DE2000 calculates colour difference between a reference and sample % using CIEDE2000 colour difference (as defined in Luo, Cui and Rigg (2000)) % % Input data can be single values or multiple values arranged in columns % LABREF can be a single value while LAB is a column % % Colour Engineering Toolbox % author: ?Phil Green % version: 1.1 % date: 17-01-2001 % modified: Zheng Yuanlin, 2010-08-23 % % set the values of parametric weighting factors KL,KC,KH if nargin>2 if length(K)>2 kL=K(1);kC=K(2);kH=K(3); end else kL=1;kC=1;kH=1; end %___________________________________________________________________ L=LABREF(:,1); a=LABREF(:,2); b=LABREF(:,3); C=(a.^2+b.^2).^0.5; Ls=LAB(:,1); as=LAB(:,2); bs=LAB(:,3); Cs=(as.^2+bs.^2).^0.5; %find G and recompute a', C' and h' Cm=(C+Cs)/2; G=0.5*(1-(Cm.^7./(Cm.^7+25^7)).^0.5); % a轴调整因子 a=(1+G).*a; as=(1+G).*as; C=(a.^2+b.^2).^0.5; h=hue_angle(a,b); Cs=(as.^2+bs.^2).^0.5; hs=hue_angle(as,bs); zcidx = (C.*Cs == 0); %% find the mean chroma and hue for each reference/sample pair Cm=(C+Cs)/2; % 彩度平均值 %% 色相角平均值 % hm=hmean_GS(h,hs,C,Cs); hm=(h+hs)/2; hm=hm-360*(abs(h-hs)>360); hm=hm+(hm<0)*360; hm(zcidx)=h(zcidx)+hs(zcidx); %% 色相角差 % Dh=angle_diff_GS(h,hs,C,Cs); Dh=hs-h; Dh = Dh - 2*180* (Dh > 180 ); Dh = Dh + 2*180* (Dh < (-180) ); Dh(zcidx)=0; %% rad=pi/180; DL=(L-Ls); % 亮度差 DC=Cs-C; % 彩度差 DH=2*((C.*Cs).^0.5).*sin(rad*(Dh)/2); % 色相差 % calculate T T=1-0.17*cos(rad*(hm-30))+0.24*cos(rad*2*hm)+0.32*cos(rad*(3*hm+6))-0.2*cos(rad*(4*hm-63)); %权重因子 SL, SC, SH SL=1+(0.015.*((L+Ls)./2-50).^2)./(20+((L+Ls)./2-50).^2).^0.5; SC=1+0.045.*Cm; SH=1+0.015.*Cm.*T; Dt=30*exp(-(((hm-275)/25).^2)); RC=2.*((Cm.^7)./(Cm.^7+25.^7)).^.5; RT=-sin(2*rad*Dt).*RC; Rot=RT.*(DC./(SC.*kC)).*(DH./(SH.*kH)); % 调整项 DE=((DL./(SL.*kL)).^2+(DC./(SC.*kC)).^2+(DH./(SH.*kH)).^2+Rot).^0.5; return
色差公式计算类:
internal class ColorDifferences { /// <CIE 1976 L*a*b*色差公式> /// </summary> /// <param name="L1"></param> /// <param name="a1"></param> /// <param name="b1"></param> /// <param name="L2"></param> /// <param name="a2"></param> /// <param name="b2"></param> /// <returns></returns> public static double delta_Eab(double L1, double a1, double b1, double L2, double a2, double b2) { double Eab = 0; // △Eab double chafang_L = 0; // (L1-L2)*(L1-L2) double chafang_a = 0; // (a1-a2)*(a1-a2) double chafang_b = 0; // (b1-b2)*(b1-b2) chafang_L = ( L1 - L2 ) * ( L1 - L2 ); //差-方 chafang_a = ( a1 - a2 ) * ( a1 - a2 ); chafang_b = ( b1 - b2 ) * ( b1 - b2 ); Eab = Math.Pow(chafang_L + chafang_a + chafang_b, 0.5); return Eab; } /// <CIE DE 2000色差公式> /// </summary> /// <param name="L1"></param> /// <param name="a1"></param> /// <param name="b1"></param> /// <param name="L2"></param> /// <param name="a2"></param> /// <param name="b2"></param> /// <returns></returns> public static double delta_E00(double L1, double a1, double b1, double L2, double a2, double b2) { //参考《现代颜色技术原理及应用》P88数据 double E00 = 0; //CIEDE2000色差E00 double LL1, LL2, aa1, aa2, bb1, bb2; //声明L' a' b' （1,2） double delta_LL, delta_CC, delta_hh, delta_HH; // 第二部的四个量 double kL, kC, kH; double RT = 0; //旋转函数RT double G = 0; //G表示CIELab 颜色空间a轴的调整因子,是彩度的函数. double mean_Cab = 0; //两个样品彩度的算术平均值 double SL, SC, SH, T; //------------------------------------------ //参考实验条件见P88 kL = 1; kC = 1; kH = 1; //------------------------------------------ mean_Cab = (CaiDu(a1, b1) + CaiDu(a2, b2)) / 2; double mean_Cab_pow7 = Math.Pow(mean_Cab, 7); //两彩度平均值的7次方 G = 0.5*(1-Math.Pow(mean_Cab_pow7 / (mean_Cab_pow7 + Math.Pow(25, 7)), 0.5)); LL1 = L1; aa1 = a1 * (1 + G); bb1 = b1; LL2 = L2; aa2 = a2 * (1 + G); bb2 = b2; double CC1, CC2; //两样本的彩度值 CC1 = CaiDu(aa1, bb1); CC2 = CaiDu(aa2, bb2); double hh1, hh2; //两样本的色调角 hh1 = SeDiaoJiao(aa1, bb1); hh2 = SeDiaoJiao(aa2, bb2); delta_LL = LL1 - LL2; delta_CC = CC1 - CC2; delta_hh = SeDiaoJiao(aa1, bb1) - SeDiaoJiao(aa2, bb2); delta_HH = 2 * Math.Sin(Math.PI*delta_hh / 360) * Math.Pow(CC1 * CC2, 0.5); //-------第三步-------------- //计算公式中的加权函数SL,SC,SH,T double mean_LL = (LL1 + LL2) / 2; double mean_CC = (CC1 + CC2) / 2; double mean_hh = (hh1 + hh2) / 2; SL = 1 + 0.015 * Math.Pow(mean_LL - 50, 2) / Math.Pow(20 + Math.Pow(mean_LL - 50, 2), 0.5); SC = 1 + 0.045 * mean_CC; T = 1 - 0.17 * Math.Cos((mean_hh - 30) * Math.PI / 180) + 0.24 * Math.Cos((2 * mean_hh) * Math.PI / 180) + 0.32 * Math.Cos((3 * mean_hh + 6) * Math.PI / 180) - 0.2 * Math.Cos((4 * mean_hh - 63) * Math.PI / 180); SH = 1 + 0.015 * mean_CC * T; //------第四步-------- //计算公式中的RT double mean_CC_pow7 = Math.Pow(mean_CC, 7); double RC = 2 * Math.Pow(mean_CC_pow7 / (mean_CC_pow7 + Math.Pow(25, 7)), 0.5); double delta_xita = 30 * Math.Exp(-Math.Pow((mean_hh - 275) / 25, 2)); //△θ 以°为单位 RT = -Math.Sin((2 * delta_xita) * Math.PI / 180) * RC; double L_item, C_item, H_item; L_item = delta_LL / (kL * SL); C_item = delta_CC / (kC * SC); H_item = delta_HH / (kH * SH); E00 = Math.Pow(L_item * L_item + C_item * C_item + H_item * H_item + RT * C_item * H_item, 0.5); return E00; } //彩度计算 private static double CaiDu(double a,double b) { double Cab = 0; Cab = Math.Pow(a * a + b * b, 0.5); return Cab; } //色调角计算 private static double SeDiaoJiao(double a, double b) { double h = 0; double hab = 0; h = (180 / Math.PI) * Math.Atan(b / a); //有正有负 if (a>0&&b>0) { hab = h; } else if (a<0&&b>0) { hab = 180 + h; } else if (a<0&&b<0) { hab = 180 + h; } else //a>0&&b<0 { hab = 360 + h; } return hab; } }
【转载自】
[图文]第3章 - - 色差、色深与白度 - 百度文库 https://wenku.baidu.com/view/e07185ccfad6195f302ba65d.html
色差公式的过去、现在、未来[译]_郑眼观色_新浪博客 http://blog.sina.com.cn/s/blog_4b892b79010007z0.html
CIEDE2000色差公式 – 色彩管理网 https://www.colortell.com/3424.html
色差公式的发展历程 – 色彩管理网 https://www.colortell.com/3233.html
新色差公式CIE2000 - 百度文库 https://wenku.baidu.com/view/74febcdaad51f01dc281f199.html?sxts=1561381510082
[转载]关于CIEDE2000的描述_GaoCheng_新浪博客 http://blog.sina.com.cn/s/blog_b15eea350102y1qh.html