画像表示装置及び画像表示方法
【課題】マルチバンドカメラから取得されたマルチバンド画像に対し、スペクトルの表示や色度値の表現等を行うためのＧＵＩ表示方法を提案し、ユーザの画像評価の負荷を軽減すること。
【解決手段】コンピュータ１１によって、２つ以上のバンドを撮影できるマルチバンドカメラ１で取得したマルチバンド画像１２の内のあるバンドの画像を表示し、その画像上で、少なくとも１画素以上が指定されたとき、スペクトル計算部１３によって、上記マルチバンド画像における上記指定された少なくとも１画素以上の信号値から、分光スペクトルを推定し、その推定した分光スペクトルから色度情報を求め、該色度情報の色度値を色度図としてスペクトル出力用画面１４に表示する。そして、この表示した色度図において、任意の部分の指定を受けると、その指定された任意の部分の拡大表示を行う。
【解決手段】コンピュータ１１によって、２つ以上のバンドを撮影できるマルチバンドカメラ１で取得したマルチバンド画像１２の内のあるバンドの画像を表示し、その画像上で、少なくとも１画素以上が指定されたとき、スペクトル計算部１３によって、上記マルチバンド画像における上記指定された少なくとも１画素以上の信号値から、分光スペクトルを推定し、その推定した分光スペクトルから色度情報を求め、該色度情報の色度値を色度図としてスペクトル出力用画面１４に表示する。そして、この表示した色度図において、任意の部分の指定を受けると、その指定された任意の部分の拡大表示を行う。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、マルチバンドカメラにより取得されたマルチバンド画像からスペクトルグラフ表示、色度値表示、及び、複数の対象点どうしの相互の相対値を表示する画像表示装置及び画像表示方法に関する。
【背景技術】
【０００２】
マルチバンドカメラとは一般に、通常の３バンドカラーのカメラのＲ，Ｇ，Ｂフィルタと比べて、より特定の波長透過特性を持つ複数のバンド数で撮影するカメラを指す。バンドを切り替える手段は、例えば、特許文献１に開示されているような波長可変フィルタであったり、特許文献２に開示されているような回転フィルタを用いる。装置、目的によってバンド数は異なるが、各バンドでの透過率特性は例えば上記特許文献１の図２に示されるようになっている。
【０００３】
撮影した画像は一般にバンド数だけ存在するので、上記特許文献１の図２では１６バンドによる１６枚の画像になり、上記特許文献２の図８に示す回転フィルタを用いた場合には１０バンドによる１０枚の画像になる。
【０００４】
マルチバンドカメラで撮像されたマルチバンド画像は、バンド毎に画像として表示することができる。
【０００５】
上記特許文献１においても、また上記特許文献２でも、スペクトルを推定する方法については述べられているが、実際にソフトウエアから操作する方法については触れられていない。
【特許文献１】特開２００１−０９９７１０号公報
【特許文献２】特開平１１−０９６３３３号公報
【発明の開示】
【発明が解決しようとする課題】
【０００６】
スペクトルの表示は一般に、撮影した画像の全体ではなく、着目している特定の場所を対象に行う。その為には、指定は任意の場所をソフトウエアから指定し、その場所だけでスペクトルを表示できた方が良い。
【０００７】
また、３８０ｎｍから７８０ｎｍを１ｎｍおきに撮影し、４０１バンドの画像を撮った場合は、１ｎｍずつの４０１バンドの画像から精度良くスペクトルを推定できるが、４０１バンドよりも少ない、即ち１ｎｍよりも粗いサンプリング波長間隔で画像を撮った場合は、サンプリング波長間の広い、少ない画像データから４０１バンドの場合と略同等のスペクトルを推定する必要がある。推定方法によっても、スペクトルが異なることが考えられるため、任意に推定方法を選ぶことができたほうが良い。
【０００８】
さらに、エリアを複数指定した場合や複数の画像からスペクトルを表示する場合は、スペクトル同士の比較も望まれる。
【０００９】
スペクトルが求められると、それを利用し測色値を計算することができる。ここでもやはり、エリアを複数定義している場合は、ある基準に対してΔＥなど色差を計算し、表示することができると、ユーザの画像に対する理解を深めることができる。さらに、基準のエリアから、画像の各ピクセルがどの位の色差（ΔＥ）になっているか、一目でわかるような表示方法が望まれる。
【００１０】
本発明は、上記の点に鑑みてなされたもので、マルチバンドカメラから取得されたマルチバンド画像に対し、スペクトルの表示や色度値の表現等を行うためのＧＵＩ表示方法を提案し、ユーザの画像評価の負荷を軽減することが可能な画像表示装置及び画像表示方法を提供することを目的とする。
【課題を解決するための手段】
【００１１】
上記の目的を達成するために、本発明の第１の態様による画像表示装置は、
２つ以上のバンドを撮影できるマルチバンドカメラで取得したマルチバンド画像を表示する画像表示装置において、
上記マルチバンド画像の内のあるバンドの画像を表示する画像表示手段と、
上記画像表示手段に表示された上記あるバンドの画像上で、少なくとも１画素以上を指定する画素指定手段と、
上記マルチバンド画像における上記画素指定手段によって指定された上記少なくとも１画素以上の信号値から、分光スペクトルを推定するスペクトル推定手段と、
上記スペクトル推定手段で推定した分光スペクトルから、色度情報を求める色度情報算出手段と、
上記スペクトル推定手段で推定した分光スペクトルを表示すると共に、上記色度情報算出手段で求めた上記色度情報の色度値を色度図として表示するスペクトル表示手段と、
を具備し、
上記スペクトル表示手段は、上記表示した色度図において、任意の部分の拡大表示機能を備えることを特徴とする。
【００１２】
また、上記の目的を達成するために、本発明の第２の態様による画像表示方法は、
２つ以上のバンドを撮影できるマルチバンドカメラで取得したマルチバンド画像を表示する画像表示方法において、
上記マルチバンド画像の内のあるバンドの画像を表示する工程と、
その表示された上記あるバンドの画像上で、少なくとも１画素以上を指定する工程と、
上記マルチバンド画像における上記指定された上記少なくとも１画素以上の信号値から、分光スペクトルを推定する工程と、
その推定した分光スペクトルから、色度情報を求める工程と、
上記推定した分光スペクトルを表示すると共に、上記求めた色度情報の色度値を色度図として表示する工程と、
を有し、
上記表示する工程は、上記表示した色度図において、任意の部分の拡大表示を行う工程を含むことを特徴とする。
【発明の効果】
【００１３】
本発明によれば、マルチバンドカメラから取得されたマルチバンド画像に対し、スペクトルの表示や色度値の表現等を行うためのＧＵＩ表示方法を提案し、ユーザの画像評価の負荷を軽減することが可能な画像表示装置及び画像表示方法を提供することができる。
【発明を実施するための最良の形態】
【００１４】
以下、本発明を実施するための最良の形態を図面を参照して説明する。
【００１５】
［第１実施形態］
まず、本発明の第１実施形態として、画像上の任意の場所を対象とした、信号レベルのスペクトルを表示する方法を説明する。
【００１６】
本実施形態に係る画像表示装置の構成を示す図１（Ａ）において、マルチバンドカメラ１は、コンピュータ１１によって制御されている。マルチバンドカメラ１は一般に、撮像素子であるＣＣＤを内蔵し、画像を撮影して電子データとして画像データをコンピュータ１１に送る。撮影された画像は、マルチバンド画像１２としてコンピュータ１１の管理するメモリ上にある。なお、画像はその都度撮影されても良いし、過去に撮影した画像をハードディスク等の記憶媒体から読み出してきてメモリ上にマルチバンド画像１２として展開されていても良い。従って、本実施形態の処理を行う時は必ずしもマルチバンドカメラ１が無くても良い。
【００１７】
コンピュータ１１は、上記マルチバンド画像１２から、スペクトル計算部１３によってスペクトルを計算し、スペクトル出力用画面１４にグラフ表示を行う。本実施形態では、コンピュータ１１はＩＢＭ社のＰＣ−ＡＴ互換機を用い、ＯＳとしてＷｉｎｄｏｗｓ（登録商標）２０００を利用する。スペクトル出力用画面１４はＷｉｎｄｏｗｓ２０００の管理するウィンドウ画面であり、その画面はＣＲＴなどのモニタに表示される。
【００１８】
マルチバンドカメラ１は１８バンドとすると、上記マルチバンド画像１２は１８枚のバンド画像でなり、コンピュータ１１はユーザの指定に従ったバンドの画像を表示する。例えば、バンド６の画像が図２のように表示する。ユーザは、このバンド画像画面１５を見ながら、スペクトルを表示したい場所を特定する。この図２の例では、コンピュータ１１のポインティングデバイスであるマウスのドラッグ＆ドロップ機能（矩形の左上、右下を指定する方法）を使って、矩形を５つ指定している。なお、エリアの数や形はユーザによって異なる形でも良く、この実施形態の範疇に入る。
【００１９】
そしてユーザが、図示していない「信号値のスペクトル表示」のメニューを選択すると、コンピュータ１１は、図３のようなＧＵＩ（グラフィカルユーザインターフェース）、即ち補間方法選択画面１６を表示して、補間の方法を選択させる。
【００２０】
このような補間方法選択画面１６で、直線補間が選択された場合には、コンピュータ１１は、スペクトル計算部１３により、それぞれの矩形内の画素値、つまり、これはマルチバンドカメラ１のＣＣＤ素子の信号値、の平均値を算出する。そして、その輝度をバンド毎にプロットして直線で結んだグラフをスペクトル出力用画面１４として、図１（Ｂ）に示すように表示する。なお、このグラフの縦軸は信号の強度である。ＣＣＤ素子の出力が１６ｂｉｔであるときは、信号強度（輝度）として、０から６５５３５までの値をとることができる。また、横軸は波長を表している。この例では、バンド１が３８０ｎｍ、バンド２が４１８ｎｍの中心波長であり、同様にバンド１８まで中心波長を表す波長での輝度をプロットしていく。そして、選択された補間方法は、この例では直線補間であったので、バンド間は直線で補間する。なお、該スペクトル出力用画面１４の右上には、エリア選択リストボックス１４１が配されており、このエリア選択リストボックス１４１で選択することで、上記バンド画像画面１５で特定された各エリアごとの信号強度のスペクトルを、適宜切り替えながら観察することができるようになっている。
【００２１】
また、ユーザが図３の補間方法選択画面１６で、スプライン補正を選択したときには、補正間隔を１０ｎｍ、５ｎｍ、１ｎｍのうちから選択して、直線補間のときと同様に、バンド１〜バンド１８でのｘ座標、ｙ座標を求めてからスプライン補正を行う。
【００２２】
以上のようにして、画像上の任意の場所について、信号レベルのスペクトルを表示することができる。
【００２３】
画像上の任意の場所について、分光スペクトルを表示する方法を、本実施形態の変形例として説明する。
【００２４】
この場合は、図４に示すように、図１（Ａ）と同様の構成に加えて、参照画像１７を利用する。また、スペクトル計算部１３が分光スペクトルを計算するものとなっている。
【００２５】
即ち、分光スペクトルは、通常の物体に対しては分光反射率、顕微鏡等のように透過物体を観察する場合は分光透過率のスペクトルとしてグラフ表示を行う。従って、反射率か透過率を計算するための基準となる画像、つまり参照画像１７が必要となる。ここで、反射率スペクトルを求める場合は、参照画像１７の対象は標準白色板を撮影したものになる。また、透過率スペクトルを求める場合の参照画像１７は、観察するものを取り除き、照明光だけの画像を用いる。
【００２６】
つまり、参照画像１７は、マルチバンド画像１２を撮影する時と同じ手順で、撮影の対象を標準白色板などに代えて撮影する。従って、撮影は、マルチバンド画像１２の撮影の他に、そのような標準白色板などの撮影を１回、余分に行うことが必要となる。しかしながら、マルチバンド画像１２の中に、参照する場所を一緒に写し込むことができれば、マルチバンド撮影する撮影対象の撮影のみで良い。
【００２７】
以下に、その方法を説明する。図５（Ａ）は、撮影対象を撮影したマルチバンド画像１２の内、ユーザに指定された或いは初期画像として設定されたバンドの画像を表示したバンド画像画面１５を示す図である。この例では、色票が撮影対象になっている。もし、この色票の中に、反射率が極めて１００％に近いものがあった場合は、これを標準白色板と同様な扱いと定義することができる。この例では、図５（Ａ）の左下の矩形がその場所であり、ユーザはマウスのドラッグ＆ドロップ操作で、その場所を指定する。或いは、コンピュータ１１で自動的にその場所を認識しても良い。この例では、指定された場所が矩形１５１で図示されている。コンピュータ１１は、その矩形１５１の信号値（輝度）の平均値を使って参照画像１７を作成する。なお、こうして作成された参照画像１７を、図５（Ｂ）に示すような参照画像画面１８に表示して、ユーザに確認を受けることができる。
【００２８】
なお、参照画像１７は、光源に対して反射率が１００％の対象物が理想であるため、信号値（輝度）の明るいところを探すためには、図６のようなＧＵＩ、即ちバンド選択画面１９を表示し、マウス位置の１つ、或いは複数の画素をエリアとして定義し、その画素の信号値（輝度）を表示すると、評価を行い易い。ここで、図６の中央左側の「Ｌｅｖｅｌ」以下はバンド番号１９１を表し、バンド１〜バンド１８が表示されている。「Ｌｅｖｅｌ」の右側の数値「（１９７４，０）」は、現在のバンド画像画面１５に表示された画像上のマウス位置を表す。その下の数字で信号値を表示している。例えば、バンド１では「１７６（１／２５６）」と表示されているが、これは、１６ｂｉｔデータとして信号値が１７６であることを示している。また括弧内の数値は、コンピュータ１１のディスプレイは０〜２５５の２５６階調で輝度を表すことが多いため、この２５６階調換算で値が幾つになるかを示している。この場合は１である。
【００２９】
上記スペクトル計算部１３では、分光反射／透過率を計算して、スペクトル出力用画面１４へ表示を行う。その方法を手順に従って説明する。対象のマルチバンド画像１２に対して、分光反射／透過率のスペクトルを表示する対象の場所を指定する。図２はその様子を表し、入力方法は前述した手順で行う。この画像には参照画像１７が対になっている。作成方法は前述したように対象画像の一部から作るか、参照板だけの画像を改めて撮影する。
【００３０】
次に反射率／透過率を求める。この例では色票を使っているため、反射率である。反射率は次の式で求めることができる。
【００３１】
対象画像の指定バンドの反射率＝（対象画像の指定バンドの指定エリアの平均値）／（参照画像の指定バンドの指定エリアの平均値）
図２の例では、バンド数は１８なので、バンド毎に値を求める。また、指定エリアが５つあるので、それぞれのエリアで反射率を求める。例えば、図２の左下のエリア１は白色板であったとして、この場所から参照画像１７を前述した方法で作成した場合は、全てのバンドにおいてエリア１の反射率は１になる。それ以外のエリアでは一般に参照画像１７よりも信号値は小さいはずなので、０≦反射率≦１．０を示すのが普通である。
【００３２】
そして、求めた反射率をスペクトル出力用画面１４にグラフとして表示する。なお、この場合も、ユーザが図７に表すようなＧＵＩ、即ち補間方法選択画面２０で補間の方法を選択できるようにする。図７に示す状態では、直線補間が選ばれている。グラフにプロットする方法は、前述の信号レベルのスペクトルをグラフ表示させた時と同様な手順になる。つまり、図１（Ｂ）のグラフで横軸に各バンドで用いたフィルタの中心波長の位置を波長換算で求め、縦軸に反射率として値をプロットする。従って、前述の場合は縦軸は信号値であったが、分光反射率では、一般に１．０を最大値とした反射率で描画を行う。直線補間では、プロットの間を前述したように直線で補間し、スペクトルグラフを完成させる。
【００３３】
また、図７の補間方法選択画面２０で、補間方法としてスプライン補正を選択したときは、補間を１０ｎｍ、５ｎｍ、１ｎｍで行うか選択させた後、スプライン曲線で近似を行えば良い。
【００３４】
なお、図５（Ａ）及び（Ｂ）を参照して、ユーザが目で確認しながら、参照画像１７を作成する方法を説明したが、それ以外に参照画像１７を作成する変形例について述べる。
【００３５】
参照板は照明光を理想的には１００％反射するようにしたものであるため、図５（Ａ）に示すバンド画像画面１５で表示されるバンド画像内に参照板に相当する場所があれば、その場所の信号輝度は他と比べて高いはずである。従って、図５（Ａ）のようなバンド画像において全ての画素の信号値（輝度）を調べて、明るい順に並べる。その中で明るい方から例えば、３０画素分を参照板位置での画素情報として、その３０画素の信号値の平均値から図５（Ｂ）の参照画像画面１８に示すような参照画像１７を作成することができる。
【００３６】
もう１つの変形例として、露出オーバの部分を除いて、参照画像１７を作成する方法を説明する。図５（Ａ）に示すバンド画像画面１５で表示されるようなバンド画像は、適度な露出時間を設定して撮影したつもりでも、部分的に非常に明るい部分が存在してしまう可能性がある。例えば、参照板では参照板自身の材質、照明位置が調整されているので、正反射部分が除かれている場合が多いが、正式な参照板が用意できなかった場合や、照明条件を適切にできなかった場合は、参照としたい画像の部分に正反射部分を含むことが考えられる。この場合は、前述した方法で、バンド画像の全ての画素の信号値（輝度）を調べて、明るい順に並べるが、その後で、一定の輝度以上のものを取り除く。例えば、信号値の最大値が２５５であった場合は、並べられた画素のうち、２４０以上の信号値は正反射部分として除外する。そして、残った画素の中から、明るい輝度の３０画素分を同様な方法で参照画像１７を作成することができる。なお、この例ではスレッシュホールドを２４０としたが、マルチバンド画像１２の入力手段、信号値の幅、撮影物、ユーザの指定で他の値にしても同様な方法が成り立つ。
【００３７】
次に、上記特許文献２に開示されているようなスペクトル推定方法で、分光スペクトルの推定を行う例について述べる。
【００３８】
図８は、その場合の構成を示す図で、マルチバンドカメラ１、コンピュータ１１、マルチバンド画像１２、及び、スペクトル出力用画面１４の構成と機能は図１（Ａ）と同じである。分光スペクトルの推定は、スペクトル計算部１３で行う。ここでの推定方法は、上記特許文献２で記述されている。また、その推定の為に必要なファイルが特性ファイル２１として用意され、その特性ファイル２１をスペクトル計算部１３が参照する。
【００３９】
次に手順を示す。画像は今までと同様に、図２のような形で表示され、スペクトルを計算したい場所を特定する。図２では５つのエリアが指定されている。そして、その後、ユーザが分光反射／透過率のスペクトル表示メニューを選択したとき、図７の補間方法選択画面２０を表示し、補間方法として、３番目の「推定」が選択されると、先に選択した５つのエリアで上記特許文献２の方法でスペクトルの推定を行う。その結果は、図１（Ｂ）のようなグラフで表示される。ただし、この場合は縦軸は反射率であり、一般的には最大値を１．０とした０〜１．０の値になる。
【００４０】
なお、以上のいずれの方法でもエリアは、図２に示すようなバンド画像画面１５を使って予め指定しておいた。しかしながら、推定処理が高速に実行され、ユーザにストレスを感じさせない程度であるならば、マウスで１画素、或いは複数の画素を選択させて、その位置のスペクトルをリアルタイムで計算して表示するようにしても良い。或いは、それらの混合で例えば、エリア１〜エリア５は予め指定しておき、エリア６は現在のマウス位置にしてスペクトルを表示するというようにしても良い。
【００４１】
また、エリアが複数ある場合のスペクトルの比較は、図１（Ｂ）のようにグラフをエリア毎に表示しエリアの選択を変えることで行っていたが、図１（Ｂ）のウィンドウをエリア毎に作成し、それらを並べてタイル表示として、比較できるようにしても良い。或いは、１つのウィンドウにエリア毎に色を変えるなどの工夫をして同時に表示しても良い。さらに、相対的に信号値や分光反射／透過率のスペクトルを比較する場合は、図９に示すように、あるエリアを基準とした相対値を表示しても良い。この図では、エリア１を基準とした信号値（輝度）の相対値を縦軸として、グラフ表示を行っている。
【００４２】
また、図１０に示すように、図４の構成と図８の構成を組み合わせ、適宜ユーザが選択できるようにしても良い。
【００４３】
以上のような方法によって、画像の中でスペクトルを求めたい場所を特定し、グラフを表示することができる。推定方法が複数である場合もそれらを選択することが可能である。また、参照画像１７を作成する場合も、指定したエリアから作成することができる。以上のような方法によって、スペクトルを表示したい場所を画像を見ながら指定し、そのエリアでのスペクトルを表示することができる。また、グラフを描画する対象のエリアは次々に変えることができるので、エリア毎のスペクトルの違いを一目にして区別することができる効果がある。
【００４４】
［第２実施形態］
次に、本発明の第２実施形態を説明する。本実施形態は、マルチバンドカメラ１で撮影したマルチバンド画像に対し、特定のエリアを指定し、その場所の色度値、色度図を表示するというものである。
【００４５】
色度値は、分光反射率／透過率のスペクトルから推定するため、構成は上記第１実施形態と同様であり、スペクトル計算部１３により上記第１実施形態で説明した方法で分光反射率／透過率のスペクトルを計算した後、後述する方法で色度値が求められる。そして結果は、スペクトル出力用画面１４と同様に同じウィンドウ、或いは別ウィンドウで画面に表示される。色度値は前述したように分光反射率／透過率のスペクトルから求められるため、その推定方法によって、図１（Ａ）、図４、図８、または図１０の構成が考えられる。
【００４６】
次に、色度値ＸＹＺ、ｘｙｚ、Ｌ＊，ａ＊，ｂ＊を表示する方法を説明する。
【００４７】
まず、測色したい部分を指定して、その場所の分光反射／透過率スペクトルを推定する。その方法は、上記第１実施形態で説明した通りである。その結果、推定された分光反射／透過率スペクトルをＲ（λ）であったとする。この形態ではλは３８０［ｎｍ］〜７８０［ｎｍ］まで、１［ｎｍ］間隔で求められたとする。この時の照明光をＳ（λ）とすると、
Ｘ＝ＫΣＳ（λ）ｘ（λ）Ｒ（λ）
Ｙ＝ＫΣＳ（λ）ｙ（λ）Ｒ（λ）
Ｚ＝ＫΣＳ（λ）ｚ（λ）Ｒ（λ）
Ｋ＝１００／（ΣＳ（λ）ｙ（λ））
で表される。ただし、ｘ（λ）、ｙ（λ）、ｚ（λ）はＸＹＺ表色系における等色関数の値であり、Σはλを３８０〜７８０とする。ここで求まったＸＹＺより、ｘｙｚは次の式で求めることができる。
【００４８】
ｘ＝Ｘ／（Ｘ＋Ｙ＋Ｚ）
ｙ＝Ｙ／（Ｘ＋Ｙ＋Ｚ）
ｚ＝Ｚ／（Ｘ＋Ｙ＋Ｚ）
さらに、レンダリング照明光をｅ（λ）として、Ｌ＊，ａ＊，ｂ＊を求めると、次の公式で表すことができる。
【数１】
【００４９】
以上の公式から求めた値を、図１１に示すようなウィンドウ２２に表示する。このウィンドウ２２の右側では、エリア１に対して、（ｘ，ｙ）が（０．３１３，０．３２９）、Ｙが９９．９９３、（Ｌ＊，ａ＊，ｂ＊）が（９９．９９７，０．０００，０．００２）であることを示している。さらに（ｘ，ｙ）の値は、ウィンドウ２２の左側の色度図２２１上にプロットされている。この例では、エリアを２つ定義しているため、それぞれがプロットされている（点２２２−１，２２２−２）。
【００５０】
色度図２２１は、ユーザの利便性を考え、該色度図２２１中の可視領域である釣鐘２２３内に、コンピュータ１１で表すことのできる最も近い色をカラーで表示している（釣鐘２２３の右下が赤、左下が青、上が緑）。一方、グレースケールで使いたいユーザは、図１２のウィンドウ２３に示すように、色度図２３１を、単に釣鐘２３２の外側を黒、内側を灰色とし、さらに、（ｘ，ｙ）が正確に読み取れるように目盛りを付けたものとして表示する。この図１１のようなウィンドウ２２と図１２のようなウィンドウ２３とは、マウスの右クリックのメニューなどの指示により切り替えるようにすると便利である。また、ここでは図示していないが、カラー表示上で目盛り付きの色度図を選択させても良い。さらに、色度図２２１，２３１の部分的な情報を細かく見たい場合は、マウスの右クリックによるメニュー選択によって、図１３に示すように色度図２２１，２３１の部分的な拡大図が表示できると便利である。
【００５１】
なお、以上の方法では、エリアは、図２に示すようなバンド画像画面１５を使って予め指定しておいた。しかしながら、推定処理が高速に実行され、ユーザにストレスを感じさせない程度であるならば、マウスで１画素、或いは複数の画素を選択させて、その位置の色度値、色度図をリアルタイムで表示するようにしても良い。或いは、それらの混合で例えば、エリア１〜エリア５は予め指定しておき、エリア６は現在のマウス位置にして色度値、色度図を表示するというようにしても良い。さらにこの経過で、分光反射／透過率スペクトルも算出しているので、そのグラフを同時に同じウィンドウか別ウィンドウで表示することもできる。
【００５２】
以上のような方法によって、マルチバンドカメラ１で撮影した画像に対し、特定のエリアを指定し、その場所の色度値、色度図を表示することができる。エリアを複数設定した場合はその違いについて、色度図上でプロット位置を同時に表示して比較することができる効果がある。
【００５３】
［第３実施形態］
次に、本発明の第３実施形態を説明する。本実施形態は、マルチバンドカメラ１で撮影した画像に対し、特定のエリアを指定し、その場所を基準に他のエリアとの色差を表示するというものである。
【００５４】
本実施形態における構成は、上記第２実施形態と同様である。
【００５５】
本実施形態では、色差を画像として表示する。まず、基準となるエリアを指定して色度値を求める。その方法は、上記第２実施形態で説明したとおりである。例えば、図５（Ａ）に示すバンド画像画面１５で基準エリアを指定し、色度値を求める。その値を（Ｌｏ＊，ａｏ＊，ｂｏ＊）とする。一方、色差ΔＥは、公式より、
ΔＥ＝（（Ｌｏ＊−Ｌ１＊）２＋（ａｏ＊−ａ１＊）２＋（ｂｏ＊−ｂ１＊）２）１／２
で表される。（Ｌ１＊，ａ１＊，ｂ１＊）は、比較対象の場所の色度値である。
【００５６】
本実施形態では、エリア１を基準とした色差を、画面の全ピクセルに対して行い、その結果を画像として表す。今、左上の画素をＩ０とし、順にＩ２、Ｉ３、として表し、右下の画素をＩｎとする。もし、画像が２０４８×２０４８ピクセルであった時はｎ＝４１９４３０３である。この時の全ピクセルに対して色差を求めるので、その時の色差をΔＥ０〜ΔＥｎとして求める。次に、この中で最大の色差の値をΔＥｍａｘ、最小の色差の値をΔＥｍｉｎとすると、ｉ番目の画素で画像に表示するグレースケール上のピクセル値Ｄｉは次の式で表すことができる。
【００５７】
Ｄｉ＝（ΔＥｉ−ΔＥｍｉｎ）／（ΔＥｍａｘ−ΔＥｍｉｎ）×２５５
ただし、０≦ｉ≦ｎであり、ウィンドウ画像上で表すことのできる階調は０〜２５５であるとする。よって、元の画像ウィンドウの画素Ｉｉに対して、ピクセル値Ｄｉで置き換えると、基準エリアに対する色差の画像表示を行うことができる。
【００５８】
図１４にその例を示す。この例では、基準エリアを予め矩形として決めるのではなく、リアルタイムでマウスポインタ２４の位置の１画素或いは周囲の数画素で構成するエリアを基準エリアとしている。図１４（Ａ）で基準エリアをマウスで指定している様子を示している。図１４（Ｂ）ではその結果として、色差画像を表示している。
【００５９】
次に、その変形例として、色差ΔＥが指定した基準値、つまり、スレッシュホールドよりも小さい、或いは大きい場合にのみ、元の画像とは異なる色をつける方法について述べる。
【００６０】
基準のエリアを決めて画像像の全画素について色差ΔＥを求めるところは、前述した方法と全く同じである。そして、画素Ｉ０〜Ｉｎに対して求まった色差をΔＥ０〜ΔＥｎであったとする。ここで、スレッシュホールドとなる基準の色差ΔＥｓよりも各画素の色差ΔＥｉが大きいか小さいかを判定する。ここでは、基準値よりも小さい場合は、その画素を赤（（Ｒ，Ｇ，Ｂ）＝（２５５，０，０））で塗りつぶし、それ以上だった場合は何もせず、元の画像（画素）を表示するものとする。その結果を図１５に示す。この図で赤に塗りつぶされている部分（ハッチングして示す部分）は色差が基準値よりも小さく、基準エリアと色が近いことが分かる。同様に、あるスレッシュホールドよりも大きいところだけ色づけを行えば、基準エリアとの色差が大きい場所を調べることができる。
【００６１】
なお、この例でも今までと同様に、推定処理が速くユーザのストレスを感じさせない程度なら、マウスで１画素、或いは複数の画素を選択させて、色差を表した画像をリアルタイムで表示しても良い。
【００６２】
以上のような方法によって、マルチバンドカメラ１で撮影した画像に対し、特定のエリアを指定し、その場所を基準に他のエリアとの色差を表示することができる。また、基準エリア以外を画素単位にして色差を画像として表示することができる。これによってユーザは色差について、その大きさ、画像上での位置を明確にしながら比較、評価することができる効果がある。
【００６３】
以上、実施形態に基づいて説明したような本発明によって、マルチバンドカメラ１のデータを扱うソフトウエアのＧＵＩで、着目するエリアを１つ、或いは複数指定して、その場所の推定されたスペクトルを適時観察することができる。また、ここでのスペクトルとはマルチスペクトルカメラのＣＣＤ素子など入力手段側で得られた信号値（輝度）でも良いし、分光反射率、或いは分光透過率でも良い。ユーザは指定したエリアのスペクトルを観察するだけではなく、エリアどうしの相対値をグラフとしてみることもできる。
【００６４】
また、分光反射／透過率スペクトルが求まると、測色値、色度値を求めて、その値や色度図を表示することができる。特に、色度値では基準となるエリア（１画素でも良い）からの相対値を色差（ΔＥ）として求めて、指定よりも小さい色差、或いは大きい色差をグレースケールとは別の色を着色して画像表示することで、撮影された画像の評価をユーザは視覚的に行うことができる。特に、スペクトルの推定、測色値、色度値、色差のどのようなデータでもターゲット、基準エリアの場所を画像上でマウスで指定することにより、リアルタイムでこれらの値、グラフ、画像を観察することで、ユーザが行うマルチバンドカメラ１で撮られた画像の評価をやり易くする効果がある。
【００６５】
以上実施形態に基づいて本発明を説明したが、本発明は上述した実施形態に限定されるものではなく、本発明の要旨の範囲内で種々の変形や応用が可能なことは勿論である。
【００６６】
（付記）
前記の具体的実施形態から、以下のような構成の発明を抽出することができる。
【００６７】
（１） ２つ以上のバンドを撮影できるマルチバンドカメラで取得したマルチバンド画像を表示する画像表示装置において、
上記マルチバンド画像の内のあるバンドの画像を表示する画像表示手段と、
上記画像表示手段に表示された上記あるバンドの画像上で、少なくとも１画素以上を指定する画素指定手段と、
上記マルチバンド画像における上記画素指定手段によって指定された上記少なくとも１画素以上の信号値から、スペクトルを推定するスペクトル推定手段と、
上記スペクトル推定手段で推定したスペクトルを表示するスペクトル表示手段と、
を具備することを特徴とする画像表示装置。
【００６８】
（対応する実施形態）
この（１）に記載の画像表示装置に関する実施形態は、第１乃至第３実施形態が対応する。ここで、それらの実施形態において、コンピュータ１１が上記画像表示手段及び画素指定手段に対応し、スペクトル計算部１３が上記スペクトル推定手段に対応し、スペクトル出力用画面１４が上記スペクトル表示手段に対応する。
（作用効果）
この（１）に記載の画像表示装置によれば、画像の中でスペクトルを求めたい場所を特定し、グラフを表示することができる。
【００６９】
（２） 上記スペクトルは、信号値の輝度スペクトルであることを特徴とする（１）に記載の画像表示装置。
【００７０】
（対応する実施形態）
この（２）に記載の画像表示装置に関する実施形態は、第１実施形態が対応する。
（作用効果）
この（２）に記載の画像表示装置によれば、マルチバンド画像の中でスペクトルを求めたい場所の信号値の輝度スペクトルのグラフを表示することができる。
【００７１】
（３） 上記スペクトルは、分光スペクトルであることを特徴とする（１）に記載の画像表示装置。
【００７２】
（対応する実施形態）
この（３）に記載の画像表示装置に関する実施形態は、第１実施形態が対応する。
（作用効果）
この（３）に記載の画像表示装置によれば、マルチバンド画像の中でスペクトルを求めたい場所の信号値の分光反射／透過率スペクトルのグラフを表示することができる。
【００７３】
（４） 上記スペクトル推定手段は、上記マルチバンド画像の各バンド画像間の輝度データを補間することによってスペクトルを求めることを特徴とする（１）に記載の画像表示装置。
【００７４】
（対応する実施形態）
この（４）に記載の画像表示装置に関する実施形態は、第１実施形態が対応する。
（作用効果）
この（４）に記載の画像表示装置によれば、マルチバンド画像の各バンド画像間の輝度データを補間してスペクトルを求めることで、切れ目のないグラフを表示することができる。
【００７５】
（５） 上記補間は、線形補間であることを特徴とする（４）に記載の画像表示装置。
【００７６】
（対応する実施形態）
この（５）に記載の画像表示装置に関する実施形態は、第１実施形態が対応する。
【００７７】
（６） 上記補間は、スプライン曲線補間であることを特徴とする（４）に記載の画像表示装置。
【００７８】
（対応する実施形態）
この（６）に記載の画像表示装置に関する実施形態は、第１実施形態が対応する。
【００７９】
（７） 上記スペクトル推定手段は、上記マルチバンド画像に対応する参照画像を用い、上記マルチバンド画像及び上記参照画像の両者に対応した上記少なくとも１画素以上の信号値から、スペクトルを推定することを特徴とする（３）に記載の画像表示装置。
【００８０】
（対応する実施形態）
この（７）に記載の画像表示装置に関する実施形態は、第１実施形態が対応する。
【００８１】
（８） 上記マルチバンド画像の少なくとも１画素以上の信号値から、上記参照画像を作成する参照画像作成手段を更に具備することを特徴とする（７）に記載の画像表示装置。
【００８２】
（対応する実施形態）
この（８）に記載の画像表示装置に関する実施形態は、第１実施形態が対応する。ここで、その実施形態におけるコンピュータ１１が上記参照画像作成手段に対応する。
（作用効果）
この（８）に記載の画像表示装置によれば、マルチバンド画像から参照画像を作成するので、別途、参照画像を撮像する必要がない。
【００８３】
（９） 上記参照画像作成手段は、
上記画像表示手段に表示された上記あるバンドの画像上で、所望のエリアの形状と位置の指定を受けて、使用する画素を選択する画素選択手段と、
上記画素選択手段によって選択された使用する画素の値から上記参照画像を作成する作成手段と、
を有することを特徴とする（８）に記載の画像表示装置。
【００８４】
（対応する実施形態）
この（９）に記載の画像表示装置に関する実施形態は、第１実施形態が対応する。ここで、その実施形態におけるコンピュータ１１が上記画素選択手段及び作成手段に対応する。
（作用効果）
この（９）に記載の画像表示装置によれば、参照画像を作成する場合も指定したエリアから作成することができる。
【００８５】
（１０） 上記作成手段は、上記画素選択手段によって選択された使用する画素の値の平均値を元に上記参照画像を作成することを特徴とする（９）に記載の画像表示装置。
【００８６】
（対応する実施形態）
この（１０）に記載の画像表示装置に関する実施形態は、第１実施形態が対応する。
【００８７】
（１１） 上記参照画像作成手段は、上記マルチバンド画像の内のあるバンドの全画素から輝度値の高い画素の場所を自動的に求めて画素値を決め、その画素値を元に上記参照画像を作成することを特徴とする（８）に記載の画像表示装置。
【００８８】
（対応する実施形態）
この（１１）に記載の画像表示装置に関する実施形態は、第１実施形態が対応する。
（作用効果）
この（１１）に記載の画像表示装置によれば、ユーザの手補煩わすことなく自動的に参照画像を作成できる。
【００８９】
（１２） 上記参照画像作成手段は、上記輝度値の高い画素の場所を求める際に、輝度の高さのスレッシュホールドを設定し、それよりも輝度値の高い画素は露出オーバとして取り除き、上記スレッシュホールド以下の輝度値を持つ画素の中から上記輝度値の高い画素の場所を求めることを特徴とする（１１）に記載の画像表示装置。
【００９０】
（対応する実施形態）
この（１２）に記載の画像表示装置に関する実施形態は、第１実施形態が対応する。
（作用効果）
この（１２）に記載の画像表示装置によれば、不適切な画素の輝度値を除くことで、より適切な参照画像を作成することができる。
【００９１】
（１３） 上記参照画像作成手段は、指定した１つ以上の画素の輝度情報を全バンドにおいて一目で確認可能に提示する手段を有することを特徴とする（８）に記載の画像表示装置。
【００９２】
（対応する実施形態）
この（１３）に記載の画像表示装置に関する実施形態は、第１実施形態が対応する。ここで、その実施形態において、コンピュータ１１、スペクトル計算部１３、及びスペクトル出力用画面１４が上記一目で確認可能に提示する手段に対応する。
（作用効果）
この（１３）に記載の画像表示装置によれば、グラフを描画する対象のエリアは次々に変えることができるので、エリア毎のスペクトルの違いを一目にして区別することができるようになる。
【００９３】
（１４） 上記スペクトル推定手段で推定した分光スペクトルから、色度情報を求める色度情報算出手段を更に具備することを特徴とする（３）に記載の画像表示装置。
【００９４】
（対応する実施形態）
この（１４）に記載の画像表示装置に関する実施形態は、第２実施形態が対応する。ここで、その実施形態において、スペクトル計算部１３が上記色度情報算出手段に対応する。
（作用効果）
この（１４）に記載の画像表示装置によれば、マルチバンドカメラで撮影した画像に対し、特定のエリアを指定し、その場所の色度情報を表示することができる。また、エリアを複数設定した場合はその違いについて、色度図上でプロット位置を同時に表示して比較することができる効果がある。
【００９５】
（１５） 上記スペクトル表示手段は、上記色度情報算出手段で求めた上記色度情報の色度値（ｘ，ｙ）を色度図として表示することを特徴とする（１４）に記載の画像表示装置。
【００９６】
（対応する実施形態）
この（１５）に記載の画像表示装置に関する実施形態は、第２実施形態が対応する。
（作用効果）
この（１５）に記載の画像表示装置によれば、マルチバンドカメラで撮影した画像に対し、特定のエリアを指定し、その場所の色度値、色度図を表示することができる。
【００９７】
（１６） 上記スペクトル表示手段は、上記表示した色度図において、任意の部分の拡大表示機能を備えることを特徴とする（１５）に記載の画像表示装置。
【００９８】
（対応する実施形態）
この（１６）に記載の画像表示装置に関する実施形態は、第２実施形態が対応する。
（作用効果）
この（１６）に記載の画像表示装置によれば、色度図の部分的な情報を細かく見ることができる。
【００９９】
（１７） 上記色度情報算出手段は、同一の画像或いは異なる画像で２つ以上の色度情報を求め、それを比較する機能を備えることを特徴とする（１４）に記載の画像表示装置。
【０１００】
（対応する実施形態）
この（１７）に記載の画像表示装置に関する実施形態は、第３実施形態が対応する。
（作用効果）
この（１７）に記載の画像表示装置によれば、マルチバンドカメラで撮影した画像に対し、特定のエリアを指定し、その場所を基準に他のエリアとの色差を表示することができる。
【０１０１】
（１８） 上記色度情報算出手段における２つ以上の色度情報の比較方法は、色度値（ｘ，ｙ及びＬ＊，ａ＊，ｂ＊）の数値を比較することであることを特徴とする（１７）に記載の画像表示装置。
【０１０２】
（対応する実施形態）
この（１８）に記載の画像表示装置に関する実施形態は、第３実施形態が対応する。
（作用効果）
この（１８）に記載の画像表示装置によれば、ユーザは、色差について、その大きさを明確にしながら比較、評価することができる。
【０１０３】
（１９） 上記色度値の比較は、着目する１つ以上のある画素の色度値を基準値として、他の画素はその基準値との色差を求めるものであり、
上記スペクトル表示手段は、その求めた色差を画像で表現する、
ことを特徴とする（１８）に記載の画像表示装置。
【０１０４】
（対応する実施形態）
この（１９）に記載の画像表示装置に関する実施形態は、第３実施形態が対応する。
（作用効果）
この（１９）に記載の画像表示装置によれば、基準エリア以外を画素単位にして色差を画像として表示することができる。これによってユーザは、色差について、その大きさ、画像上での位置を明確にしながら比較、評価することができる。
【０１０５】
（２０） 上記画素指定手段は、
上記画像表示手段に表示された上記あるバンドの画像上の任意の位置にポインタを移動させる移動手段を有し、
該移動手段による上記ポインタの移動に合わせて、そのポインタが指し示す位置の画素を自動的に指定する、
ことを特徴とする（１）、（１２）、（１３）、（１４）、（１５）、（１７）、（１８）、及び（１９）の何れかに記載の画像表示装置。
【０１０６】
（対応する実施形態）
この（２０）に記載の画像表示装置に関する実施形態は、第１乃至第３実施形態が対応する。ここで、それらの実施形態において、コンピュータ１１が上記移動手段に対応する。
（作用効果）
この（２０）に記載の画像表示装置によれば、ポインタの移動に伴ってリアルタイムに、その位置のスペクトルを表示できる。
【０１０７】
（２１） ２つ以上のバンドを撮影できるマルチバンドカメラで取得したマルチバンド画像を表示する画像表示装置において、
上記マルチバンド画像の内のあるバンドの画像を表示する画像表示手段と、
上記画像表示手段に表示された上記あるバンドの画像から少なくとも２つ以上の指定画素を指定する画素指定手段と、
上記画素指定手段によって指定された指定画素各々のスペクトルを求めるスペクトル算出手段と、
上記画素指定手段によって指定された指定画素の内、基準となる１つの指定画素の上記スペクトル算出手段で求められたスペクトルを選択し、その選択されたスペクトルと他の指定画素の上記スペクトル算出手段で求められたスペクトルとの相対値を表示する表示手段と、
を具備することを特徴とする画像表示装置。
【０１０８】
（対応する実施形態）
この（２１）に記載の画像表示装置に関する実施形態は、第１実施形態が対応する。ここで、その実施形態において、コンピュータ１１が上記画像表示手段及び画素指定手段に対応し、スペクトル計算部１３が上記スペクトル算出手段に対応し、スペクトル出力用画面１４が上記表示手段に対応する。
（作用効果）
この（２１）に記載の画像表示装置によれば、画像の中でスペクトルを求めたいエリアを複数指定したとき、エリア毎のスペクトルの違いを一目にして区別することができる。
【０１０９】
（２２） ２つ以上のバンドを撮影できるマルチバンドカメラで取得したマルチバンド画像を表示する画像表示方法において、
上記マルチバンド画像の内のあるバンドの画像を表示する工程と、
その表示された上記あるバンドの画像上で、少なくとも１画素以上を指定する工程と、
上記マルチバンド画像における上記指定された上記少なくとも１画素以上の信号値から、スペクトルを推定する工程と、
その推定したスペクトルを表示する工程と、
を有することを特徴とする画像表示方法。
【０１１０】
（対応する実施形態）
この（２２）に記載の画像表示方法に関する実施形態は、第１乃至第３実施形態が対応する。
（作用効果）
この（２２）に記載の画像表示方法によれば、画像の中でスペクトルを求めたい場所を特定し、グラフを表示することができる。
【０１１１】
（２３） 上記スペクトルは、信号値の輝度スペクトルであることを特徴とする（２２）に記載の画像表示方法。
【０１１２】
（対応する実施形態）
この（２３）に記載の画像表示方法に関する実施形態は、第１実施形態が対応する。
（作用効果）
この（２３）に記載の画像表示方法によれば、マルチバンド画像の中でスペクトルを求めたい場所の信号値の輝度スペクトルのグラフを表示することができる。
【０１１３】
（２４） 上記スペクトルは、分光スペクトルであることを特徴とする（２２）に記載の画像表示方法。
【０１１４】
（対応する実施形態）
この（２４）に記載の画像表示方法に関する実施形態は、第１実施形態が対応する。
（作用効果）
この（２４）に記載の画像表示方法によれば、マルチバンド画像の中でスペクトルを求めたい場所の信号値の分光反射／透過率スペクトルのグラフを表示することができる。
【０１１５】
（２５） 上記スペクトルを推定する工程は、上記マルチバンド画像の各バンド画像間の輝度データを補間することによってスペクトルを求めることを特徴とする（２２）に記載の画像表示方法。
【０１１６】
（対応する実施形態）
この（２５）に記載の画像表示方法に関する実施形態は、第１実施形態が対応する。
（作用効果）
この（２５）に記載の画像表示方法によれば、マルチバンド画像の各バンド画像間の輝度データを補間してスペクトルを求めることで、切れ目のないグラフを表示することができる。
【０１１７】
（２６） 上記補間は、線形補間であることを特徴とする（２５）に記載の画像表示方法。
【０１１８】
（対応する実施形態）
この（２６）に記載の画像表示方法に関する実施形態は、第１実施形態が対応する。
【０１１９】
（２７） 上記補間は、スプライン曲線補間であることを特徴とする（２５）に記載の画像表示方法。
【０１２０】
（対応する実施形態）
この（２７）に記載の画像表示方法に関する実施形態は、第１実施形態が対応する。
【０１２１】
（２８） 上記スペクトルを推定する工程は、上記マルチバンド画像に対応する参照画像を用い、上記マルチバンド画像及び上記参照画像の両者に対応した上記少なくとも１画素以上の信号値から、スペクトルを推定することを特徴とする（２４）に記載の画像表示方法。
【０１２２】
（対応する実施形態）
この（２８）に記載の画像表示方法に関する実施形態は、第１実施形態が対応する。
【０１２３】
（２９） 上記マルチバンド画像の少なくとも１画素以上の信号値から、上記参照画像を作成する工程を更に具備することを特徴とする（２８）に記載の画像表示方法。
【０１２４】
（対応する実施形態）
この（２９）に記載の画像表示方法に関する実施形態は、第１実施形態が対応する。
（作用効果）
この（２９）に記載の画像表示方法によれば、マルチバンド画像から参照画像を作成するので、別途、参照画像を撮像する必要がない。
【０１２５】
（３０） 上記参照画像を作成する工程は、
上記表示された上記あるバンドの画像上で、所望のエリアの形状と位置の指定を受けて、使用する画素を選択する工程と、
その選択された使用する画素の値から上記参照画像を作成する工程と、
を有することを特徴とする（２９）に記載の画像表示方法。
【０１２６】
（対応する実施形態）
この（３０）に記載の画像表示方法に関する実施形態は、第１実施形態が対応する。
（作用効果）
この（３０）に記載の画像表示方法によれば、参照画像を作成する場合も指定したエリアから作成することができる。
【０１２７】
（３１） 上記参照画像を作成する工程は、上記選択された使用する画素の値の平均値を元に上記参照画像を作成することを特徴とする（３０）に記載の画像表示方法。
【０１２８】
（対応する実施形態）
この（３１）に記載の画像表示方法に関する実施形態は、第１実施形態が対応する。
【０１２９】
（３２） 上記参照画像を作成する工程は、上記マルチバンド画像の内のあるバンドの全画素から輝度値の高い画素の場所を自動的に求めて画素値を決め、その画素値を元に上記参照画像を作成することを特徴とする（２９）に記載の画像表示方法。
【０１３０】
（対応する実施形態）
この（３２）に記載の画像表示方法に関する実施形態は、第１実施形態が対応する。
（作用効果）
この（３２）に記載の画像表示方法によれば、ユーザの手補煩わすことなく自動的に参照画像を作成できる。
【０１３１】
（３３） 上記参照画像を作成する工程は、上記輝度値の高い画素の場所を求める際に、輝度の高さのスレッシュホールドを設定し、それよりも輝度値の高い画素は露出オーバとして取り除き、上記スレッシュホールド以下の輝度値を持つ画素の中から上記輝度値の高い画素の場所を求めることを特徴とする（３２）に記載の画像表示方法。
【０１３２】
（対応する実施形態）
この（３３）に記載の画像表示方法に関する実施形態は、第１実施形態が対応する。
（作用効果）
この（３３）に記載の画像表示方法によれば、不適切な画素の輝度値を除くことで、より適切な参照画像を作成することができる。
【０１３３】
（３４） 上記参照画像を作成する工程は、指定した１つ以上の画素の輝度情報を全バンドにおいて一目で確認可能に提示する工程を有することを特徴とする（２９）に記載の画像表示方法。
【０１３４】
（対応する実施形態）
この（３４）に記載の画像表示方法に関する実施形態は、第１実施形態が対応する。
（作用効果）
この（３４）に記載の画像表示方法によれば、グラフを描画する対象のエリアは次々に変えることができるので、エリア毎のスペクトルの違いを一目にして区別することができるようになる。
【０１３５】
（３５） 上記スペクトルを推定する工程で推定した分光スペクトルから、色度情報を求める工程を更に有することを特徴とする（２４）に記載の画像表示方法。
【０１３６】
（対応する実施形態）
この（３５）に記載の画像表示方法に関する実施形態は、第２実施形態が対応する。
（作用効果）
この（３５）に記載の画像表示方法によれば、マルチバンドカメラで撮影した画像に対し、特定のエリアを指定し、その場所の色度情報を表示することができる。また、エリアを複数設定した場合はその違いについて、色度図上でプロット位置を同時に表示して比較することができる効果がある。
【０１３７】
（３６） 上記色度情報を求める工程で求めた上記色度情報の色度値（ｘ，ｙ）を色度図として表示することを特徴とする（３５）に記載の画像表示方法。
【０１３８】
（対応する実施形態）
この（３６）に記載の画像表示方法に関する実施形態は、第２実施形態が対応する。
（作用効果）
この（３６）に記載の画像表示方法によれば、マルチバンドカメラで撮影した画像に対し、特定のエリアを指定し、その場所の色度値、色度図を表示することができる。
【０１３９】
（３７） 上記表示した色度図において、任意の部分を拡大表示することを特徴とする（３６）に記載の画像表示方法。
【０１４０】
（対応する実施形態）
この（３７）に記載の画像表示方法に関する実施形態は、第２実施形態が対応する。
（作用効果）
この（３７）に記載の画像表示方法によれば、色度図の部分的な情報を細かく見ることができる。
【０１４１】
（３８） 上記色度情報を求める工程は、同一の画像或いは異なる画像で２つ以上の色度情報を求め、それを比較することを特徴とする（３５）に記載の画像表示方法。
【０１４２】
（対応する実施形態）
この（３８）に記載の画像表示方法に関する実施形態は、第３実施形態が対応する。
（作用効果）
この（３８）に記載の画像表示方法によれば、マルチバンドカメラで撮影した画像に対し、特定のエリアを指定し、その場所を基準に他のエリアとの色差を表示することができる。
【０１４３】
（３９） 上記色度情報を求める工程における２つ以上の色度情報の比較方法は、色度値（ｘ，ｙ及びＬ＊，ａ＊，ｂ＊）の数値を比較することであることを特徴とする（３８）に記載の画像表示方法。
【０１４４】
（対応する実施形態）
この（３９）に記載の画像表示方法に関する実施形態は、第３実施形態が対応する。
（作用効果）
この（３９）に記載の画像表示方法によれば、ユーザは、色差について、その大きさを明確にしながら比較、評価することができる。
【０１４５】
（４０） 上記色度値の比較は、着目する１つ以上のある画素の色度値を基準値として、他の画素はその基準値との色差を求めるものであり、
その求めた色差を画像で表現する、
ことを特徴とする（３９）に記載の画像表示方法。
【０１４６】
（対応する実施形態）
この（４０）に記載の画像表示方法に関する実施形態は、第３実施形態が対応する。
（作用効果）
この（４０）に記載の画像表示方法によれば、基準エリア以外を画素単位にして色差を画像として表示することができる。これによってユーザは、色差について、その大きさ、画像上での位置を明確にしながら比較、評価することができる。
【０１４７】
（４１） 上記少なくとも１画素以上を指定する工程は、
上記表示された上記あるバンドの画像上の任意の位置にポインタを移動させ、
そのポインタの移動に合わせて、そのポインタが指し示す位置の画素を自動的に指定する、
ことを特徴とする（２２）、（３３）、（３４）、（３５）、（３６）、（３８）、（３９）、及び（４０）の何れかに記載の画像表示方法。
【０１４８】
（対応する実施形態）
この（４１）に記載の画像表示方法に関する実施形態は、第１乃至第３実施形態が対応する。
（作用効果）
この（４１）に記載の画像表示方法によれば、ポインタの移動に伴ってリアルタイムに、その位置のスペクトルを表示できる。
【０１４９】
（４２） ２つ以上のバンドを撮影できるマルチバンドカメラで取得したマルチバンド画像を表示する画像表示方法において、
上記マルチバンド画像の内のあるバンドの画像を表示する工程と、
その表示された上記あるバンドの画像から少なくとも２つ以上の指定画素を指定する工程と、
その指定された指定画素各々のスペクトルを求める工程と、
上記指定された指定画素の内、基準となる１つの指定画素の上記求められたスペクトルを選択し、その選択されたスペクトルと他の指定画素の上記スペクトルとの相対値（比）を表示する工程と、
を有することを特徴とする画像表示方法。
【０１５０】
（対応する実施形態）
この（４２）に記載の画像表示方法に関する実施形態は、第１実施形態が対応する。
（作用効果）
この（４２）に記載の画像表示方法によれば、画像の中でスペクトルを求めたいエリアを複数指定したとき、エリア毎のスペクトルの違いを一目にして区別することができる。
【０１５１】
（４３） ２つ以上のバンドを撮影できるマルチバンドカメラで取得したマルチバンド画像を表示する画像表示装置において、
上記マルチバンド画像の内のあるバンドの画像を表示する画像表示手段と、
上記画像表示手段に表示された上記あるバンドの画像上で、少なくとも１画素以上を指定する画素指定手段と、
上記マルチバンド画像における上記画素指定手段によって指定された上記少なくとも１画素以上の信号値から、分光スペクトルを推定するスペクトル推定手段と、
上記スペクトル推定手段で推定した分光スペクトルから、色度情報を求める色度情報算出手段と、
上記スペクトル推定手段で推定した分光スペクトルを表示すると共に、上記色度情報算出手段で求めた上記色度情報の色度値を色度図として表示するスペクトル表示手段と、
を具備し、
上記スペクトル表示手段は、上記表示した色度図において、任意の部分の拡大表示機能を備えることを特徴とする画像表示装置。
【０１５２】
（対応する実施形態）
この（４３）に記載の画像表示装置に関する実施形態は、第２実施形態が対応する。ここで、その実施形態において、コンピュータ１１が上記画像表示手段及び画素指定手段に対応し、スペクトル計算部１３が上記スペクトル推定手段及び上記色度情報算出手段に対応し、コンピュータ、スペクトル計算部１３及びスペクトル出力用画面１４が上記スペクトル表示手段に対応する。
（作用効果）
この（４３）に記載の画像表示装置によれば、画像の中でスペクトルを求めたい場所を特定し、グラフを表示することができる。また、マルチバンド画像の中でスペクトルを求めたい場所の信号値の分光反射／透過率スペクトルのグラフを表示することができる。更に、マルチバンドカメラで撮影した画像に対し、特定のエリアを指定し、その場所の色度情報を表示することができる。そして、エリアを複数設定した場合はその違いについて、色度図上でプロット位置を同時に表示して比較することができる効果がある。更に、マルチバンドカメラで撮影した画像に対し、特定のエリアを指定し、その場所の色度値、色度図を表示することができる。そして、色度図において、任意の部分の拡大表示が行えるので、色度図の部分的な情報を細かく見ることができる。
【０１５３】
（４４） ２つ以上のバンドを撮影できるマルチバンドカメラで取得したマルチバンド画像を表示する画像表示方法において、
上記マルチバンド画像の内のあるバンドの画像を表示する工程と、
その表示された上記あるバンドの画像上で、少なくとも１画素以上を指定する工程と、
上記マルチバンド画像における上記指定された上記少なくとも１画素以上の信号値から、分光スペクトルを推定する工程と、
その推定した分光スペクトルから、色度情報を求める工程と、
上記推定した分光スペクトルを表示すると共に、上記求めた色度情報の色度値を色度図として表示する工程と、
を有し、
上記表示する工程は、上記表示した色度図において、任意の部分の拡大表示を行う工程を含むことを特徴とする画像表示方法。
【０１５４】
（対応する実施形態）
この（４４）に記載の画像表示方法に関する実施形態は、第２実施形態が対応する。
（作用効果）
この（４４）に記載の画像表示方法によれば、画像の中でスペクトルを求めたい場所を特定し、グラフを表示することができる。また、マルチバンド画像の中でスペクトルを求めたい場所の信号値の分光反射／透過率スペクトルのグラフを表示することができる。更に、マルチバンドカメラで撮影した画像に対し、特定のエリアを指定し、その場所の色度情報を表示することができる。そして、エリアを複数設定した場合はその違いについて、色度図上でプロット位置を同時に表示して比較することができる効果がある。更に、マルチバンドカメラで撮影した画像に対し、特定のエリアを指定し、その場所の色度値、色度図を表示することができる。そして、色度図において、任意の部分の拡大表示が行えるので、色度図の部分的な情報を細かく見ることができる。
【図面の簡単な説明】
【０１５５】
【図１】（Ａ）は本発明の第１実施形態に係る画像表示装置の構成を示す図であり、（Ｂ）はスペクトル出力用画面を示す図である。
【図２】バンド画像画面を示す図である。
【図３】補間方法選択画面を示す図である。
【図４】第１実施形態に係る画像表示装置の変形例の構成を示す図である。
【図５】（Ａ）は色票を撮影した場合のバンド画像を表示したバンド画像画面を示す図であり、（Ｂ）は作成した参照画像を表示した参照画像画面を示す図である。
【図６】バンド選択画面を示す図である。
【図７】変形例における補間方法選択画面を示す図である。
【図８】第１実施形態に係る画像表示装置の別の変形例の構成を示す図である。
【図９】エリア１を基準とした信号値（輝度）の相対値を表示するグラフを示す図である。
【図１０】第１実施形態に係る画像表示装置の更に別の変形例の構成を示す図である。
【図１１】本発明の第２実施形態に係る画像表示装置における表示例を示す図である。
【図１２】第２実施形態に係る画像表示装置における別の表示例を示す図である。
【図１３】第２実施形態に係る画像表示装置における更に別の表示例を示す図である。
【図１４】（Ａ）は本発明の第３実施形態に係る画像表示装置において基準エリアをマウスで指定している様子を示す図であり、（Ｂ）はその結果として表示される色差画像を示す図である。
【図１５】第３実施形態に係る画像表示装置における別の表示例を示す図である。
【符号の説明】
【０１５６】
１…マルチバンドカメラ、 １１…コンピュータ、 １２…マルチバンド画像、 １３…スペクトル計算部、 １４…スペクトル出力用画面、 １５…バンド画像画面、 １６…補間方法選択画面、 １７…参照画像、 １８…参照画像画面、 １９…バンド選択画面、 ２０…補間方法選択画面、 ２１…特性ファイル、 ２２，２３…ウィンドウ、 ２４…マウスポインタ、 １４１…エリア選択リストボックス、 １５１…矩形、 １９１…バンド番号、 ２２１…色度図、 ２２２−１，２２２−２…点、 ２２３，２３２…釣鐘、 ２３１…色度図。
【技術分野】
【０００１】
本発明は、マルチバンドカメラにより取得されたマルチバンド画像からスペクトルグラフ表示、色度値表示、及び、複数の対象点どうしの相互の相対値を表示する画像表示装置及び画像表示方法に関する。
【背景技術】
【０００２】
マルチバンドカメラとは一般に、通常の３バンドカラーのカメラのＲ，Ｇ，Ｂフィルタと比べて、より特定の波長透過特性を持つ複数のバンド数で撮影するカメラを指す。バンドを切り替える手段は、例えば、特許文献１に開示されているような波長可変フィルタであったり、特許文献２に開示されているような回転フィルタを用いる。装置、目的によってバンド数は異なるが、各バンドでの透過率特性は例えば上記特許文献１の図２に示されるようになっている。
【０００３】
撮影した画像は一般にバンド数だけ存在するので、上記特許文献１の図２では１６バンドによる１６枚の画像になり、上記特許文献２の図８に示す回転フィルタを用いた場合には１０バンドによる１０枚の画像になる。
【０００４】
マルチバンドカメラで撮像されたマルチバンド画像は、バンド毎に画像として表示することができる。
【０００５】
上記特許文献１においても、また上記特許文献２でも、スペクトルを推定する方法については述べられているが、実際にソフトウエアから操作する方法については触れられていない。
【特許文献１】特開２００１−０９９７１０号公報
【特許文献２】特開平１１−０９６３３３号公報
【発明の開示】
【発明が解決しようとする課題】
【０００６】
スペクトルの表示は一般に、撮影した画像の全体ではなく、着目している特定の場所を対象に行う。その為には、指定は任意の場所をソフトウエアから指定し、その場所だけでスペクトルを表示できた方が良い。
【０００７】
また、３８０ｎｍから７８０ｎｍを１ｎｍおきに撮影し、４０１バンドの画像を撮った場合は、１ｎｍずつの４０１バンドの画像から精度良くスペクトルを推定できるが、４０１バンドよりも少ない、即ち１ｎｍよりも粗いサンプリング波長間隔で画像を撮った場合は、サンプリング波長間の広い、少ない画像データから４０１バンドの場合と略同等のスペクトルを推定する必要がある。推定方法によっても、スペクトルが異なることが考えられるため、任意に推定方法を選ぶことができたほうが良い。
【０００８】
さらに、エリアを複数指定した場合や複数の画像からスペクトルを表示する場合は、スペクトル同士の比較も望まれる。
【０００９】
スペクトルが求められると、それを利用し測色値を計算することができる。ここでもやはり、エリアを複数定義している場合は、ある基準に対してΔＥなど色差を計算し、表示することができると、ユーザの画像に対する理解を深めることができる。さらに、基準のエリアから、画像の各ピクセルがどの位の色差（ΔＥ）になっているか、一目でわかるような表示方法が望まれる。
【００１０】
本発明は、上記の点に鑑みてなされたもので、マルチバンドカメラから取得されたマルチバンド画像に対し、スペクトルの表示や色度値の表現等を行うためのＧＵＩ表示方法を提案し、ユーザの画像評価の負荷を軽減することが可能な画像表示装置及び画像表示方法を提供することを目的とする。
【課題を解決するための手段】
【００１１】
上記の目的を達成するために、本発明の第１の態様による画像表示装置は、
２つ以上のバンドを撮影できるマルチバンドカメラで取得したマルチバンド画像を表示する画像表示装置において、
上記マルチバンド画像の内のあるバンドの画像を表示する画像表示手段と、
上記画像表示手段に表示された上記あるバンドの画像上で、少なくとも１画素以上を指定する画素指定手段と、
上記マルチバンド画像における上記画素指定手段によって指定された上記少なくとも１画素以上の信号値から、分光スペクトルを推定するスペクトル推定手段と、
上記スペクトル推定手段で推定した分光スペクトルから、色度情報を求める色度情報算出手段と、
上記スペクトル推定手段で推定した分光スペクトルを表示すると共に、上記色度情報算出手段で求めた上記色度情報の色度値を色度図として表示するスペクトル表示手段と、
を具備し、
上記スペクトル表示手段は、上記表示した色度図において、任意の部分の拡大表示機能を備えることを特徴とする。
【００１２】
また、上記の目的を達成するために、本発明の第２の態様による画像表示方法は、
２つ以上のバンドを撮影できるマルチバンドカメラで取得したマルチバンド画像を表示する画像表示方法において、
上記マルチバンド画像の内のあるバンドの画像を表示する工程と、
その表示された上記あるバンドの画像上で、少なくとも１画素以上を指定する工程と、
上記マルチバンド画像における上記指定された上記少なくとも１画素以上の信号値から、分光スペクトルを推定する工程と、
その推定した分光スペクトルから、色度情報を求める工程と、
上記推定した分光スペクトルを表示すると共に、上記求めた色度情報の色度値を色度図として表示する工程と、
を有し、
上記表示する工程は、上記表示した色度図において、任意の部分の拡大表示を行う工程を含むことを特徴とする。
【発明の効果】
【００１３】
本発明によれば、マルチバンドカメラから取得されたマルチバンド画像に対し、スペクトルの表示や色度値の表現等を行うためのＧＵＩ表示方法を提案し、ユーザの画像評価の負荷を軽減することが可能な画像表示装置及び画像表示方法を提供することができる。
【発明を実施するための最良の形態】
【００１４】
以下、本発明を実施するための最良の形態を図面を参照して説明する。
【００１５】
［第１実施形態］
まず、本発明の第１実施形態として、画像上の任意の場所を対象とした、信号レベルのスペクトルを表示する方法を説明する。
【００１６】
本実施形態に係る画像表示装置の構成を示す図１（Ａ）において、マルチバンドカメラ１は、コンピュータ１１によって制御されている。マルチバンドカメラ１は一般に、撮像素子であるＣＣＤを内蔵し、画像を撮影して電子データとして画像データをコンピュータ１１に送る。撮影された画像は、マルチバンド画像１２としてコンピュータ１１の管理するメモリ上にある。なお、画像はその都度撮影されても良いし、過去に撮影した画像をハードディスク等の記憶媒体から読み出してきてメモリ上にマルチバンド画像１２として展開されていても良い。従って、本実施形態の処理を行う時は必ずしもマルチバンドカメラ１が無くても良い。
【００１７】
コンピュータ１１は、上記マルチバンド画像１２から、スペクトル計算部１３によってスペクトルを計算し、スペクトル出力用画面１４にグラフ表示を行う。本実施形態では、コンピュータ１１はＩＢＭ社のＰＣ−ＡＴ互換機を用い、ＯＳとしてＷｉｎｄｏｗｓ（登録商標）２０００を利用する。スペクトル出力用画面１４はＷｉｎｄｏｗｓ２０００の管理するウィンドウ画面であり、その画面はＣＲＴなどのモニタに表示される。
【００１８】
マルチバンドカメラ１は１８バンドとすると、上記マルチバンド画像１２は１８枚のバンド画像でなり、コンピュータ１１はユーザの指定に従ったバンドの画像を表示する。例えば、バンド６の画像が図２のように表示する。ユーザは、このバンド画像画面１５を見ながら、スペクトルを表示したい場所を特定する。この図２の例では、コンピュータ１１のポインティングデバイスであるマウスのドラッグ＆ドロップ機能（矩形の左上、右下を指定する方法）を使って、矩形を５つ指定している。なお、エリアの数や形はユーザによって異なる形でも良く、この実施形態の範疇に入る。
【００１９】
そしてユーザが、図示していない「信号値のスペクトル表示」のメニューを選択すると、コンピュータ１１は、図３のようなＧＵＩ（グラフィカルユーザインターフェース）、即ち補間方法選択画面１６を表示して、補間の方法を選択させる。
【００２０】
このような補間方法選択画面１６で、直線補間が選択された場合には、コンピュータ１１は、スペクトル計算部１３により、それぞれの矩形内の画素値、つまり、これはマルチバンドカメラ１のＣＣＤ素子の信号値、の平均値を算出する。そして、その輝度をバンド毎にプロットして直線で結んだグラフをスペクトル出力用画面１４として、図１（Ｂ）に示すように表示する。なお、このグラフの縦軸は信号の強度である。ＣＣＤ素子の出力が１６ｂｉｔであるときは、信号強度（輝度）として、０から６５５３５までの値をとることができる。また、横軸は波長を表している。この例では、バンド１が３８０ｎｍ、バンド２が４１８ｎｍの中心波長であり、同様にバンド１８まで中心波長を表す波長での輝度をプロットしていく。そして、選択された補間方法は、この例では直線補間であったので、バンド間は直線で補間する。なお、該スペクトル出力用画面１４の右上には、エリア選択リストボックス１４１が配されており、このエリア選択リストボックス１４１で選択することで、上記バンド画像画面１５で特定された各エリアごとの信号強度のスペクトルを、適宜切り替えながら観察することができるようになっている。
【００２１】
また、ユーザが図３の補間方法選択画面１６で、スプライン補正を選択したときには、補正間隔を１０ｎｍ、５ｎｍ、１ｎｍのうちから選択して、直線補間のときと同様に、バンド１〜バンド１８でのｘ座標、ｙ座標を求めてからスプライン補正を行う。
【００２２】
以上のようにして、画像上の任意の場所について、信号レベルのスペクトルを表示することができる。
【００２３】
画像上の任意の場所について、分光スペクトルを表示する方法を、本実施形態の変形例として説明する。
【００２４】
この場合は、図４に示すように、図１（Ａ）と同様の構成に加えて、参照画像１７を利用する。また、スペクトル計算部１３が分光スペクトルを計算するものとなっている。
【００２５】
即ち、分光スペクトルは、通常の物体に対しては分光反射率、顕微鏡等のように透過物体を観察する場合は分光透過率のスペクトルとしてグラフ表示を行う。従って、反射率か透過率を計算するための基準となる画像、つまり参照画像１７が必要となる。ここで、反射率スペクトルを求める場合は、参照画像１７の対象は標準白色板を撮影したものになる。また、透過率スペクトルを求める場合の参照画像１７は、観察するものを取り除き、照明光だけの画像を用いる。
【００２６】
つまり、参照画像１７は、マルチバンド画像１２を撮影する時と同じ手順で、撮影の対象を標準白色板などに代えて撮影する。従って、撮影は、マルチバンド画像１２の撮影の他に、そのような標準白色板などの撮影を１回、余分に行うことが必要となる。しかしながら、マルチバンド画像１２の中に、参照する場所を一緒に写し込むことができれば、マルチバンド撮影する撮影対象の撮影のみで良い。
【００２７】
以下に、その方法を説明する。図５（Ａ）は、撮影対象を撮影したマルチバンド画像１２の内、ユーザに指定された或いは初期画像として設定されたバンドの画像を表示したバンド画像画面１５を示す図である。この例では、色票が撮影対象になっている。もし、この色票の中に、反射率が極めて１００％に近いものがあった場合は、これを標準白色板と同様な扱いと定義することができる。この例では、図５（Ａ）の左下の矩形がその場所であり、ユーザはマウスのドラッグ＆ドロップ操作で、その場所を指定する。或いは、コンピュータ１１で自動的にその場所を認識しても良い。この例では、指定された場所が矩形１５１で図示されている。コンピュータ１１は、その矩形１５１の信号値（輝度）の平均値を使って参照画像１７を作成する。なお、こうして作成された参照画像１７を、図５（Ｂ）に示すような参照画像画面１８に表示して、ユーザに確認を受けることができる。
【００２８】
なお、参照画像１７は、光源に対して反射率が１００％の対象物が理想であるため、信号値（輝度）の明るいところを探すためには、図６のようなＧＵＩ、即ちバンド選択画面１９を表示し、マウス位置の１つ、或いは複数の画素をエリアとして定義し、その画素の信号値（輝度）を表示すると、評価を行い易い。ここで、図６の中央左側の「Ｌｅｖｅｌ」以下はバンド番号１９１を表し、バンド１〜バンド１８が表示されている。「Ｌｅｖｅｌ」の右側の数値「（１９７４，０）」は、現在のバンド画像画面１５に表示された画像上のマウス位置を表す。その下の数字で信号値を表示している。例えば、バンド１では「１７６（１／２５６）」と表示されているが、これは、１６ｂｉｔデータとして信号値が１７６であることを示している。また括弧内の数値は、コンピュータ１１のディスプレイは０〜２５５の２５６階調で輝度を表すことが多いため、この２５６階調換算で値が幾つになるかを示している。この場合は１である。
【００２９】
上記スペクトル計算部１３では、分光反射／透過率を計算して、スペクトル出力用画面１４へ表示を行う。その方法を手順に従って説明する。対象のマルチバンド画像１２に対して、分光反射／透過率のスペクトルを表示する対象の場所を指定する。図２はその様子を表し、入力方法は前述した手順で行う。この画像には参照画像１７が対になっている。作成方法は前述したように対象画像の一部から作るか、参照板だけの画像を改めて撮影する。
【００３０】
次に反射率／透過率を求める。この例では色票を使っているため、反射率である。反射率は次の式で求めることができる。
【００３１】
対象画像の指定バンドの反射率＝（対象画像の指定バンドの指定エリアの平均値）／（参照画像の指定バンドの指定エリアの平均値）
図２の例では、バンド数は１８なので、バンド毎に値を求める。また、指定エリアが５つあるので、それぞれのエリアで反射率を求める。例えば、図２の左下のエリア１は白色板であったとして、この場所から参照画像１７を前述した方法で作成した場合は、全てのバンドにおいてエリア１の反射率は１になる。それ以外のエリアでは一般に参照画像１７よりも信号値は小さいはずなので、０≦反射率≦１．０を示すのが普通である。
【００３２】
そして、求めた反射率をスペクトル出力用画面１４にグラフとして表示する。なお、この場合も、ユーザが図７に表すようなＧＵＩ、即ち補間方法選択画面２０で補間の方法を選択できるようにする。図７に示す状態では、直線補間が選ばれている。グラフにプロットする方法は、前述の信号レベルのスペクトルをグラフ表示させた時と同様な手順になる。つまり、図１（Ｂ）のグラフで横軸に各バンドで用いたフィルタの中心波長の位置を波長換算で求め、縦軸に反射率として値をプロットする。従って、前述の場合は縦軸は信号値であったが、分光反射率では、一般に１．０を最大値とした反射率で描画を行う。直線補間では、プロットの間を前述したように直線で補間し、スペクトルグラフを完成させる。
【００３３】
また、図７の補間方法選択画面２０で、補間方法としてスプライン補正を選択したときは、補間を１０ｎｍ、５ｎｍ、１ｎｍで行うか選択させた後、スプライン曲線で近似を行えば良い。
【００３４】
なお、図５（Ａ）及び（Ｂ）を参照して、ユーザが目で確認しながら、参照画像１７を作成する方法を説明したが、それ以外に参照画像１７を作成する変形例について述べる。
【００３５】
参照板は照明光を理想的には１００％反射するようにしたものであるため、図５（Ａ）に示すバンド画像画面１５で表示されるバンド画像内に参照板に相当する場所があれば、その場所の信号輝度は他と比べて高いはずである。従って、図５（Ａ）のようなバンド画像において全ての画素の信号値（輝度）を調べて、明るい順に並べる。その中で明るい方から例えば、３０画素分を参照板位置での画素情報として、その３０画素の信号値の平均値から図５（Ｂ）の参照画像画面１８に示すような参照画像１７を作成することができる。
【００３６】
もう１つの変形例として、露出オーバの部分を除いて、参照画像１７を作成する方法を説明する。図５（Ａ）に示すバンド画像画面１５で表示されるようなバンド画像は、適度な露出時間を設定して撮影したつもりでも、部分的に非常に明るい部分が存在してしまう可能性がある。例えば、参照板では参照板自身の材質、照明位置が調整されているので、正反射部分が除かれている場合が多いが、正式な参照板が用意できなかった場合や、照明条件を適切にできなかった場合は、参照としたい画像の部分に正反射部分を含むことが考えられる。この場合は、前述した方法で、バンド画像の全ての画素の信号値（輝度）を調べて、明るい順に並べるが、その後で、一定の輝度以上のものを取り除く。例えば、信号値の最大値が２５５であった場合は、並べられた画素のうち、２４０以上の信号値は正反射部分として除外する。そして、残った画素の中から、明るい輝度の３０画素分を同様な方法で参照画像１７を作成することができる。なお、この例ではスレッシュホールドを２４０としたが、マルチバンド画像１２の入力手段、信号値の幅、撮影物、ユーザの指定で他の値にしても同様な方法が成り立つ。
【００３７】
次に、上記特許文献２に開示されているようなスペクトル推定方法で、分光スペクトルの推定を行う例について述べる。
【００３８】
図８は、その場合の構成を示す図で、マルチバンドカメラ１、コンピュータ１１、マルチバンド画像１２、及び、スペクトル出力用画面１４の構成と機能は図１（Ａ）と同じである。分光スペクトルの推定は、スペクトル計算部１３で行う。ここでの推定方法は、上記特許文献２で記述されている。また、その推定の為に必要なファイルが特性ファイル２１として用意され、その特性ファイル２１をスペクトル計算部１３が参照する。
【００３９】
次に手順を示す。画像は今までと同様に、図２のような形で表示され、スペクトルを計算したい場所を特定する。図２では５つのエリアが指定されている。そして、その後、ユーザが分光反射／透過率のスペクトル表示メニューを選択したとき、図７の補間方法選択画面２０を表示し、補間方法として、３番目の「推定」が選択されると、先に選択した５つのエリアで上記特許文献２の方法でスペクトルの推定を行う。その結果は、図１（Ｂ）のようなグラフで表示される。ただし、この場合は縦軸は反射率であり、一般的には最大値を１．０とした０〜１．０の値になる。
【００４０】
なお、以上のいずれの方法でもエリアは、図２に示すようなバンド画像画面１５を使って予め指定しておいた。しかしながら、推定処理が高速に実行され、ユーザにストレスを感じさせない程度であるならば、マウスで１画素、或いは複数の画素を選択させて、その位置のスペクトルをリアルタイムで計算して表示するようにしても良い。或いは、それらの混合で例えば、エリア１〜エリア５は予め指定しておき、エリア６は現在のマウス位置にしてスペクトルを表示するというようにしても良い。
【００４１】
また、エリアが複数ある場合のスペクトルの比較は、図１（Ｂ）のようにグラフをエリア毎に表示しエリアの選択を変えることで行っていたが、図１（Ｂ）のウィンドウをエリア毎に作成し、それらを並べてタイル表示として、比較できるようにしても良い。或いは、１つのウィンドウにエリア毎に色を変えるなどの工夫をして同時に表示しても良い。さらに、相対的に信号値や分光反射／透過率のスペクトルを比較する場合は、図９に示すように、あるエリアを基準とした相対値を表示しても良い。この図では、エリア１を基準とした信号値（輝度）の相対値を縦軸として、グラフ表示を行っている。
【００４２】
また、図１０に示すように、図４の構成と図８の構成を組み合わせ、適宜ユーザが選択できるようにしても良い。
【００４３】
以上のような方法によって、画像の中でスペクトルを求めたい場所を特定し、グラフを表示することができる。推定方法が複数である場合もそれらを選択することが可能である。また、参照画像１７を作成する場合も、指定したエリアから作成することができる。以上のような方法によって、スペクトルを表示したい場所を画像を見ながら指定し、そのエリアでのスペクトルを表示することができる。また、グラフを描画する対象のエリアは次々に変えることができるので、エリア毎のスペクトルの違いを一目にして区別することができる効果がある。
【００４４】
［第２実施形態］
次に、本発明の第２実施形態を説明する。本実施形態は、マルチバンドカメラ１で撮影したマルチバンド画像に対し、特定のエリアを指定し、その場所の色度値、色度図を表示するというものである。
【００４５】
色度値は、分光反射率／透過率のスペクトルから推定するため、構成は上記第１実施形態と同様であり、スペクトル計算部１３により上記第１実施形態で説明した方法で分光反射率／透過率のスペクトルを計算した後、後述する方法で色度値が求められる。そして結果は、スペクトル出力用画面１４と同様に同じウィンドウ、或いは別ウィンドウで画面に表示される。色度値は前述したように分光反射率／透過率のスペクトルから求められるため、その推定方法によって、図１（Ａ）、図４、図８、または図１０の構成が考えられる。
【００４６】
次に、色度値ＸＹＺ、ｘｙｚ、Ｌ＊，ａ＊，ｂ＊を表示する方法を説明する。
【００４７】
まず、測色したい部分を指定して、その場所の分光反射／透過率スペクトルを推定する。その方法は、上記第１実施形態で説明した通りである。その結果、推定された分光反射／透過率スペクトルをＲ（λ）であったとする。この形態ではλは３８０［ｎｍ］〜７８０［ｎｍ］まで、１［ｎｍ］間隔で求められたとする。この時の照明光をＳ（λ）とすると、
Ｘ＝ＫΣＳ（λ）ｘ（λ）Ｒ（λ）
Ｙ＝ＫΣＳ（λ）ｙ（λ）Ｒ（λ）
Ｚ＝ＫΣＳ（λ）ｚ（λ）Ｒ（λ）
Ｋ＝１００／（ΣＳ（λ）ｙ（λ））
で表される。ただし、ｘ（λ）、ｙ（λ）、ｚ（λ）はＸＹＺ表色系における等色関数の値であり、Σはλを３８０〜７８０とする。ここで求まったＸＹＺより、ｘｙｚは次の式で求めることができる。
【００４８】
ｘ＝Ｘ／（Ｘ＋Ｙ＋Ｚ）
ｙ＝Ｙ／（Ｘ＋Ｙ＋Ｚ）
ｚ＝Ｚ／（Ｘ＋Ｙ＋Ｚ）
さらに、レンダリング照明光をｅ（λ）として、Ｌ＊，ａ＊，ｂ＊を求めると、次の公式で表すことができる。
【数１】
【００４９】
以上の公式から求めた値を、図１１に示すようなウィンドウ２２に表示する。このウィンドウ２２の右側では、エリア１に対して、（ｘ，ｙ）が（０．３１３，０．３２９）、Ｙが９９．９９３、（Ｌ＊，ａ＊，ｂ＊）が（９９．９９７，０．０００，０．００２）であることを示している。さらに（ｘ，ｙ）の値は、ウィンドウ２２の左側の色度図２２１上にプロットされている。この例では、エリアを２つ定義しているため、それぞれがプロットされている（点２２２−１，２２２−２）。
【００５０】
色度図２２１は、ユーザの利便性を考え、該色度図２２１中の可視領域である釣鐘２２３内に、コンピュータ１１で表すことのできる最も近い色をカラーで表示している（釣鐘２２３の右下が赤、左下が青、上が緑）。一方、グレースケールで使いたいユーザは、図１２のウィンドウ２３に示すように、色度図２３１を、単に釣鐘２３２の外側を黒、内側を灰色とし、さらに、（ｘ，ｙ）が正確に読み取れるように目盛りを付けたものとして表示する。この図１１のようなウィンドウ２２と図１２のようなウィンドウ２３とは、マウスの右クリックのメニューなどの指示により切り替えるようにすると便利である。また、ここでは図示していないが、カラー表示上で目盛り付きの色度図を選択させても良い。さらに、色度図２２１，２３１の部分的な情報を細かく見たい場合は、マウスの右クリックによるメニュー選択によって、図１３に示すように色度図２２１，２３１の部分的な拡大図が表示できると便利である。
【００５１】
なお、以上の方法では、エリアは、図２に示すようなバンド画像画面１５を使って予め指定しておいた。しかしながら、推定処理が高速に実行され、ユーザにストレスを感じさせない程度であるならば、マウスで１画素、或いは複数の画素を選択させて、その位置の色度値、色度図をリアルタイムで表示するようにしても良い。或いは、それらの混合で例えば、エリア１〜エリア５は予め指定しておき、エリア６は現在のマウス位置にして色度値、色度図を表示するというようにしても良い。さらにこの経過で、分光反射／透過率スペクトルも算出しているので、そのグラフを同時に同じウィンドウか別ウィンドウで表示することもできる。
【００５２】
以上のような方法によって、マルチバンドカメラ１で撮影した画像に対し、特定のエリアを指定し、その場所の色度値、色度図を表示することができる。エリアを複数設定した場合はその違いについて、色度図上でプロット位置を同時に表示して比較することができる効果がある。
【００５３】
［第３実施形態］
次に、本発明の第３実施形態を説明する。本実施形態は、マルチバンドカメラ１で撮影した画像に対し、特定のエリアを指定し、その場所を基準に他のエリアとの色差を表示するというものである。
【００５４】
本実施形態における構成は、上記第２実施形態と同様である。
【００５５】
本実施形態では、色差を画像として表示する。まず、基準となるエリアを指定して色度値を求める。その方法は、上記第２実施形態で説明したとおりである。例えば、図５（Ａ）に示すバンド画像画面１５で基準エリアを指定し、色度値を求める。その値を（Ｌｏ＊，ａｏ＊，ｂｏ＊）とする。一方、色差ΔＥは、公式より、
ΔＥ＝（（Ｌｏ＊−Ｌ１＊）２＋（ａｏ＊−ａ１＊）２＋（ｂｏ＊−ｂ１＊）２）１／２
で表される。（Ｌ１＊，ａ１＊，ｂ１＊）は、比較対象の場所の色度値である。
【００５６】
本実施形態では、エリア１を基準とした色差を、画面の全ピクセルに対して行い、その結果を画像として表す。今、左上の画素をＩ０とし、順にＩ２、Ｉ３、として表し、右下の画素をＩｎとする。もし、画像が２０４８×２０４８ピクセルであった時はｎ＝４１９４３０３である。この時の全ピクセルに対して色差を求めるので、その時の色差をΔＥ０〜ΔＥｎとして求める。次に、この中で最大の色差の値をΔＥｍａｘ、最小の色差の値をΔＥｍｉｎとすると、ｉ番目の画素で画像に表示するグレースケール上のピクセル値Ｄｉは次の式で表すことができる。
【００５７】
Ｄｉ＝（ΔＥｉ−ΔＥｍｉｎ）／（ΔＥｍａｘ−ΔＥｍｉｎ）×２５５
ただし、０≦ｉ≦ｎであり、ウィンドウ画像上で表すことのできる階調は０〜２５５であるとする。よって、元の画像ウィンドウの画素Ｉｉに対して、ピクセル値Ｄｉで置き換えると、基準エリアに対する色差の画像表示を行うことができる。
【００５８】
図１４にその例を示す。この例では、基準エリアを予め矩形として決めるのではなく、リアルタイムでマウスポインタ２４の位置の１画素或いは周囲の数画素で構成するエリアを基準エリアとしている。図１４（Ａ）で基準エリアをマウスで指定している様子を示している。図１４（Ｂ）ではその結果として、色差画像を表示している。
【００５９】
次に、その変形例として、色差ΔＥが指定した基準値、つまり、スレッシュホールドよりも小さい、或いは大きい場合にのみ、元の画像とは異なる色をつける方法について述べる。
【００６０】
基準のエリアを決めて画像像の全画素について色差ΔＥを求めるところは、前述した方法と全く同じである。そして、画素Ｉ０〜Ｉｎに対して求まった色差をΔＥ０〜ΔＥｎであったとする。ここで、スレッシュホールドとなる基準の色差ΔＥｓよりも各画素の色差ΔＥｉが大きいか小さいかを判定する。ここでは、基準値よりも小さい場合は、その画素を赤（（Ｒ，Ｇ，Ｂ）＝（２５５，０，０））で塗りつぶし、それ以上だった場合は何もせず、元の画像（画素）を表示するものとする。その結果を図１５に示す。この図で赤に塗りつぶされている部分（ハッチングして示す部分）は色差が基準値よりも小さく、基準エリアと色が近いことが分かる。同様に、あるスレッシュホールドよりも大きいところだけ色づけを行えば、基準エリアとの色差が大きい場所を調べることができる。
【００６１】
なお、この例でも今までと同様に、推定処理が速くユーザのストレスを感じさせない程度なら、マウスで１画素、或いは複数の画素を選択させて、色差を表した画像をリアルタイムで表示しても良い。
【００６２】
以上のような方法によって、マルチバンドカメラ１で撮影した画像に対し、特定のエリアを指定し、その場所を基準に他のエリアとの色差を表示することができる。また、基準エリア以外を画素単位にして色差を画像として表示することができる。これによってユーザは色差について、その大きさ、画像上での位置を明確にしながら比較、評価することができる効果がある。
【００６３】
以上、実施形態に基づいて説明したような本発明によって、マルチバンドカメラ１のデータを扱うソフトウエアのＧＵＩで、着目するエリアを１つ、或いは複数指定して、その場所の推定されたスペクトルを適時観察することができる。また、ここでのスペクトルとはマルチスペクトルカメラのＣＣＤ素子など入力手段側で得られた信号値（輝度）でも良いし、分光反射率、或いは分光透過率でも良い。ユーザは指定したエリアのスペクトルを観察するだけではなく、エリアどうしの相対値をグラフとしてみることもできる。
【００６４】
また、分光反射／透過率スペクトルが求まると、測色値、色度値を求めて、その値や色度図を表示することができる。特に、色度値では基準となるエリア（１画素でも良い）からの相対値を色差（ΔＥ）として求めて、指定よりも小さい色差、或いは大きい色差をグレースケールとは別の色を着色して画像表示することで、撮影された画像の評価をユーザは視覚的に行うことができる。特に、スペクトルの推定、測色値、色度値、色差のどのようなデータでもターゲット、基準エリアの場所を画像上でマウスで指定することにより、リアルタイムでこれらの値、グラフ、画像を観察することで、ユーザが行うマルチバンドカメラ１で撮られた画像の評価をやり易くする効果がある。
【００６５】
以上実施形態に基づいて本発明を説明したが、本発明は上述した実施形態に限定されるものではなく、本発明の要旨の範囲内で種々の変形や応用が可能なことは勿論である。
【００６６】
（付記）
前記の具体的実施形態から、以下のような構成の発明を抽出することができる。
【００６７】
（１） ２つ以上のバンドを撮影できるマルチバンドカメラで取得したマルチバンド画像を表示する画像表示装置において、
上記マルチバンド画像の内のあるバンドの画像を表示する画像表示手段と、
上記画像表示手段に表示された上記あるバンドの画像上で、少なくとも１画素以上を指定する画素指定手段と、
上記マルチバンド画像における上記画素指定手段によって指定された上記少なくとも１画素以上の信号値から、スペクトルを推定するスペクトル推定手段と、
上記スペクトル推定手段で推定したスペクトルを表示するスペクトル表示手段と、
を具備することを特徴とする画像表示装置。
【００６８】
（対応する実施形態）
この（１）に記載の画像表示装置に関する実施形態は、第１乃至第３実施形態が対応する。ここで、それらの実施形態において、コンピュータ１１が上記画像表示手段及び画素指定手段に対応し、スペクトル計算部１３が上記スペクトル推定手段に対応し、スペクトル出力用画面１４が上記スペクトル表示手段に対応する。
（作用効果）
この（１）に記載の画像表示装置によれば、画像の中でスペクトルを求めたい場所を特定し、グラフを表示することができる。
【００６９】
（２） 上記スペクトルは、信号値の輝度スペクトルであることを特徴とする（１）に記載の画像表示装置。
【００７０】
（対応する実施形態）
この（２）に記載の画像表示装置に関する実施形態は、第１実施形態が対応する。
（作用効果）
この（２）に記載の画像表示装置によれば、マルチバンド画像の中でスペクトルを求めたい場所の信号値の輝度スペクトルのグラフを表示することができる。
【００７１】
（３） 上記スペクトルは、分光スペクトルであることを特徴とする（１）に記載の画像表示装置。
【００７２】
（対応する実施形態）
この（３）に記載の画像表示装置に関する実施形態は、第１実施形態が対応する。
（作用効果）
この（３）に記載の画像表示装置によれば、マルチバンド画像の中でスペクトルを求めたい場所の信号値の分光反射／透過率スペクトルのグラフを表示することができる。
【００７３】
（４） 上記スペクトル推定手段は、上記マルチバンド画像の各バンド画像間の輝度データを補間することによってスペクトルを求めることを特徴とする（１）に記載の画像表示装置。
【００７４】
（対応する実施形態）
この（４）に記載の画像表示装置に関する実施形態は、第１実施形態が対応する。
（作用効果）
この（４）に記載の画像表示装置によれば、マルチバンド画像の各バンド画像間の輝度データを補間してスペクトルを求めることで、切れ目のないグラフを表示することができる。
【００７５】
（５） 上記補間は、線形補間であることを特徴とする（４）に記載の画像表示装置。
【００７６】
（対応する実施形態）
この（５）に記載の画像表示装置に関する実施形態は、第１実施形態が対応する。
【００７７】
（６） 上記補間は、スプライン曲線補間であることを特徴とする（４）に記載の画像表示装置。
【００７８】
（対応する実施形態）
この（６）に記載の画像表示装置に関する実施形態は、第１実施形態が対応する。
【００７９】
（７） 上記スペクトル推定手段は、上記マルチバンド画像に対応する参照画像を用い、上記マルチバンド画像及び上記参照画像の両者に対応した上記少なくとも１画素以上の信号値から、スペクトルを推定することを特徴とする（３）に記載の画像表示装置。
【００８０】
（対応する実施形態）
この（７）に記載の画像表示装置に関する実施形態は、第１実施形態が対応する。
【００８１】
（８） 上記マルチバンド画像の少なくとも１画素以上の信号値から、上記参照画像を作成する参照画像作成手段を更に具備することを特徴とする（７）に記載の画像表示装置。
【００８２】
（対応する実施形態）
この（８）に記載の画像表示装置に関する実施形態は、第１実施形態が対応する。ここで、その実施形態におけるコンピュータ１１が上記参照画像作成手段に対応する。
（作用効果）
この（８）に記載の画像表示装置によれば、マルチバンド画像から参照画像を作成するので、別途、参照画像を撮像する必要がない。
【００８３】
（９） 上記参照画像作成手段は、
上記画像表示手段に表示された上記あるバンドの画像上で、所望のエリアの形状と位置の指定を受けて、使用する画素を選択する画素選択手段と、
上記画素選択手段によって選択された使用する画素の値から上記参照画像を作成する作成手段と、
を有することを特徴とする（８）に記載の画像表示装置。
【００８４】
（対応する実施形態）
この（９）に記載の画像表示装置に関する実施形態は、第１実施形態が対応する。ここで、その実施形態におけるコンピュータ１１が上記画素選択手段及び作成手段に対応する。
（作用効果）
この（９）に記載の画像表示装置によれば、参照画像を作成する場合も指定したエリアから作成することができる。
【００８５】
（１０） 上記作成手段は、上記画素選択手段によって選択された使用する画素の値の平均値を元に上記参照画像を作成することを特徴とする（９）に記載の画像表示装置。
【００８６】
（対応する実施形態）
この（１０）に記載の画像表示装置に関する実施形態は、第１実施形態が対応する。
【００８７】
（１１） 上記参照画像作成手段は、上記マルチバンド画像の内のあるバンドの全画素から輝度値の高い画素の場所を自動的に求めて画素値を決め、その画素値を元に上記参照画像を作成することを特徴とする（８）に記載の画像表示装置。
【００８８】
（対応する実施形態）
この（１１）に記載の画像表示装置に関する実施形態は、第１実施形態が対応する。
（作用効果）
この（１１）に記載の画像表示装置によれば、ユーザの手補煩わすことなく自動的に参照画像を作成できる。
【００８９】
（１２） 上記参照画像作成手段は、上記輝度値の高い画素の場所を求める際に、輝度の高さのスレッシュホールドを設定し、それよりも輝度値の高い画素は露出オーバとして取り除き、上記スレッシュホールド以下の輝度値を持つ画素の中から上記輝度値の高い画素の場所を求めることを特徴とする（１１）に記載の画像表示装置。
【００９０】
（対応する実施形態）
この（１２）に記載の画像表示装置に関する実施形態は、第１実施形態が対応する。
（作用効果）
この（１２）に記載の画像表示装置によれば、不適切な画素の輝度値を除くことで、より適切な参照画像を作成することができる。
【００９１】
（１３） 上記参照画像作成手段は、指定した１つ以上の画素の輝度情報を全バンドにおいて一目で確認可能に提示する手段を有することを特徴とする（８）に記載の画像表示装置。
【００９２】
（対応する実施形態）
この（１３）に記載の画像表示装置に関する実施形態は、第１実施形態が対応する。ここで、その実施形態において、コンピュータ１１、スペクトル計算部１３、及びスペクトル出力用画面１４が上記一目で確認可能に提示する手段に対応する。
（作用効果）
この（１３）に記載の画像表示装置によれば、グラフを描画する対象のエリアは次々に変えることができるので、エリア毎のスペクトルの違いを一目にして区別することができるようになる。
【００９３】
（１４） 上記スペクトル推定手段で推定した分光スペクトルから、色度情報を求める色度情報算出手段を更に具備することを特徴とする（３）に記載の画像表示装置。
【００９４】
（対応する実施形態）
この（１４）に記載の画像表示装置に関する実施形態は、第２実施形態が対応する。ここで、その実施形態において、スペクトル計算部１３が上記色度情報算出手段に対応する。
（作用効果）
この（１４）に記載の画像表示装置によれば、マルチバンドカメラで撮影した画像に対し、特定のエリアを指定し、その場所の色度情報を表示することができる。また、エリアを複数設定した場合はその違いについて、色度図上でプロット位置を同時に表示して比較することができる効果がある。
【００９５】
（１５） 上記スペクトル表示手段は、上記色度情報算出手段で求めた上記色度情報の色度値（ｘ，ｙ）を色度図として表示することを特徴とする（１４）に記載の画像表示装置。
【００９６】
（対応する実施形態）
この（１５）に記載の画像表示装置に関する実施形態は、第２実施形態が対応する。
（作用効果）
この（１５）に記載の画像表示装置によれば、マルチバンドカメラで撮影した画像に対し、特定のエリアを指定し、その場所の色度値、色度図を表示することができる。
【００９７】
（１６） 上記スペクトル表示手段は、上記表示した色度図において、任意の部分の拡大表示機能を備えることを特徴とする（１５）に記載の画像表示装置。
【００９８】
（対応する実施形態）
この（１６）に記載の画像表示装置に関する実施形態は、第２実施形態が対応する。
（作用効果）
この（１６）に記載の画像表示装置によれば、色度図の部分的な情報を細かく見ることができる。
【００９９】
（１７） 上記色度情報算出手段は、同一の画像或いは異なる画像で２つ以上の色度情報を求め、それを比較する機能を備えることを特徴とする（１４）に記載の画像表示装置。
【０１００】
（対応する実施形態）
この（１７）に記載の画像表示装置に関する実施形態は、第３実施形態が対応する。
（作用効果）
この（１７）に記載の画像表示装置によれば、マルチバンドカメラで撮影した画像に対し、特定のエリアを指定し、その場所を基準に他のエリアとの色差を表示することができる。
【０１０１】
（１８） 上記色度情報算出手段における２つ以上の色度情報の比較方法は、色度値（ｘ，ｙ及びＬ＊，ａ＊，ｂ＊）の数値を比較することであることを特徴とする（１７）に記載の画像表示装置。
【０１０２】
（対応する実施形態）
この（１８）に記載の画像表示装置に関する実施形態は、第３実施形態が対応する。
（作用効果）
この（１８）に記載の画像表示装置によれば、ユーザは、色差について、その大きさを明確にしながら比較、評価することができる。
【０１０３】
（１９） 上記色度値の比較は、着目する１つ以上のある画素の色度値を基準値として、他の画素はその基準値との色差を求めるものであり、
上記スペクトル表示手段は、その求めた色差を画像で表現する、
ことを特徴とする（１８）に記載の画像表示装置。
【０１０４】
（対応する実施形態）
この（１９）に記載の画像表示装置に関する実施形態は、第３実施形態が対応する。
（作用効果）
この（１９）に記載の画像表示装置によれば、基準エリア以外を画素単位にして色差を画像として表示することができる。これによってユーザは、色差について、その大きさ、画像上での位置を明確にしながら比較、評価することができる。
【０１０５】
（２０） 上記画素指定手段は、
上記画像表示手段に表示された上記あるバンドの画像上の任意の位置にポインタを移動させる移動手段を有し、
該移動手段による上記ポインタの移動に合わせて、そのポインタが指し示す位置の画素を自動的に指定する、
ことを特徴とする（１）、（１２）、（１３）、（１４）、（１５）、（１７）、（１８）、及び（１９）の何れかに記載の画像表示装置。
【０１０６】
（対応する実施形態）
この（２０）に記載の画像表示装置に関する実施形態は、第１乃至第３実施形態が対応する。ここで、それらの実施形態において、コンピュータ１１が上記移動手段に対応する。
（作用効果）
この（２０）に記載の画像表示装置によれば、ポインタの移動に伴ってリアルタイムに、その位置のスペクトルを表示できる。
【０１０７】
（２１） ２つ以上のバンドを撮影できるマルチバンドカメラで取得したマルチバンド画像を表示する画像表示装置において、
上記マルチバンド画像の内のあるバンドの画像を表示する画像表示手段と、
上記画像表示手段に表示された上記あるバンドの画像から少なくとも２つ以上の指定画素を指定する画素指定手段と、
上記画素指定手段によって指定された指定画素各々のスペクトルを求めるスペクトル算出手段と、
上記画素指定手段によって指定された指定画素の内、基準となる１つの指定画素の上記スペクトル算出手段で求められたスペクトルを選択し、その選択されたスペクトルと他の指定画素の上記スペクトル算出手段で求められたスペクトルとの相対値を表示する表示手段と、
を具備することを特徴とする画像表示装置。
【０１０８】
（対応する実施形態）
この（２１）に記載の画像表示装置に関する実施形態は、第１実施形態が対応する。ここで、その実施形態において、コンピュータ１１が上記画像表示手段及び画素指定手段に対応し、スペクトル計算部１３が上記スペクトル算出手段に対応し、スペクトル出力用画面１４が上記表示手段に対応する。
（作用効果）
この（２１）に記載の画像表示装置によれば、画像の中でスペクトルを求めたいエリアを複数指定したとき、エリア毎のスペクトルの違いを一目にして区別することができる。
【０１０９】
（２２） ２つ以上のバンドを撮影できるマルチバンドカメラで取得したマルチバンド画像を表示する画像表示方法において、
上記マルチバンド画像の内のあるバンドの画像を表示する工程と、
その表示された上記あるバンドの画像上で、少なくとも１画素以上を指定する工程と、
上記マルチバンド画像における上記指定された上記少なくとも１画素以上の信号値から、スペクトルを推定する工程と、
その推定したスペクトルを表示する工程と、
を有することを特徴とする画像表示方法。
【０１１０】
（対応する実施形態）
この（２２）に記載の画像表示方法に関する実施形態は、第１乃至第３実施形態が対応する。
（作用効果）
この（２２）に記載の画像表示方法によれば、画像の中でスペクトルを求めたい場所を特定し、グラフを表示することができる。
【０１１１】
（２３） 上記スペクトルは、信号値の輝度スペクトルであることを特徴とする（２２）に記載の画像表示方法。
【０１１２】
（対応する実施形態）
この（２３）に記載の画像表示方法に関する実施形態は、第１実施形態が対応する。
（作用効果）
この（２３）に記載の画像表示方法によれば、マルチバンド画像の中でスペクトルを求めたい場所の信号値の輝度スペクトルのグラフを表示することができる。
【０１１３】
（２４） 上記スペクトルは、分光スペクトルであることを特徴とする（２２）に記載の画像表示方法。
【０１１４】
（対応する実施形態）
この（２４）に記載の画像表示方法に関する実施形態は、第１実施形態が対応する。
（作用効果）
この（２４）に記載の画像表示方法によれば、マルチバンド画像の中でスペクトルを求めたい場所の信号値の分光反射／透過率スペクトルのグラフを表示することができる。
【０１１５】
（２５） 上記スペクトルを推定する工程は、上記マルチバンド画像の各バンド画像間の輝度データを補間することによってスペクトルを求めることを特徴とする（２２）に記載の画像表示方法。
【０１１６】
（対応する実施形態）
この（２５）に記載の画像表示方法に関する実施形態は、第１実施形態が対応する。
（作用効果）
この（２５）に記載の画像表示方法によれば、マルチバンド画像の各バンド画像間の輝度データを補間してスペクトルを求めることで、切れ目のないグラフを表示することができる。
【０１１７】
（２６） 上記補間は、線形補間であることを特徴とする（２５）に記載の画像表示方法。
【０１１８】
（対応する実施形態）
この（２６）に記載の画像表示方法に関する実施形態は、第１実施形態が対応する。
【０１１９】
（２７） 上記補間は、スプライン曲線補間であることを特徴とする（２５）に記載の画像表示方法。
【０１２０】
（対応する実施形態）
この（２７）に記載の画像表示方法に関する実施形態は、第１実施形態が対応する。
【０１２１】
（２８） 上記スペクトルを推定する工程は、上記マルチバンド画像に対応する参照画像を用い、上記マルチバンド画像及び上記参照画像の両者に対応した上記少なくとも１画素以上の信号値から、スペクトルを推定することを特徴とする（２４）に記載の画像表示方法。
【０１２２】
（対応する実施形態）
この（２８）に記載の画像表示方法に関する実施形態は、第１実施形態が対応する。
【０１２３】
（２９） 上記マルチバンド画像の少なくとも１画素以上の信号値から、上記参照画像を作成する工程を更に具備することを特徴とする（２８）に記載の画像表示方法。
【０１２４】
（対応する実施形態）
この（２９）に記載の画像表示方法に関する実施形態は、第１実施形態が対応する。
（作用効果）
この（２９）に記載の画像表示方法によれば、マルチバンド画像から参照画像を作成するので、別途、参照画像を撮像する必要がない。
【０１２５】
（３０） 上記参照画像を作成する工程は、
上記表示された上記あるバンドの画像上で、所望のエリアの形状と位置の指定を受けて、使用する画素を選択する工程と、
その選択された使用する画素の値から上記参照画像を作成する工程と、
を有することを特徴とする（２９）に記載の画像表示方法。
【０１２６】
（対応する実施形態）
この（３０）に記載の画像表示方法に関する実施形態は、第１実施形態が対応する。
（作用効果）
この（３０）に記載の画像表示方法によれば、参照画像を作成する場合も指定したエリアから作成することができる。
【０１２７】
（３１） 上記参照画像を作成する工程は、上記選択された使用する画素の値の平均値を元に上記参照画像を作成することを特徴とする（３０）に記載の画像表示方法。
【０１２８】
（対応する実施形態）
この（３１）に記載の画像表示方法に関する実施形態は、第１実施形態が対応する。
【０１２９】
（３２） 上記参照画像を作成する工程は、上記マルチバンド画像の内のあるバンドの全画素から輝度値の高い画素の場所を自動的に求めて画素値を決め、その画素値を元に上記参照画像を作成することを特徴とする（２９）に記載の画像表示方法。
【０１３０】
（対応する実施形態）
この（３２）に記載の画像表示方法に関する実施形態は、第１実施形態が対応する。
（作用効果）
この（３２）に記載の画像表示方法によれば、ユーザの手補煩わすことなく自動的に参照画像を作成できる。
【０１３１】
（３３） 上記参照画像を作成する工程は、上記輝度値の高い画素の場所を求める際に、輝度の高さのスレッシュホールドを設定し、それよりも輝度値の高い画素は露出オーバとして取り除き、上記スレッシュホールド以下の輝度値を持つ画素の中から上記輝度値の高い画素の場所を求めることを特徴とする（３２）に記載の画像表示方法。
【０１３２】
（対応する実施形態）
この（３３）に記載の画像表示方法に関する実施形態は、第１実施形態が対応する。
（作用効果）
この（３３）に記載の画像表示方法によれば、不適切な画素の輝度値を除くことで、より適切な参照画像を作成することができる。
【０１３３】
（３４） 上記参照画像を作成する工程は、指定した１つ以上の画素の輝度情報を全バンドにおいて一目で確認可能に提示する工程を有することを特徴とする（２９）に記載の画像表示方法。
【０１３４】
（対応する実施形態）
この（３４）に記載の画像表示方法に関する実施形態は、第１実施形態が対応する。
（作用効果）
この（３４）に記載の画像表示方法によれば、グラフを描画する対象のエリアは次々に変えることができるので、エリア毎のスペクトルの違いを一目にして区別することができるようになる。
【０１３５】
（３５） 上記スペクトルを推定する工程で推定した分光スペクトルから、色度情報を求める工程を更に有することを特徴とする（２４）に記載の画像表示方法。
【０１３６】
（対応する実施形態）
この（３５）に記載の画像表示方法に関する実施形態は、第２実施形態が対応する。
（作用効果）
この（３５）に記載の画像表示方法によれば、マルチバンドカメラで撮影した画像に対し、特定のエリアを指定し、その場所の色度情報を表示することができる。また、エリアを複数設定した場合はその違いについて、色度図上でプロット位置を同時に表示して比較することができる効果がある。
【０１３７】
（３６） 上記色度情報を求める工程で求めた上記色度情報の色度値（ｘ，ｙ）を色度図として表示することを特徴とする（３５）に記載の画像表示方法。
【０１３８】
（対応する実施形態）
この（３６）に記載の画像表示方法に関する実施形態は、第２実施形態が対応する。
（作用効果）
この（３６）に記載の画像表示方法によれば、マルチバンドカメラで撮影した画像に対し、特定のエリアを指定し、その場所の色度値、色度図を表示することができる。
【０１３９】
（３７） 上記表示した色度図において、任意の部分を拡大表示することを特徴とする（３６）に記載の画像表示方法。
【０１４０】
（対応する実施形態）
この（３７）に記載の画像表示方法に関する実施形態は、第２実施形態が対応する。
（作用効果）
この（３７）に記載の画像表示方法によれば、色度図の部分的な情報を細かく見ることができる。
【０１４１】
（３８） 上記色度情報を求める工程は、同一の画像或いは異なる画像で２つ以上の色度情報を求め、それを比較することを特徴とする（３５）に記載の画像表示方法。
【０１４２】
（対応する実施形態）
この（３８）に記載の画像表示方法に関する実施形態は、第３実施形態が対応する。
（作用効果）
この（３８）に記載の画像表示方法によれば、マルチバンドカメラで撮影した画像に対し、特定のエリアを指定し、その場所を基準に他のエリアとの色差を表示することができる。
【０１４３】
（３９） 上記色度情報を求める工程における２つ以上の色度情報の比較方法は、色度値（ｘ，ｙ及びＬ＊，ａ＊，ｂ＊）の数値を比較することであることを特徴とする（３８）に記載の画像表示方法。
【０１４４】
（対応する実施形態）
この（３９）に記載の画像表示方法に関する実施形態は、第３実施形態が対応する。
（作用効果）
この（３９）に記載の画像表示方法によれば、ユーザは、色差について、その大きさを明確にしながら比較、評価することができる。
【０１４５】
（４０） 上記色度値の比較は、着目する１つ以上のある画素の色度値を基準値として、他の画素はその基準値との色差を求めるものであり、
その求めた色差を画像で表現する、
ことを特徴とする（３９）に記載の画像表示方法。
【０１４６】
（対応する実施形態）
この（４０）に記載の画像表示方法に関する実施形態は、第３実施形態が対応する。
（作用効果）
この（４０）に記載の画像表示方法によれば、基準エリア以外を画素単位にして色差を画像として表示することができる。これによってユーザは、色差について、その大きさ、画像上での位置を明確にしながら比較、評価することができる。
【０１４７】
（４１） 上記少なくとも１画素以上を指定する工程は、
上記表示された上記あるバンドの画像上の任意の位置にポインタを移動させ、
そのポインタの移動に合わせて、そのポインタが指し示す位置の画素を自動的に指定する、
ことを特徴とする（２２）、（３３）、（３４）、（３５）、（３６）、（３８）、（３９）、及び（４０）の何れかに記載の画像表示方法。
【０１４８】
（対応する実施形態）
この（４１）に記載の画像表示方法に関する実施形態は、第１乃至第３実施形態が対応する。
（作用効果）
この（４１）に記載の画像表示方法によれば、ポインタの移動に伴ってリアルタイムに、その位置のスペクトルを表示できる。
【０１４９】
（４２） ２つ以上のバンドを撮影できるマルチバンドカメラで取得したマルチバンド画像を表示する画像表示方法において、
上記マルチバンド画像の内のあるバンドの画像を表示する工程と、
その表示された上記あるバンドの画像から少なくとも２つ以上の指定画素を指定する工程と、
その指定された指定画素各々のスペクトルを求める工程と、
上記指定された指定画素の内、基準となる１つの指定画素の上記求められたスペクトルを選択し、その選択されたスペクトルと他の指定画素の上記スペクトルとの相対値（比）を表示する工程と、
を有することを特徴とする画像表示方法。
【０１５０】
（対応する実施形態）
この（４２）に記載の画像表示方法に関する実施形態は、第１実施形態が対応する。
（作用効果）
この（４２）に記載の画像表示方法によれば、画像の中でスペクトルを求めたいエリアを複数指定したとき、エリア毎のスペクトルの違いを一目にして区別することができる。
【０１５１】
（４３） ２つ以上のバンドを撮影できるマルチバンドカメラで取得したマルチバンド画像を表示する画像表示装置において、
上記マルチバンド画像の内のあるバンドの画像を表示する画像表示手段と、
上記画像表示手段に表示された上記あるバンドの画像上で、少なくとも１画素以上を指定する画素指定手段と、
上記マルチバンド画像における上記画素指定手段によって指定された上記少なくとも１画素以上の信号値から、分光スペクトルを推定するスペクトル推定手段と、
上記スペクトル推定手段で推定した分光スペクトルから、色度情報を求める色度情報算出手段と、
上記スペクトル推定手段で推定した分光スペクトルを表示すると共に、上記色度情報算出手段で求めた上記色度情報の色度値を色度図として表示するスペクトル表示手段と、
を具備し、
上記スペクトル表示手段は、上記表示した色度図において、任意の部分の拡大表示機能を備えることを特徴とする画像表示装置。
【０１５２】
（対応する実施形態）
この（４３）に記載の画像表示装置に関する実施形態は、第２実施形態が対応する。ここで、その実施形態において、コンピュータ１１が上記画像表示手段及び画素指定手段に対応し、スペクトル計算部１３が上記スペクトル推定手段及び上記色度情報算出手段に対応し、コンピュータ、スペクトル計算部１３及びスペクトル出力用画面１４が上記スペクトル表示手段に対応する。
（作用効果）
この（４３）に記載の画像表示装置によれば、画像の中でスペクトルを求めたい場所を特定し、グラフを表示することができる。また、マルチバンド画像の中でスペクトルを求めたい場所の信号値の分光反射／透過率スペクトルのグラフを表示することができる。更に、マルチバンドカメラで撮影した画像に対し、特定のエリアを指定し、その場所の色度情報を表示することができる。そして、エリアを複数設定した場合はその違いについて、色度図上でプロット位置を同時に表示して比較することができる効果がある。更に、マルチバンドカメラで撮影した画像に対し、特定のエリアを指定し、その場所の色度値、色度図を表示することができる。そして、色度図において、任意の部分の拡大表示が行えるので、色度図の部分的な情報を細かく見ることができる。
【０１５３】
（４４） ２つ以上のバンドを撮影できるマルチバンドカメラで取得したマルチバンド画像を表示する画像表示方法において、
上記マルチバンド画像の内のあるバンドの画像を表示する工程と、
その表示された上記あるバンドの画像上で、少なくとも１画素以上を指定する工程と、
上記マルチバンド画像における上記指定された上記少なくとも１画素以上の信号値から、分光スペクトルを推定する工程と、
その推定した分光スペクトルから、色度情報を求める工程と、
上記推定した分光スペクトルを表示すると共に、上記求めた色度情報の色度値を色度図として表示する工程と、
を有し、
上記表示する工程は、上記表示した色度図において、任意の部分の拡大表示を行う工程を含むことを特徴とする画像表示方法。
【０１５４】
（対応する実施形態）
この（４４）に記載の画像表示方法に関する実施形態は、第２実施形態が対応する。
（作用効果）
この（４４）に記載の画像表示方法によれば、画像の中でスペクトルを求めたい場所を特定し、グラフを表示することができる。また、マルチバンド画像の中でスペクトルを求めたい場所の信号値の分光反射／透過率スペクトルのグラフを表示することができる。更に、マルチバンドカメラで撮影した画像に対し、特定のエリアを指定し、その場所の色度情報を表示することができる。そして、エリアを複数設定した場合はその違いについて、色度図上でプロット位置を同時に表示して比較することができる効果がある。更に、マルチバンドカメラで撮影した画像に対し、特定のエリアを指定し、その場所の色度値、色度図を表示することができる。そして、色度図において、任意の部分の拡大表示が行えるので、色度図の部分的な情報を細かく見ることができる。
【図面の簡単な説明】
【０１５５】
【図１】（Ａ）は本発明の第１実施形態に係る画像表示装置の構成を示す図であり、（Ｂ）はスペクトル出力用画面を示す図である。
【図２】バンド画像画面を示す図である。
【図３】補間方法選択画面を示す図である。
【図４】第１実施形態に係る画像表示装置の変形例の構成を示す図である。
【図５】（Ａ）は色票を撮影した場合のバンド画像を表示したバンド画像画面を示す図であり、（Ｂ）は作成した参照画像を表示した参照画像画面を示す図である。
【図６】バンド選択画面を示す図である。
【図７】変形例における補間方法選択画面を示す図である。
【図８】第１実施形態に係る画像表示装置の別の変形例の構成を示す図である。
【図９】エリア１を基準とした信号値（輝度）の相対値を表示するグラフを示す図である。
【図１０】第１実施形態に係る画像表示装置の更に別の変形例の構成を示す図である。
【図１１】本発明の第２実施形態に係る画像表示装置における表示例を示す図である。
【図１２】第２実施形態に係る画像表示装置における別の表示例を示す図である。
【図１３】第２実施形態に係る画像表示装置における更に別の表示例を示す図である。
【図１４】（Ａ）は本発明の第３実施形態に係る画像表示装置において基準エリアをマウスで指定している様子を示す図であり、（Ｂ）はその結果として表示される色差画像を示す図である。
【図１５】第３実施形態に係る画像表示装置における別の表示例を示す図である。
【符号の説明】
【０１５６】
１…マルチバンドカメラ、 １１…コンピュータ、 １２…マルチバンド画像、 １３…スペクトル計算部、 １４…スペクトル出力用画面、 １５…バンド画像画面、 １６…補間方法選択画面、 １７…参照画像、 １８…参照画像画面、 １９…バンド選択画面、 ２０…補間方法選択画面、 ２１…特性ファイル、 ２２，２３…ウィンドウ、 ２４…マウスポインタ、 １４１…エリア選択リストボックス、 １５１…矩形、 １９１…バンド番号、 ２２１…色度図、 ２２２−１，２２２−２…点、 ２２３，２３２…釣鐘、 ２３１…色度図。
【特許請求の範囲】
【請求項１】
２つ以上のバンドを撮影できるマルチバンドカメラで取得したマルチバンド画像を表示する画像表示装置において、
上記マルチバンド画像の内のあるバンドの画像を表示する画像表示手段と、
上記画像表示手段に表示された上記あるバンドの画像上で、少なくとも１画素以上を指定する画素指定手段と、
上記マルチバンド画像における上記画素指定手段によって指定された上記少なくとも１画素以上の信号値から、分光スペクトルを推定するスペクトル推定手段と、
上記スペクトル推定手段で推定した分光スペクトルから、色度情報を求める色度情報算出手段と、
上記スペクトル推定手段で推定した分光スペクトルを表示すると共に、上記色度情報算出手段で求めた上記色度情報の色度値を色度図として表示するスペクトル表示手段と、
を具備し、
上記スペクトル表示手段は、上記表示した色度図において、任意の部分の拡大表示機能を備えることを特徴とする画像表示装置。
【請求項２】
２つ以上のバンドを撮影できるマルチバンドカメラで取得したマルチバンド画像を表示する画像表示方法において、
上記マルチバンド画像の内のあるバンドの画像を表示する工程と、
その表示された上記あるバンドの画像上で、少なくとも１画素以上を指定する工程と、
上記マルチバンド画像における上記指定された上記少なくとも１画素以上の信号値から、分光スペクトルを推定する工程と、
その推定した分光スペクトルから、色度情報を求める工程と、
上記推定した分光スペクトルを表示すると共に、上記求めた色度情報の色度値を色度図として表示する工程と、
を有し、
上記表示する工程は、上記表示した色度図において、任意の部分の拡大表示を行う工程を含むことを特徴とする画像表示方法。
【請求項１】
２つ以上のバンドを撮影できるマルチバンドカメラで取得したマルチバンド画像を表示する画像表示装置において、
上記マルチバンド画像の内のあるバンドの画像を表示する画像表示手段と、
上記画像表示手段に表示された上記あるバンドの画像上で、少なくとも１画素以上を指定する画素指定手段と、
上記マルチバンド画像における上記画素指定手段によって指定された上記少なくとも１画素以上の信号値から、分光スペクトルを推定するスペクトル推定手段と、
上記スペクトル推定手段で推定した分光スペクトルから、色度情報を求める色度情報算出手段と、
上記スペクトル推定手段で推定した分光スペクトルを表示すると共に、上記色度情報算出手段で求めた上記色度情報の色度値を色度図として表示するスペクトル表示手段と、
を具備し、
上記スペクトル表示手段は、上記表示した色度図において、任意の部分の拡大表示機能を備えることを特徴とする画像表示装置。
【請求項２】
２つ以上のバンドを撮影できるマルチバンドカメラで取得したマルチバンド画像を表示する画像表示方法において、
上記マルチバンド画像の内のあるバンドの画像を表示する工程と、
その表示された上記あるバンドの画像上で、少なくとも１画素以上を指定する工程と、
上記マルチバンド画像における上記指定された上記少なくとも１画素以上の信号値から、分光スペクトルを推定する工程と、
その推定した分光スペクトルから、色度情報を求める工程と、
上記推定した分光スペクトルを表示すると共に、上記求めた色度情報の色度値を色度図として表示する工程と、
を有し、
上記表示する工程は、上記表示した色度図において、任意の部分の拡大表示を行う工程を含むことを特徴とする画像表示方法。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図１２】
【図１３】
【図１４】
【図１５】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図１２】
【図１３】
【図１４】
【図１５】
【公開番号】特開２００８−１５１７８１（Ｐ２００８−１５１７８１Ａ）
【公開日】平成２０年７月３日（２００８．７．３）
【国際特許分類】
【出願番号】特願２００７−３１２８１１（Ｐ２００７−３１２８１１）
【出願日】平成１９年１２月３日（２００７．１２．３）
【分割の表示】特願２００３−３４８６２５（Ｐ２００３−３４８６２５）の分割
【原出願日】平成１５年１０月７日（２００３．１０．７）
【出願人】（００００００３７６）オリンパス株式会社 (11,466)
【Ｆターム（参考）】
【公開日】平成２０年７月３日（２００８．７．３）
【国際特許分類】
【出願日】平成１９年１２月３日（２００７．１２．３）
【分割の表示】特願２００３−３４８６２５（Ｐ２００３−３４８６２５）の分割
【原出願日】平成１５年１０月７日（２００３．１０．７）
【出願人】（００００００３７６）オリンパス株式会社 (11,466)
【Ｆターム（参考）】
[ Back to top ]