"入力信号値に対応した3次元の格子点を備え、各格子点には入力信号値に対応する出力信号値が記憶されている、色変換用の3次元ルックアップテーブルの作成方法であって、前記3次元ルックアップテーブルに、1入力信号値に対応した3次元の格子点を備え、各格子点には入力信号値に対する表示器の理想とする表示色が記憶されている3次元理想特性を求める過程と、2入力信号値に対応した3次元の格子点を備え、各格子点には入力信号値に対する表示器の測定された実際の表示色が記憶されている3次元実測特性を求める過程と、33次元ルックアップテーブルの格子点群のうち、出力信号値を求めようとしている目標点に係る理想特性値が、3次元実測特性に係る実測特性値で取り囲まれているかどうかを判定する過程と、4目標点に係る理想特性値が、3次元実測特性に係る実測特性値で取り囲まれている場合には、目標点を取り囲んでいる3次元実測特性の各格子点の実測特性値を用いて、目標点の出力信号値を3次元補間により求める過程と、5前記4の処理で求められた各格子点の出力信号値からノイズ成分を除去する過程を経て求められた、入力信号値に対応する出力信号値を記憶させることを特徴とする3次元ルックアップテーブルの作成方法。しかしながら、このような構成を有する従来例の場合には、次のような問題がある。すなわち、従来手法によると、3次元実測特性を求める過程でノイズ成分を含む場合があり、その結果、3次元ルックアップテーブルの精度が悪くなるという問題点がある。なお、従来の四面体補間法による3次元ルックアップテーブルの出力信号値決定手法によると、4つの実測点で囲われない目標点があった場合に、その目標点に係る出力信号値を次のようにして求めている。以下、図10を参照して説明する。図10において、符号T1~T4は表示色が実測された実測点である。また、符号T0は、4つの実測点で囲われない目標点である。従来手法によると、目標点T0から近い2つの実測点T1とT2とを結ぶ線分に、目標点T0から垂線を下し、その交点T0 ́を求める。そして、実測点T1と交点T0 ́との距離と、格子点T2と交点T0 ́との距離に応じて、2つの実測点T1、T2の表示色を比例配分し、それを目標点T0の出力信号値と推定している。このように、従来手法によると、4つの実測点で囲われない目標点があった場合に、その目標点に係る出力信号値としては、近傍2点の実測点の表示色を直線補間して決定しているので、四面体補間法により求められた他の格子点の出力信号値との整合性が悪く、結果として、3次元ルックアップテーブルの精度が悪いという問題点がある。本発明は、このような事情に鑑みてなされたものであって、3次元ルックアップテーブルの出力信号値の精度を高めることを目的とする。本発明は、このような目的を達成するために、次のような構成をとる。すなわち、請求項1に記載の発明は、入力信号値に対応した3次元の格子点を備え、各格子点には入力信号値に対応する出力信号値が記憶されている、色変換用の3次元ルックアップテーブルの作成方法であって、前記3次元ルックアップテーブルに、1入力信号値に対応した3次元の格子点を備え、各格子点には入力信号値に対する表示器の理想とする表示色が記憶されている3次元理想特性を求める過程と、2入力信号値に対応した3次元の格子点を備え、各格子点には入力信号値に対する表示器の測定された実際の表示色が記憶されている3次元実測特性を求める過程と、33次元ルックアップテーブルの格子点群のうち、出力信号値を求めようとしている目標点に係る理想特性値が、3次元実測特性に係る実測特性値で取り囲まれているかどうかを判定する過程と、4目標点に係る理想特性値が、3次元実測特性に係る実測特性値で取り囲まれている場合には、目標点を取り囲んでいる3次元実測特性の各格子点の実測特性値を用いて、目標点の出力信号値を3次元補間により求める過程と、5前記4の処理で求められた各格子点の出力信号値からノイズ成分を除去する過程を経て求められた、入力信号値に対応する出力信号値を記憶させることを特徴とする。請求項1記載の発明によると、3次元実測特性を求める過程で混入したノイズ成分を除去するので、3次元ルックアップテーブルの出力信号値の精度を一層高めることができる。本発明によれば、3次元実測特性を求める過程で含まれるノイズ成分を除去するので、3次元ルックアップテーブルの精度を高めることができる。以下、図面を参照して本発明の第1実施例を説明する。図1は、第1実施例に係る表示装置の概略構成を示したブロック図である。本実施例に係る表示装置は液晶表示装置であって、入力側から順に配置された入力側ルックアップテーブルLUT1と、3次元ルックアップテーブルLUT2と、出力側ルックアップテーブルLUT3と、液晶表示器4とを備えている。出力側ルックアップテーブル3は、液晶表示器4のγ特性をキャンセルして、そのγ特性をほぼ直線化γ=1するためのものであって、入力信号値R2,G2,B2に応じて個別に設けられた1次元のルックアップテーブルである。一方、入力側ルックアップテーブル1は、出力側ルックアップテーブル3で液晶表示器4のγ特性がキャンセルされていることを前提として、表示装置のγ特性を任意に設定するためのものであって、入力信号値R0,G0,B0に応じて個別に設けられた1次元のルックアップテーブルである。3次元ルックアップテーブル2は色変換用に用いられるもので、入力信号値R1,G1,B1に対応した3次元の格子点を備え、各格子点には入力信号値R1,G1,B1に対応する出力信号値R2,G2,B2が記憶されている。本実施例の特徴である3次元ルックアップテーブル2の出力信号値パラメータの決定手法は後に詳しく説明する。本実施例装置によると、液晶表示器4に応じたビット数、例えばそれぞれが8ビット又は10ビットの入力信号値R0,G0,B0が入力側ルックアップテーブル1に与えられる。入力信号値R0,G0,B0は、入力側ルックアップテーブル1に予め記憶されている任意のγ特性によって変換されて信号値R1,G1,B1として出力される。この信号値R1,G1,B1が3次元ルックアップテーブル2に与えられる。3次元ルックアップテーブル2は、予め設定された色変換特性でもって入力信号値R1,G1,B1を変換して、信号値R2,G2,B2を出力する。この信号値R2,G2,B2は、出力側ルックアップテーブル3で、液晶表示器4のγ特性とは逆の特性で変換される。その出力信号値R3,G3,B3が液晶表示器4に与えられる。以上のように本実施例装置によると、出力側ルックアップテーブル3で液晶表示器4のγ特性がキャンセルされるので、入力側ルックアップテーブル1に設定された任意のγ特性で表示を行うことができる。また、3次元ルックアップテーブル2によって任意の色変換を行うことができるので、液晶表示器4の固有の色特性に影響されないように色再現をすることが可能である。なお、入出力のルックアップテーブル1、3の機能を3次元ルックアップテーブル2にもたせることもできるが、そうすると3次元ルックアップテーブル2のデータ量が膨大になる。本実施例では、入出力のルックアップテーブル1、3を3次元ルックアップテーブル2とは個別に設けているので、3次元ルックアップテーブル2のデータ量を低く抑えることができる。次に、3次元ルックアップテーブル2の出力信号値の決定手法を図2のフローチャートを参照して説明する。ステップS1:理想特性の算出まず、入力信号値R,G,Bに対して、どのような色再現をしたいのかを表した表示装置の理想特性を算出する。理想特性は、入力信号値に対応した3次元の格子点を備え、各格子点には入力信号値に対する表示器の理想とする表示色が記憶されている。本実施例では、表示色をXYZ表色系の三刺激値X,Y,Zで表しているが、Lab系で表してもよい。以下に、理想特性を算出する具体例を説明する。入力信号値R,G,BとXYZ表色系の赤、緑、青の三刺激値X,Y,Zと間には、次式1で示されるような関係がある。そこで、まず、1式中の変換行列Mを算出する。具体的には、赤の入力信号R1=255、G1=0、B1=0を与えたときの、XYZ表色系の三刺激値X,Y,Zを示す第1関係式を得る。三刺激値X,Y,Zの値は、ICCプロファイルに記述されたデバイス特性表示器4の特性から知ることができる。なお、実際に測定した表示器4の特性から三刺激値X,Y,Zの値を決めてもよいし、あるいは、sRGBやNTSCといった規格で定められている特性から三刺激値X,Y,Zの値を決めてもよい。赤の入力信号を与えたときと同様にして、緑の入力信号R=0、G=255、B=0を与えたときの第2関係式、青の入力信号R=0、G=0、B=255を与えたときの第3関係式を求める。前記第1~第3の関係式に基づいて、変換行例Mの各要素m00m01m22を決定する。変換行列Mが算出されたら、その変換行列Mと、入力側ルックアップテーブル1に設定したγ特性輝度特性とを用いて、任意の入力信号R0,G0,B0で特定される格子点の三刺激値X,Y,Z、すなわち、理想とする表示色を算出する。以下に具体的な算出手順を説明する。まず、次式2で示されるような任意の入力信号Viを設定する。次に、入力側ルックアップテーブル1に設定されたγ特性を考慮した信号R0i ́,G0i ́,B0i ́を求める。R、G、Bのそれぞれのγ特性をγr、γg、γbとすると、信号R0i ́,G0i ́,B0i ́は、次式3で求められる。γ特性を考慮した信号R0i ́,G0i ́,B0i ́と、先に求めた変換行列Mとを用いて、理想特性の各格子点の三刺激値X,Y,Zを次式4により求める。上記の演算を全ての任意入力格子点について実行することにより、表示装置の理想特性を得る。ステップS2:実測特性の取得表示装置に任意入力信号前記2式参照を与え、そのときに液晶表示器4で表示される色を測定して、実測された三刺激値Xi,Yi,Ziを得る。全ての格子点について測定して得られる3次元実測特性は、入力信号値に対応した3次元の格子点を備え、各格子点には入力信号値に対する表示器の測定された実際の表示色本実施例では、三刺激値Xi,Yi,Ziが記憶されている。ステップS3、S4:四面体補間の可否の判定3次元ルックアップテーブル3の格子点群のうち、出力信号値を求めようとしている目標点に係る理想特性値が、3次元実測特性に係る実測測定値で取り囲まれているかどうかを判定する。以下その手順を詳しく説明する。まず、3次元ルックアップテーブル3の格子点群のうち、出力信号値を求めようとしている目標点に係る理想特性値の色空間内の座標をP0以下、目標点P0と呼ぶとする。そして、3次元実測特性の各格子点に対応する実測点この実測点も色空間内の座標で表されるのうちで、目標点P0に最も近い1つの実測点P1を見出す。以下、図3を参照して説明する。見出した実測点P1を1つの頂点として、目標点P0を内包する6面体について、P1を頂点とする三角錐、すなわち四面体P1,P2,P3,P4を仮定する図3中の符号Aで示した四面体。そして、平面の方程式を使って、目標点P0が四面体P1,P2,P3,P4に囲まれているかどうかを判定する。具体的には次のような手順で判定する。まず、四面体P1,P2,P3,P4の全ての平面について方程式を求める。そして、目標点P0と各平面の位置関係を調べ、四つの平面に対して目標点P0が全て正方向に位置していると判定されれば、目標点P0は四面体P1,P2,P3,P4に囲まれていると判断する。平面の方程式は次のようにして求めることができる。図4を参照する。例えば、実測点P1,P2,P3を例に採る。この平面の方程式を、axbyczd=0で表すと、係数a、b、cは、点P1からP2へ向かうベクトルV1と、点P1からP3へ向かうベクトルV2との外積で得られる式5参照。ただし、d=-1×a×V1xb×V1yc×V1zである。点P4が平面P1,P2,P3の正の方向s≧0となるようにする。ここで、s=a×P4xb×P4yc×P4zdである。求められた平面の方程式に対して点P4が負の方向に位置した場合はa,b,c,dに-1を掛ける。以上のような判定を行って、目標点P0が四面体P1,P2,P3,P4に囲まれていないと判定された場合は、実測点P1を平面P2,P3,P4を挟んで反対側の実測点P1 ́に移動し、目標点P0が四面体P1 ́,P2,P3,P4図3中に符号Bで示すに囲まれているかどうかを同様に判定する。目標点P0が四面体P1 ́,P2,P3,P4に囲まれていな場合、目標点P0が平面P1 ́,P2,P3の外側負の方向にあれば、目標点P0は、四面体P1 ́,P2,P3,P4 ́図3中に符号Cで示すに囲まれていることになる。一方、目標点P0が平面P1 ́,P2,P4の外側にあれば、目標点P0は、四面体P1 ́,P2,P3 ́,P4図3中に符号Dで示すに囲まれていることになる。また、目標点P0が平面P1 ́,P3,P4の外側にあれば、目標点P0は、四面体P1 ́,P2 ́,P3,P4図3中に符号Eで示すに囲まれていることになる。ステップS5四面体補間によるパラメータ算出以上のようにして、目標点P0を取り囲む四面体の存在が確認されると、次の手順で目標点P0に対応する出力信号値R2,G2,B2を四面体補間により求める。以下、図5を参照して説明する。いま、目標点P0が四面体P1,P2,P3,P4に囲まれているとする。この四面体P1,P2,P3,P4を、更に目標点P0を頂点とする4つの四面体v1、v2、v3、v4に分割する。そして、これらの四面体v1~v4の各体積を用い、次式6に従って目標点の出力信号値R2,G2,B2を算出する。RGBP0=RGBP1×V1RGBP2×V2RGBP3×V3RGBP4×V4/V6式6において、RGBP0は、3次元ルックアップテーブル2の目標点P0に対応する出力信号値R2,G2,B2である。また、RGBP1~RGBP4は、3次元実測特性の実測点P1~P4に対応する入力信号値R,G,Bである。さらに、V1~V4は、四面体v1~v4の体積である。また、Vは、四面体P1,P2,P3,P4の体積である。本実施例では、目標点に係る理想特性値が、3次元実測特性に係る実測特性値で取り囲まれている場合に、四面体補間によって目標点の出力信号値を求めているが、3次元補間としては、四面体補間の他に、スプライン補間、三角柱プリズム補間、六面体立方体補間、ピラミッド補間、体積補間などを用いることができる。図2のフローチャートに戻る。ステップS4において、目標点P0を取り囲む四面体が存在しないと判定されると、ステップS6に進む。ステップS6:出力信号値の推定目標点に係る理想特性値が、3次元実測特性に係る実測特性値で取り囲まれていない場合には、3次元ルックアップテーブルの格子点群において、目標点の近傍にある第1格子点の出力信号値と、この第1格子点を通って目標点近傍に向かう仮想ライン上にある格子点のうち、少なくとも第1格子点と隣接している第2格子点の出力信号値とを用いて、目標点の出力信号値の候補を求める。そして、この処理を目標点の近傍にある複数個の格子点について繰り返し行って出力信号値の候補をそれぞれ求め、これらの出力信号値の候補の平均値を目標点の出力信号値とする。以下、図6及び図7を参照して具体的に説明する。図6において、符号Aで示した領域の格子点群は、四面体補間によって出力信号値が求められた3次元ルックアップテーブル2の格子点群である。一方、符号Bで示した領域の格子点群は、これらの格子点に係る理想特性値が、3次元実測特性に係る実測特性値で取り囲まれていないので、四面体補間で出力信号値を求めることができなかった格子点群である。いま、領域B内の目標点P0の出力信号値を求める場合を説明する。まず、目標点P0の近傍にある第1格子点を抽出する。3次元格子を考えた場合に、目標点を取り囲む近傍格子は26個ある。このうちステップS5の四面体補間で出力信号値が求められた格子点を第1格子点とする。説明の簡単のために、図7に示すように、目標点P0の近傍に3つの第1格子点P2、P3、P4があったとする。以下、第1格子点P2に絞って説明する。この第1格子点P2を通って目標点近傍に向かう仮想ライン上にあって、第1格子点と隣接している第2格子点P6を抽出する。そして、第2格子点P6から第1格子点P2に向かうベクトルVP2-P6を想定する。このベクトルVP2-P6を平行移動してP2を基点とするベクトルVP2-P6P2を想定し、そのベクトルV2×P2-P6の先端に目標点の候補P01があると推定する。そうすると目標点の第1候補P01の出力信号値RGBP01は、格子点P2およびP6の各出力信号値RGBP2、RGBP6を使って、次式7で求めることができる。RGBP01=2×RGBP2-RGBP67同様にして、第1格子点P3と、これに隣接する第2格子点P8とから、目標点の第2候補P02の出力信号値RGBP02は、次式8で求めることができる。RGBP02=2×RGBP3-RGBP88また、第1格子点P4と、これに隣接する第2格子点P10とから、目標点の第3候補P03の出力信号値RGBP03は、次式9で求めることができる。RGBP03=2×RGBP4-RGBP109図6及び図7で示した符号P0 ́は、目標点の第1~第3候補P01~P03の集合を示している。以上のようにして求められた目標点の候補点P01、P02、P03は、それぞれ格子点P6からP2に向かって生じている出力信号値の変化を候補点P01に反映させており、同様に、格子点P8からP3に向かって生じている出力信号値の変化を候補点P02に、格子点P10からP4に向かって生じている出力信号値の変化を候補点P03に、それぞれ反映させている。したがて、最終の目標点の出力信号値RGBP0としては、次式10に示すように、これらの候補点P01、P02、P03、の出力信号値の平均値を採れば良い。RGBP0=RGBP01RGBP02RGBP03RGBP0N/N10ここで、Nは、第1格子点の個数である。本実施例では、目標点の出力信号の候補を、第1格子点の出力信号値と、仮想ライン上で第1格子点と隣接する第2格子点の出力信号値とを使って求めているが、仮想ライン上で第2格子点に隣接する第3格子点の出力信号値を更に考慮して求めるようにしてもよい。図2のフローチャートに戻る。ステップS6で目標点の出力信号値が推定されると、次のステップS7に進み、3次元ルックアップテーブル2の全ての格子点について出力信号値を算出したかどうかを判断し、全ての格子点について出力信号値が求まるまで、ステップS3~S6を繰り返し実行する。ステップS8:ノイズ成分の除去全ての格子点について出力信号値が求まると、ステップS8に進んで出力信号値群の中に含まれるノイズを除去する。すなわち、3次元ルックアップテーブル2の格子点の出力信号値を求めるために、上述したステップS2で説明したように、液晶表示器4に表示される色を実測している。このような表示色の実測過程において、実測された3次元実測特性にノイズ成分が混入することがある。このノイズ成分が3次元ルックアップテーブル2の格子点の出力信号値に混入する結果、格子点の出力信号値の変化が不連続なものになる。格子点の出力信号値は、本来、連続性があると予測されるので、上記のノイズ成分をラプラシアンフィルタやメディアンフィルタを使って除去する。以上のようにして3次元ルックアップテーブル2の全ての格子点の出力信号値を決定することができる。本実施例によれば、3次元ルックアップテーブルの格子点群のうち、出力信号値を求めようとしている目標点に係る理想特性値が、3次元実測特性に係る実測特性値で取り囲まれていない場合には、目標点の近傍にある第1格子点の出力信号値と、この第1格子点を通って目標点近傍に向かう仮想ライン上にあって、第1格子点と隣接している第2格子点の出力信号値とを用いて、目標点の出力信号値の候補を求め、上記の処理を目標点の近傍にある複数個の格子点について繰り返し行って出力信号値の候補をそれぞれ求め、これらの出力信号値の候補の平均値を目標点の出力信号値としているので、3次元ルックアップテーブルの各格子点の出力信号値を精度良く求めることができる。第2実施例は、表示装置のカラー較正用プログラムに関する。図8は、表示装置のカラー較正用プログラムを実行するコンピュータと、カラー較正される対象である表示装置の概略構成を示したブロック図である。表示装置10は、第1実施例で説明したものと同様の構成であるので、ここでの説明は省略する。コンピュータ20は、中央演算処理装置CPU21と、データの読み出し専用メモリROM22と、データの一次記憶用のランダムアクセスメモリRAM23を備え、これらの要素がバス24を介して接続されている。また、コンパクトディスクCDROM25を駆動するCD-ROMドライバ26がコントローラ27を介してバス24に接続されている。さらに、ハードディスク28を駆動するハードディスクドライバ29がコントローラ30を介してバス24に接続されている。操作機器としてのキーボード31やマウス32がコントローラ33を介してバス24に接続されている。そして、表示装置10も、コントローラ34を介してバス24に接続されている。以下、図9に示した表示装置のカラー較正用プログラムのフローチャートを参照して第2実施例の動作を説明する。本実施例に係る表示装置のカラー較正用プログラムは、例えば記憶媒体としてのCD-ROM25などで提供されてコンピュータ20に取り込まれ、そのプログラムはハードディスク28に格納される。なお、カラー較正用プログラムは、インターネットなどを用いた通信で提供されてもよい。カラー較正時、CPU21は、ハードディスク28内のカラー較正用プログラムを読み出してその処理を実行する。ステップT1:理想特性の取得このステップT1では、入力信号値に対応した3次元の格子点を備え、各格子点には入力信号値に対する表示器の理想とする表示色が記憶されている3次元理想特性を取得する。具体的には、第1実施例のステップS1と同様である。ただし、変換行例Mを算出する際の基準となるRGB値や三刺激値X,Y,Zはオペレータがキーボード31などを操作して入力する。これらの値が入力されるとCPU21が1式に基づいて変換行例Mを算出する。さらに、オペレータがキーボード31などを操作して、入力側ルックアップテーブル31のγ特性を入力すると、CPU21は、先に算出した変換行例Mと、そのγ特性とを用いて、4式に従って理想特性の各格子点の三刺激値X,Y,Zを算出する。なお、理想特性の取得を上述した演算によらずに、別途作成した理想特性をコンピュータ20に直接に取り込むようにしてもよい。ステップT2:実測特性の取得このステップT2では、入力信号値に対応した3次元の格子点を備え、各格子点には入力信号値に対する表示器の測定された実際の表示色が記憶されている3次元実測特性を取得する。具体的には、第1実施例のステップS2と同様である。ただし、実測された三刺激値Xi,Yi,Ziは、オペレータがキーボード31などを操作して入力する。ステップT3、T4:四面体補間の可否の判定ここでは、CPU21が、3次元ルックアップテーブルの格子点群のうち、出力信号値を求めようとしている目標点に係る理想特性値が、3次元実測特性に係る実測特性値で取り囲まれているかどうかを判定する。具体的には、第1実施例のステップS3、S4と同様であるので、ここでの説明は省略する。ステップT5:四面体補間によるパラメータ算出ここでは、目標点に係る理想特性値が、3次元実測特性に係る実測特性値で取り囲まれている場合には、目標点を取り囲んでいる3次元実測特性の各格子点の実測特性値を用いて、目標点の出力信号値を3次元補間例えば四面体補間により求める。具体的には、第1実施例のステップS5と同様であるので、ここでの説明は省略する。ステップT6:パラメータの推定ここでは、目標点に係る理想特性値が、3次元実測特性に係る実測特性値で取り囲まれていない場合には、3次元ルックアップテーブルの格子点群において、目標点の近傍にある第1格子点の出力信号値と、この第1格子点を通って目標点近傍に向かう仮想ライン上にある格子点のうち、少なくとも第1格子点と隣接している第2格子点の出力信号値とを用いて、目標点の出力信号値の候補を求める。そして、この処理を目標点の近傍にある複数個の格子点について繰り返し行って出力信号値の候補をそれぞれ求め、これらの出力信号値の候補の平均値を目標点の出力信号値とする。具体的には、第1実施例のステップS6と同様であるので、ここでの説明は省略する。ステップT7:ステップT6で目標点の出力信号値が推定されると、次のステップT7に進み、3次元ルックアップテーブル2の全ての格子点について出力信号値パラメータを算出したかどうかを判断し、全ての格子点について出力信号値が求まるまで、ステップT3~T6を繰り返し実行する。ステップT8:ノイズ成分の除去3次元ルックアップテーブルの全ての格子点について出力信号値が求まると、出力信号値群に含まれるノイズ成分をラプラシアンフィルタやメディアンフィルタを使って除去する。ステップT9:出力信号値パラメータの更新出力信号値群からノイズ成分を除去した後、これらのデータをコントラーラ34を介して表示装置10に送り、3次元ルックアップテーブル2の出力信号値を更新する。なお、入力側ルックアップテーブル1のγ特性や、出力側ルックアップテーブル3のγ特性も、コンピュータ20を使って更新するようにしてもよい。以上のように、第2実施例によれば、3次元ルックアップテーブル2のパラメータを精度良く設定することができるので、表示装置のカラー再現性を良好に維持することができる。なお、上述の実施例で説明したノイズ成分の除去過程は、3次元実測特性を利用して3次元ルックアップテーブルの出力信号値を得ている場合に好適に用いることができ、3次元補間の手法などは特に限定されない。"