(57)【要約】  【課題】光学的フレアの影響を低減して意図した通りの色再現を行い得る表示システム等を提供する。  【解決手段】フレア算出装置１３は、所定のテストカラー画像データをプロジェクタ１へ出力してスクリーン上に画像を投影し、投影された画像をテスト画像撮影カメラ４により撮像させてその撮影画像データを取得し、取得した撮影画像データと元のテストカラー画像データとに基づいて、この表示システムにおける光学的フレアの影響を表現する表示特性データを算出する。フレア補正装置１２は、画像データ記憶装置１１からカラー画像データを取得して、フレア算出装置１３により算出された表示特性データを用いて補正し、補正後のカラー画像データを画像データ記憶装置１１へ再び出力して記憶させる。【選択図】図２
2005189542.tif 000002
【特許請求の範囲】
【請求項１】
カラー画像を表示するためのカラー画像表示装置と、
カラー画像データを補正することにより、上記カラー画像表示装置へ出力するための補正カラー画像データを生成する画像補正装置と、
を具備し、
上記画像補正装置は、上記カラー画像表示装置へ出力する複数のテストカラー画像データと、これら複数のテストカラー画像データに各対応して該カラー画像表示装置に表示されるテストカラー画像の表示色の空間分布と、の関係に基づいて、上記カラー画像データから、上記カラー画像表示装置の光学的フレアを補正する補正カラー画像データを算出するものであることを特徴とする表示システム。
【請求項２】
上記画像補正装置は、上記テストカラー画像データに対応して上記カラー画像表示装置に表示されるテストカラー画像の表示色の空間分布を測定するためのテストカラー画像測定手段を有して構成されたものであることを特徴とする請求項１に記載の表示システム。
【請求項３】
上記テストカラー画像測定手段は、輝度計、色彩計、分光放射輝度計、モノクロカメラ、カラーカメラ、マルチバンドカメラ、の内の少なくとも１つを含むものであることを特徴とする請求項２に記載の表示システム。
【請求項４】
上記画像補正装置は、上記テストカラー画像データと、該テストカラー画像データに対応して上記カラー画像表示装置に表示されるテストカラー画像の表示色の空間分布と、に基づいて該カラー画像表示装置の表示特性データを算出する表示特性算出手段を有して構成され、上記表示特性算出手段により算出した表示特性データに基づいて、上記補正カラー画像データを算出するものであることを特徴とする請求項１から請求項３の何れか一に記載の表示システム。
【請求項５】
上記画像補正装置は、上記表示特性データに基づいて上記カラー画像データに係るフレア分布データを算出し、該カラー画像データからこのフレア分布データを減算することにより、上記補正カラー画像データを算出するフレア算出手段をさらに有して構成されたものであることを特徴とする請求項４に記載の表示システム。
【請求項６】
複数画素でなるカラー画像データの各画素データを成分とするベクトルをＰ、このベクトルＰの成分数をＮ（Ｎは自然数）、Ｎ×Ｎの単位行列をＥ、上記カラー画像表示装置の表示特性のＮ×Ｎの行列による表現をＭ、任意の定数をＫ、とすると、上記フレア分布データを表すベクトルＦは、次の数式、
2005189542.tif 000021
により与えられることを特徴とする請求項５に記載の表示システム。
【請求項７】
コンピュータに、
カラー画像表示装置へ複数のテストカラー画像データを出力して複数のテストカラー画像を表示させる第１の手順、
上記第１の手順により上記カラー画像表示装置に表示された複数のテストカラー画像の表示色の空間分布をそれぞれ取得する第２の手順、
上記複数のテストカラー画像データと、これら複数のテストカラー画像データに各対応して上記第２の手順により取得したテストカラー画像の表示色の空間分布と、に基づいて、カラー画像データから、上記カラー画像表示装置の光学的フレアを補正する補正カラー画像データを算出する第３の手順、
上記第３の手順により算出された補正カラー画像データを上記カラー画像表示装置へ出力して表示させる第４の手順、
を実行させるための表示プログラム。
【請求項８】
カラー画像表示装置へ複数のテストカラー画像データを出力して複数のテストカラー画像を表示させる第１の手順と、
上記第１の手順により上記カラー画像表示装置に表示された複数のテストカラー画像の表示色の空間分布をそれぞれ取得する第２の手順と、
上記複数のテストカラー画像データと、これら複数のテストカラー画像データに各対応して上記第２の手順により取得したテストカラー画像の表示色の空間分布と、に基づいて、カラー画像データから、上記カラー画像表示装置の光学的フレアを補正する補正カラー画像データを算出する第３の手順と、
上記第３の手順により算出された補正カラー画像データを上記カラー画像表示装置へ出力して表示させる第４の手順と、
を含むことを特徴とする表示方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、光学的フレアの影響を補正して画像の表示を行う表示システム、表示プログラム、表示方法に関する。
【背景技術】
【０００２】
近年、電子商取引や電子美術館・博物館などの実現を目的として、被写体の画像をディスプレイ上に正確な色で再現する技術の研究が活発に行われている。
【０００３】
これらの研究は、画像入力装置や画像表示装置の色特性を測定して、これらの装置の色特性情報に基づいて色信号の補正を行うものとなっている。ここで用いられる装置の色特性情報のフォーマットを標準化することは、色再現システムを普及させる上で重要であり、International Color Consortium（ＩＣＣ）では装置の色特性情報をデバイスプロファイルとして規定している（http://www.color.org）。
【０００４】
上述したようなＩＣＣのデバイスプロファイルや現行のカラー画像システムでは、カラー画像機器または画像データに対して、画像の空間座標には依存しない情報として色特性を規定し、この色特性情報に基づいて色再現を行っている。
【発明の開示】
【発明が解決しようとする課題】
【０００５】
上述したようなＩＣＣのデバイスプロファイルや現行のカラー画像システムでは、表示装置における光学的なフレアのような、画像内のある位置の色が他の位置の色に及ぼす影響についてまでを考慮することができるものとはなっていない。そのために、光学的なフレアの影響が生じるような表示装置では、十分に正確な色再現を行うことができない場合があった。
【０００６】
本発明は上記事情に鑑みてなされたものであり、光学的フレアの影響を低減して意図した通りの色再現を行い得る表示システム、表示プログラム、表示方法を提供することを目的としている。
【課題を解決するための手段】
【０００７】
上記の目的を達成するために、第１の発明による表示システムは、カラー画像を表示するためのカラー画像表示装置と、カラー画像データを補正することにより上記カラー画像表示装置へ出力するための補正カラー画像データを生成する画像補正装置と、を具備し、上記画像補正装置は、上記カラー画像表示装置へ出力する複数のテストカラー画像データと、これら複数のテストカラー画像データに各対応して該カラー画像表示装置に表示されるテストカラー画像の表示色の空間分布と、の関係に基づいて、上記カラー画像データから、上記カラー画像表示装置の光学的フレアを補正する補正カラー画像データを算出するものである。
【０００８】
また、第２の発明による表示システムは、上記第１の発明による表示システムにおいて、上記画像補正装置が、上記テストカラー画像データに対応して上記カラー画像表示装置に表示されるテストカラー画像の表示色の空間分布を測定するためのテストカラー画像測定手段を有して構成されたものである。
【０００９】
さらに、第３の発明による表示システムは、上記第２の発明による表示システムにおいて、上記テストカラー画像測定手段が、輝度計、色彩計、分光放射輝度計、モノクロカメラ、カラーカメラ、マルチバンドカメラ、の内の少なくとも１つを含むものである。
【００１０】
第４の発明による表示システムは、上記第１から第３の発明による表示システムにおいて、上記画像補正装置が、上記テストカラー画像データと、該テストカラー画像データに対応して上記カラー画像表示装置に表示されるテストカラー画像の表示色の空間分布と、に基づいて該カラー画像表示装置の表示特性データを算出する表示特性算出手段を有して構成され、上記表示特性算出手段により算出した表示特性データに基づいて、上記補正カラー画像データを算出するものである。
【００１１】
第５の発明による表示システムは、上記第４の発明による表示システムにおいて、上記画像補正装置が、上記表示特性データに基づいて上記カラー画像データに係るフレア分布データを算出し、該カラー画像データからこのフレア分布データを減算することにより、上記補正カラー画像データを算出するフレア算出手段をさらに有して構成されたものである。
【００１２】
第６の発明による表示システムは、上記第５の発明による表示システムにおいて、複数画素でなるカラー画像データの各画素データを成分とするベクトルをＰ、このベクトルＰの成分数をＮ（Ｎは自然数）、Ｎ×Ｎの単位行列をＥ、上記カラー画像表示装置の表示特性のＮ×Ｎの行列による表現をＭ、任意の定数をＫ、とすると、上記フレア分布データを表すベクトルＦは、次の数式、
2005189542.tif 000003
により与えられる。
【００１３】
第７の発明による表示プログラムは、コンピュータに、カラー画像表示装置へ複数のテストカラー画像データを出力して複数のテストカラー画像を表示させる第１の手順、上記第１の手順により上記カラー画像表示装置に表示された複数のテストカラー画像の表示色の空間分布をそれぞれ取得する第２の手順、上記複数のテストカラー画像データとこれら複数のテストカラー画像データに各対応して上記第２の手順により取得したテストカラー画像の表示色の空間分布とに基づいてカラー画像データから上記カラー画像表示装置の光学的フレアを補正する補正カラー画像データを算出する第３の手順、上記第３の手順により算出された補正カラー画像データを上記カラー画像表示装置へ出力して表示させる第４の手順、を実行させるためのプログラムである。
【００１４】
第８の発明による表示方法は、カラー画像表示装置へ複数のテストカラー画像データを出力して複数のテストカラー画像を表示させる第１の手順と、上記第１の手順により上記カラー画像表示装置に表示された複数のテストカラー画像の表示色の空間分布をそれぞれ取得する第２の手順と、上記複数のテストカラー画像データとこれら複数のテストカラー画像データに各対応して上記第２の手順により取得したテストカラー画像の表示色の空間分布とに基づいてカラー画像データから上記カラー画像表示装置の光学的フレアを補正する補正カラー画像データを算出する第３の手順と、上記第３の手順により算出された補正カラー画像データを上記カラー画像表示装置へ出力して表示させる第４の手順と、を含む方法である。
【発明の効果】
【００１５】
本発明の表示システム、表示プログラム、表示方法によれば、光学的フレアの影響を低減して意図した通りの色再現を行うことが可能となる。
【発明を実施するための最良の形態】
【００１６】
まず、実施例を具体的に説明する前に、本発明に用いられる原理について詳細に説明する。この原理は、入力画像データが、表示装置の光学的フレア（以下では適宜、単にフレアと述べる。）などの影響を受けて異なる表示画像データとなる場合において、どのような補正入力画像データを入力すれば本来の入力画像データと同一の表示画像データが得られるかを求めるものである。
【００１７】
［原理］
表示装置へ入力される画素数Ｎの入力画像データを（ｐ1，ｐ2，…，ｐN）tとする。ここに、右肩の添え字ｔは転置を表している。また、この入力画像データに対して、実際に表示される画像の光分布を画素数Ｎの画像データとして離散化して表した表示画像データを（ｇ1，ｇ2，…，ｇN）tとする。なお、この表示画像データは、例えば、デジタルカメラにより表示装置に表示されている画像を撮影することにより、取得することができる。この表示画像データは、表示装置のフレア等の影響を受けるために、一般には入力画像データと一致せず、他の画素位置の信号により表示される光の一部が重畳されたものとなる。さらに、表示装置への入力画像信号が０である場合にも、表示画像データは一般には０にはならない。このときの表示画像データをバイアス（ｏ1，ｏ2，…，ｏN）tとする。これらの影響を考慮して、入力画像データと表示画像データとの関係をモデル化したものが数式１である。  
［数１］
2005189542.tif 000004
【００１８】
この数式１の行列式における各要素の文字に対応する大文字を用いて、該数式１を簡単に表記したものが次の数式２である。  
［数２］
2005189542.tif 000005
なお、この数式２を含む以下の数式等においては、行列やベクトルを肉太文字で表しているが、表記上の点から、文章中では通常の太さの文字で代用することにする。
【００１９】
上記数式１あるいは数式２は、ある画素位置ｎ（ｎ＝１〜Ｎ）の表示光が、該画素位置ｎの周囲の他の画素位置に広がって、該他の画素位置の表示光に加算される現象を、全ての画素位置ｎについてモデル化したものである。ここで、数式２に示したような行列Ｍ（数式１においては、要素ｍij（ｉ＝１〜Ｎ，ｊ＝１〜Ｎ）でなる行列）を、表示装置の表示特性と呼ぶ。
【００２０】
上記表示画像データＧは、本来は、上記入力画像データＰと一致することが望ましいが、上述したように、表示画像はフレア等の影響を受けるために、一般には、入力画像データと表示画像データとは異なる。
【００２１】
そこで、入力画像データを補正し、補正した入力画像データを用いて取得した表示画像データが、本来の入力画像データと一致するか、もしくは近似するように、該補正した入力画像データを算出することを考える。この補正入力画像データをＰ’とすると、該補正入力画像データＰ’に対応する表示画像データである補正表示画像データＧ’は、次の数式３に示すようになる。  
［数３］
2005189542.tif 000006
【００２２】
この数式３に示した補正表示画像データＧ’が、本来の入力画像データＰと一致する条件式は、次の数式４に示すようになる。  
［数４］
2005189542.tif 000007
【００２３】
この数式４を満たすためには、上記数式３を用いることにより、次の数式５に示すような補正入力画像データＰ’を用いれば良い。  
［数５］
2005189542.tif 000008
ここに、表示特性Ｍは既知であるとする。
【００２４】
この数式５に示すような表示特性Ｍは、上述したようにＮ×Ｎの行列であり、例えば１２８０×１０２４の画素でなる表示装置の場合には、Ｎ＝１２８０×１０２４＝１３１０７２０となるために、一般にデータサイズは非常に大きくなるといえる。
【００２５】
これに対して、表示特性Ｍが特別な特性をもつ場合には、演算をより簡単に行うことができる可能性がある。例えば、表示装置におけるフレア等による光の広がり方が、画素位置にほぼ依存せず、同じ分布で近似することができる場合を考える。このときの表示特性をＭ’（ｍ’1，ｍ’2，…，ｍ’N）とすると、上記数式１または数式２は、この表示特性Ｍ’と入力画像データＰとのコンボリューションを用いて、次の数式６に示すように表される。  
［数６］
2005189542.tif 000009
ここに、記号「＊」はコンボリューション演算を表している。
【００２６】
この数式６における入力画像データＰを補正入力画像データＰ’とし、さらに上記数式４に示した条件を用いることにより、次の数式７が得られる。  
［数７］
2005189542.tif 000010
【００２７】
表示特性Ｍ’が既知である場合には、この数式７から補正入力画像データＰ’を算出することが可能である。すなわち、２つの画像のコンボリューション画像（ここではＰ−Ｏ）から、一方の画像（ここではＭ’）が既知である場合に、他方の画像（ここではＰ’）を算出するデコンボリューションに関する技術は公知である。例えば、文献１（Rosenfeld and A. C. Kak, Digital Picture Processing, Academic Press １９７６（長尾  真  監訳  ディジタル画像処理  近代科学社  １９７８））の第７章に記載されているような手法を用いることができる。
【００２８】
上記数式５に示したような補正方法は、要素数の多い行列Ｍの逆行列の演算が含まれ、また、上記数式７や文献１に記載されたような補正方法は、コンボリューションの逆演算などが含まれるために、演算が複雑であり、処理系の負荷が大きくなったり処理時間を要したりすることが考えられる。
【００２９】
これに対して、より簡便な補正方法として、フレア量を入力画像データから差し引くことで補正入力画像データを算出することを考える。今、表示画像データＧが、数式８に示すような入力画像データＰとバイアスＯとフレア等の影響により生じるフレア分布データであるフレア成分Ｆとの和により表されるものとしてモデル化する。  
［数８］
2005189542.tif 000011
【００３０】
この数式８におけるフレア成分Ｆは、さらに数式２を用いることにより、次の数式９に示すように表される。  
［数９］
2005189542.tif 000012
ここに、ＥはＮ×Ｎの単位行列を表している。
【００３１】
この数式９で与えられるＦを用いて、Ｐ−Ｆ−Ｍ-1Ｏを補正入力画像データＰ’として入力した場合に、補正表示画像データＧ’は、数式３から、
［数１０］
2005189542.tif 000013
となる。
【００３２】
ここで、Ｍの非対角成分が１よりも小さい場合、補正表示画像データＧ’のフレア補正誤差である数式１０の第２項の−（Ｍ−Ｅ）2Ｐは、補正前の表示画像データＧのフレアによる影響である数式９の（Ｍ−Ｅ）Ｐよりも小さい値となり、フレアの影響が改善されることが分かる。
【００３３】
さらにこの結果を考慮して、次の数式１１により与えられるＦ
［数１１］
2005189542.tif 000014
を用いてＰ−Ｆ−Ｍ-1Ｏを補正入力画像データＰ’として入力した場合、補正表示画像データＧ’は、数式３から、次の数式１２に示すようになる。  
［数１２］
2005189542.tif 000015
【００３４】
この数式１２の第２項の（Ｍ−Ｅ）3Ｐは、フレア補正誤差であり、次数が３次となっていて、より一層フレアの影響が小さくなっていることが分かる。
【００３５】
これを、同様にしてＫ次の項まで行った場合のフレアＦは、次の数式１３に示すように与えられる。  
［数１３］
2005189542.tif 000016
【００３６】
従って、この数式１３に示すようなフレアＦを補正するための補正入力画像データＰ’は、次の数式１４に示すようになる。  
［数１４］
2005189542.tif 000017
【００３７】
この数式１４におけるＫを大きな値にするほど、フレアの影響をより小さく抑制する補正を行うことができる。あるいは、計算量や計算時間などを考慮してＫを適宜の値とすることにより、処理系の負荷を抑制しながらフレアを良好に取り除く補正を行うことができる。
【００３８】
このような方法によるフレア補正に対しても、フレア等による光の広がり方が画素位置に依存することなく同じ分布であると近似できる場合には、数式９のＦはコンボリューションを用いて、次の数式１５に示すように置き換えることができる。  
［数１５］
2005189542.tif 000018
【００３９】
同様に数式１３は、
［数１６］
2005189542.tif 000019
と置き換えることができる。ここに、Ｅ’は、画像の中心位置に対応する成分が１、それ以外の成分が０となる列ベクトルである。
【００４０】
従って、数式１４に対応する補正入力画像データＰ’は、次の数式１７に示すようになる。  
［数１７］
2005189542.tif 000020
ここに、Ｏ~（「~」を「Ｏ」の上側に表記することができないために、このような表記法をとる。）は、ＯをＭ’でデコンボリューションしたものである。
【００４１】
ここまでの説明は、画像データを１チャンネルのデータとして扱ったものであったが、画像がカラー画像である場合は、画像データは一般に３チャンネルのデータとして扱われる。従って、この場合には、ＲＧＢ各チャンネル毎に上述したような表示特性Ｍまたは表示特性Ｍ’を算出して、算出された表示特性に基づきフレア補正を行うことになる。
【００４２】
さらに、カラー画像が４原色以上の原色数の画像データとして表される場合にも、同様に、各チャンネルの画像に対して上述したようなフレア補正処理を行うことにより、カラー画像として本来意図した表示を達成することが可能となる。
【００４３】
また、表示装置やデジタルカメラ等のカラー画像機器の信号値は、対応する輝度と非線形な関係にある場合がある。このような場合には、各信号の非線形性を補正してから上述したような処理を行う必要があるが、こうした階調特性の補正は公知の技術であるために、ここでは説明を省略している。すなわち、上述した説明は、信号の非線形性を補正した後の、リニアな空間における原理を説明したものとなっている。
【００４４】
以下、図面を参照して本発明の実施例を説明する。
【実施例１】
【００４５】
図１から図１０は本発明の実施例１を示したものであり、図１は表示システムの構成の概要を示す図である。
【００４６】
この表示システムは、画像を投影するためのカラー画像表示装置たるプロジェクタ１と、このプロジェクタ１により投影される補正画像を生成するための画像補正装置２と、上記プロジェクタ１から画像が投影されるカラー画像表示装置たるスクリーン３と、このスクリーン３に表示される画像の全領域を撮影し得るように配置された例えばデジタルカメラ等でなるテストカラー画像測定手段でありカラーカメラたるテスト画像撮影カメラ４と、を有して構成されている。
【００４７】
上記テスト画像撮影カメラ４は、広い意味での画像補正装置に含まれるものであり、撮影レンズの光学的特性に起因する画像の面内ムラや、あるいは撮像素子の面内感度ムラに起因する画像の面内ムラなどを補正することができる回路等を備えていて、出力する例えばＲＧＢのデジタル画像データは、これらの補正が行われたものであるとする。さらに、このテスト画像撮影カメラ４は、入力光強度に対して線形な応答で信号を出力するものとする。
【００４８】
次に、このような表示システムの動作について説明する。
【００４９】
まず、表示システムにおいて、光学的フレアを補正するのに必要な表示特性データを取得するときの動作は、以下のようになっている。
【００５０】
画像補正装置２は、予め内部に記憶している所定のテストカラー画像データをプロジェクタ１へ出力する。
【００５１】
プロジェクタ１は、該画像補正装置２から供給されたテストカラー画像データを、スクリーン３に投影する。
【００５２】
画像補正装置２は、このスクリーン３に表示されている上記テストカラー画像データに対応した表示色分布の画像を、テスト画像撮影カメラ４を制御することにより撮影させて転送させ、取り込む。
【００５３】
そして、該画像補正装置２は、テスト画像撮影カメラ４から取得したカラー画像データと、上記プロジェクタ１へ供給した元のテストカラー画像データと、を用いて、カラー画像データを補正するために用いる表示特性データを算出する。
【００５４】
次に、表示特性データ取得後の表示システムにより、通常の画像を投影するときの動作は、以下のようになっている。なお、この通常の画像を投影するときには、上記テスト画像撮影カメラ４は不要であるために、例えば取り外されている。
【００５５】
また、画像補正装置２の内部には、輝度に対して線形となるように変換されたカラー画像データが、予め記憶されているものとする。
【００５６】
画像補正装置２は、予め内部に記憶しているカラー画像データを、上述したように算出した表示特性データを用いて補正してから、補正カラー画像データとして内部に記憶しておく。
【００５７】
その後に、操作者により表示するカラー画像データが選択されたときに、該選択に対応する補正カラー画像データをプロジェクタ１へ出力する。
【００５８】
プロジェクタ１は、画像補正装置２から供給された補正カラー画像データに基づいて、スクリーン３へ画像を投影して表示する。
【００５９】
これにより、スクリーン３には、フレア等の影響を補正し得るような画像が表示され、観察者は、画像の表示者が本来意図した画像を観察することが可能となる。
【００６０】
なお、本実施例では、プロジェクタ１へ入力される画像データと、テスト画像撮影カメラ４から出力される画像データと、画像補正装置２で処理される画像データとは、全て、横１２８０画素×縦１０２４画素で、かつ各画素がＲＧＢの３色でなる、３チャンネルの画像データであるとする。
【００６１】
次に、図２は、画像補正装置２の構成を示すブロック図である。
【００６２】
この画像補正装置２は、予め記憶している所定のテストカラー画像データを上記プロジェクタ１へ出力して上記テスト画像撮影カメラ４から該テストカラー画像データに係る撮影されたカラー画像データ（撮影画像データ）を取得し、取得した撮影画像データと元のテストカラー画像データとに基づいて表示特性データを算出するフレア算出装置１３と、表示用のカラー画像データを記憶するものであってこのカラー画像データを後述するフレア補正装置１２により補正して得られる補正カラー画像データも記憶する画像データ記憶装置１１と、この画像データ記憶装置１１からカラー画像データを取得して上記フレア算出装置１３により算出された表示特性データを用いて補正し、補正後のカラー画像データ（補正カラー画像データ）を該画像データ記憶装置１１へ再び出力して記憶させるフレア補正装置１２と、を有して構成されている。
【００６３】
次に、このような画像補正装置２の動作について説明する。
【００６４】
まず、表示特性データを取得するときには、以下のように動作する。
【００６５】
上記フレア算出装置１３は、上記テストカラー画像データを上記プロジェクタ１へ出力して、テストカラー画像をスクリーン３上に表示させる。フレア算出装置１３は、これと同期して、テスト画像撮影カメラ４を制御することにより、スクリーン３の表示画像を撮影させる。そして、フレア算出装置１３は、撮影したカラー画像データを転送させることにより取得して、取得したカラー画像データと、元のテストカラー画像データと、に基づいて、表示特性データを算出し記憶する。ここに、フレア算出装置１３が算出する表示特性データは、２種類あり、１つが上述した数式１または数式２で定義されるＮ×Ｎの行列Ｍ、他の１つが数式６に用いられているＭ’である。ただし、画像データがＲＧＢ３チャンネルであることに対応して、各表記はＲＧＢ３チャンネル分のデータであるものとする。なお、画像データの画素数は１２８０×１０２４＝１３１０７２０であるために、この例ではＮは１３１０７２０となる。フレア算出装置１３が表示特性データＭ，Ｍ’を算出する処理については、後で図３を参照してさらに詳細に説明する。
【００６６】
続いて、カラー画像データを補正するときには、以下のように動作する。
【００６７】
フレア補正装置１２は、画像データ記憶装置１１に記憶されているカラー画像データを読み出すとともに、フレア補正方法に応じて、フレア算出装置１３から上述した２種類の表示特性データＭまたはＭ’の何れかを入力して、読み出した表示特性データに基づきフレア補正を行い、補正カラー画像データを算出する。
【００６８】
なお、本実施例の説明では、カラー画像データが上記原理の説明における入力画像データＰに対応し、補正カラー画像データが該原理の説明における補正入力画像データＰ’に対応している。上述したように、カラー画像データおよび補正カラー画像データは、各画素がＲＧＢ３チャンネルの画像データであるために、「Ｐ」と「Ｐ’」は、ぞれぞれＲＧＢ３チャンネルのデータを表すものとする。
【００６９】
該フレア補正装置１２は、次に説明するような第１〜第４の補正モジュールを有して構成されており、上記フレア算出装置１３から入力した表示特性データＭまたはＭ’を、これらの第１〜第４の補正モジュールにおいて用いて、補正カラー画像データを算出するようになっている。
【００７０】
まず、第１の補正モジュールは、表示特性データＭを入力して、数式５により補正カラー画像データＰ’を算出する。なお、バイアスＯは、全ての成分が０となるようなテストカラー画像をプロジェクタ１からスクリーン３に投影して、このスクリーン３の画像を上記テスト画像撮影カメラ４を用いて撮影する等により予め測定されて、フレア補正装置１２内に記憶されているものとする。
【００７１】
次に、第２の補正モジュールは、表示特性データＭ’を入力して、数式７に基づくデコンボリューション等の演算を行うことにより、補正カラー画像データＰ’を算出する。
【００７２】
さらに、第３の補正モジュールは、フレア算出手段であって、表示特性データＭを入力して、数式１４により補正カラー画像データＰ’を算出する。なお、定数Ｋは、画像補正装置２の操作者が任意に設定することができるようになっている。
【００７３】
そして、第４の補正モジュールは、フレア算出手段であって、表示特性データＭ’を入力して、数式１７により補正カラー画像データＰ’を算出する。なお、ここでも、数式１７における数式１６部分の項数（この項数は、上記定数Ｋに対応している。）は、画像補正装置２の操作者が任意に設定することができるようになっている。
【００７４】
このようにして、フレア補正装置１２により算出された補正カラー画像データは、該フレア補正装置１２から画像データ記憶装置１１へ出力されて、該画像データ記憶装置１１により記憶される。
【００７５】
そして、補正カラー画像データを投影して観察するときには、以下のように行う。
【００７６】
操作者が、画像補正装置２を操作することにより、該画像補正装置２に記憶されている所望の補正カラー画像データを選択する。すると、選択に対応した補正カラー画像データが画像データ記憶装置１１から読み出されてプロジェクタ１へ出力される。プロジェクタ１が、入力した補正カラー画像データに基づく画像をスクリーン３へ投影することにより、該スクリーン３上に表示されるカラー画像を光学的フレアの影響が低減された状態で観察することが可能となる。
【００７７】
図３は、フレア算出装置１３の構成を示すブロック図である。
【００７８】
このフレア算出装置１３は、テストカラー画像データや後述するような幾何補正パターン画像データを記憶しており必要に応じて上記プロジェクタ１と後述する撮影画像入力装置２２と後述する補正データ算出装置２３とへ出力するテスト画像出力装置２１と、上記テスト画像撮影カメラ４を制御することにより撮影されたカラー画像データを該テスト画像撮影カメラ４から入力するものであり上記幾何補正パターン画像データに基づいて後述するように幾何補正に必要な座標変換テーブルを算出して該テスト画像撮影カメラ４から入力したカラー画像データを幾何補正してから出力する撮影画像入力装置２２と、上記テスト画像出力装置２１から取得した元のテストカラー画像データと上記撮影画像入力装置２２を介して取得した撮影され幾何補正されたカラー画像データとに基づき上記表示特性Ｍと表示特性Ｍ’とを算出する表示特性算出手段たる補正データ算出装置２３と、この補正データ算出装置２３により算出された表示特性Ｍ，Ｍ’を記憶し必要に応じて上記フレア補正装置１２へ出力する補正データ記憶装置２４と、を有して構成されている。
【００７９】
次に、このようなフレア算出装置１３の動作について説明する。
【００８０】
テスト画像出力装置２１は、表示特性測定用のテストカラー画像データをプロジェクタ１へ出力するとともに、テストカラー画像データを出力した旨の信号を撮影画像入力装置２２へ送信する。
【００８１】
また、テスト画像出力装置２１は、プロジェクタ１へ出力したテストカラー画像データに関する情報を、補正データ算出装置２３へも出力する。
【００８２】
撮影画像入力装置２２は、テスト画像出力装置２１からの上記信号を受けると、テスト画像撮影カメラ４を制御することにより、プロジェクタ１からスクリーン３に投影されているテストカラー画像を撮影させる。こうしてテスト画像撮影カメラ４により撮影されたカラー画像は、撮影画像データとして撮影画像入力装置２２へ転送される。撮影画像入力装置２２は、取得した撮影画像データを、補正データ算出装置２３へ出力する。
【００８３】
補正データ算出装置２３は、テスト画像出力装置２１から入力した元のテストカラー画像データに関する情報と、撮影画像入力装置２２から入力した撮影画像データと、に基づいて、表示特性データを算出する処理を行う。
【００８４】
すなわち、この補正データ算出装置２３は、２種類の表示特性データＭ，Ｍ’に応じた２種類の表示特性データ算出モジュールを有して構成されていて、第１の表示特性データ算出モジュールは表示特性データＭを算出し、第２の表示特性データ算出モジュールは表示特性データＭ’を算出するものである。このとき、何れの表示特性データ算出モジュールを用いるかは、画像補正装置２の操作者が選択できるように構成されている。
【００８５】
次に、図４はテスト画像出力装置２１が出力する幾何補正パターン画像データを示す図、図５は十字パターンの中心位置の座標情報が記録されたテキストデータを示す図である。
【００８６】
テスト画像出力装置２１は、テストカラー画像データを出力するに先立って、まず、例えば図４に示すような幾何補正パターン画像データをプロジェクタ１へ出力する。
【００８７】
このときにテスト画像出力装置２１が出力する幾何補正パターン画像データは、図４に示すように、白色の背景に、黒色の十字パターンが、等間隔で縦４×横５に配列された画像データとなっている。
【００８８】
各十字パターンの中心位置の座標情報（幾何補正パターンデータ）は、図５に示すような形式により、テキストデータとしてテスト画像出力装置２１から撮影画像入力装置２２へ出力されるようになっている。
【００８９】
この図５に示す例では、最も左上に位置する十字パターンの中心位置を座標１、その右隣りに位置する十字パターンの中心位置を座標２、などとして、最も右下に位置する座標２０の十字パターンの中心位置までの画素位置が記載されている。このときの座標系は、例えば、最も左上角の画素を座標（０，０）、最も右下角の画素を座標（１２７９，１０２３）として、各画素の座標を表すような座標系が採用されている。
【００９０】
撮影画像入力装置２２は、後述するように、このような座標情報と、テスト画像撮影カメラ４から取得した幾何補正パターン画像の撮影画像データと、に基づいて、テストカラー画像データの空間座標と、テスト画像撮影カメラ４により撮影される画像の空間座標と、の対応関係を与える座標変換テーブルを作成する。
【００９１】
幾何補正に係る座標変換テーブルの作成が終了すると、テスト画像出力装置２１は、次に、テストカラー画像データをプロジェクタ１へ出力する。
【００９２】
図６はテスト画像出力装置２１が出力するテストカラー画像データにおける分割される領域の様子を示す図、図７は分割された領域の座標情報が記録されたテキストデータを示す図である。
【００９３】
このテストカラー画像データは、図６に示すように、１２８０画素×１０２４画素でなる全領域を縦４×横５に均等に分割した領域（２５６画素×２５６画素の領域）の内の、１つの領域のみでＲＧＢの内の何れか１色を例えば最大輝度で表示するような画像データとなっている。そして、分割された領域の全てに関して、ＲＧＢでなる３色の表示をそれぞれ行うために、６０種類のテストカラー画像データが用意されて順次表示されるようになっている。
【００９４】
このように合計２０個の領域に分割した理由は、各画素毎に処理を行う場合には、ＲＧＢの各３色について１画素単位で時系列的に順次発光させることになり、データを取得するのに要する時間が長くなり過ぎてしまうためであり、また、１画素単位で発光させると、そのフレアの影響を他の画素位置で検出しようとしても光量が足りず、また各画素の最大輝度にムラがあるなどのデータの安定性に欠ける可能性があるからである。こうして、ブロック単位で処理を行うことにより、十分な光量の安定したデータに基づいて、短時間に処理を行うことが可能となる。
【００９５】
このテストカラー画像データの領域の座標情報（パターンデータ）は、図７に示すような形式により、テキストデータとしてテスト画像出力装置２１から補正データ算出装置２３へ出力されるようになっている。
【００９６】
この図７に示す例では、上記図５において用いたのと同様の座標系が用いられており、左上角に位置する領域をパターン１、その右隣りに位置する領域をパターン２、などとして、右下角に位置するパターン２０の領域までの画素位置が記載されている。
【００９７】
より具体的には、（０，０，２５６，２５６）と記載されたパターン１は、領域内の左上角の画素位置が（０，０）で、この画素位置から（２５６，２５６）だけの大きさを有する領域であることを示している。従って、例えば（１０２４，７６８，２５６，２５６）と記載されたパターン２０は、領域内の左上角の画素位置が（１０２４，７６８）で、この画素位置から（２５６，２５６）だけの大きさを有する領域となっていることが分かる。
【００９８】
補正データ算出装置２３は、後述するように、このような座標情報と、テスト画像撮影カメラ４から取得したテストカラー画像の撮影画像データと、に基づいて、表示特性データＭ、または表示特性データＭ’を算出するようになっている。
【００９９】
次に、図８は、撮影画像入力装置２２の構成を示すブロック図である。
【０１００】
この撮影画像入力装置２２は、テスト画像出力装置２１からの信号に基づいて上記テスト画像撮影カメラ４に撮影を行わせるように制御するカメラ制御装置３１と、上記テスト画像撮影カメラ４により撮影された画像データを入力して記憶する撮影画像記憶装置３２と、この撮影画像記憶装置３２に記憶された幾何補正パターン画像に係る撮影画像と上記テスト画像出力装置２１からの幾何補正パターン画像に対応する座標情報とに基づいて幾何補正テーブルを算出する幾何補正データ算出装置３３と、この幾何補正データ算出装置３３により算出された幾何補正テーブルに基づいて上記撮影画像記憶装置３２に記憶されたテストカラー画像に係る画像データを幾何補正し上記補正データ算出装置２３へ出力する幾何補正装置３４と、を有して構成されている。
【０１０１】
続いて、このような撮影画像入力装置２２の動作について説明する。
【０１０２】
カメラ制御装置３１は、上記テスト画像出力装置２１から画像データをプロジェクタ１へ出力した旨の信号を入力すると、テスト画像撮影カメラ４を制御し撮影を行わせるためのコマンドを、該テスト画像撮影カメラ４へ出力する。
【０１０３】
撮影画像記憶装置３２は、テスト画像撮影カメラ４から伝送された画像データを入力して記憶する。該撮影画像記憶装置３２は、撮影画像データが幾何補正パターン画像に係るものである場合には該撮影画像データを幾何補正データ算出装置３３へ、また、テストカラー画像データに係るものである場合には該撮影画像データを幾何補正装置３４へ、それぞれ出力する。
【０１０４】
幾何補正データ算出装置３３は、撮影画像記憶装置３２から幾何補正パターン画像に係る撮影画像を入力するとともに、テスト画像出力装置２１から幾何補正パターン画像に対応する座標情報を入力して、幾何補正テーブルを算出する処理を行う。
【０１０５】
この幾何補正テーブルは、テスト画像撮影カメラ４から入力された画像データの座標を、テスト画像出力装置２１から出力される画像データの座標へ変換するためのテーブルデータであり、具体的には次のように算出される。
【０１０６】
まず、撮影画像記憶装置３２から入力された幾何補正パターン画像に係る撮影画像から、十字パターンを検出してその中心座標を求める。次に、得られた２０組の十字パターンの中心座標と、上記テスト画像出力装置２１から入力した幾何補正パターン画像に対応する座標と、の対応関係に基づき、幾何補正テーブルを算出する。
【０１０７】
十字パターンを検出する技術や、２０組のサンプル座標の対応関係から幾何補正テーブルを算出する技術は、多くの公知のものがあり、これらを適宜利用することができるが、ここではその詳細については説明を省略する。
【０１０８】
このようにして幾何補正データ算出装置３３により算出された幾何補正テーブルは、幾何補正装置３４へ出力される。
【０１０９】
幾何補正装置３４は、幾何補正データ算出装置３３から上述したように予め算出された幾何補正テーブルを入力するとともに、撮影画像記憶装置３２からテストカラー画像データに係る撮影画像を入力し、該幾何補正テーブルを参照してテストカラー画像データに係る撮影画像の座標変換を行い、変換後の画像データを補正データ算出装置２３へ出力する。
【０１１０】
補正データ算出装置２３は、テスト画像出力装置２１から入力したテスト画像の座標情報と、撮影画像入力装置２２から入力した幾何補正後のテストカラー画像の撮影画像と、に基づいて、表示特性データＭと表示特性データＭ’との少なくとも一方を算出し、補正データ記憶装置２４へ出力する。
【０１１１】
このような補正データ算出装置２３の処理について、図９および図１０を参照して説明する。図９はテストカラー画像に係る撮影画像に対して設定されるサンプル領域の様子を示す図、図１０は発光領域と発光領域以外のサンプル領域との様子を示す図である。
【０１１２】
補正データ算出装置２３は、表示特性データを求めるために、幾何補正後のテストカラー画像の撮影画像から所定のサンプル領域における信号値の取得を行う。
【０１１３】
このサンプル領域は、例えば図９に示すように設定されるようになっている。すなわち、サンプル領域を９画素×９画素の領域として設定して、このサンプル領域を、上記図５に示した２０個のテスト画像の発光領域における各中心座標部分に位置するように、等間隔に縦４×横５で配列して、サンプル領域Ｓ１〜Ｓ２０とする。
【０１１４】
そして、図１０に示すように、テストカラー画像の発光領域以外のサンプル領域（図１０に示す例では、左上角が発光領域となっているために、サンプル領域Ｓ２〜Ｓ２０）について信号値を求めて、各サンプル領域内における各画素の信号値の和（サンプル領域が９画素×９画素である場合には、８１画素分の信号値の和）をそれぞれ算出して平均値をとり、これらを各サンプル領域の中心座標に対するフレア信号とする。こうして、発光領域以外の１９個のサンプル領域の各中心座標のみについてのフレア信号分布をまず求める。なお、このように複数画素分のデータを加算して平均値をとっているのは、フレアの影響による光量がそれほど大きくないためと、データを平均化して信頼性を向上するためである。フレアには高周波成分はないと想定することができるために、このような処理が可能となっている。そして、サンプル領域のみの信号値に基づき処理を行うことにより、処理時間を短縮することが可能となる。
【０１１５】
その後、これら１９個のフレア信号に基づいて、補間処理を行うことにより、他の全ての画素位置におけるフレア信号を求める。このとき、図１０に示す例のように、発光領域が画像の四隅である場合には、近隣画素位置のフレア信号から外挿を行うことにより、フレア信号を求めることになる。
【０１１６】
こうして、１つのテストカラー画像に対して全画素のフレア信号を求め、このような処理を、図５に示したような２０個のテストカラー画像の全てについて行う。なお、ここでは３チャンネルのＲＧＢカラー画像を例に挙げているために、合計で６０個のテストカラー画像について上述した処理を行うことになる。
【０１１７】
１つのチャンネル毎に２０個ずつ得られたフレア信号分布は、上記図６に示したような２０個の発光領域の中心画素のみが発光したとき（すなわち、中心画素のみが発光画素となったとき）のフレア信号分布と見なされる。実際の発光は、２５６画素×２５６画素でなる発光領域全体の発光により行われるために、６５５３６で除算して１つの発光画素当たりのフレア信号分布に変換される。その他の発光画素位置に対応するフレア信号分布については、近隣の発光画素位置に対応するフレア信号分布を用いて、補間処理を行うことにより算出する。
【０１１８】
こうして、全ての画素が発光位置となったときのフレア信号分布が算出されるが、全領域が上述したように１２８０画素×１０２４画素で構成される場合には、全部で１３１０７２０個のフレア信号分布が算出されることになる。
【０１１９】
こうして、１３１０７２０画素各々のフレア信号を含んで構成されるフレア信号分布が、１３１０７２０画素の発光位置に一対一に対応するように１３１０７２０個作成されて、これが１３１０７２０行１３１０７２０列の行列式でなる表示特性データＭとなる。このような表示特性データは、上述したように、３つのチャンネルのそれぞれについて生成される。ここに、表示特性データＭは、行列式の各要素ｍijにおけるｊが発光画素の座標に対応し、ｉがフレア信号を取得する画素の座標に対応している。
【０１２０】
また、補正データ算出装置２３の第２の表示特性データ算出モジュールにより算出され、フレア補正装置１２の第２の補正モジュールまたは第４の補正モジュールにおいて用いられる表示特性データＭ’は、２０個の発光領域に各対応する２０個のフレア分布を、発光領域の中心座標が画像の中心座標となるようにそれぞれ座標移動した後に、これら２０個の平均値を求めることにより算出する。
【０１２１】
フレア補正装置１２は、このようにして算出した表示特性データＭまたはＭ’を用いて、上述したようにカラー画像データの補正を行い、補正カラー画像データとして画像データ記憶装置１１へ出力する。
【０１２２】
なお、通常の画像表示装置と同様に、得られた補正カラー画像データに対して、プロジェクタの階調特性を考慮した階調補正を行うが、この階調補正は色再現処理に関連する公知の技術であるために、ここではその説明を省略する。
【０１２３】
また、本実施例では、カラー表示装置としてプロジェクタを例に挙げて説明したが、本発明はこれに限定されるものではなく、例えば、ＣＲＴ、液晶パネル等の任意の画像表示装置に適用することができる。
【０１２４】
さらに、テストカラー画像データに対応する表示色の空間分布を取得する手段として、上述ではＲＧＢ方式のデジタルカメラで構成されるテスト画像撮影カメラ（カラーカメラ）を用いたが、モノクロカメラや４バンド以上のバンド数を有するマルチバンドカメラを用いることも可能である。あるいは、図９に示した例のような、測定サンプル数が比較的少ない場合には、表示色の空間分布を取得する手段として、カメラに代えて、分光放射輝度計、輝度計、色彩計等のスポット測定を行う測定器を用いることもできる。この場合には、測定の精度を上げ得ると期待することができる。
【０１２５】
そして、上述では、プロジェクタが投影する画像データやテスト画像撮影カメラにより取得する画像データが、全て横１２８０画素×縦１０２４画素である場合を例に示したが、異なる画素数の構成であっても良いのは勿論、表示する画素数と撮像する画素数とが異なっていても構わない。一般に、表示する画素数と撮像する画素数とは、任意の組み合わせを用いることが可能である。この場合には、補正カラー画像データのサイズに対応して表示特性データの算出を行うことになる。
【０１２６】
加えて、幾何補正パターンにおける十字パターンの数、テストカラー画像における発光領域の数、フレア信号測定のためのサンプル領域の数は、上述の例では全て２０としたが、これに限らず、夫々を独立に任意の数に設定することができる。また、測定精度や測定時間等を考慮して、画像補正装置の操作者が所望に設定することができるように構成しても構わない。
【０１２７】
なお、上述では、画像データを予めフレア補正した後に記憶しておき、スクリーンに投影する際には、すでに補正されている画像データを用いるようにしているが、十分な処理速度を確保することができる場合には、映像ソースから入力される画像データをリアルタイムでフレア補正して表示させるようにしても構わない。
【０１２８】
そして、上述では、ハードウェア的に処理を行う表示システムを例に挙げて説明したが、これに限らず、モニタ等の表示装置やデジタルカメラ等の測定装置が接続されたコンピュータにおいて、表示プログラムを実行することにより同等の機能を実現するようにしても良いし、こうした構成のシステム等に適用される表示方法であっても構わない。
【０１２９】
このような実施例１によれば、任意の画素位置の表示色が、それ以外の画素位置からの光により受ける影響を良好に低減することができ、高い色再現性をもったカラー画像を表示し得る表示システムを実現することが可能となる。
【０１３０】
また、テストカラー画像データに対応する表示色の空間分布を測定するテストカラー画像測定手段を備えたために、カラー画像表示装置の表示特性を正確かつ簡便に測定することが可能となる。そして、これにより、カラー画像表示装置の経時的変化にも対応することができる。
【０１３１】
特に、このテストカラー画像測定手段として、デジタルカメラ等のカラーカメラを用いることにより、より簡便に表示色の空間分布を取得することが可能となる。
【０１３２】
一方、テストカラー画像測定手段として、輝度計、色彩計、分光放射輝度計などを用いる場合には、より正確に表示特性を測定することが可能になる。あるいは、テストカラー画像測定手段としてモノクロカメラを用いることにより、低コストな装置構成とすることができる。さらには、テストカラー画像測定手段としてマルチバンドカメラを用いることにより、高精度に表示特性を取得することができるとともに、空間的な測定も高精度に行うことができる。
【０１３３】
そして、テストカラー画像データと、このテストカラー画像データに対応する表示色の空間分布と、に基づいて、カラー画像表示装置の表示特性データを算出して用いることにより、フレアのモデルに基づいた正確なフレア補正を行うことが可能となる。
【０１３４】
さらに、フレア算出手段によりフレア分布データを算出し、算出されたフレア分布データに基づいて補正カラー画像データを算出するようにしたために、補正画像データの算出を簡便に行うことが可能となる。
【０１３５】
また、フレア分布データを表すベクトルＦとして数式１３の表現を用いたために、この数式１３における定数Ｋを適切な値に設定することにより、フレア分布データの計算負荷と精度とを考慮した最適なフレア補正を行うことが可能となる。同様に、フレア分布データを表すベクトルＦとして数式１６の表現を用いる場合にも、右辺の項数を適切な値に設定することにより、計算負荷と精度とを考慮した最適なフレア補正を行うことが可能となる。
【０１３６】
こうして、光学的フレアの影響を低減して意図した通りの色再現を行い得る表示システムとなる。
【０１３７】
なお、本発明は上述した実施例に限定されるものではなく、発明の主旨を逸脱しない範囲内において種々の変形や応用が可能であることは勿論である。
【産業上の利用可能性】
【０１３８】
本発明は、表示システム、表示プログラム、表示方法一般に広く利用することができ、特に、光学的フレアの影響によりある画素の光が他の画素の色再現に影響を与える可能性のある表示システム、表示プログラム、表示方法に好適に用いることができる。
【図面の簡単な説明】
【０１３９】
【図１】本発明の実施例１における表示システムの構成の概要を示す図。
【図２】上記実施例１における画像補正装置の構成を示すブロック図。
【図３】上記実施例１におけるフレア算出装置の構成を示すブロック図。
【図４】上記実施例１において、テスト画像出力装置が出力する幾何補正パターン画像データを示す図。
【図５】上記実施例１において、十字パターンの中心位置の座標情報が記録されたテキストデータを示す図。
【図６】上記実施例１において、テスト画像出力装置が出力するテストカラー画像データにおける分割される領域の様子を示す図。
【図７】上記実施例１において、分割された領域の座標情報が記録されたテキストデータを示す図。
【図８】上記実施例１における撮影画像入力装置の構成を示すブロック図。
【図９】上記実施例１において、テストカラー画像に係る撮影画像に対して設定されるサンプル領域の様子を示す図。
【図１０】上記実施例１において、発光領域と発光領域以外のサンプル領域との様子を示す図。
【符号の説明】
【０１４０】
１…プロジェクタ（カラー画像表示装置）
２…画像補正装置
３…スクリーン（カラー画像表示装置）
４…テスト画像撮影カメラ（テストカラー画像測定手段、カラーカメラ）
１１…画像データ記憶装置
１２…フレア補正装置（フレア算出手段を含む）
１３…フレア算出装置
２１…テスト画像出力装置
２２…撮影画像入力装置
２３…補正データ算出装置（表示特性算出手段）
２４…補正データ記憶装置
３１…カメラ制御装置
３２…撮影画像記憶装置
３３…幾何補正データ算出装置
３４…幾何補正装置
代理人    弁理士      伊    藤      進
【図１】
2005189542.tif 000022
【図２】
2005189542.tif 000023
【図３】
2005189542.tif 000024
【図４】
2005189542.tif 000025
【図５】
2005189542.tif 000026
【図６】
2005189542.tif 000027
【図７】
2005189542.tif 000028
【図８】
2005189542.tif 000029
【図９】
2005189542.tif 000030
【図１０】
2005189542.tif 000031