色差のソースを表示

'''色差'''（しきさ、{{lang-en-short|color difference, color distance}}）とは、色彩学において2つの[[色]]の間に定義される指標の一つである。色差が大きいほど区別しやすく、色差が小さいほど区別しにくくなる。 従来は官能評価することしかできなかった「色の差」の概念が、色差を用いることにより定量的に検討できるようになった。重要な色の判定等において、これら特性の定量化は極めて重要である。通常はデバイス非依存の色空間における[[ユークリッド距離]]で定義される。

== ユークリッド距離 ==
ほとんどの色差の定義は色空間における距離であり、距離の決定方法はユークリッド距離である。ある色空間がRGB (Red, Green, Blue) の組み合わせから構成され、そこで色差を計算する場合、R, G, Bからなる線形3次元空間で構成される色空間で計算する。
: <math>\mathrm{distance} = \sqrt{ (R_2-R_1)^2+(G_2-G_1)^2 + (B_2-B_1)^2 }</math>
単純化のため、しばしば平方根を取り除いて計算が行われる:
: <math>\mathrm{distance}^2 = { (R_2-R_1)^2+(G_2-G_1)^2 + (B_2-B_1)^2 }</math>
この式は、単一の色と色同士を比較し、近いかどうかを判断する目的で利用できる。しかし自乗された色差を加算していく場合、色差の[[分散 (確率論)|分散値]]となってしまう。

RGB値をヒトの視覚に一致させる試みは、これまで数多く行われてきて、各成分は等しく重み付けされた（red 30%, green 59%, blue 11%）ものが考案された。しかしながら, これらの試みでは色の決定において良い結果は得られず、明るさの影響を受け、ヒトの視覚がこれらの色に対して少ない許容値を持つことがわかった。より近い近似式は係数2, 4, および 3を用いて: <math>\sqrt{ 2 \times \Delta R^2 + 4 \times \Delta G^2 + 3 \times \Delta B^2 }</math> <br>より簡単な近似式は、<ref>http://www.compuphase.com/cmetric.htm</ref>
: <math>\bar r = {{C_{1,R}  + C_{2,R} } \over 2}</math>
: <math>\Delta R = C_{1,R}  - C_{2,R}</math>
: <math>\Delta G = C_{1,G}  - C_{2,G}</math>
: <math>\Delta B = C_{1,B}  - C_{2,B}</math>
: <math>\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 }</math>
他にもいくつかの色差に関する数式があり、例えばHSVのように色相を円状にしたり、円筒座標や三角錐座標に色を配置したりしている。しかしながらそれらはRGBの派生であり、単純なユークリッド距離に近い傾向を表すヒトの視覚特性を考慮していない。

== LAB Delta E ==
[[国際照明委員会]] (CIE) は、色差を{{Math|Δ''E''<sup>*</sup><sub>''ab''</sub>}} (または {{Math|Δ''E*''}}, {{Math|''dE''*}}, dE, あるいは "delta E") と定義している。[[デルタ]]は[[ギリシャ文字]]で距離を表すことが多く、'''E'''はドイツ語の''Empfindung (感覚)''を意味する。[[ヘルマン・フォン・ヘルムホルツ]]と[[エヴァルト・ヘリング]]の業績に因んで名付けられた。<ref name="google">{{cite book|title=Color Vision: Perspectives from Different Disciplines|url=https://books.google.com/books?id=DrduOSrOFegC&pg=PA188&lpg=PA188|date=1998|publisher=Walter de Gruyter|page=188|last1=Backhaus|last2=Kliegl|last3=Werner|first1=W.|first2=R.|first3=J. S.|accessdate=2014-12-02|isbn=9783110154313}}</ref><ref name="google2">{{cite book|last=Valberg|first=A.|title=Light Vision Color|url=https://books.google.com/books?id=OoESifAi9ZsC|date=2005|publisher=Wiley|page=278|accessdate=2014-12-02|isbn=9780470849026}}</ref>

様々な研究成果に基いて、複数の[[丁度可知差異]]と呼ばれる{{Math|Δ''E''}}の値が提案されている。この丁度可知差異の値は、経験論に基づかずに"1.0"という値が広く知られているが、近年の研究では、Mahy et al. (1994) の研究によると{{Math|2.3 Δ''E''}}が丁度可知差異として確認されている。しかし、元となる[[Lab色空間|CIE Lab色空間]]における知覚的非均等性があるため、より優れたCIE 1994および2000の計算式がCIEにより定義された。<ref>Real World Color Management, Second Edition (Bruce Fraser)</ref> これらの知覚的非均等性を考慮することは、ヒトの目が違う色同士よりも同じ色同士に対してより敏感に色覚することからも、重要である。丁度可知差異の概念が意味を成すためにも、この非均等性を考慮にいれるべきである。さもないと、ある２つの色の間で、目には知覚不可能でもスペクトル上は明確な差異がある場合において{{Math|Δ''E''}}の値が意味を成さなくなってしまう。<ref>[https://archive.is/20121203082912/http://www.aim-dtp.net/aim/evaluation/cie_de/index.htm Evaluation of the CIE Color Difference Formulas]</ref>

=== CIE76 (delta E76) ===
CIE 1976色差計算式は、初めて定義された色差式であり、CIELAB座標におけるユークリッド距離として計算される。CIELAB色空間は設計当初の意図よりも、特に彩度の大きい領域で知覚的に均等ではないことが明らかになってきたため、この計算式はCIE 1994 および 2000の計算式に置き換えられている。 すなわち、彩度が大きい領域では、色差が大きすぎるように計算されてしまう。

[[Lab色空間|L*a*b*]]色空間において、２つの色 <math>({L^*_1},{a^*_1},{b^*_1})</math> および <math>({L^*_2},{a^*_2},{b^*_2})</math>を用いて:
:<math>\Delta E_{ab}^* = \sqrt{ (L^*_2-L^*_1)^2+(a^*_2-a^*_1)^2 + (b^*_2-b^*_1)^2 }</math>
<math>\Delta E_{ab}^* \approx 2.3</math> が丁度可知差異に相当する。<ref>{{cite book|last=Sharma|first=Gaurav|title=Digital Color Imaging Handbook|url=https://books.google.com/books?id=OxlBqY67rl0C&pg=PA31&vq=1.42&dq=jnd+gaurav+sharma|edition=1.7.2|year=2003|publisher=[[CRC Press]]|isbn=0-8493-0900-X}}</ref>

=== CIE94 (delta E94) ===
CIE 1976計算式は、L*a*b*色空間はそのまま使用し、知覚的非均等性に対応するよう拡張されている。自動車用塗料の許容値から派生された、特定用途の重み付け係数を導入している。<ref>{{cite web|url=http://www.colorwiki.com/wiki/Delta_E:_The_Color_Difference|title=Delta E: The Color Difference|accessdate=2009-04-16|publisher=Colorwiki.com}}</ref>

Δ''E'' (1994)は、L*a*b*色空間から計算される、[[彩度|L*C*h*色空間]]における明度、彩度、色相の差異から算出される。 基準色<ref>Called such because the operator is not [[:en:Commutative commutative]]. </ref> を<math>(L^*_1,a^*_1,b^*_1)</math>、比較色を <math>(L^*_2,a^*_2,b^*_2)</math>とすると、色差は:<ref>{{cite web|url=http://www.brucelindbloom.com/Eqn_DeltaE_CIE94.html|title=Delta E (CIE 1994)|accessdate=2011-03-23|publisher=Brucelindbloom.com|first=Bruce Justin|last=Lindbloom}}</ref><ref>{{cite web|url=http://www.colorpro.com/info/software/heggie.html|title=Colour Difference Software by David Heggie|accessdate=2009-04-16|date=1995-12-19|publisher=Colorpro.com}}</ref><ref>{{cite report|title=Colorimetry - Part 4: CIE 1976 L*a*b* Colour Space|series=Draft Standard|accessdate=2011-03-23|publisher=CIE|id=CIE DS 014-4.3/E:2007|date=2007}}</ref>
:<math>\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 }</math>
ここで:
: <math>\Delta L^* = L^*_1 - L^*_2</math>
: <math>C^*_1 = \sqrt{ {a^*_1}^2 + {b^*_1}^2 }</math>
: <math>C^*_2 = \sqrt{ {a^*_2}^2 + {b^*_2}^2 }</math>
: <math>\Delta C^*_{ab} = C^*_1 - C^*_2</math>
: <math>\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 }</math>
: <math>\Delta a^* = a^*_1 - a^*_2</math>
: <math>\Delta b^* = b^*_1 - b^*_2</math>
: <math>S_L = 1</math>
: <math>S_C = 1+K_1 C^*_1</math>
: <math>S_H = 1+K_2 C^*_1</math>
''k<sub>C</sub>'' および ''k<sub>H</sub>'' は、通常の同じ値をとり、重み係数 ''k<sub>L</sub>''、''K''<sub>1</sub> および ''K''<sub>2</sub> は、適用用途により異なる:

{| class="wikitable" border="1"
|-
!
! グラフィック・アーツ
! テキスタイル
|-
| <math>k_L</math>
| 1
| 2
|-
| <math>K_1</math>
| 0.045
| 0.048
|-
| <math>K_2</math>
| 0.015
| 0.014
|}

幾何学的には、<math>\Delta H^*_{ab}</math>は弦の長さに等しく対応する。
<ref>{{cite book|last=Klein|first=Georg A.|title=Industrial Color Physics|page=147|isbn=978-1-4419-1196-4}}</ref>

=== CIEDE2000 (delta E00) ===
CIE 1994の定義は、知覚的均等性を十分に確保できていなかったため、CIEは定義を見直し、5箇所の訂正を追加した:<ref>{{cite journal|last=Sharma|first=Gaurav|year=2005|title=The CIEDE2000 color-difference formula: Implementation notes, supplementary test data, and mathematical observations|url=http://www.ece.rochester.edu/~gsharma/ciede2000/ciede2000noteCRNA.pdf|journal=Color Research & Applications|volume=30|issue=1|pages=21–30|publisher=[[Wiley Interscience]]|ref=CITEREFSharma2005|first2=Wencheng|last2=Wu|first3=Edul N.|last3=Dalal|doi=10.1002/col.20070}}</ref><ref>{{cite web|url=http://www.brucelindbloom.com/Eqn_DeltaE_CIE2000.html|title=Delta E (CIE 2000)|accessdate=2009-04-16|publisher=Brucelindbloom.com|first=Bruce Justin|last=Lindbloom}}</ref>
* 青色領域 (色相角度275°付近) における問題に対処するための色相回転項 (R<sub>T</sub>) :<ref>[http://www.brucelindbloom.com/MunsellCalcHelp.html#BluePurple The "Blue Turns Purple" Problem], Bruce Lindbloom</ref>
* 中立色の補正項(L*C*h*の違いにおける最大値)
* 明度補正項 (S<sub>L</sub>)
* 彩度補正項 (S<sub>C</sub>)
* 色相補正項 (S<sub>H</sub>)
:<math>\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} }</math>

:: <small>'''注:''' 以下の式はラジアンでなく度を用いる; 特に''R<sub>T</sub>''への影響が大きい。</small>
:: <small>''k<sub>L</sub>''、 ''k<sub>C </sub>''および ''k<sub>H</sub>'' は、通常同一値。</small>

:<math>\Delta L^\prime = L^*_2 - L^*_1</math>

:<math>\bar{L} = \frac{L^*_1 + L^*_2}{2} \quad \bar{C} = \frac{C^*_1 + C^*_2}{2}</math>

:<math>
  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)
</math>

:<math>
  \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
</math>

:<math>
  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
</math>

:: <small>'''注:''' 逆正接 (tan<sup>&#x2212;1</sup>) は通常コンピューター上では、[[atan2]]用いて、−π から π ラジアンの範囲での<code>atan2(b,a&#x27;)</code> により計算できる。一方で色相は0から360度の範囲であるため、変換が必要である。</small><small>''a&#x27;'' および ''b'' がゼロの場合、逆正接は不定値となる (すなわち対応する ''C&#x27;'' もゼロとなる)。この場合、色相角はゼロを代入する。{{Harvnb|Sharma|2005|loc=eqn. 7}}参照のこと。</small>

: <math>
  \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}
</math>
:: <small>'''注:''' ''C&#x2032;''<sub>1</sub> あるいは ''C&#x2032;''<sub>2</sub> がゼロの場合、&#x394;h&#x2032; は無関係となり、ゼロが代入される。 </small><small>{{Harvnb|Sharma|2005|loc=eqn. 10}}参照のこと。</small>
<div class="cx-template-editor-source-container" lang="en" dir="ltr" style="display: none;"><div class="cx-template-editor-source"><div class="cx-template-editor-title">Harvnb</div></div></div>

: <math>
  \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}
</math>
:: <small>'''注:''' ''C&#x2032;''<sub>1</sub> あるいは ''C&#x2032;''<sub>2</sub> がゼロの場合、H&#x2032; は ''h&#x2032;''<sub>1</sub>+''h&#x2032;''<sub>2</sub> となる (2で除算しない; 本質的には、ひとつの角度が不定の場合、もう一方の角度を使って平均を取る。不定の角度がゼロであることを前提としている)。 </small><small>{{Harvnb|Sharma|2005|loc=eqn. 7 and p. 23}}参照のこと。インターネット上のほとんどの実装に、平均色相角の計算エラーがある、と述べている。</small>

:<math>
  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)
</math>

:<math>
  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
</math>

:<math>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]</math>

=== 許容値 ===
[[ファイル:CIExy1931_MacAdam.png|右|サムネイル|276x276ピクセル|[[CIE 1931 色空間|CIE 1931色空間]]におけるマクアダム図。表示されている楕円は最大10倍も実際のサイズと異なる。]]
'''許容値'''は、「基準色に対してどのような色が見分けがつかないほど近いか？」を定義する。もしも2色間の距離が'''知覚的均等'''であれば、答えは単純に「丁度可知差異以内のすべての色」である。 このことは色域内において知覚的均等であることを前提としている。 さもないと、許容値は参照値を含む関数となってしまい、複雑化してしまう。

例えば[[CIE 1931 色空間]]は、L* (明度)　が一定の場合の許容値を表す線が[[デビッド・マクアダム#マクアダムの楕円|マクアダムの楕円]]として定義されている。図で示すとおり、許容値を表す楕円のサイズは、場所によりバラバラであることがわかる。 このことにより、[[CIE Luv色空間|CIELUV]] と [[Lab色空間|CIELAB]]色空間の定義へとつながった。

より一般的に、明度を変動させると、許容値は[[楕円体]]となる。上記のように距離の表現における重み付け係数を増加させると、各々の係数の軸において楕円体の体積を増やす方向に働く。<ref name="xrite">{{Cite web|url=http://www.xrite.com/documents/literature/en/L10-024_Color_Tolerance_en.pdf|title=A guide to Understanding Color Tolerancing|accessdate=2014-12-02|author=Susan Hughes|date=14 January 1998}}</ref>

== 参照 ==
* [[Lab色空間|CIELAB]]

== 脚注 ==
{{脚注ヘルプ}}
<!--=== 注釈 ===
{{Reflist|group="注釈"}}-->
=== 出典 ===
{{reflist|30em}}

== 外部リンク ==
* [http://www.brucelindbloom.com/ColorDifferenceCalc.html Bruce Lindbloomの色差計算]
* [http://www.ece.rochester.edu/~gsharma/ciede2000/ CIE dE2000色差式]: Gaurav Sharma氏によるMATLABおよびExcelでの実装


{{DEFAULTSORT:しきさ}}
[[Category:色]]
[[Category:色覚]]