(57)【要約】【課題】  色度座標系における微小色領域の色差を計算できるようにすると共に、人間の目の知覚に合った目標色と再現色の色差をより高精度に評価できるようにする。【解決手段】  所定の主観評価によって得られた境界を示す制御色３２及び中心色を示す目標色３１を色度座標系の色度座標値に変換する３Ｄ−ＬＵＴ２０と、ここで変換された色度座標系で当該目標色３１の原点を中心色とし、当該中心色の周囲に分布する制御色３２の色度座標値を結んで色弁別用の楕円を作成する画像処理手段１８と、ここで作成された楕円のパラメータａ，ｂ，θを算出すると共に、目標色３１と再現色３３との色相差Ｈｕｅを算出し、楕円のパラメータａ，ｂ，θ及び色相差Ｈｕｅに基づいて重み係数を求め、かつ、目標色３１の色度座標値と再現色３３の色度座標値との色差を求め、重み係数を目標色３１と再現色３３との色差に演算するＤＳＰ１９とを備えるものである。【選択図】    図１
2006237887.tif 000002
【特許請求の範囲】
【請求項１】
色度座標系における目標色と再現色との色差を評価する色評価支援装置において、
灰色背景に目標色及び制御色を配置して複数の被験者により色の差を目視観察させ、前記制御色を変化させたとき、前記目標色の周囲について、当該目標色及び制御色とが異なると判断される境界を求める処理を主観評価としたとき、
前記主観評価によって得られた前記境界を示す制御色及び中心色を示す目標色を色度座標系の色度座標値に変換する色変換手段と、
前記色変換手段により変換された色度座標系で当該目標色の原点を中心色とし、当該中心色の周囲に分布する制御色の色度座標値を結んで色弁別用の楕円を作成する画像処理手段と、
前記画像処理手段によって作成された前記楕円のパラメータを算出すると共に、前記目標色と再現色との色相差を算出し、前記楕円のパラメータ及び色相差に基づいて前記重み係数を求め、かつ、前記目標色の色度座標値と再現色の色度座標値との色差を求め、前記重み係数を前記目標色と再現色との色差に演算する演算手段とを備えることを特徴とする色評価支援装置。
【請求項２】
前記演算手段は、
前記画像処理手段によって作成された前記楕円の外郭線を色弁別境界線としたとき、前記中心色の原点から前記色弁別境界線に至る距離を算出すると共に当該距離の逆数を算出して重み係数を求めることを特徴とする請求項１に記載の色差評価支援装置。
【請求項３】
前記色度座標系が三次元色空間座標系である場合であって、
前記演算手段は、
前記三次元色空間座標系の各格子点毎に楕円パラメータを演算することを特徴とする請求項１に記載の色差評価支援装置。
【請求項４】
前記演算手段は、
前記各格子点毎に演算された楕円パラメータを補間処理して格子間に位置する楕円のパラメータを算出することを特徴とする請求項１に記載の色評価支援装置。
【請求項５】
前記演算手段は、
前記三次元色空間座標系において、４つの格子点を結んだ三角錐内に前記目標位置が位置するように、当該目標位置から最も近い４つの格子点を抽出し、
抽出された４つの前記格子点の楕円のパラメータを補間処理して目標位置の楕円パラメータを算出することを特徴とする請求項４に記載の色評価支援装置。
【請求項６】
色度座標系における目標色と再現色との色差を評価する色評価方法において、  灰色背景に目標色及び制御色を配置して複数の被験者により色の差を目視観察させ、前記制御色を変化させたとき、前記目標色の周囲について、当該目標色及び制御色とが異なると判断される境界を求める処理を主観評価としたとき、
前記主観評価によって得られた前記境界を示す制御色及び中心色を示す目標色を色度座標系の色度座標値に変換する工程と、
変換された前記色度座標系で当該目標色の原点を中心色とし、当該中心色の周囲に分布する制御色の色度座標値を結んで色弁別用の楕円を作成する工程と、
作成された前記楕円のパラメータを算出すると共に、前記目標色と再現色との色相差を算出し、前記楕円のパラメータ及び色相差に基づいて前記重み係数を求め、かつ、前記目標色の色度座標値と再現色の色度座標値との色差を求め、前記重み係数を前記目標色と再現色との色差に演算する工程とを有することを特徴とする色評価方法。
【請求項７】
作成された前記楕円の外郭線を色弁別境界線としたとき、前記中心色の原点から前記色弁別境界線に至る距離を算出すると共に当該距離の逆数を算出して重み係数を求める工程を有することを特徴とする請求項６に記載の色差評価方法。
【請求項８】
前記色度座標系が三次元色空間座標系である場合であって、
前記三次元色空間座標系の各格子点毎に楕円パラメータを演算する工程を有することを特徴とする請求項６に記載の色差評価方法。
【請求項９】
前記各格子点毎に演算された楕円パラメータを補間処理して格子間に位置する楕円のパラメータを算出する工程を有することを特徴とする請求項６に記載の色評価方法。
【請求項１０】
前記三次元色空間座標系において、４つの格子点を結んだ三角錐内に前記目標位置が位置するように、当該目標位置から最も近い４つの格子点を抽出し、
抽出された４つの前記格子点の楕円のパラメータを補間処理して目標位置の楕円パラメータを算出する工程を有することを特徴とする請求項９に記載の色評価方法。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、カラープリンタやカラー複写機、これらのカラー複合機等における色再現性に関して、再現色を目標色に色合わせ調整するカラーマネージメントシステム等に適用して好適な色評価支援装置及び色評価方法に関するものである。
【背景技術】
【０００２】
従来、カラープリンタやカラー複写機、これらのカラー複合機等における色再現性に関して、再現色を目標色に色合わせ調整する場合に、ＪＩＳＺ８７２３に基づく、Ｄ５０や、Ｄ６５といった標準の光又は補助標準の光が用いられ、製品の色再現スペックとして、主に標準光Ｄ５０に色合わせ処理を行っている。しかしながら、同一の光源下で色合わせ処理を行ったとしても、プリントサンプルを観察する角度を変えたり、他の光源に変えたりすると、ターゲット色票との色が異なって見えてしまう場合がある。
【０００３】
この色が異なって見えるという度合いを色差により評価する手法としては、非特許文献１であるＪＩＳＺ８７１９を利用することが提案されている。色が異なって見える原因は、ターゲット色票の分光特性と、プリントサンプルの分光特性が異なるためであって、２つの分光特性において、ある照明光と分光被視感度とが掛け合わされて積分された値を比較したとき、その積分値が等しいと色が一致するが、観察する角度が変わった場合や、光源が変わった場合に、同様に計算された積分値が異なってしまうためである。
【０００４】
この種の色差評価方法においては、色差計算式が使用されるが、例えば、国際照明学会（以下ＣＩＥという）Ｌａｂ１９７６色差式を用いた場合に、高クロマ領域での色差が大きな値となる欠点を有している。この欠点を修正した非特許文献２には、クロマ量の違いが人間の目で感じる大きさに近づくように色差を修正している。色差は、色の知覚的な相違を定量的に表したものであり、Ｌ*ａ*ｂ*色度図で距離として求められる。距離として求めた色差は、知覚的には均等にならない。
【０００５】
非特許文献３では、色度座標上で色差を求める場合に色領域毎に色弁別域を定めている。この色弁別域において、所定の色（中心色）に対して多数の方向で加法混色による等色実験を行った結果、一定方向の色に対して、毎回とも同じ色にならず、視覚が識別できる範囲で変動した。この多数の方向に対して、この視覚変動の標準偏差をｘｙ色度図にプロットすると、中心色の色度座標の周囲に楕円が得られるというものである。
【０００６】
一般に、ある２つのプリントサンプルが、ある光源下では同じ色に見える（Ｌａ*ｂ*値が等しくなること）が、観察する位置を変えたり、他の光源下では異なった色に見えること（Ｌ*ａ*ｂ*値が異なること）を照明光条件等色、又は、光源に依存して色が変化するのでその現象を光源依存と呼ばれることが多い。この現象を２つのプリントサンプルと光源の分光特性からＬａｂ色空間におけるＬ*，ａ*，ｂ*値を計算する流れと共に説明する。
【０００７】
ここに１つの光源Ｐ１（λ）と１つのプリントサンプル（有色試料）Ｒ１（λ）が存在するとき、人間の目に入るときの分光特性は、Ｐ１（λ）×Ｒ１（λ）で表わされる。このとき、λは３８０〜７３０ｎｍまで可視域の波長である。人間の目の錐体細胞の分光特性は、赤、緑、青の色それぞれがｘ（λ），ｙ（λ），ｚ（λ）で示されるので、各錐体へのエネルギーは、Ｐ１（λ）×Ｒ１（λ）×ｘ（λ）、Ｐ１（λ）×Ｒ１（λ）×ｙ（λ）、Ｐ１（λ）×Ｒ１（λ）×ｚ（λ）となる。
【０００８】
一方、測色に関して、三刺激値Ｘ’Ｙ’Ｚ’は、上述した波長３８０〜７３０ｎｍまでの積分値で表され、その値は（１）式により求められる。
【０００９】
【数１】
2006237887.tif 000003
【００１０】
ここで係数ｋは、Ｙ’値の最大値を１．０に正規化するための補正係数である。
（１）式により、三刺激値Ｘ’Ｙ’Ｚ’が求まると、Ｌａｂ色空間におけるＬ＊，ａ＊，ｂ＊値（単にＬａｂ値という）は、公式（ＪＩＳ２０６８）通りに計算することができる。
【００１１】
これを２つの試料Ｒ２，Ｒ２に対し、まず、ある光源Ｐ１で測色してＬａｂ値を求める。試料Ｒ１とＲ２の分光特性は異なるが、Ｌ*ａ*ｂ*色空間座標系における色度ａ*ｂ*値は同じである。つまり、２つの試料Ｒ２，Ｒ２は条件等色している。これは、試料Ｒ１とＲ２の分光特性が異なっても、色度ａ*ｂ*値が同じになるのは、三刺激値Ｘ’Ｙ’Ｚ’が積分値であるため、分光特性が異なっても積分値が同じになるためである。
【００１２】
次に、同じ試料Ｒ１，Ｒ２に対して、別の光源Ｐ２で測色してＬａｂ値を求める。この測色の結果、試料Ｒ１の色度は、例えば、光源Ｐ１の時とほぼ同じになる。これは試料Ｒ１の分光特性がフラットであるためである。逆に試料Ｒ２は色度が大きく変わって見える。Ｐ１光源下で同じ色に見えた試料Ｒ１，Ｒ２は、Ｐ２光源下では異なった色に見えている。つまり、上述した条件等色が崩れており、光源に依存した現象といえる。このような現象が存在するため、プリンタの色合わせ処理等を行う場合に、ある特定の光源で色合わせを行っても、他の光源では色が異なるという現象が発生するのである。
【００１３】
【非特許文献１】ＪＩＳＺ８７１９「条件等色指数−照明光条件等色の評価方法」
【非特許文献２】”Industrial Colour-Difference Evaluation，”CIE Technical Report 116（1995）
【非特許文献３】D.L.MacAdam ”Visual Sensitivities to Color Differences in daylight、”Journal of the Optical Society of America，32.247-274(1942)
【発明の開示】
【発明が解決しようとする課題】
【００１４】
ところで、従来方式のカラーマネージメントシステムによれば、カラープリンタやカラー複写機、これらのカラー複合機等における色再現性に関して、再現色を目標色に色合わせ調整する場合に、次のような問題がある。
【００１５】
ｉ．非特許文献１によるＪＩＳＺ８７１９に示されるように、特定の光源下で色合わせ処理を行っても、見る角度によって、また、他の光源下で色が変わってしまうことは仕方がないことであるが、その色が変わる程度を、微小な色領域の色差で表すことが困難になっている。
【００１６】
ii．非特許文献１によるＪＩＳＺ８７１９には、標準光Ｄ５０や、Ｄ６５等の光源を使用して、ターゲット色票とプリントサンプルとの色差を求める手法が示されているが、更に高精度な色再現性が必要とされる場合、ＪＩＳＺ８７１９による色差を求める方法が完全に人間の目の特性を計算できていないのが現状である。このことで、ＪＩＳＺ８７１９の色差による計算では十分な色再現性が必ずしも保証されない。
【００１７】
iii．非特許文献２では、ＣＩＥＬａｂ１９７６色差式を用いた場合に、高クロマ領域での色差が大きな値となるので、その欠点を修正して、クロマ量の違いが人間の目で感じる大きさに近づくように色差を修正している。非特許文献３では、色度座標上で色差を求める場合に色領域毎に色弁別域を定めている。この色弁別域において、中心色の色度座標の周囲に楕円が得られているが、この楕円を折り込んで色差を計算することについては何ら開示されておらず、昨今のカラー画像形成装置において、高精度な色再現を要求される中で、新たな色差計算式による手法が望まれている。
【００１８】
iv．因みに、ＣＩＥ１９９４色差式によれば、ＣＩＥＬａｂ１９７６色差式に比べて改良されているものの、色領域毎に微小色領域を設定していないばかりか、色弁別楕円を組み込んだ色差計算処理をしていないので、高精度な色再現性に対処できていないのが現状である。
【００１９】
なお、現在、提案されている色差式には、上述した非特許文献２及び３の他に
・  CIE  Lab1976
・  CIE  94    CIE-116
・  CIEDE2000  CIE-142：Technical Report，Improvement to Industrial Colour-difference evaluation,（2001）
・  CMC（1:1）  F.J.J.Clarket，et al：Modification to the JPC 79 colour-difference formula，J.Soc.Dyers Col．,100，pp.128-132（1984）
・MLAB  H.Komatsubara：Evaluation of color difference, AIC Color 85, 44(1985)
・HUNTER
・NC-IIIC K.Ikeda et al：Improvement of uniformity in colour space in term of color specification and colour difference evaluation, CIE proceedings 22nd, pp.49-50（1991)
等が挙げられるが、いずれの色差式も色差という視覚特性を計算式により示そうとしており、色差という視覚特性は改良されているが、数式では、視覚の非線形な特性を完全に吸収できずに、条件によっては(非線形な部分では）、見た目の色差とは一致しない場合が多い。
【００２０】
そこで、この発明は上述した課題を解決したものであって、色度座標系における微小色領域の色差を計算できるようにすると共に、人間の目の知覚に合った目標色と再現色の色差をより高精度に評価できるようにした色評価支援装置及び色評価方法を提供することを目的とする。
【課題を解決するための手段】
【００２１】
上記課題を解決するために、本発明に係る色評価支援装置は、色度座標系における目標色と再現色との色差を評価する色評価支援装置において、灰色背景に目標色及び制御色を配置して複数の被験者により色の差を目視観察させ、制御色を変化させたとき、目標色の周囲について、当該目標色及び制御色とが異なると判断される境界を求める処理を主観評価としたとき、この主観評価によって得られた境界を示す制御色及び中心色を示す目標色を色度座標系の色度座標値に変換する色変換手段と、この色変換手段により変換された色度座標系で当該目標色の原点を中心色とし、当該中心色の周囲に分布する制御色の色度座標値を結んで色弁別用の楕円を作成する画像処理手段と、この画像処理手段によって作成された楕円のパラメータを算出すると共に、目標色と再現色との色相差を算出し、楕円のパラメータ及び色相差に基づいて重み係数を求め、かつ、目標色の色度座標値と再現色の色度座標値との色差を求め、重み係数を目標色と再現色との色差に演算する演算手段とを備えることを特徴とするものである。
【００２２】
本発明に係る色差評価支援装置によれば、色度座標系における目標色と再現色との色差を評価する場合に、まず、主観評価では、灰色背景に目標色及び制御色を配置して複数の被験者により色の差が目視観察され、制御色を変化させたとき、目標色の周囲について、当該目標色及び制御色とが異なると判断される境界が求められる。これを前提にして、色変換手段は、主観評価によって得られた境界を示す制御色及び中心色を示す目標色を色度座標系の色度座標値に変換する。画像処理手段は、色変換手段により変換された色度座標系で当該目標色の原点を中心色とし、当該中心色の周囲に分布する制御色の色度座標値を結んで色弁別用の楕円を作成する。
【００２３】
演算手段は、画像処理手段によって作成された楕円のパラメータを算出すると共に、目標色と再現色との色相差を算出し、楕円のパラメータ及び色相差に基づいて重み係数を求める。例えば、画像処理手段によって作成された楕円の外郭線を色弁別境界線としたとき、中心色の原点から色弁別境界線に至る距離を算出すると共に当該距離の逆数を算出して重み係数を求める。更に、演算手段は、目標色の色度座標値と再現色の色度座標値との色差を求め、重み係数を目標色と再現色との色差に演算するようになされる。
【００２４】
従って、色度座標系における微小色領域の色差を計算できるようになるので、従来方式に比べて人間の目の知覚に合った目標色と再現色の色差をより高精度に評価できるようになる。
【００２５】
本発明に係る色評価方法は、色度座標系における目標色と再現色との色差を評価する方法において、灰色背景に目標色及び制御色を配置して複数の被験者により色の差を目視観察させ、制御色を変化させたとき、目標色の周囲について、当該目標色及び制御色とが異なると判断される境界を求める処理を主観評価としたとき、主観評価によって得られた境界を示す制御色及び中心色を示す目標色を色度座標系の色度座標値に変換する工程と、変換された色度座標系で当該目標色の原点を中心色とし、当該中心色の周囲に分布する制御色の色度座標値を結んで色弁別用の楕円を作成する工程と、作成された楕円のパラメータを算出すると共に、目標色と再現色との色相差を算出し、楕円のパラメータ及び色相差に基づいて重み係数を求め、かつ、目標色の色度座標値と再現色の色度座標値との色差を求め、重み係数を目標色と再現色との色差に演算する工程とを有することを特徴とするものである。
【００２６】
本発明に係る色差評価方法によれば、色度座標系における目標色と再現色との色差を評価する場合に、色度座標系における微小色領域の色差を計算できるようになるので、従来方式に比べて人間の目の知覚に合った目標色と再現色の色差をより高精度に評価できるようになる。これにより、再現色を目標色に色合わせ調整するカラーマネージメントシステム等に当該方法を十分に応用できるようになる。
【発明の効果】
【００２７】
本発明に係る色差評価支援装置及び色評価方法によれば、色度座標系における目標色と再現色との色差を評価する場合に、重み係数を目標色と再現色との色差に演算する演算手段を備え、当該演算手段は、色度座標系で目標色の原点を中心色とし、当該中心色の周囲に分布する制御色の色度座標値を結んで作成された色弁別用の楕円のパラメータを算出すると共に、当該目標色と再現色との色相差を算出し、楕円のパラメータ及び色相差に基づいて重み係数を求め、かつ、目標色の色度座標値と再現色の色度座標値との色差を求めるものである。
【００２８】
この構成によって、色度座標系における微小色領域の色差を計算できるようになるので、従来方式に比べて人間の目の知覚に合った目標色と再現色の色差をより高精度に評価できるようになる。従って、再現色を目標色に色合わせ処理するカラーマネージメントシステム等に当該装置を十分に応用できるようになる。
【発明を実施するための最良の形態】
【００２９】
以下、図面を参照しながら、この発明の実施例に係る色評価支援装置及び色評価方法について説明をする。  
図１は、本発明に係る実施例としての色評価支援装置１００の構成例を示す概念図である。  
図１に示す色評価支援装置１００は、主観評価によって得られる色度座標系における目標色と再現色との色差を評価する装置である。ここに主観評価とは、灰（グレー）色背景に目標色及び制御色を配置して複数の被験者により色の差を目視観察させ、制御色を変化させたとき、目標色の周囲について、当該目標色及び制御色とが異なると判断される境界を求める処理をいう。色評価支援装置１００には、デスクトップ型のパーソナルコンピュータ（以下パソコンという）や、ノート型のパソコンが使用される。
【００３０】
色評価支援装置１００は、システムバス１１、ＲＯＭ（Read Only Memory）１２、ワーク用のＲＡＭ（Random Access Memory）１３、入出力装置（以下Ｉ／Ｏインターフェース１４という）、ＣＰＵ（Central Processing Unit；中央処理ユニット）１５、入力ツール１６、モニタ１７、画像処理手段１８、デジタル信号処理装置（以下ＤＳＰ１９という）、３次元ルックアップテーブル（以下３Ｄ−ＬＵＴ２０という）及びハードディスク（以下ＨＤＤ２１という）を有して構成される。
【００３１】
システムバス１１にはＲＯＭ１２が接続され、当該装置全体を制御するためのシステム起動用のプログラムデータＤＰが格納される。ＲＡＭ１３には、システム立ち上げ時のプログラムデータＤＰや、色評価支援時の制御コマンド等を一時記憶するように使用される。ＣＰＵ１５は、電源がオンされると、ＲＯＭ１２からプログラムデータＤＰを読み出してＲＡＭ１３に展開してシステムを起動し、当該装置全体を制御するようになされる。
【００３２】
システムバス１１にはＩ／Ｏインターフェース１４が接続される。Ｉ／Ｏインターフェース１４には、入力ツール１６、モニタ１７及び測色装置２００が接続される。入力ツール１６は、キーボード１６ａ及びマウス１６ｂを備えて構成される。
【００３３】
キーボード１６ａやマウス１６ｂ等は、ＣＰＵ１５に対して、グレー色背景３０上の目標色３１及び制御色３２から楕円のパラメータや、再現色３３と目標色３１との距離差（以下色差ΔＥ*abという）を算出するような指示を入力する際に使用される。例えば、キーボード１６ａは、ユーザにより操作され、色差算出時において、その指示をＣＰＵ１５に操作データＤ３１として出力するように操作される。同様にして、マウス１６ｂはＣＰＵ１５に操作データＤ３２を出力するように操作される。以下Ｄ３１及びＤ３２を一括して操作データＤ３という。
【００３４】
モニタ１７には、Ｌ*ａ*ｂ*色座標系（三次元色空間座標系）を示す３次元の格子映像や、Ｘ軸を色座標ａ*とし、Ｙ軸を色座標ｂ*とする二次元色座標で楕円映像等を表示データＤ２に基づいて表示するようになされる。モニタ１７には液晶表示パネルやブラウン管ディスプレイが使用される。その他に、プラズマディスプレイ（ＰＤＰ）やエレクトロルミネセンス（ＥＬ）素子を利用した平面表示装置が使用される。
【００３５】
この例で、Ｉ／Ｏインターフェース１４には測色装置２００が接続され、所定の光源下に配置されたグレー色背景３０上の目標色３１及び制御色３２を測色して測色データＤ１を当該色評価支援装置内の画像処理手段１８に出力する。測色装置２００には分光放射輝度計が使用される。画像処理手段１８は測色装置２００から測色データＤ１を入力して画像処理を実行する。例えば、画像処理手段１８は赤（Ｒ）、緑（Ｇ）、青（Ｂ）の３色の測色データＤ１から制御色及び目標色に係る輝度や色度情報を分離する。
【００３６】
画像処理手段１８には、システムバス１１を通じて、例えば、３Ｄ−ＬＵＴ２０が接続される。３Ｄ−ＬＵＴ２０は、色変換手段を構成し、主観評価によって得られた境界を示す制御色及び中心色を示す目標色に係る色度情報を、例えば、Ｌ*ａ*ｂ*色座標系の色度座標値（Ｌ*ａ*ｂ*値；以下Ｌａｂ値ともいう）に変換するようになされる。３Ｄ−ＬＵＴ２０には、ＥＥＰＲＯＭ等の不揮発メモリが使用される。
【００３７】
色変換手段は、Ｌ*ａ*ｂ*色座標系のルックテーブルの他に、Ｘ，Ｙ，Ｚ空間座標系、Ｙｘｙ空間座標系、Ｌ*ｕ*ｖ*空間座標系でもよく、また、色変換手段には、三次元色空間座標系のルックアップテーブルに限られることはなく、二次元座標系や、一次元系の色変換テーブル機能を有するものであってもよい。二次元座標系には、ａ＊ｂ＊色座標テーブル、ｘ，ｙ色座標テーブル、ｕ＊，ｖ＊色座標テーブルなどが使用できる。１次元系の場合には、Ｌ＊系テーブル、Ｙ系テーブルが使用できる。このように、ルックアップテーブルは、一次元以上であればよい。
【００３８】
画像処理手段１８は、３Ｄ−ＬＵＴ２０により変換された色度座標系で当該目標色の原点を中心色とし、当該中心色の周囲に分布する制御色の色度座標値を結んで色弁別用の楕円を作成する。画像処理手段１８には、例えば、システムバス１１を通じてＤＳＰ１９が接続される。ＤＳＰ１９は演算手段の一例であり、ＣＰＵ１５の制御を受けてグレー色背景３０上の目標色３１及び制御色３２に関して、画像処理手段１８によって作成された楕円のパラメータａ，ｂ，θを算出すると共に、目標色と再現色との色相差Ｈｕｅを算出する。演算途中のＬａｂ値は、ＲＡＭ１３やＨＤＤ２１等に格納される。
【００３９】
ＤＳＰ１９は、色評価時等において、楕円のパラメータａ，ｂ，θ及び色相差Ｈｕｅに基づいて重み係数を求め、かつ、目標色の色度座標値と再現色の色度座標値との色差を求め、重み係数を目標色と再現色との色差に演算するようになされる。例えば、ＤＳＰ１９は、画像処理手段１８によって作成された楕円の外郭線を色弁別境界線としたとき、中心色の原点から色弁別境界線に至る距離を算出すると共に当該距離の逆数を算出して重み係数を求める。
【００４０】
ＤＳＰ１９は、色度座標系がＬ*ａ*ｂ*色座標系である場合であって、例えば、Ｌ*ａ*ｂ*色座標系の各格子点毎に楕円パラメータａ，ｂ，θを演算する。この例で、ＤＳＰ１９は、各格子点毎に演算された楕円パラメータａ，ｂ，θを補間処理して格子間に位置する楕円のパラメータａ，ｂ，θを算出する。このとき、ＤＳＰ１９は、Ｌ*ａ*ｂ*色座標系において、４つの格子点を結んだ三角錐内に目標位置が位置するように、当該目標位置から最も近い４つの格子点を抽出し、ここに抽出された４つの格子点の楕円のパラメータａ，ｂ，θを補間処理して目標位置の楕円パラメータａ，ｂ，θを算出する。
【００４１】
ＣＰＵ１５は、キーボード１６ａ及びマウス１６ｂ等の操作指示に従って、例えば、ＤＳＰ１９の演算結果を受信し、再現色３３と目標色３１との色差ΔＥ*ab算出時の重み係数を決定するように操作される。
【００４２】
図２は、グレー色背景３０上の目標色３１及び制御色３２から成る主観対象例を示す構成図である。  
図２に示す主観対象例によれば、色度座標系における目標色と再現色との色差を評価する場合に、グレー色背景３０上に目標色３１及び制御色３２が配置されて構成される。目標色３１は、色度座標系において、楕円のパラメータａ，ｂ，θを求めたい色である。制御色３２は、自由に変えられる色である。グレー色背景３０は画用紙等のアートボードでよい。目標色３１及び制御色３２は、モニタや、色票（色紙）などで制御できればよい。
【００４３】
このグレー色背景３０上に配置された目標色３１及び制御色３２は、複数の被験者により色の差が目視観察される。この主観評価例で、制御色３２を変化させたとき、目標色３１の周囲について、当該目標色３１及び制御色３２とが異なると判断される境界を求めるようになされる。
【００４４】
図３は、色合わせ処理時の目標色３１及び再現色３３との比較例を示す図である。図３に示す比較例によれば、色評価の一例となる色合わせ処理時に、目標色３１と再現色３３とが比較される。この色合わせ処理では、目標色３１と再現色３３とが等しくなるように、例えば、再現色３３を出力したカラー画像形成系のイエロー（Ｙ）色、マゼンタ（Ｍ）色、シアン（Ｃ）色、ブラック（ＢＫ）色画像データが補正される。目標色３１は、例えば、カラー原稿の色サンプルである。
【００４５】
この種の色合わせ処理では、カラー原稿をスキャナで読み込んでＲＧＢ画像データを取得し、これをＹ色、Ｍ色、Ｃ色、ＢＫ色画像データに色変換し、この色変換後のＹ色、Ｍ色、Ｃ色、ＢＫ色画像データに基づいてカラー画像形成系からカラーサンプル画像が出力される。カラーサンプル画像には目標色３１を再現した再現色３３が形成される。このカラーサンプル画像を評価して色合わせ調整を行う場合に、例えば、Ｌ*ａ*ｂ*色座標系において、再現色３３と目標色３１との距離差（以下色差ΔＥ*abという）を算出し、等色か否かを判別するようになされる。
【００４６】
この例では、上述の判別の際に、３Ｄ−ＬＵＴ２０により変換されたＬ*ａ*ｂ*色座標系で当該目標色３１の原点を中心色とし、当該中心色の周囲に分布する制御色３２の色度座標値を結んで色弁別用の楕円を作成する。この楕円のパラメータａ，ｂ，θを算出すると共に、目標色３１と再現色３３との色相差Ｈｕｅを算出する。この楕円のパラメータａ，ｂ，θ及び色相差Ｈｕｅに基づいて重み係数を求め、かつ、目標色３１の色度座標値と再現色３３の色度座標値との色差を求め、重み係数を目標色３１と再現色３３との色差に演算するようになされる。
【００４７】
図４Ａ，Ｂ及び図５は、Ｌ*ａ*ｂ*色座標系における楕円の作成例（その１、２）を示す図である。この例では、その説明を分かりやすくするため、紙面と直交する軸をＬ＊軸としたとき、当該紙面に投影されたａ*ｂ*色座標系における弁別楕円（以下単に楕円という）の作成例について説明する。
【００４８】
図４Ａにおいて、実線に示す横軸はＸ座標軸に相当する、ａ*ｂ*色座標系における色度ａ*軸である。同様に示す縦軸はＹ座標軸に相当する、ａ*ｂ*色座標系における色度ｂ*軸である。また、波線に示す横軸は、ある格子点を原点（０）とする、ａ*ｂ*色座標系における色度ａ*軸である。同様に示す縦軸は、そのａ*ｂ*色座標系における色度ｂ*軸である。格子点の原点には中心色である目標色３１が設定される。
【００４９】
このａ*ｂ*色座標系において、弁別楕円を求めるには、まず、図４Ａに示す制御色３２を色相０°の方向に変化させる。次に被験者は、図２に示した目標色３１と制御色３２を見比べて、制御色３２を変化させて、これを目視評価し、制御色３２と目標色３１とが異なると判断する境界を何度も繰り返して求める。この境界は、複数の被験者の平均値を採る。この点を色相０°での境界点と設定する。
この境界点は、色相０°における目標色３１と異なると判断された制御色３２の色度座標値（ａ*ｂ*）によって示される。
【００５０】
この境界点を楕円外郭線の起点にして、色相を３０°に設定する。この設定は、目標色３１の周囲に円を描き、これに沿って制御色３２を３０°だけ方向をずらして観察する。この状態で、前述のように、被験者は、目標色３１と制御色３２を見比べて、制御色３２を変化させて、これを目視評価し、制御色３２と目標色３１とが異なると判断する境界を何度も繰り返して求める。この境界は、複数の被験者の平均値を採る。この点を色相３０°での境界点と設定する。この境界点は、色相３０°における目標色３１と異なると判断された制御色３２の色度座標値（ａ*ｂ*）によって示される。上述した処理を３０°刻みで１２回、すなわち、目標色（中心色）の原点０を中心にして３６０°だけ繰り返すと、図４Ｂに示すようなほぼ楕円を想像させるプロットが得られる。
【００５１】
次に、図４Ｂにおいて、図４Ａで得られたプロットに適合する楕円を求める。例えば、１２個の境界点を結ぶと、折れ線楕円が形成されるが、この折れ線楕円を既存の図形修正機能を利用して数値表現可能な楕円に修復する。図形修正機能は、楕円パラメータａ，ｂ，θが求まれば、どんなフィッティング処理手法でもよく、何ら限定されない。これにより、図５に示すａ*ｂ*色座標系において、ある格子点を原点（０）とする目標色３１と制御色３２とを判別するための楕円を作成することができる。
【００５２】
その後、図５に示す楕円のパラメータａ，ｂ，θを求める。パラメータａは、楕円の短軸方向の距離Ｘであり、ａ*ｂ*色座標系によって示される。パラメータｂは、楕円の長軸方向の距離Ｙであり、ａ*ｂ*色座標系によって示される。回転角θは、ａ*ｂ*色座標系のａ*色座標軸と、楕円の短軸とが成す角度であって、ａ*色座標軸を基準にした楕円の回転角度である。
【００５３】
このような目標色３１の楕円パラメータａ，ｂ，θは、３Ｄ−ＬＵＴ２０の格子点分だけ求めるようになされる。この例では、主観評価処理を繰り返すことで、目標色３１の楕円パラメータａ，ｂ，θを求めて３Ｄ−ＬＵＴ２０に格納する。この例では、Ｌ*ａ*ｂ*色座標系における楕円パラメータａ，ｂ，θを求める方法について説明したが、Ｌ*ａ*ｂ*色座標系以外、また２次元色座標系、１次元系であっても、同様にして主観評価処理は行えばよい。
【００５４】
図６は、Ｌ*ａ*ｂ*色座標系における楕円パラメータａ，ｂ，θの補間例を示す図である。  
この例で、目標色３１のＬ*ａ*ｂ*値が得られても、３Ｄ−ＬＵＴ２０の格子上にそのＬ*ａ*ｂ*値に相当する楕円パラメータａ，ｂ，θが存在しない場合がある。例えば、格子間にターゲットの目標色３１が設定され、そのＬ*ａ*ｂ*値に相当する楕円パラメータａ，ｂ，θが存在しない場合である。この場合は、補間処理を利用して、図６に示すように、楕円パラメータａ，ｂ，θを求めればよい。
【００５５】
ここで、格子間に設定されたターゲット目標色３１のＬａｂ値をＬ*tａ*tｂ*tとし、その目標色３１の楕円パラメータａ，ｂ，θをａx，ｂx，θxとする。次に、目標色３１に最も近くて、当該ターゲット目標色３１を含有する格子点４つを検索する。この４格子点で、これらの基準点となる第０番のＬａｂ値をＬ*0ａ*0ｂ*0とし、その楕円パラメータａ，ｂ，θをａ0，ｂ0，θ0とする。同様にして、第１のＬａｂ値をＬ*1ａ*1ｂ*1とし、楕円パラメータａ，ｂ，θをａ1，ｂ1，θ1とする。第２のＬａｂ値をＬ*2ａ*2ｂ*2とし、楕円パラメータａ，ｂ，θをａ2，ｂ2  ，θ2とする。第３のＬａｂ値をＬ*3ａ*3ｂ*3とし、楕円パラメータａ，ｂ，θをａ3  ，ｂ3，θ3とする。
【００５６】
この例で、ターゲット目標色３１を取り囲む４つの格子点の目標色３１のＬ*ａ*ｂ*値と、ターゲット目標色３１のＬ*ａ*ｂ*値との間には、定数ｋ１，ｋ２，ｋ３としたとき、（２）式の関係が成り立つ。すなわち、
【００５７】
【数２】
2006237887.tif 000004
【００５８】
により与えられる。
この（２）式を計算して、定数ｋ１，ｋ２，ｋ３を求める。この関係は、楕円パラメータａ，ｂ，θにも成り立つので、ターゲット目標色３１のＬ*tａ*tｂ*t値に対する楕円パラメータａx，ｂx，θxは、式（３）、すなわち、
【００５９】
【数３】
2006237887.tif 000005
【００６０】
により与えられる（三角錐補間処理方法）。
これにより、主観評価処理で得られなかった、格子間の目標色３１のＬ*tａ*tｂ*t値に対する楕円パラメータａx，ｂx，θxを求めることができる。この例では、Ｌａｂ値を例にして説明したが、他のＸ，Ｙ，Ｚ空間座標系、Ｙｘｙ空間座標系、Ｌ*ｕ*ｖ*空間座標系を利用する場合でも同様な補間処理方法を採ることができる。補間処理は、三角錐補間処理方法に限定されることはなく、他の補間処理方法でもよい。
【００６１】
続いて、本発明に係る色評価方法について、色評価支援装置１００の動作例を説明する。図７は、色評価支援装置１００における動作例を示すフローチャートである。
この例では、図４Ａに示したように、予めグレー色背景３１に目標色３１及び制御色３２を配置して複数の被験者により色の差を目視観察させ、制御色３２を変化させたとき、目標色３１の周囲について、当該目標色３１及び制御色３２とが異なると判断される境界を求める処理を主観評価としたとき、この主観評価によって得られた境界を示す制御色３２及び中心色を示す目標色３１をＬ*ａ*ｂ*色座標系のＬ*ａ*ｂ*値に変換する工程と、図４Ｂ及び図５に示したように、そのＬ*ａ*ｂ*色座標系で当該目標色３１の原点を中心色とし、当該中心色の周囲に分布する制御色３２のＬ*ａ*ｂ*値を結んで色弁別用の楕円を作成する工程と、ここで作成された楕円のパラメータａ，ｂ，θが算出され、この楕円パラメータａ，ｂ，θの値がＬ*ａ*ｂ*色座標系の各格子点毎に３Ｄ−ＬＵＴ２０に格納される工程が既に完了している場合であって、図３に示したＬ*ａ*ｂ*色座標系における目標色３１と再現色３３との色差を評価する場合を例に挙げる。
【００６２】
これらを動作条件にして、図７に示すフローチャートのステップＡ１で初期設定を実行する。ユーザは、キーボード１６ａ又はマウス１６ｂを使用して、目標色３１と再現色３３との色差計算の実行要否や、格子間の楕円パラメータａ，ｂ，θに関して補間処理の実行要否等を設定するようになされる。キーボード１６ａ又はマウス１６ｂを使用して得られた操作データＤ３は、Ｉ／Ｏインターフェース１４を通じてＣＰＵ１５に出力される。
【００６３】
次に、ステップＡ２に移行してＣＰＵ１５は、操作データＤ３に基づいてＬ*ａ*ｂ*色座標系に目標色３１と再現色３３とを設定して色差評価を開始する。このとき、比較対象となる図３に示した目標色３１と再現色３３とを測色装置２００を使用して色評価支援装置１００に取り込む。もちろん、これに限られることはなく、既に、色評価支援装置１００内に格納されている目標色３１と再現色３３に係る測色データＤ１を使用してもよい。そして、ステップＡ３で、ＣＰＵ１５は目標色３１と再現色３３との色差計算の実行要否に基づいて制御を分岐する。この例で、目標色３１と再現色３３との色差計算を実行要が設定されている場合は、ステップＡ４に移行して、ＣＰＵ１５は、目標色３１と制御色３２との色弁別境界を示す楕円パラメータａ，ｂ，θであるＸ座標軸ａ、Ｙ座標軸ｂ、回転角度θを目標色３１のＬ*ａ*ｂ*値に基づいて３Ｄ−ＬＵＴ２０から選択して読み出すようにＤＳＰ１９を制御する。この制御で目標色３１のＬ*ａ*ｂ*値と３Ｄ−ＬＵＴ２０からの楕円パラメータａ，ｂ，θとが決定される。
【００６４】
その後、ステップＡ５に移行してＣＰＵ１５は、補間処理の実行要否に基づいて制御を分岐する。補間処理の実行要が設定されている場合は、ステップＡ６に移行してＣＰＵ１５は補間処理を実行するようにＤＳＰ１９を制御する。このとき、ＤＳＰ１９は、各格子点毎に演算された楕円パラメータａ，ｂ，θを補間処理して格子間に位置する楕円のパラメータａ，ｂ，θを算出する。
【００６５】
この例で、ターゲット目標色３１を取り囲む４つの格子点の目標色３１のＬ*ａ*ｂ*値と、ターゲット目標色３１のＬ*ａ*ｂ*値との間には、先に説明した（２）式の関係が成り立つので、この（２）式を計算して、定数ｋ１，ｋ２，ｋ３が求められる。この関係は、楕円パラメータａ，ｂ，θにも成り立つので、（３）式に基づいて、ターゲット目標色３１のＬ*tａ*tｂ*t値に対する楕円パラメータａx，ｂx，θxが求められる。これにより、主観評価処理で得られなかった、格子間の目標色３１のＬ*tａ*tｂ*t値に対する楕円パラメータａx，ｂx，θxを求めることができる。
【００６６】
次に、ステップＡ７でＤＳＰ１９は、従来方式の色差計算方法に基づいて目標色３１と再現色３３との色差ΔＥoldを算出する。ΔＥoldは、目標色３１の色度座標値と再現色３３の色度座標値との色差（距離＃１）である。ここで、目標色３１のＬａｂ値をＬａｂ-togetとし、再現色３３のＬａｂ値をＬａｂ-reproductionとする。まず、従来方式と同様にして、色差ΔＥoldを（４）式、すなわち、
【００６７】
【数４】
2006237887.tif 000006
【００６８】
により求める。この（４）式に係る演算はＤＳＰ１９が実行する。この演算により、目標色３１と再現色３３との色差ΔＥoldを算出することができる。なお、補間処理の実行要が設定されていない場合（実行否の場合）は、ステップＡ７で最も近似した格子点のパラメータａ，ｂ，θを使用してΔＥoldを求めるようになる。
【００６９】
その後、ステップＡ８に移行してＤＳＰ１９は重み係数αを計算する。このとき、ＤＳＰ１９は、目標色３１のａ*ｂ*値と再現色３３のａ*ｂ*値の各差分Δ*ａ、Δ*ｂから色相Ｈueを算出する。例えば、ＤＳＰ１９は、目標色３１と再現色３３との色度差（短軸）Δａを（５）式、色度差（長軸）Δｂを（６）式、色相差（回転角）Ｈueを（７）式により求める。
Δａ＝ａ-reproduction−ａ-toget                    ・・・・・（５）
Δｂ＝ｂ-reproduction−ｂ-toget                    ・・・・・（６）
ΔＨue＝arctan（Δｂ／Δａ）                         ・・・・・（７）
【００７０】
ＤＳＰ１９は、色度差がいずれもΔａ＞０及びΔｂ＞０となる場合であって、色相差Ｈueが、０°＜ΔＨue＜９０°となる場合に分けて計算を実行する。
その後、ＤＳＰ１９は、楕円パラメータａ，ｂ，θ及び色相Ｈueに基づいて、
Ｌ*ａ*ｂ*色座標系における目標色３１の中心色の原点から色弁別境界線に至るまでの距離＃２を算出し、当該距離＃２に基づいて重み係数αを求める。重み係数αは、パラメータａ，ｂ，θから成る楕円の角度ΔＨueでの原点から距離＃２の逆数である。つまり、（８）式、すなわち、
重み係数α＝１／原点からの距離＃２                   ・・・・・（８）
によって与えられる。
【００７１】
ここで、回転角ΔＨue’における楕円の外郭線上で新たな評価点の色差ΔＥnewを求める場合、まず、（９）式より楕円の回転角の計算を実行する。
ΔＨue’＝ΔＨue−θ                         ・・・・・（９）
【００７２】
次に、楕円の公式により、原点から見た楕円の外郭線上の新たな評価点のａ*軸成分の距離Ｘ及びｂ*軸成分の距離Ｙの値をそれぞれ求める。この例で、ａ*軸のパラメータを色度値ａ*とし、回転角をΔＨue’としたとき、ａ*軸成分の距離Ｘは、（１０）式、すなわち、
距離Ｘ＝パラメータａ*cos（ΔＨue’×π／１８０）    ・・・・・（１０）
により与えられる。
【００７３】
また、ｂ*軸のパラメータを色度値ｂ*とし、回転角をΔＨue’としたとき、ｂ*軸成分の距離Ｙは、（１１）式、すなわち、
距離Ｙ＝パラメータｂ*sin（ΔＨue’×π／１８０）    ・・・・・（１１）
により与えられる。
【００７４】
これらの関係により、原点から新たな評価点に至る距離＃２は、（１２）式、すなわち、
【００７５】
【数５】
2006237887.tif 000007
【００７６】
により与えられる。これら（５）〜（１２）式に係る演算はＤＳＰ１９が実行する。これらの演算により重み係数αを求めることができる。
【００７７】
そして、ステップＡ９に移行して、ＤＳＰ１９は、目標色３１と再現色３３との色差ΔＥoldと、重み係数αを使用して新たな色差ΔＥnewを算出する。このとき、ＤＳＰ１９は、目標色３１と再現色３３との色差（距離＃１）に重み係数αを演算する。例えば、（１２）式を計算した値の逆数を重み係数αとしたとき、新たな色差ΔＥnewは、（１３）式、すなわち、
ΔＥnew＝ΔＥold×重み係数α                      ・・・・・（１３）
により与えられる。この（１３）式に係る演算をＤＳＰ１９が実行する。この（１３）式を以下で新色差式ともいう。これにより、距離＃１に、距離＃２の逆数を掛け合わせてものを目標色３１と再現色３３との色差ΔＥnewとすることができる。
【００７８】
その後、ステップＡ１０でＣＰＵ１５は、目標色３１と再現色３３との色差ΔＥnewをモニタ１７に表示したり、また、これらの色評価結果をＨＤＤ２１に保存するようになされる。その後、ステップＡ１１に移行して色評価処理を終了するか否かを判断する。例えば、ＣＰＵ１５は電源オフ情報を検出して色評価処理を終了する。電源オフ情報が検出されない場合は、ステップＡ１に戻って上述した処理を繰り返すようになされる。
【００７９】
このように、実施例としての色差評価支援装置及び色評価方法によれば、Ｌ*ａ*ｂ*色座標系における目標色３１と再現色３３との色差ΔＥnewを評価する場合に、３Ｄ−ＬＵＴ２０は、主観評価処理によって得られた色弁別境界を示す制御色３２及び中心色を示す目標色３１をＬ*ａ*ｂ*色座標系のＬ*ａ*ｂ*値に変換する。画像処理手段１８は、３Ｄ−ＬＵＴ２０により変換されたＬ*ａ*ｂ*色座標系で当該目標色３１の原点を中心色とし、当該中心色の周囲に分布する制御色３２のＬ*ａ*ｂ*値を結んで色弁別用の楕円を作成する。
【００８０】
ＤＳＰ１９は、画像処理手段１８によって作成された楕円のパラメータａ，ｂ，θを算出すると共に、目標色３１と再現色３３との色相差Ｈueを算出し、楕円のパラメータａ，ｂ，θ及び色相差Ｈueに基づいて重み係数αを求める。この例では、画像処理手段１８によって作成された楕円の中心色の原点から色弁別境界線に至る距離＃２を算出すると共に当該距離＃２の逆数を算出して重み係数αを求める。更に、ＤＳＰ１９は、目標色３１のＬ*ａ*ｂ*値と再現色３３のＬ*ａ*ｂ*値との色差ΔＥoldを求め、重み係数αを目標色３１と再現色３３との色差ΔＥoldに演算するようになされる。
【００８１】
従って、Ｌ*ａ*ｂ*色座標系における微小色領域の色差ΔＥnewを計算できるようになるので、従来方式に比べて人間の目の知覚に合った目標色３１と再現色３３の色差ΔＥnewをより高精度に評価できるようになる。これにより、再現色３３を目標色３１に色合わせ処理するカラーマネージメントシステム等に当該色差評価支援装置１００を十分に応用できるようになる。
【００８２】
この実施例で重み係数αは、距離＃２の値が大きくなるにつれて、小さくなるものであれば、どんな値でもよい。色差ΔＥnewは、重み係数αを色差ΔＥoldに乗算するものの他に、上記条件が満たされれば、色差ΔＥnewは、重み係数αと色差ΔＥoldとを加算（足し算）するものであってもよい。
【００８３】
なお、既に文献等に記載されているＬａｂ座標系での色差ΔＥoldを用いる場合には、主観評価処理を行わずに、従来方式の研究結果データをそのまま使用して演算してもよい。また、３Ｄ−ＬＵＴ２０から読み出す楕円パラメータａ，ｂ，θについては、目標色３１のＬ*ａ*ｂ*値に基づいて選択する場合について説明したが、これに限られることはなく、再現色３３のＬ*ａ*ｂ*値に基づいて選択してもよい。また、目標色３１と再現色３３の平均Ｌ*ａ*ｂ*値に基づいて楕円パラメータａ，ｂ，θを選択してもよい。
【００８４】
図８〜図１１は、新色差式に基づく色差計算例（その１〜４）を示すグラフである。図８Ａ〜図１１Ａは、Ｌ*ａ*ｂ*色座標系における再現色の円及び目標色の楕円を示している。図８Ａ〜図１１Ａにおいて、横軸は色度ａ*軸であり、縦軸は、色度ｂ*軸である。
【００８５】
図８Ｂ〜図１１Ｂは、色差又は距離対色相［°］の関係例を示すグラフである。図８Ｂ〜図１１Ｂにおいて、横軸は色相［°］であり、０°〜３６０°（３０°刻みの等分目盛り）である。縦軸は、色差又は距離であり、０〜４（０．５刻みの等分目盛り）である。これらのグラフは、表示データＤ２に基づいてモニタ１７に表示される。
【００８６】
［色差計算例１］
図８Ａに示す第１の色差計算例では、図３に示した目標色３１のＬ*ａ*ｂ*値が（Ｌ*，ａ*，ｂ*）＝（５０，０，０）とき、再現色３３のＬ*ａ*ｂ*値が（Ｌ*，ａ*，ｂ*）＝（５０，ａx，ｂx）となるときの、パラメータａx，ｂxを求める場合を例に挙げる。
【００８７】
この目標色３１のＬ*ａ*ｂ*値に対応する楕円パラメータａ，ｂ，θはａ＝０．５，ｂ＝２．０（ａ：ｂ＝１：４），θ＝２０°である場合を例に採る。なお、再現色３３の色度値ａ*，ｂ*は、図８Ａに示す原点０，０から距離「１」に相当する円とする。この場合、従来方式では、色差は常に「１」である。
【００８８】
本発明方式によれば、図８Ａに示す弁別楕円が採用され、新たな色差式として先に説明した（１３）式が採用されるので、弁別楕円の外郭線上が同じ色差になる。
【００８９】
図８Ｂに示す波線は、旧色差であり、色差は常に「１」である。細い実線は、旧色差であって、従来方式の計算式（ΔＥold）により得られる距離＃１であり、太い実線は、新色差であって、上述の条件で、すなわち、本発明方式の計算式（ΔＥnew）により得られる距離＃１×αである。図８Ｂに示すように、色相により距離＃１が変化すると、それに反比例して新色差が変化していることが分かる。
【００９０】
この例で新色差は、色相（回転角θ）＝２０°及び２００°で最大となり、色相＝１１０°及び２８０°で最小となっている。色相＝４５°、１７０°及び２２５°で色差＝０°を跨いでいることが分かる。このように、Ｌ*ａ*ｂ*色座標系における微小色領域の色差ΔＥnewを計算できるようになるので、従来方式に比べて人間の目の知覚に合った目標色３１と再現色３３の色差ΔＥnewをより高精度に評価できるようになる。
【００９１】
［色差計算例２］
図９Ａに示す第２の色差計算例では、第１の色差計算例と同様にして、再現色３３のＬ*ａ*ｂ*値が（Ｌ*，ａ*，ｂ*）＝（５０，ａx，ｂx）となるときの、パラメータａx，ｂxを求める場合であって、目標色３１のＬ*ａ*ｂ*値に対応する楕円パラメータａ，ｂ，θがａ＝０．８１６，ｂ＝１．４１４（ａ：ｂ＝１：２），θ＝１０°である場合を例に採る。なお、再現色３３の色度値ａ*，ｂ*は、図９Ａに示す原点０，０から距離「１」に相当する円とする。この場合、従来方式では、色差は常に「１」である。
【００９２】
本発明方式によれば、図９Ａに示す弁別楕円が採用され、新たな色差式として（１３）式が採用されるので、弁別楕円の外郭線上が同じ色差になる。  
図９Ｂに示す波線は、旧色差であり、色差は常に「１」である。細い実線は、旧色差であって、従来方式の計算式（ΔＥold）により得られる距離＃１であり、太い実線は、新色差であって、本発明方式の計算式（ΔＥnew）により得られる距離＃１×αである。図９Ｂに示すように、色相により距離＃１が変化すると、それに反比例して新色差が変化していることが分かる。
【００９３】
この例で新色差は、色相（回転角θ）＝１０°及び１９０°で最大となり、色相＝１００°及び２８０°で最小となっている。色相＝４５°、１５５°、２２５°及び３４０°で色差＝０°を跨いでいることが分かる。このように、第１の色差計算例と同様にして、Ｌ*ａ*ｂ*色座標系における微小色領域の色差ΔＥnewを計算できるようになるので、従来方式に比べて人間の目の知覚に合った目標色３１と再現色３３の色差ΔＥnewをより高精度に評価できるようになる。
【００９４】
［色差計算例３］
図１０Ａに示す第３の色差計算例では、第１及び第２の色差計算例と同様にして、再現色３３のＬ*ａ*ｂ*値が（Ｌ*，ａ*，ｂ*）＝（５０，ａx，ｂx）となるときの、パラメータａx，ｂxを求める場合であって、目標色３１のＬ*ａ*ｂ*値に対応する楕円パラメータａ，ｂ，θがａ＝０．８１６，ｂ＝１．２２５（ａ：ｂ＝２：３），θ＝１０°である場合を例に採る。なお、再現色３３の色度値ａ*，ｂ*は、図１０Ａに示す原点０，０から距離「１」に相当する円とする。この場合、従来方式では、色差は常に１である。
【００９５】
本発明方式によれば、図１０Ａに示す弁別楕円が採用され、新たな色差式として（１３）式が採用されるので、弁別楕円の外郭線上が同じ色差になる。  
図１０Ｂに示す波線は、旧色差であり、色差は常に「１」である。細い実線は、旧色差であって、従来方式の計算式（ΔＥold）により得られる距離＃１であり、太い実線は、新色差であって、本発明方式の計算式（ΔＥnew）により得られる距離＃１×αである。図１０Ｂに示すように、色相により距離＃１が変化すると、それに反比例して新色差が変化していることが分かる。
【００９６】
この例で新色差は、色相（回転角θ）＝１０°及び１９０°で最大となり、色相＝１００°及び２８５°で最小となっている。色相＝４５°、１５０°、２３０°及び３３０°で色差＝０°を跨いでいることが分かる。このように、第１及び第２の色差計算例と同様にして、Ｌ*ａ*ｂ*色座標系における微小色領域の色差ΔＥnewを計算できるようになるので、従来方式に比べて人間の目の知覚に合った目標色３１と再現色３３の色差ΔＥnewをより高精度に評価できるようになる。なお、第１〜第３の色差計算例から判るように、新色差と旧色差との差分が大きくなるほど、楕円パラメータａ，ｂ，θの比が大きいことも明確になった（第３の色差計算例の楕円に比べて、第１の色差計算例の楕円が扁平している）。
【００９７】
［色差計算例４］
図１１Ａに示す第４の色差計算例では、第１〜第３の色差計算例と同様にして、再現色３３のＬ*ａ*ｂ*値が（Ｌ*，ａ*，ｂ*）＝（５０，ａx，ｂx）となるときの、パラメータａx，ｂxを求める場合であって、目標色３１のＬ*ａ*ｂ*値に対応する楕円パラメータａ，ｂ，θがａ＝０．５，ｂ＝２．０（ａ：ｂ＝１：４），θ＝２０°である場合を例に採る。なお、再現色３３の色度値ａ*，ｂ*は、図１１Ａに示す原点０，０から距離「２」に相当する円とする。この場合、従来方式では、色差は常に２である。
【００９８】
本発明方式によれば、図１１Ａに示す弁別楕円が採用され、新たな色差式として（１３）式が採用されるので、弁別楕円の外郭線上が同じ色差になる。  
図１１Ｂに示す波線は、旧色差であり、色差は常に「２」である。細い実線は、旧色差であって、従来方式の計算式（ΔＥold）により得られる距離＃１であり、太い実線は、新色差であって、本発明方式の計算式（ΔＥnew）により得られる距離＃１×αである。図１１Ｂに示すように、色相により距離＃１が変化すると、それに反比例して新色差が変化していることが分かる。
【００９９】
この例で新色差は、色相＝１５°及び１９５°で最大となり、色相＝１０５°及び２８５°で最小となっている。色相＝４０°、１６０°、２２０°及び３５０°で色差＝０°を跨いでいることが分かる。このように、第１〜第３の色差計算例と同様にして、Ｌ*ａ*ｂ*色座標系における微小色領域の色差ΔＥnewを計算できるようになるので、従来方式に比べて人間の目の知覚に合った目標色３１と再現色３３の色差ΔＥnewをより高精度に評価できるようになる。
【産業上の利用可能性】
【０１００】
この発明は、カラープリンタやカラー複写機、これらのカラー複合機等における再現色を目標色に色合わせ調整するカラーマネージメントシステムに適用して極めて好適である。
【図面の簡単な説明】
【０１０１】
【図１】本発明に係る実施例としての色評価支援装置１００の構成例を示す概念図である。
【図２】グレー色背景上の目標色及び制御色から成る主観対象例を示す構成図である。
【図３】色合わせ処理時の目標色及び再現色との比較例を示す図である。
【図４】（Ａ），（Ｂ）は、Ｌ*ａ*ｂ*色座標系における楕円の作成例（その１）を示す図である。
【図５】Ｌ*ａ*ｂ*色座標系における楕円の作成例（その２）を示す図である。
【図６】Ｌ*ａ*ｂ*色座標系における楕円パラメータａ，ｂ，θの補間例を示す図である。
【図７】色評価支援装置１００における動作例を示すフローチャートである。
【図８】（Ａ）及び（Ｂ）は、新色差式に基づく色差計算例（その１）を示すグラフである。
【図９】（Ａ）及び（Ｂ）は、新色差式に基づく色差計算例（その２）を示すグラフである。
【図１０】（Ａ）及び（Ｂ）は、新色差式に基づく色差計算例（その３）を示すグラフである。
【図１１】（Ａ）及び（Ｂ）は、新色差式に基づく色差計算例（その４）を示すグラフである。
【符号の説明】
【０１０２】
１１  システムバス
１２  ＲＯＭ
１３  ＲＡＭ
１５  ＣＰＵ（制御手段）
１６  入力ツール
１７  モニタ
１８  画像処理手段
１９  ＤＳＰ（制御手段）
２０  ３Ｄ−ＬＵＴ（不揮発メモリ）
３０  グレー色背景
３１  目標色
３２  制御色
３３  再現色
１００  色評価支援装置
２００  測色装置
【図１】
2006237887.tif 000008
【図２】
2006237887.tif 000009
【図３】
2006237887.tif 000010
【図４】
2006237887.tif 000011
【図５】
2006237887.tif 000012
【図６】
2006237887.tif 000013
【図７】
2006237887.tif 000014
【図８】
2006237887.tif 000015
【図９】
2006237887.tif 000016
【図１０】
2006237887.tif 000017
【図１１】
2006237887.tif 000018