颜色差异

颜色差异（Color difference），亦称颜色距离，是色彩学上的一个关注点。它量化了一个概念。在未量化之前，人们只能用形容词来大概描述这个概念，这使得对颜色要求严格的工作者们很不方便。颜色差异可以通过色彩空间内的欧氏距离简单计算得出，也可以使用国际照明委员会较为复杂、均匀的人类知觉公式计算。

## 欧氏距离

很多日常使用的“颜色差异”，是直接通过在一个“设备无关”的色彩空间里，进行欧氏距离的计算得到的。给定一个RGB（红绿蓝）的色彩空间，最简单的差异计算方式就是在这个三维空间里求两个点间的距离：

$$\mathrm{\Delta_{RGB}} = \sqrt{ (R_2-R_1)^2+(G_2-G_1)^2 + (B_2-B_1)^2 }$$

如果还要简单一点，经常可以去掉平方根：

$$\mathrm{\Delta_{RGB}}^2 = { (R_2-R_1)^2+(G_2-G_1)^2 + (B_2-B_1)^2 }$$

这样求得的距离可以用来比较两个颜色中哪一个更接近某个给定的颜色。如果把这些平方距离加在一起，会得到一组距离间的方差。

有不少人尝试将RGB三值加上权重，希望可以让得到的结果更加符合人类感官。一种做法是使用2、4、3：

$\sqrt{ 2 \times \Delta R^2 + 4 \times \Delta G^2 + 3 \times \Delta B^2 }$,

对于0~255的数值，有一种近似是：

$$\bar r = {{C_{1,R} + C_{2,R} } \over 2}$$

$$\Delta R = C_{1,R} - C_{2,R}$$

$$\Delta G = C_{1,G} - C_{2,G}$$

$$\Delta B = C_{1,B} - C_{2,B}$$

$$\Delta C = \sqrt {\left( {2 + {{\bar r} \over {256}}} \right) \times \Delta R^2 + 4 \times \Delta G^2 + \left( {2 + {{255 - \bar r} \over {256}}} \right) \times \Delta B^2 }$$

## Delta E

国际照明委员会（CIE）称他们的度量标准为“1=ΔE*ab”（也作“1=ΔE*”，另有“1=dE*”、“dE”、“Delta E”）。其中，Δ 是一个常用来表示“差异”的希腊字母，E表示Empfindung，也就是德文中的“感觉”。富有影响的赫尔曼·冯·亥姆霍兹和埃瓦尔德·赫林曾使用过它。

两个颜色之间的最小可觉差(JND)值，也就是有多大的1=ΔE*差距才能刚好被察觉，目前还没有定论。尽管并没有实践经验支持，'1.0'这个值经常被提到并用作JND。可是，1994年，Mahy et al研究并评估出了一个2.31=ΔE的JND值。然而，支撑该理论的CIELAB色彩空间中，知觉的非均匀特性（换言之，有些颜色在改变时，人眼对其十分敏感，而有的就分辨得不太清楚）导致这种理论不太站得住脚。这使得CIE用了接下来的几年时间修缮他们的定义，最终产生了更完善的（据CIE）1994年和2000年的公式。这些非均匀特性很重要，因为人眼对某些特定的颜色更敏感。一个好用的度量方法，应该把这些因素都考虑在内，才能使“刚刚好能被察觉到的差异”这种提法有意义。不然的话，可能出现这样的情况：对于两组分别含有一对不同的颜色、同样适用某个ΔE值的颜色组，其中一组颜色的差异能被人眼觉察出来，另一组却不能。

L、a和b三个值（见下方公式）一般是介于(-1,1)之间的双精度数值。Lab图像在保存时，L被存储为无符号（即“不分正负”）8位元整数，a、b则被存储为有符号的8位元整数，以保持文件大小。

### CIE76

1976年的这个色彩差异公式是首个用较为均匀的Lab空间计算色彩差异值的公式。Lab比RGB等空间在感官上均匀一点，所以得到的结果会更好。不过后来人们发现Lab色彩空间，尤其是在饱和度较高的区域里，并没有设计时预想的那么“感官上均匀”，所以几次更新了1=ΔE公式。CIE76公式会高估较为饱和的颜色之间的差异。

应用$({L^*_1},{a^*_1},{b^*_1})$和$({L^*_2},{a^*_2},{b^*_2})$两个L*a*b*色彩空间的颜色：

$$\Delta E_{ab}^* = \sqrt{ (L^*_2-L^*_1)^2+(a^*_2-a^*_1)^2 + (b^*_2-b^*_1)^2 }$$

本公式对应的最小可觉差：$\Delta E_{ab}^* \approx 2.3$

### CIE94

在汽车喷漆行业里，人们测量了大量有关人眼容差情况的数据。随着这些“应用特定”的数据的发布，人们将1976年的定义进行了进一步的发展，以更好地应对感知非均匀特性。该公式仍然采用的是Lab色彩空间。

1994版ΔE是在L*C*h*色彩空间下定义的。它的亮度、色品（或称“浓度”）、色调是从L*a*b*坐标系里计算来的。对于一个基准色$(L^*_1,a^*_1,b^*_1)$和另一个颜色$(L^*_2,a^*_2,b^*_2)$，它们的差异是：

$$\Delta E_{94}^* = \sqrt{ \left(\frac{\Delta L^*}{k_L S_L}\right)^2 + \left(\frac{\Delta C^*_{ab}}{k_C S_C}\right)^2 + \left(\frac{\Delta H^*_{ab}}{k_H S_H}\right)^2 }$$

其中：

$$\Delta L^* = L^*_1 - L^*_2$$

$$C^*_1 = \sqrt{ {a^*_1}^2 + {b^*_1}^2 }$$

$$C^*_2 = \sqrt{ {a^*_2}^2 + {b^*_2}^2 }$$

$$\Delta C^*_{ab} = C^*_1 - C^*_2$$

$$\Delta H^*_{ab} = \sqrt{ {\Delta E^*_{ab}}^2 - {\Delta L^*}^2 - {\Delta C^*_{ab}}^2 } = \sqrt{ {\Delta a^*}^2 + {\Delta b^*}^2 - {\Delta C^*_{ab}}^2 }$$

$$\Delta a^* = a^*_1 - a^*_2$$

$$\Delta b^* = b^*_1 - b^*_2$$

$$S_L = 1$$

$$S_C = 1+K_1 C^*_1$$

$$S_H = 1+K_2 C^*_1$$

另外，其中的kC和kH一般都取1，系数（权值）kL、K1和K2则依应用场合不同而取不同的值。如：

图形艺术｜纹理、纺织品

$k_L$｜1｜2

$K_1$｜0.045｜0.048

$K_2$｜0.015｜0.014

从几何学上讲，$\Delta H^*_{ab}$这个数量即是这两种颜色的相同色度圈的弦长的算术平均值。

### CIEDE2000

鉴于1994年的公式并没有充分解决感知非均匀特性的问题，CIE再次修缮了定义，并加入了5个修订系数：

- 色调旋转项(RT)，用来应对常出问题的蓝色区域（色相角度275°左右）。
- 中性色调补偿（对应L*C*h差异）
- 亮度补偿（SL）
- 色度补偿（SC）
- 色调补偿（SH）

$$\Delta E_{00}^* = \sqrt{ \left(\frac{\Delta L'}{k_L S_L}\right)^2 + \left(\frac{\Delta C'}{k_C S_C}\right)^2 + \left(\frac{\Delta H'}{k_H S_H}\right)^2 + R_T \frac{\Delta C'}{k_C S_C}\frac{\Delta H'}{k_H S_H} }$$

注：下面的公式请使用角度而非弧度，尤其在RT上需要注意。

kL、kC、kH一般取1。

$$\Delta L^\prime = L^*_2 - L^*_1$$

$$\bar{L} = \frac{L^*_1 + L^*_2}{2} \quad \bar{C} = \frac{C^*_1 + C^*_2}{2}$$

$$<pre> a_1^\prime = a_1^* + \frac{a_1^*}{2} \left( 1 - \sqrt{\frac{\bar{C}^7}{\bar{C}^7 + 25^7}} \right) \quad a_2^\prime = a_2^* + \frac{a_2^*}{2} \left( 1 - \sqrt{\frac{\bar{C}^7}{\bar{C}^7 + 25^7}} \right) </pre>$$

$$<pre> \bar{C}^\prime = \frac{C_1^\prime + C_2^\prime}{2} \mbox{ and } \Delta{C'}=C'_2-C'_1 \quad \mbox{where } C_1^\prime = \sqrt{a_1^{'^2} + b_1^{*^2}} \quad C_2^\prime = \sqrt{a_2^{'^2} + b_2^{*^2}} \quad </pre>$$

$$<pre> h_1^\prime=\text{atan2} (b_1^*, a_1^\prime) \mod 360^\circ, \quad h_2^\prime=\text{atan2} (b_2^*, a_2^\prime) \mod 360^\circ </pre>$$

Note: The inverse tangent (tan−1) can be computed using a common library routine `atan2(b, a&prime;)` which usually has a range from −π to π radians; color specifications are given in 0 to 360 degrees, so some adjustment is needed. The inverse tangent is indeterminate if both a′ and b are zero (which also means that the corresponding C′ is zero); in that case, set the hue angle to zero. See .

$$<pre> \Delta h' = \begin{cases} h_2^\prime - h_1^\prime & \left| h_1^\prime - h_2^\prime \right| \leq 180^\circ \\ h_2^\prime - h_1^\prime + 360^\circ & \left| h_1^\prime - h_2^\prime \right| > 180^\circ, h_2^\prime \leq h_1^\prime \\ h_2^\prime - h_1^\prime - 360^\circ & \left| h_1^\prime - h_2^\prime \right| > 180^\circ, h_2^\prime > h_1^\prime \end{cases} </pre>$$

Note: When either C′1 or C′2 is zero, then Δh′ is irrelevant and may be set to zero. See .

$$<pre> \Delta H^\prime = 2 \sqrt{C_1^\prime C_2^\prime} \sin (\Delta h^\prime/2), \quad \bar{H}^\prime=\begin{cases} (h_1^\prime + h_2^\prime + 360^\circ)/2 & \left| h_1^\prime - h_2^\prime \right| > 180^\circ \\ (h_1^\prime + h_2^\prime)/2 & \left| h_1^\prime - h_2^\prime \right| \leq 180^\circ \end{cases} </pre>$$

Note: When either C′1 or C′2 is zero, then H′ is h′1+h′2 (no divide by 2; essentially, if one angle is indeterminate, then use the other angle as the average; relies on indeterminate angle being set to zero). See stating most implementations on the internet at the time had "an error in the computation of average hue".

$$<pre> T = 1 - 0.17 \cos ( \bar{H}^\prime - 30^\circ ) + 0.24 \cos (2\bar{H}^\prime) + 0.32 \cos (3\bar{H}^\prime + 6^\circ ) - 0.20 \cos (4\bar{H}^\prime - 63^\circ) </pre>$$

$$<pre> S_L = 1 + \frac{0.015 \left( \bar{L} - 50 \right)^2}{\sqrt{20 + {\left(\bar{L} - 50 \right)}^2} } \quad S_C = 1+0.045 \bar{C}^\prime \quad S_H = 1+0.015 \bar{C}^\prime T </pre>$$

$$R_T = -2 \sqrt{\frac{\bar{C}'^7}{\bar{C}'^7+25^7}} \sin \left[ 60^\circ \cdot \exp \left( -\left[ \frac{\bar{H}'-275^\circ}{25^\circ} \right]^2 \right) \right]$$

## 应用

### 差异分类

在讨论人类感知时，一般使用ΔE计算色彩差异。如前所述，一般认为在ΔE系统中，1.0是刚好注意得到的量。取决于产品的受众，有时候会改用2.5作为可以接受的限制。例如一般情况下，人们会说“叶子和草地是绿的”，但仔细看就能发现叶子和草地都有从黄绿到蓝绿等不同的差异。

这类差异在很多产品的品质控制中非常重要：在生产和运输中经常会有一些来自内部或外部的随机扰动，造成实际上显示的颜色和之前校准出的结果有一些差异。在最理想的情况下，要做到“一般的观察者”不能察觉这类变化。胶印过程中通常允许最高2到4个单位的ΔE色彩差异；而家用打印机造成的色差可能更大。

## 外部链接

- Michel Leonard 的 CIE ΔE 2000 实现，涵盖 20 多种编程语言。
- Bruce Lindbloom's color difference calculator . Uses all metrics defined herein.
- The CIEDE2000 Color-Difference Formula , by Gaurav Sharma. Implementations in MATLAB and Excel.
- Color Similarity Tool , by Dragan Bajcic, Implementations in PHP and Javascript.
- python-colormath . Implementation in Python.
- Precise Color Communication Konica Minolta Sensing.