色差

色彩科学において、色差（しきさ）あるいは 色の距離 （いろのきょり）は、2つの色の間に定義される指標の一つである。 従来は官能評価することしかできなかった「色の差」の概念が、定量的に検討できるようになった。重要な色の判定等において、これら特性の定量化は極めて重要である。通常はデバイス非依存の色空間におけるユークリッド距離で定義される。 ほとんどの色差の定義は色空間における距離であり、距離の決定方法はユークリッド距離である。ある色空間がRGB (Red, Green, Blue) の組み合わせから構成され、そこで色差を計算する場合、R, G, Bからなる線形3次元空間で構成される色空間で計算する。 d i s t a n c e = ( R 2 − R 1 ) 2 + ( G 2 − G 1 ) 2 + ( B 2 − B 1 ) 2 {\displaystyle \mathrm {distance} ={\sqrt {(R_{2}-R_{1})^{2}+(G_{2}-G_{1})^{2}+(B_{2}-B_{1})^{2}}}} 単純化のため、しばしば平方根を取り除いて計算が行われる: d i s t a n c e 2 = ( R 2 − R 1 ) 2 + ( G 2 − G 1 ) 2 + ( B 2 − B 1 ) 2 {\displaystyle \mathrm {distance} ^{2}={(R_{2}-R_{1})^{2}+(G_{2}-G_{1})^{2}+(B_{2}-B_{1})^{2}}} この式は、単一の色と色同士を比較し、近いかどうかを判断する目的で利用できる。しかし自乗された色差を加算していく場合、色差の分散値となってしまう。 RGB値をヒトの視覚に一致させる試みは、これまで数多く行われてきて、各成分は等しく重み付けされた (red 30%, green 59%, blue 11%)ものが考案された。しかしながら, これらの試みでは色の決定において良い結果は得られず、明るさの影響を受け、ヒトの視覚がこれらの色に対して少ない許容値を持つことがわかった。より近い近似式は係数2, 4, および 3を用いて: 2 × Δ R 2 + 4 × Δ G 2 + 3 × Δ B 2 {\displaystyle {\sqrt {2\times \Delta R^{2}+4\times \Delta G^{2}+3\times \Delta B^{2}}}}   より簡単な近似式は、 r ¯ = C 1 , R + C 2 , R 2 {\displaystyle {\bar {r}}={{C_{1,R}+C_{2,R}} \over 2}} Δ R = C 1 , R − C 2 , R {\displaystyle \Delta R=C_{1,R}-C_{2,R}} Δ G = C 1 , G − C 2 , G {\displaystyle \Delta G=C_{1,G}-C_{2,G}} Δ B = C 1 , B − C 2 , B {\displaystyle \Delta B=C_{1,B}-C_{2,B}} Δ C = ( 2 + r ¯ 256 ) × Δ R 2 + 4 × Δ G 2 + ( 2 + 255 − r ¯ 256 ) × Δ B 2 {\displaystyle \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}}}} 他にもいくつかの色差に関する数式があり、例えばHSVのように色相を円状にしたり、円筒座標や三角錐座標に色を配置したりしている。しかしながらそれらはRGBの派生であり、単純なユークリッド距離に近い傾向を表すヒトの視覚特性を考慮していない。 国際照明委員会 (CIE) は、色差をΔE*ab (または ΔE*, dE*, dE, あるいは "delta E") と定義している。デルタはギリシャ文字で距離を表すことが多く、Eはドイツ語のEmpfindung (感覚)を意味する。ヘルマン・フォン・ヘルムホルツとエヴァルト・ヘリングの業績に因んで名付けられた。 様々な研究成果に基いて、複数の丁度可知差異と呼ばれるΔEの値が提案されている。この丁度可知差異の値は、経験論に基づかずに"1.0"という値が広く知られているが、近年の研究では、Mahy et al. (1994) の研究によると2.3 ΔEが丁度可知差異として確認されている。しかし、元となるCIE Lab色空間における知覚的非均一性があるため、より優れたCIE 1994および2000の計算式がCIEにより定義された。 これらの知覚的非均一性を考慮することは、ヒトの目が違う色同士よりも同じ色同士に対してより敏感に色覚することからも、重要である。丁度可知差異の概念が意味を成すためにも、この非均一性を考慮にいれるべきである。さもないと、ある２つの色の間で、目には知覚不可能でもスペクトル上は明確な差異がある場合においてΔEの値が意味を成さなくなってしまう。 CIE 1976色差計算式は、初めて定義された色差式であり、CIELAB座標におけるユークリッド距離として計算される。CIELAB色空間は設計当初の意図よりも、特に彩度の大きい領域で知覚的に非均一であることが明らかになってきたため、この計算式はCIE 1994 および 2000の計算式に置き換えられている。 すなわち、彩度の大きいの領域では、色差が大きすぎるように計算されてしまう。 L*a*b*色空間において、２つの色  ( L 1 ∗ , a 1 ∗ , b 1 ∗ ) {\displaystyle ({L_{1}^{*}},{a_{1}^{*}},{b_{1}^{*}})} および  ( L 2 ∗ , a 2 ∗ , b 2 ∗ ) {\displaystyle ({L_{2}^{*}},{a_{2}^{*}},{b_{2}^{*}})} を用いて: Δ E a b ∗ = ( L 2 ∗ − L 1 ∗ ) 2 + ( a 2 ∗ − a 1 ∗ ) 2 + ( b 2 ∗ − b 1 ∗ ) 2 {\displaystyle \Delta E_{ab}^{*}={\sqrt {(L_{2}^{*}-L_{1}^{*})^{2}+(a_{2}^{*}-a_{1}^{*})^{2}+(b_{2}^{*}-b_{1}^{*})^{2}}}} Δ E a b ∗ ≈ 2.3 {\displaystyle \Delta E_{ab}^{*}\approx 2.3}  が丁度可知差異に相当する。 CIE 1976計算式は、L*a*b*色空間はそのまま使用し、知覚的非均一性に対応するよう拡張されている。自動車用塗料の許容値から派生された、特定用途の重み付け係数を導入している。 ΔE (1994)は、L*a*b*色空間から計算される、L*C*h*色空間における明度、彩度、色相の差異から算出される。 基準色 を ( L 1 ∗ , a 1 ∗ , b 1 ∗ ) {\displaystyle (L_{1}^{*},a_{1}^{*},b_{1}^{*})} 、比較色を  ( L 2 ∗ , a 2 ∗ , b 2 ∗ ) {\displaystyle (L_{2}^{*},a_{2}^{*},b_{2}^{*})} とすると、色差は: Δ E 94 ∗ = ( Δ L ∗ k L S L ) 2 + ( Δ C a b ∗ k C S C ) 2 + ( Δ H a b ∗ k H S H ) 2 {\displaystyle \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}}}} ここで: Δ L ∗ = L 1 ∗ − L 2 ∗ {\displaystyle \Delta L^{*}=L_{1}^{*}-L_{2}^{*}} C 1 ∗ = a 1 ∗ 2 + b 1 ∗ 2 {\displaystyle C_{1}^{*}={\sqrt {{a_{1}^{*}}^{2}+{b_{1}^{*}}^{2}}}} C 2 ∗ = a 2 ∗ 2 + b 2 ∗ 2 {\displaystyle C_{2}^{*}={\sqrt {{a_{2}^{*}}^{2}+{b_{2}^{*}}^{2}}}} Δ C a b ∗ = C 1 ∗ − C 2 ∗ {\displaystyle \Delta C_{ab}^{*}=C_{1}^{*}-C_{2}^{*}} Δ H a b ∗ = Δ E a b ∗ 2 − Δ L ∗ 2 − Δ C a b ∗ 2 = Δ a ∗ 2 + Δ b ∗ 2 − Δ C a b ∗ 2 {\displaystyle \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}}}} Δ a ∗ = a 1 ∗ − a 2 ∗ {\displaystyle \Delta a^{*}=a_{1}^{*}-a_{2}^{*}} Δ b ∗ = b 1 ∗ − b 2 ∗ {\displaystyle \Delta b^{*}=b_{1}^{*}-b_{2}^{*}} S L = 1 {\displaystyle S_{L}=1} S C = 1 + K 1 C 1 ∗ {\displaystyle S_{C}=1+K_{1}C_{1}^{*}} S H = 1 + K 2 C 1 ∗ {\displaystyle S_{H}=1+K_{2}C_{1}^{*}} kC および kH は、通常の同じ値をとり、重み係数 kL、K1 および K2 は、適用用途により異なる: 幾何学的には、 Δ H a b ∗ {\displaystyle \Delta H_{ab}^{*}} は弦の長さに等しく対応する。 CIE 1994の定義は、知覚的均一性を十分に確保できていなかったため、CIEは定義を見直し、5箇所の訂正を追加した: 青色領域 (色相角度275°付近) における問題に対処するための色相回転項 (RT) : 中立色の補正項(L*C*h*の違いにおける最大値) 明度補正項 (SL) 彩度補正項 (SC) 色相補正項 (SH) Δ E 00 ∗ = ( Δ L ′ k L S L ) 2 + ( Δ C ′ k C S C ) 2 + ( Δ H ′ k H S H ) 2 + R T Δ C ′ k C S C Δ H ′ k H S H {\displaystyle \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 は、通常同一値。 Δ L ′ = L 2 ∗ − L 1 ∗ {\displaystyle \Delta L^{\prime }=L_{2}^{*}-L_{1}^{*}} L ¯ = L 1 ∗ + L 2 ∗ 2 C ¯ = C 1 ∗ + C 2 ∗ 2 {\displaystyle {\bar {L}}={\frac {L_{1}^{*}+L_{2}^{*}}{2}}\quad {\bar {C}}={\frac {C_{1}^{*}+C_{2}^{*}}{2}}} a 1 ′ = a 1 ∗ + a 1 ∗ 2 ( 1 − C ¯ 7 C ¯ 7 + 25 7 ) a 2 ′ = a 2 ∗ + a 2 ∗ 2 ( 1 − C ¯ 7 C ¯ 7 + 25 7 ) {\displaystyle 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)} C ¯ ′ = C 1 ′ + C 2 ′ 2  and  Δ C ′ = C 2 ′ − C 1 ′ where  C 1 ′ = a 1 ′ 2 + b 1 ∗ 2 C 2 ′ = a 2 ′ 2 + b 2 ∗ 2 {\displaystyle {\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 ′ = atan2 ( b 1 ∗ , a 1 ′ ) mod 360 ∘ , h 2 ′ = atan2 ( b 2 ∗ , a 2 ′ ) mod 360 ∘ {\displaystyle 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 }} 注: 逆正接 (tan−1) は通常コンピューター上では、−π から π ラジアンの範囲でのatan2(b,a') を用いて計算できる。一方で色相は0から360度の範囲であるため、変換が必要である。a' および b がゼロの場合、逆正接は不定値となる (すなわち対応する C' もゼロとなる)。この場合、色相角はゼロを代入する。Sharma 2005, eqn. 7参照のこと。 Δ h ′ = { h 2 ′ − h 1 ′ | h 1 ′ − h 2 ′ | ≤ 180 ∘ h 2 ′ − h 1 ′ + 360 ∘ | h 1 ′ − h 2 ′ | > 180 ∘ , h 2 ′ ≤ h 1 ′ h 2 ′ − h 1 ′ − 360 ∘ | h 1 ′ − h 2 ′ | > 180 ∘ , h 2 ′ > h 1 ′ {\displaystyle \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}}} 注: C′1 あるいは C′2 がゼロの場合、Δh′ は無関係となり、ゼロが代入される。 Sharma 2005, eqn. 10参照のこと。 Harvnb Δ H ′ = 2 C 1 ′ C 2 ′ sin ⁡ ( Δ h ′ / 2 ) , H ¯ ′ = { ( h 1 ′ + h 2 ′ ) / 2 | h 1 ′ − h 2 ′ | ≤ 180 ∘ ( h 1 ′ + h 2 ′ + 360 ∘ ) / 2 | h 1 ′ − h 2 ′ | > 180 ∘ , h 1 ′ + h 2 ′ < 360 ∘ ( h 1 ′ + h 2 ′ − 360 ∘ ) / 2 | h 1 ′ − h 2 ′ | > 180 ∘ , h 1 ′ + h 2 ′ ≥ 360 ∘ {\displaystyle \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 })/2&\left|h_{1}^{\prime }-h_{2}^{\prime }\right|\leq 180^{\circ }\\(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 }<360^{\circ }\\(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 }\geq 360^{\circ }\end{cases}}} 注: C′1 あるいは C′2 がゼロの場合、H′ は h′1+h′2 となる (2で除算しない; 本質的には、ひとつの角度が不定の場合、もう一方の角度を使って平均を取る。不定の角度がゼロであることを前提としている)。 Sharma 2005, eqn. 7 and p. 23参照のこと。インターネット上のほとんどの実装に、平均色相角の計算エラーがある、と述べている。 T = 1 − 0.17 cos ⁡ ( H ¯ ′ − 30 ∘ ) + 0.24 cos ⁡ ( 2 H ¯ ′ ) + 0.32 cos ⁡ ( 3 H ¯ ′ + 6 ∘ ) − 0.20 cos ⁡ ( 4 H ¯ ′ − 63 ∘ ) {\displaystyle 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 + 0.015 ( L ¯ − 50 ) 2 20 + ( L ¯ − 50 ) 2 S C = 1 + 0.045 C ¯ ′ S H = 1 + 0.015 C ¯ ′ T {\displaystyle 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 C ¯ ′ 7 C ¯ ′ 7 + 25 7 sin ⁡ [ 60 ∘ ⋅ exp ⁡ ( − [ H ¯ ′ − 275 ∘ 25 ∘ ] 2 ) ] {\displaystyle 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]} 許容値は、「基準色に対してどのような色が見分けがつかないほど近いか？」を定義する。もしも2色間の距離が知覚的均一であれば、答えは単純に「丁度可知差異以内のすべての色」である。 このことは色域内において知覚的均一であることを前提としている。 さもないと、許容値は参照値を含む関数となってしまい、複雑化してしまう。 例えばCIE 1931 色空間は、L* (明度)　が一定の場合の許容値を表す線がマクアダムの楕円として定義されている。図で示すとおり、許容値を表す楕円のサイズは、場所によりバラバラであることがわかる。 このことにより、CIELUV と CIELAB色空間の定義へとつながった。 より一般的に、明度を変動させると、許容値は楕円体となる。上記のように距離の表現における重み付け係数を増加させると、各々の係数の軸のおいて楕円体の体積を増やす方向に働く。 CIELAB ^ http://www.compuphase.com/cmetric.htm ^ Backhaus, W.; Kliegl, R.; Werner, J. S. (1998). Color Vision: Perspectives from Different Disciplines. Walter de Gruyter. p. 188. ISBN 9783110154313. https://books.google.com/books?id=DrduOSrOFegC&pg=PA188&lpg=PA188 2014年12月2日閲覧。.  ^ Valberg, A. (2005). Light Vision Color. Wiley. p. 278. ISBN 9780470849026. https://books.google.com/books?id=OoESifAi9ZsC 2014年12月2日閲覧。.  ^ Real World Color Management, Second Edition (Bruce Fraser) ^ Evaluation of the CIE Color Difference Formulas ^ Sharma, Gaurav (2003). Digital Color Imaging Handbook (1.7.2 ed.). CRC Press. ISBN 0-8493-0900-X. https://books.google.com/books?id=OxlBqY67rl0C&pg=PA31&vq=1.42&dq=jnd+gaurav+sharma.  ^ “Delta E: The Color Difference”. Colorwiki.com. 2009年4月16日閲覧。 ^ Called such because the operator is not en:Commutative commutative. ^ Lindbloom, Bruce Justin. “Delta E (CIE 1994)”. Brucelindbloom.com. 2011年3月23日閲覧。 ^ “Colour Difference Software by David Heggie”. Colorpro.com (1995年12月19日). 2009年4月16日閲覧。 ^ ﻿Colorimetry - Part 4: CIE 1976 L*a*b* Colour Space﻿ (Report). Draft Standard. CIE. (2007). CIE DS 014-4.3/E:2007.  ^ Klein, Georg A.. Industrial Color Physics. p. 147. ISBN 978-1-4419-1196-4.  ^ Sharma, Gaurav; Wu, Wencheng; Dalal, Edul N. (2005年). “The CIEDE2000 color-difference formula: Implementation notes, supplementary test data, and mathematical observations”. Color Research & Applications (Wiley Interscience) 30 (1): 21–30. doi:10.1002/col.20070. http://www.ece.rochester.edu/~gsharma/ciede2000/ciede2000noteCRNA.pdf.  ^ Lindbloom, Bruce Justin. “Delta E (CIE 2000)”. Brucelindbloom.com. 2009年4月16日閲覧。 ^ The "Blue Turns Purple" Problem, Bruce Lindbloom ^ Susan Hughes (1998年1月14日). “A guide to Understanding Color Tolerancing”. 2014年12月2日閲覧。 Bruce Lindbloomの色差計算 CIE dE2000色差式: Gaurav Sharma氏によるMATLABおよびExcelでの実装