(57)【特許請求の範囲】
【請求項１】
入力信号値に対応した３次元の格子点を備え、各格子点には入力信号値に対応する出力信号値が記憶されている、色変換用の3次元ルックアップテーブルの作成方法であって、
前記3次元ルックアップテーブルに、
（１）入力信号値に対応した3次元の格子点を備え、各格子点には入力信号値に対する表示器の理想とする表示色が記憶されている３次元理想特性を求める過程と、
（２）入力信号値に対応した3次元の格子点を備え、各格子点には入力信号値に対する表示器の測定された実際の表示色が記憶されている３次元実測特性を求める過程と、
（３）3次元ルックアップテーブルの格子点群のうち、出力信号値を求めようとしている目標点に係る理想特性値が、３次元実測特性に係る実測特性値で取り囲まれているかどうかを判定する過程と、
（４）目標点に係る理想特性値が、３次元実測特性に係る実測特性値で取り囲まれている場合には、目標点を取り囲んでいる３次元実測特性の各格子点の実測特性値を用いて、目標点の出力信号値を３次元補間により求める過程と、
（５）前記（４）の処理で求められた各格子点の出力信号値からノイズ成分を除去する過程
を経て求められた、入力信号値に対応する出力信号値を記憶させることを特徴とする３次元ルックアップテーブルの作成方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、３次元ルックアップテーブルの作成方法に係り、特に、表示装置に備えられる色変換用の３次元ルックアップテーブルを作成する技術に関する。
【背景技術】
【０００２】
従来、プリンタやスキャナー等のように、色空間がシアン（C）、マゼンタ（M）、黄（Y）、黒（K)である機器では、３次元ルックアップテーブルを用いて色域変換を行ってきた。これに対して、液晶表示装置などの表示装置は、その性質上、色再現性が高いため、何も操作せずに表示しているか、もしくは、色再現理論に従った方法（例えば、マトリクス方式）を用いてきた。しかし、色再現理論に従った方法では、表示装置の特性の良し悪しに性能が左右されてしまう。近時、表示装置に対してもより精確な色のコントロールが求められるようになってきた。色空間は、ＲＧＢ、Ｌ＊ａ＊ｂ＊、ＸＹＺ等の3次元表現であるので、精確性を追求すると、3次元ルックアップテーブルによる制御が必要になる。
【０００３】
具体的には、液晶表示装置などの表示装置は、装置固有の色特性（「色域」ともいう）を有しているので、例えば同じ赤（ｒ)、緑（ｇ）、青（ｂ)の入力信号を与えても、装置が異なると画面に映し出される色は同じにならないことが知られている。このような表示装置の色特定の違いをなくすために、色変換用の３次元ルックアップテーブルが用いられる。周知のように、色変換用の３次元ルックアップテーブルは、入力信号値（ｒ,ｇ,ｂ）と出力信号値（R,G,B）とが対応付けられた検索表であるが、模式的には、入力信号値（ｒ,ｇ,ｂ）は3次元の立方体格子を構成し、この立方体格子の各格子点に出力信号値（R,G,B）が設定されたものである。以下、３次元ルックアップテーブルの各格子点に設定される出力信号値（R,G,B）を３次元ルックアップテーブルの「パラメータ」ともいう。
【０００４】
いま、３次元ルックアップテーブルがそれぞれ８ビット（各色で２５６階調）であるとすると、表示可能な色（パラメータの個数）は約１，６７７万色になる。これら全ての色変換データを持つとなると、３次元ルックアップテーブルのデータ量が膨大になる。そこで、入力信号値（ｒ,ｇ,ｂ）の各色を例えば９段階のように、大まかにサンプリングしている。この場合、表示可能な色は７２９色になる。そして、他の色は補間によって求めるようにしている。
【０００５】
３次元ルックアップテーブルの出力信号値（パラメータ）は次のようにして決められている。
まず、表示装置の理想特性を決める。この理想特性は、入力信号値（ｒ,ｇ,ｂ）で表される立方体格子の各格子点に、目標となる表示色をＸＹＺ、あるいはＬａｂの表色系で表した値をそれぞれ設定したものである。
【０００６】
次に、表示装置の実測特性を求める。この実測特性は、入力信号値（ｒ,ｇ,ｂ）で表される立方体格子の各格子点に、実測された表示色をＸＹＺ、あるいはＬａｂの表色系で表した値をそれぞれ設定したものである。
【０００７】
理想特性の各格子点（以下、「目標点（Ｔ０）」ともいう）を取り囲む４つの実測点（Ｔ１、Ｔ２、Ｔ３、Ｔ４）の各データに基づいて、目標点（Ｔ０）の表色を与える出力信号値（R,G,B））を求める。いわゆる四面体補間法により出力信号値を決定する。このようにして３次元ルックアップテーブルの全ての出力信号値を求める（例えば、特許文献１参照）。
【特許文献１】特開２０００−２５１０５４号公報（段落００８６〜０１０１、図６）
【発明の開示】
【発明が解決しようとする課題】
【０００８】
しかしながら、このような構成を有する従来例の場合には、次のような問題がある。
すなわち、従来手法によると、３次元実測特性を求める過程でノイズ成分を含む場合があり、その結果、３次元ルックアップテーブルの精度が悪くなるという問題点がある。
【０００９】
なお、従来の四面体補間法による３次元ルックアップテーブルの出力信号値決定手法によると、４つの実測点で囲われない目標点があった場合に、その目標点に係る出力信号値を次のようにして求めている。以下、図１０を参照して説明する。
【００１０】
図１０において、符号Ｔ１〜Ｔ４は表示色が実測された実測点である。また、符号Ｔ０は、４つの実測点で囲われない目標点である。従来手法によると、目標点Ｔ０から近い２つの実測点Ｔ１とＴ２とを結ぶ線分に、目標点Ｔ０から垂線を下し、その交点Ｔ０´を求める。そして、実測点Ｔ１と交点Ｔ０´との距離と、格子点Ｔ２と交点Ｔ０´との距離に応じて、２つの実測点Ｔ１、Ｔ２の表示色を比例配分し、それを目標点Ｔ０の出力信号値と推定している。
【００１１】
このように、従来手法によると、４つの実測点で囲われない目標点があった場合に、その目標点に係る出力信号値としては、近傍２点の実測点の表示色を直線補間して決定しているので、四面体補間法により求められた他の格子点の出力信号値との整合性が悪く、結果として、３次元ルックアップテーブルの精度が悪いという問題点がある。
【００１２】
本発明は、このような事情に鑑みてなされたものであって、３次元ルックアップテーブルの出力信号値の精度を高めることを目的とする。
【課題を解決するための手段】
【００１３】
本発明は、このような目的を達成するために、次のような構成をとる。
すなわち、請求項1に記載の発明は、入力信号値に対応した３次元の格子点を備え、各格子点には入力信号値に対応する出力信号値が記憶されている、色変換用の3次元ルックアップテーブルの作成方法であって、
前記3次元ルックアップテーブルに、
（１）入力信号値に対応した3次元の格子点を備え、各格子点には入力信号値に対する表示器の理想とする表示色が記憶されている３次元理想特性を求める過程と、
（２）入力信号値に対応した3次元の格子点を備え、各格子点には入力信号値に対する表示器の測定された実際の表示色が記憶されている３次元実測特性を求める過程と、
（３）3次元ルックアップテーブルの格子点群のうち、出力信号値を求めようとしている目標点に係る理想特性値が、３次元実測特性に係る実測特性値で取り囲まれているかどうかを判定する過程と、
（４）目標点に係る理想特性値が、３次元実測特性に係る実測特性値で取り囲まれている場合には、目標点を取り囲んでいる３次元実測特性の各格子点の実測特性値を用いて、目標点の出力信号値を３次元補間により求める過程と、
（５）前記（４）の処理で求められた各格子点の出力信号値からノイズ成分を除去する過程
を経て求められた、入力信号値に対応する出力信号値を記憶させることを特徴とする。
【００１４】
請求項１記載の発明によると、３次元実測特性を求める過程で混入したノイズ成分を除去するので、３次元ルックアップテーブルの出力信号値の精度を一層高めることができる。
【発明の効果】
【００１５】
本発明によれば、３次元実測特性を求める過程で含まれるノイズ成分を除去するので、３次元ルックアップテーブルの精度を高めることができる。
【実施例１】
【００１６】
以下、図面を参照して本発明の第1実施例を説明する。
図１は、第1実施例に係る表示装置の概略構成を示したブロック図である。
本実施例に係る表示装置は液晶表示装置であって、入力側から順に配置された入力側ルックアップテーブル（ＬＵＴ）１と、３次元ルックアップテーブル（ＬＵＴ）２と、出力側ルックアップテーブル（ＬＵＴ）３と、液晶表示器４とを備えている。
【００１７】
出力側ルックアップテーブル３は、液晶表示器４のγ特性をキャンセルして、そのγ特性をほぼ直線化（γ＝１）するためのものであって、入力信号値（R２，G２，B２）に応じて個別に設けられた１次元のルックアップテーブルである。一方、入力側ルックアップテーブル１は、出力側ルックアップテーブル３で液晶表示器４のγ特性がキャンセルされていることを前提として、表示装置のγ特性を任意に設定するためのものであって、入力信号値（R０，G０，B０）に応じて個別に設けられた１次元のルックアップテーブルである。
【００１８】
3次元ルックアップテーブル２は色変換用に用いられるもので、入力信号値（R１，G１，B１）に対応した３次元の格子点を備え、各格子点には入力信号値（R１，G１，B１）に対応する出力信号値（R２，G２，B２）が記憶されている。本実施例の特徴である3次元ルックアップテーブル２の出力信号値（パラメータ）の決定手法は後に詳しく説明する。
【００１９】
本実施例装置によると、液晶表示器４に応じたビット数、例えばそれぞれが８ビット又は１０ビットの入力信号値（R０，G０，B０）が入力側ルックアップテーブル１に与えられる。入力信号値（R０，G０，B０）は、入力側ルックアップテーブル１に予め記憶されている任意のγ特性によって変換されて信号値（R１，G１，B１）として出力される。この信号値（R１，G１，B１）が3次元ルックアップテーブル２に与えられる。3次元ルックアップテーブル２は、予め設定された色変換特性でもって入力信号値（R１，G１，B１）を変換して、信号値（R２，G２，B２）を出力する。この信号値（R２，G２，B２）は、出力側ルックアップテーブル３で、液晶表示器４のγ特性とは逆の特性で変換される。その出力信号値（R３，G３，B３）が液晶表示器４に与えられる。
【００２０】
以上のように本実施例装置によると、出力側ルックアップテーブル３で液晶表示器４のγ特性がキャンセルされるので、入力側ルックアップテーブル１に設定された任意のγ特性で表示を行うことができる。また、3次元ルックアップテーブル２によって任意の色変換を行うことができるので、液晶表示器４の固有の色特性に影響されないように色再現をすることが可能である。なお、入出力のルックアップテーブル１、３の機能を3次元ルックアップテーブル２にもたせることもできるが、そうすると3次元ルックアップテーブル２のデータ量が膨大になる。本実施例では、入出力のルックアップテーブル１、３を3次元ルックアップテーブル２とは個別に設けているので、3次元ルックアップテーブル２のデータ量を低く抑えることができる。
【００２１】
次に、3次元ルックアップテーブル２の出力信号値の決定手法を図２のフローチャートを参照して説明する。
【００２２】
ステップＳ１：理想特性の算出
まず、入力信号値（R，G，B）に対して、どのような色再現をしたいのかを表した表示装置の理想特性を算出する。理想特性は、入力信号値に対応した3次元の格子点を備え、各格子点には入力信号値に対する表示器の理想とする表示色が記憶されている。本実施例では、表示色をＸＹＺ表色系の三刺激値（Ｘ，Ｙ，Ｚ）で表しているが、Ｌａｂ系で表してもよい。以下に、理想特性を算出する具体例を説明する。
【００２３】
入力信号値（R，G，B）とＸＹＺ表色系の赤、緑、青の三刺激値（Ｘ，Ｙ，Ｚ）と間には、次式（１）で示されるような関係がある。
【数１】
0005410479.tif 000002
【００２４】
そこで、まず、（１）式中の変換行列Ｍを算出する。具体的には、「赤」の入力信号（Ｒ１＝２５５、Ｇ１＝０、Ｂ１＝０）を与えたときの、ＸＹＺ表色系の三刺激値（Ｘ，Ｙ，Ｚ）を示す第1関係式を得る。三刺激値（Ｘ，Ｙ，Ｚ）の値は、ＩＣＣプロファイルに記述されたデバイス特性（表示器４の特性）から知ることができる。なお、実際に測定した表示器４の特性から三刺激値（Ｘ，Ｙ，Ｚ）の値を決めてもよいし、あるいは、ｓＲＧＢやＮＴＳＣといった規格で定められている特性から三刺激値（Ｘ，Ｙ，Ｚ）の値を決めてもよい。「赤」の入力信号を与えたときと同様にして、「緑」の入力信号（Ｒ＝０、Ｇ＝２５５、Ｂ＝０）を与えたときの第2関係式、「青」の入力信号（Ｒ＝０、Ｇ＝０、Ｂ＝２５５）を与えたときの第３関係式を求める。前記第１〜第３の関係式に基づいて、変換行例Ｍの各要素（m00,m01,・・・m22）を決定する。
【００２５】
変換行列Ｍが算出されたら、その変換行列Mと、入力側ルックアップテーブル１に設定したγ特性（輝度特性）とを用いて、任意の入力信号（R０，G０，B０）で特定される格子点の三刺激値（Ｘ，Ｙ，Ｚ）、すなわち、理想とする表示色を算出する。以下に具体的な算出手順を説明する。
【００２６】
まず、次式（２）で示されるような任意の入力信号Ｖiを設定する。
【数２】
0005410479.tif 000003
【００２７】
次に、入力側ルックアップテーブル１に設定されたγ特性を考慮した信号R０ｉ´，
G０ｉ´，B０ｉ´を求める。Ｒ、Ｇ、Ｂのそれぞれのγ特性をγｒ、γｇ、γｂとすると、信号R０ｉ´，G０ｉ´，B０ｉ´は、次式（３）で求められる。
【数３】
0005410479.tif 000004
【００２８】
γ特性を考慮した信号R０ｉ´，G０ｉ´，B０ｉ´と、先に求めた変換行列Ｍとを用いて、理想特性の各格子点の三刺激値（Ｘ，Ｙ，Ｚ）を次式（４）により求める。上記の演算を全ての任意入力（格子点）について実行することにより、表示装置の理想特性を得る。
【数４】
0005410479.tif 000005
【００２９】
ステップＳ２：実測特性の取得
表示装置に任意入力信号（前記（２）式参照）を与え、そのときに液晶表示器４で表示される色を測定して、実測された三刺激値（Ｘi，Ｙi，Ｚi）を得る。全ての格子点について測定して得られる３次元実測特性は、入力信号値に対応した3次元の格子点を備え、各格子点には入力信号値に対する表示器の測定された実際の表示色（本実施例では、三刺激値（Ｘi，Ｙi，Ｚi））が記憶されている。
【００３０】
ステップＳ３、Ｓ４：四面体補間の可否の判定
３次元ルックアップテーブル３の格子点群のうち、出力信号値を求めようとしている目標点に係る理想特性値が、３次元実測特性に係る実測測定値で取り囲まれているかどうかを判定する。以下その手順を詳しく説明する。
【００３１】
まず、３次元ルックアップテーブル３の格子点群のうち、出力信号値を求めようとしている目標点に係る理想特性値の色空間内の座標をＰ０（以下、「目標点Ｐ０」と呼ぶ）とする。そして、３次元実測特性の各格子点に対応する実測点（この実測点も色空間内の座標で表される）のうちで、目標点Ｐ０に最も近い１つの実測点Ｐ１を見出す。
【００３２】
以下、図３を参照して説明する。見出した実測点Ｐ１を１つの頂点として、目標点Ｐ０を内包する６面体について、Ｐ１を頂点とする三角錐、すなわち四面体（Ｐ１，Ｐ２，Ｐ３，Ｐ４）を仮定する（図３中の符号Ａで示した四面体）。そして、平面の方程式を使って、目標点Ｐ０が四面体（Ｐ１，Ｐ２，Ｐ３，Ｐ４）に囲まれているかどうかを判定する。具体的には次のような手順で判定する。
【００３３】
まず、四面体（Ｐ１，Ｐ２，Ｐ３，Ｐ４）の全ての平面について方程式を求める。そして、目標点Ｐ０と各平面の位置関係を調べ、四つの平面に対して目標点Ｐ０が全て正方向に位置していると判定されれば、目標点Ｐ０は四面体（Ｐ１，Ｐ２，Ｐ３，Ｐ４）に囲まれていると判断する。
【００３４】
平面の方程式は次のようにして求めることができる。図４を参照する。例えば、実測点（Ｐ１，Ｐ２，Ｐ３）を例に採る。この平面の方程式を、
ａｘ＋ｂｙ＋ｃｚ＋ｄ＝０
で表すと、係数ａ、ｂ、ｃは、点Ｐ１からＰ２へ向かうベクトルＶ１と、点Ｐ１からＰ３へ向かうベクトルＶ２との外積で得られる（式（５）参照）。
【数５】
0005410479.tif 000006
【００３５】
ただし、ｄ＝−１×（ａ×Ｖ１ｘ+ｂ×Ｖ１ｙ+ｃ×Ｖ１ｚ）である。
点Ｐ４が平面（Ｐ１，Ｐ２，Ｐ３）の正の方向（ｓ≧０）となるようにする。
ここで、ｓ＝ａ×Ｐ４ｘ+ｂ×Ｐ４ｙ+ｃ×Ｐ４ｚ+ｄ  である。
求められた平面の方程式に対して点Ｐ４が負の方向に位置した場合は（ａ，ｂ，ｃ，ｄ）に「−１」を掛ける。
【００３６】
以上のような判定を行って、目標点Ｐ０が四面体（Ｐ１，Ｐ２，Ｐ３，Ｐ４）に囲まれていないと判定された場合は、実測点Ｐ１を平面（Ｐ２，Ｐ３，Ｐ４）を挟んで反対側の実測点Ｐ１´に移動し、目標点Ｐ０が四面体（Ｐ１´，Ｐ２，Ｐ３，Ｐ４）（図３中に符号Ｂで示す）に囲まれているかどうかを同様に判定する。
【００３７】
目標点Ｐ０が四面体（Ｐ１´，Ｐ２，Ｐ３，Ｐ４）に囲まれていな場合、目標点Ｐ０が平面（Ｐ１´，Ｐ２，Ｐ３）の外側（負の方向）にあれば、目標点Ｐ０は、四面体（Ｐ１´，Ｐ２，Ｐ３，Ｐ４´）（図３中に符号Ｃで示す）に囲まれていることになる。一方、目標点Ｐ０が平面（Ｐ１´，Ｐ２，Ｐ４）の外側にあれば、目標点Ｐ０は、四面体（Ｐ１´，Ｐ２，Ｐ３´，Ｐ４）（図３中に符号Ｄで示す）に囲まれていることになる。また、目標点Ｐ０が平面（Ｐ１´，Ｐ３，Ｐ４）の外側にあれば、目標点Ｐ０は、四面体（Ｐ１´，Ｐ２´，Ｐ３，Ｐ４）（図３中に符号Ｅで示す）に囲まれていることになる。
【００３８】
ステップＳ５:四面体補間によるパラメータ算出
以上のようにして、目標点Ｐ０を取り囲む四面体の存在が確認されると、次の手順で目標点Ｐ０に対応する出力信号値（R２，G２，B２）を四面体補間により求める。以下、図５を参照して説明する。いま、目標点Ｐ０が四面体（Ｐ１，Ｐ２，Ｐ３，Ｐ４）に囲まれているとする。この四面体（Ｐ１，Ｐ２，Ｐ３，Ｐ４）を、更に目標点Ｐ０を頂点とする４つの四面体ｖ１、ｖ２、ｖ３、ｖ４に分割する。そして、これらの四面体ｖ１〜ｖ４の各体積を用い、次式（６）に従って目標点の出力信号値（R２，G２，B２）を算出する。
ＲＧＢ（Ｐ０）＝（ＲＧＢ（Ｐ１）×Ｖ１+ＲＧＢ（Ｐ２）×Ｖ２+ＲＧＢ（Ｐ３）×Ｖ３+ＲＧＢ（Ｐ４）×Ｖ４）／Ｖ                                        ・・・（６）
式（６）において、ＲＧＢ（Ｐ０）は、３次元ルックアップテーブル２の目標点Ｐ０に対応する出力信号値（R２，G２，B２）である。また、ＲＧＢ（Ｐ１）〜ＲＧＢ（Ｐ４）は、３次元実測特性の実測点Ｐ１〜Ｐ４に対応する入力信号値（R，G，B）である。さらに、Ｖ１〜Ｖ４は、四面体ｖ１〜ｖ４の体積である。また、Ｖは、四面体（Ｐ１，Ｐ２，Ｐ３，Ｐ４）の体積である。
【００３９】
本実施例では、目標点に係る理想特性値が、３次元実測特性に係る実測特性値で取り囲まれている場合に、四面体補間によって目標点の出力信号値を求めているが、３次元補間としては、四面体補間の他に、スプライン補間、三角柱（プリズム）補間、六面体（立方体）補間、ピラミッド補間、体積補間などを用いることができる。
【００４０】
図２のフローチャートに戻る。ステップＳ４において、目標点Ｐ０を取り囲む四面体が存在しないと判定されると、ステップＳ６に進む。
【００４１】
ステップＳ６：出力信号値の推定
目標点に係る理想特性値が、３次元実測特性に係る実測特性値で取り囲まれていない場合には、３次元ルックアップテーブルの格子点群において、目標点の近傍にある第１格子点の出力信号値と、この第１格子点を通って目標点近傍に向かう仮想ライン上にある格子点のうち、少なくとも第１格子点と隣接している第２格子点の出力信号値とを用いて、目標点の出力信号値の候補を求める。そして、この処理を目標点の近傍にある複数個の格子点について繰り返し行って出力信号値の候補をそれぞれ求め、これらの出力信号値の候補の平均値を目標点の出力信号値とする。
以下、図６及び図７を参照して具体的に説明する。
【００４２】
図６において、符号Ａで示した領域の格子点群は、四面体補間によって出力信号値が求められた３次元ルックアップテーブル２の格子点群である。一方、符号Ｂで示した領域の格子点群は、これらの格子点に係る理想特性値が、３次元実測特性に係る実測特性値で取り囲まれていないので、四面体補間で出力信号値を求めることができなかった格子点群である。いま、領域Ｂ内の目標点Ｐ０の出力信号値を求める場合を説明する。まず、目標点Ｐ０の近傍にある第１格子点を抽出する。３次元格子を考えた場合に、目標点を取り囲む近傍格子は２６個ある。このうちステップＳ５の四面体補間で出力信号値が求められた格子点を第１格子点とする。
【００４３】
説明の簡単のために、図７に示すように、目標点Ｐ０の近傍に３つの第１格子点Ｐ２、Ｐ３、Ｐ４があったとする。以下、第１格子点Ｐ２に絞って説明する。この第１格子点Ｐ２を通って目標点近傍に向かう仮想ライン上にあって、第１格子点と隣接している第２格子点Ｐ６を抽出する。そして、第２格子点Ｐ６から第１格子点Ｐ２に向かうベクトルＶ（Ｐ２−Ｐ６）を想定する。このベクトルＶ（Ｐ２−Ｐ６）を平行移動してＰ２を基点とするベクトルＶ（Ｐ２−Ｐ６＋Ｐ２）を想定し、そのベクトルＶ（２×Ｐ２−Ｐ６）の先端に目標点の候補Ｐ０１があると推定する。そうすると目標点の第１候補Ｐ０１の出力信号値ＲＧＢ（Ｐ０１）は、格子点Ｐ２およびＰ６の各出力信号値ＲＧＢ（Ｐ２）、ＲＧＢ（Ｐ６）を使って、次式（７）で求めることができる。
ＲＧＢ（Ｐ０１）＝２×ＲＧＢ（Ｐ２）−ＲＧＢ（Ｐ６）    ・・・（７）
【００４４】
同様にして、第１格子点Ｐ３と、これに隣接する第２格子点Ｐ８とから、目標点の第２候補Ｐ０２の出力信号値ＲＧＢ（Ｐ０２）は、次式（８）で求めることができる。
ＲＧＢ（Ｐ０２）＝２×ＲＧＢ（Ｐ３）−ＲＧＢ（Ｐ８）    ・・・（８）
また、第１格子点Ｐ４と、これに隣接する第２格子点Ｐ１０とから、目標点の第３候補Ｐ０３の出力信号値ＲＧＢ（Ｐ０３）は、次式（９）で求めることができる。
-  ＲＧＢ（Ｐ０３）＝２×ＲＧＢ（Ｐ４）−ＲＧＢ（Ｐ１０）    ・・・（９）
図６及び図７で示した符号Ｐ０´は、目標点の第１〜第３候補Ｐ０１〜Ｐ０３の集合を示している。
【００４５】
以上のようにして求められた目標点の候補点Ｐ０１、Ｐ０２、Ｐ０３は、それぞれ格子点Ｐ６からＰ２に向かって生じている出力信号値の変化を候補点Ｐ０１に反映させており、同様に、格子点Ｐ８からＰ３に向かって生じている出力信号値の変化を候補点Ｐ０２に、格子点Ｐ１０からＰ４に向かって生じている出力信号値の変化を候補点Ｐ０３に、それぞれ反映させている。したがて、最終の目標点の出力信号値ＲＧＢ（Ｐ０）としては、次式（１０）に示すように、これらの候補点Ｐ０１、Ｐ０２、Ｐ０３、・・・の出力信号値の平均値を採れば良い。
ＲＧＢ（Ｐ０）＝（ＲＧＢ（Ｐ０１）＋ＲＧＢ（Ｐ０２）＋ＲＧＢ（Ｐ０３）＋・・・
＋ＲＧＢ（Ｐ０Ｎ））／Ｎ                    ・・・（１０）
ここで、Ｎは、第１格子点の個数である。
【００４６】
本実施例では、目標点の出力信号の候補を、第１格子点の出力信号値と、仮想ライン上で第１格子点と隣接する第２格子点の出力信号値とを使って求めているが、仮想ライン上で第２格子点に隣接する第３格子点の出力信号値を更に考慮して求めるようにしてもよい。
【００４７】
図２のフローチャートに戻る。ステップＳ６で目標点の出力信号値が推定されると、次のステップＳ７に進み、３次元ルックアップテーブル２の全ての格子点について出力信号値を算出したかどうかを判断し、全ての格子点について出力信号値が求まるまで、ステップＳ３〜Ｓ６を繰り返し実行する。
【００４８】
ステップＳ８：ノイズ成分の除去
全ての格子点について出力信号値が求まると、ステップＳ８に進んで出力信号値群の中に含まれるノイズを除去する。すなわち、３次元ルックアップテーブル２の格子点の出力信号値を求めるために、上述したステップＳ２で説明したように、液晶表示器４に表示される色を実測している。このような表示色の実測過程において、実測された３次元実測特性にノイズ成分が混入することがある。このノイズ成分が３次元ルックアップテーブル２の格子点の出力信号値に混入する結果、格子点の出力信号値の変化が不連続なものになる。格子点の出力信号値は、本来、連続性があると予測されるので、上記のノイズ成分をラプラシアンフィルタやメディアンフィルタを使って除去する。
【００４９】
以上のようにして３次元ルックアップテーブル２の全ての格子点の出力信号値を決定することができる。本実施例によれば、3次元ルックアップテーブルの格子点群のうち、出力信号値を求めようとしている目標点に係る理想特性値が、３次元実測特性に係る実測特性値で取り囲まれていない場合には、目標点の近傍にある第１格子点の出力信号値と、この第１格子点を通って目標点近傍に向かう仮想ライン上にあって、第１格子点と隣接している第２格子点の出力信号値とを用いて、目標点の出力信号値の候補を求め、上記の処理を目標点の近傍にある複数個の格子点について繰り返し行って出力信号値の候補をそれぞれ求め、これらの出力信号値の候補の平均値を目標点の出力信号値としているので、３次元ルックアップテーブルの各格子点の出力信号値を精度良く求めることができる。
【実施例２】
【００５０】
第２実施例は、表示装置のカラー較正用プログラムに関する。図８は、表示装置のカラー較正用プログラムを実行するコンピュータと、カラー較正される対象である表示装置の概略構成を示したブロック図である。
【００５１】
表示装置１０は、第１実施例で説明したものと同様の構成であるので、ここでの説明は省略する。コンピュータ２０は、中央演算処理装置（ＣＰＵ）２１と、データの読み出し専用メモリ（ＲＯＭ）２２と、データの一次記憶用のランダム・アクセス・メモリ（ＲＡＭ）２３を備え、これらの要素がバス２４を介して接続されている。また、コンパクトディスク（ＣＤ）ＲＯＭ２５を駆動するＣＤ−ＲＯＭドライバ２６がコントローラ２７を介してバス２４に接続されている。さらに、ハードディスク２８を駆動するハードディスクドライバ２９がコントローラ３０を介してバス２４に接続されている。操作機器としてのキーボード３１やマウス３２がコントローラ３３を介してバス２４に接続されている。そして、表示装置１０も、コントローラ３４を介してバス２４に接続されている。
【００５２】
以下、図９に示した表示装置のカラー較正用プログラムのフローチャートを参照して第２実施例の動作を説明する。
【００５３】
本実施例に係る表示装置のカラー較正用プログラムは、例えば記憶媒体としてのＣＤ−ＲＯＭ２５などで提供されてコンピュータ２０に取り込まれ、そのプログラムはハードディスク２８に格納される。なお、カラー較正用プログラムは、インターネットなどを用いた通信で提供されてもよい。カラー較正時、ＣＰＵ２１は、ハードディスク２８内のカラー較正用プログラムを読み出してその処理を実行する。
【００５４】
ステップＴ１：理想特性の取得
このステップＴ１では、入力信号値に対応した3次元の格子点を備え、各格子点には入力信号値に対する表示器の理想とする表示色が記憶されている３次元理想特性を取得する。具体的には、第１実施例のステップＳ１と同様である。ただし、変換行例Ｍを算出する際の基準となるＲＧＢ値や三刺激値（Ｘ，Ｙ，Ｚ）はオペレータがキーボード３１などを操作して入力する。これらの値が入力されるとＣＰＵ２１が（１）式に基づいて変換行例Ｍを算出する。さらに、オペレータがキーボード３１などを操作して、入力側ルックアップテーブル３１のγ特性を入力すると、ＣＰＵ２１は、先に算出した変換行例Ｍと、そのγ特性とを用いて、（４）式に従って理想特性の各格子点の三刺激値（Ｘ，Ｙ，Ｚ）を算出する。なお、理想特性の取得を上述した演算によらずに、別途作成した理想特性をコンピュータ２０に直接に取り込むようにしてもよい。
【００５５】
ステップＴ２：実測特性の取得
このステップＴ２では、入力信号値に対応した3次元の格子点を備え、各格子点には入力信号値に対する表示器の測定された実際の表示色が記憶されている３次元実測特性を取得する。具体的には、第１実施例のステップＳ２と同様である。ただし、実測された三刺激値（Ｘi，Ｙi，Ｚi）は、オペレータがキーボード３１などを操作して入力する。
【００５６】
ステップＴ３、Ｔ４：四面体補間の可否の判定
ここでは、ＣＰＵ２１が、3次元ルックアップテーブルの格子点群のうち、出力信号値を求めようとしている目標点に係る理想特性値が、３次元実測特性に係る実測特性値で取り囲まれているかどうかを判定する。具体的には、第1実施例のステップＳ３、Ｓ４と同様であるので、ここでの説明は省略する。
【００５７】
ステップＴ５：四面体補間によるパラメータ算出
ここでは、目標点に係る理想特性値が、３次元実測特性に係る実測特性値で取り囲まれている場合には、目標点を取り囲んでいる３次元実測特性の各格子点の実測特性値を用いて、目標点の出力信号値を３次元補間（例えば四面体補間）により求める。具体的には、第1実施例のステップＳ５と同様であるので、ここでの説明は省略する。
【００５８】
ステップＴ６：パラメータの推定
ここでは、目標点に係る理想特性値が、３次元実測特性に係る実測特性値で取り囲まれていない場合には、３次元ルックアップテーブルの格子点群において、目標点の近傍にある第１格子点の出力信号値と、この第１格子点を通って目標点近傍に向かう仮想ライン上にある格子点のうち、少なくとも第１格子点と隣接している第２格子点の出力信号値とを用いて、目標点の出力信号値の候補を求める。そして、この処理を目標点の近傍にある複数個の格子点について繰り返し行って出力信号値の候補をそれぞれ求め、これらの出力信号値の候補の平均値を目標点の出力信号値とする。具体的には、第1実施例のステップＳ６と同様であるので、ここでの説明は省略する。
【００５９】
ステップＴ７：
ステップＴ６で目標点の出力信号値が推定されると、次のステップＴ７に進み、３次元ルックアップテーブル２の全ての格子点について出力信号値（パラメータ）を算出したかどうかを判断し、全ての格子点について出力信号値が求まるまで、ステップＴ３〜Ｔ６を繰り返し実行する。
【００６０】
ステップＴ８：ノイズ成分の除去
3次元ルックアップテーブルの全ての格子点について出力信号値が求まると、出力信号値群に含まれるノイズ成分をラプラシアンフィルタやメディアンフィルタを使って除去する。
【００６１】
ステップＴ９：出力信号値（パラメータ）の更新
出力信号値群からノイズ成分を除去した後、これらのデータをコントラーラ３４を介して表示装置１０に送り、3次元ルックアップテーブル２の出力信号値を更新する。なお、入力側ルックアップテーブル１のγ特性や、出力側ルックアップテーブル３のγ特性も、コンピュータ２０を使って更新するようにしてもよい。
【００６２】
以上のように、第２実施例によれば、３次元ルックアップテーブル２のパラメータを精度良く設定することができるので、表示装置のカラー再現性を良好に維持することができる。
【００６３】
なお、上述の実施例で説明したノイズ成分の除去過程は、３次元実測特性を利用して３次元ルックアップテーブルの出力信号値を得ている場合に好適に用いることができ、3次元補間の手法などは特に限定されない。
【図面の簡単な説明】
【００６４】
【図１】第1実施例に係る表示装置の概略構成を示したブロック図である。
【図２】3次元ルックアップテーブル２の出力信号値（パラメータ）の決定手法を示したフローチャートである。
【図３】四面体補間の説明に供する図である。
【図４】四面体補間の説明に供する図である。
【図５】四面体補間の説明に供する図である。
【図６】四面体補間できない目標点の出力信号値を推定するための処理の説明に供する図である。
【図７】四面体補間できない目標点の出力信号値を推定するための処理の説明に供する図である。
【図８】第２実施例に係る表示装置とコンピュータの概略構成を示したブロック図である。
【図９】表示装置のカラー較正用プログラムの手順を示したフローチャートである。
【図１０】従来の３次元ルックアップテーブルの出力信号値の決定手法の説明に供する図である。
【符号の説明】
【００６５】
１  …  入力側ルックアップテーブル
２  …  ３次元ルックアップテーブル
３  …  出力側ルックアップテーブル
１０  …  表示装置
２０  …  コンピュータ
【図１】
0005410479.tif 000007
【図２】
0005410479.tif 000008
【図３】
0005410479.tif 000009
【図４】
0005410479.tif 000010
【図５】
0005410479.tif 000011
【図６】
0005410479.tif 000012
【図７】
0005410479.tif 000013
【図８】
0005410479.tif 000014
【図９】
0005410479.tif 000015
【図１０】
0005410479.tif 000016