颜色差异

**顏色差異**（英語：**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）稱他們的度量標準為「Δ_E_\*ab」（也作「Δ_E\*_」，另有「dE\*」、「dE」、「Delta E」）。其中，Δ 是一個常用來表示「差異」的希臘字母，**E** 表示 Empfindung，也就是德文中的「感覺」。富有影響的赫爾曼・馮・亥姆霍茲和埃瓦爾德・赫林曾使用過它。

兩個顏色之間的最小可覺差 (JND) 值，也就是有多大的 Δ_E\*_差距才能剛好被察覺，目前還沒有定論。儘管並沒有實踐經驗支持，'1.0' 這個值經常被提到並用作 JND。可是，1994 年，Mahy et al 研究並評估出了一個 2.3Δ_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 色彩空間，尤其是在飽和度較高的區域裡，並沒有設計時預想的那麼「感官上均勻」，所以幾次更新了 Δ_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_、_K_1 和 _K_2 則依應用場合不同而取不同的值。如：

從幾何學上講， $\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}}$

　
: $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)$

　
: ${\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 $

　
: $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 }$

　
: **Note:** The inverse tangent (tan−1) can be computed using a common library routine `atan2(b, a′)` 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 Sharma 2005，eqn. 7.

　
: $\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}}$

　
: **Note:** When either _C′_1 or _C′_2 is zero, then Δh′ is irrelevant and may be set to zero. See Sharma 2005，eqn. 10.

　
: $\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}}$

　
: **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 Sharma 2005，eqn. 7 and p. 23 stating most implementations on the internet at the time had "an error in the computation of average hue".

　
: $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 })$

　
: $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$

　
: $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_ 色彩差異；而家用印表機造成的色差可能更大。

## 腳註

## 延伸閱讀

* Robertson, Alan R. Historical development of CIE recommended color difference equations. Color Research & Application. 1990, **15** (3): 167–170. doi:10.1002/col.5080150308.
* Melgosa, M.; Quesada, J. J. and Hita, E. Uniformity of some recent color metrics tested with an accurate color-difference tolerance dataset. Applied Optics. December 1994, **33** (34): 8069–8077. PMID 20963027. doi:10.1364/AO.33.008069.
* McDonald, Roderick; Hill, MacDonald, Nobbs, Rigg, Sinclair, Smith. Roderick McDonald , 編. Colour Physics for Industry 2E. Society of Dyers and Colourists. 1997. ISBN 0-901956-70-8.

## 外部連結

* 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.