(57)【要約】【課題】本発明は、分光反射率、分光透過率、分光放射輝度及び分光光度等の分光色を簡易に補正する色補正装置を提供する。【解決手段】分光色の画像データを色補正する際に、分光色を元の次元よりも低次元の色空間に変換し、低次元色空間内で色補正を行い、変換後の低次元色から適切な次元の分光色を生成するようにしたものである。また、色補正が鮮やかさを向上する場合には、分光色の波長形状を波長幅の狭い形状に変換するようにしたものであり、分光色の画像データの色補正を簡易に行うことができる。【選択図】　図１
2004045189.tif 000002
【特許請求の範囲】
【請求項１】
画像信号が入力され、前記画像信号を分光的に表現した第１の分光色に変換する分光画像推定部と、前記第１の分光色を前記第１の分光色の次元よりも低次元の色空間である低次元色空間に変換する空間変換部と、前記低次元色空間から色補正を行う領域である抽出領域を抽出する変更色領域設定部と、前記抽出領域の前記低次元色空間を補正する色補正部と、前記補正された低次元色空間を第２の分光色に変換し、出力する色空間逆変換部とを含む色補正装置。
【請求項２】
抽出領域が、予め決められた値以上の色空間の画素である請求項１記載の色補正装置。
【請求項３】
色空間逆変換部が、補正された低次元色空間に、第２の分光色を算出するための変化係数を追加し、前記変化係数を変動させ、第１の分光色に近似した分光色に変換し、出力する請求項１又は２記載の色補正装置。
【請求項４】
変化係数が、分割された波長領域毎に設定され、前記分割波長領域の境界では、線形補間される係数である請求項３記載の色補正装置。
【請求項５】
色空間逆変換部が、補正された低次元色空間から三刺激値を算出し、前記三刺激値に変化係数を乗算し、前記変化係数を変動させ、第２の分光色を算出し、出力する請求項１又は２記載の色補正装置。
【請求項６】
低次元色空間が、色相、明度及び彩度から構成されるＨＬＳ空間である請求項１ないし５のいずれか記載の色補正装置。
【請求項７】
低次元色空間が、明度及び知覚色度から構成されるＬ＊ａ＊ｂ＊空間である請求項１ないし５のいずれか記載の色補正装置。
【請求項８】
色補正が、抽出領域の彩度の値と、予め決められた閾値との差分を求め、前記差分に応じて、前記彩度の値に補正値を加える補正である請求項１ないし７のいずれか記載の色補正装置。
【請求項９】
画像信号が入力され、前記画像信号が分光的に表現した分光色に変換される分光画像推定部と、前記分光色から色補正を行う領域である抽出領域を抽出する変更色領域設定部と、前記抽出領域の前記分光色を補正する色補正部とを含む色補正装置。
【請求項１０】
分光色の補正が、抽出領域の波長幅の縮小である請求項９記載の色補正装置。
【請求項１１】
抽出領域が、分光色の２次差分値から極大値及び極小値を算出し、最も大きい前記極大値の両側に存在する２つの前記極小値に囲まれた波長である請求項９又は１０記載の色補正装置。
【請求項１２】
画像信号が入力され、前記画像信号が分光的に表現した第１の分光色に変換される第１のステップと、前記第１の分光色を前記第１の分光色の次元よりも低次元の色空間に変換する第２のステップと、予め決められた値以上の色空間の画素が、前記低次元色空間から色補正を行う領域として抽出される第３のステップと、前記抽出領域の前記低次元色空間を補正する第４のステップと、前記補正された低次元色空間を第２の分光色に変換し、出力する第５のステップとを含む色補正方法。
【請求項１３】
第５のステップが、分割された波長領域毎に設定され、前記分割波長領域の境界では、線形補間される変化係数を、補正された低次元色空間に加えた第３の分光色に変換し、前記変化係数を変動させ、第１の分光色に近似した第２の分光色に変換し、出力する請求項１２記載の色補正方法。
【請求項１４】
第５のステップが、補正された低次元色空間を、三刺激値に変換し、前記三刺激値から算出される第２の分光色を出力する請求項１２記載の色補正方法。
【請求項１５】
低次元色空間が、色相、明度及び彩度から構成されるＨＬＳ空間である請求項１２ないし１４のいずれか記載の色補正方法。
【請求項１６】
低次元色空間が、明度及び知覚色度から構成されるＬ＊ａ＊ｂ＊空間である請求項１２ないし１４のいずれか記載の色補正方法。
【請求項１７】
画像信号が入力され、前記画像信号が分光的に表現した分光色に変換される第１のステップと、前記分光色の２次差分値から算出される極大値の両側に存在する２つの極小値に囲まれた波長が、前記分光色から色補正を行う領域として抽出される第２のステップと、前記抽出領域の波長幅の短縮により、前記抽出領域の前記分光色を補正する第３のステップとを含む色補正方法。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、分光反射率、分光透過率、分光放射輝度及び分光光度から成る画像などの情報を補正する色補正装置及び色補正方法に関するものである。
【０００２】
【従来の技術】
カラー画像を処理する手段として、スキャナ、デジタルカメラ、プリンタ及びディスプレイ等が利用されている。これらの手段間でカラー画像データを交換する手法に、スキャナ、ディジタルスチルカメラ、ディジタルムービー、マルチスペクトルカメラ又はハイパースペクトルカメラなどの画像入力手段から入力されたカラー画像データを、画像入力手段及びディスプレイ又はプリンタなどの画像出力手段に依存しない独立な色空間に変換してから、画像出力手段に出力するカラー画像データに変換するカラーマネジメント装置がある。
【０００３】
このため、画像入力手段及び画像出力手段の組み合わせ毎に、色変換処理方法を決める必要がない。
【０００４】
また、画像入力手段から入力されたカラー画像データを、画像入力手段及び画像出力手段だけでなく、カラー画像データに含まれる撮影時に撮影対象を照らした照明にも依存しない、独立な色空間に変換しておけば、撮影対象の撮影時の照明とは異なる照明下の画像を、画像出力手段に出力するができる。
【０００５】
ここで、画像入力手段及び画像出力手段に依存しない独立な色空間としては、国際標準機関である国際照明委員会（Ｉｎｔｅｒｎａｔｉｏｎａｌ　Ｃｏｍｍｉｓｓｉｏｎ　ｏｎ　Ｉｌｌｕｍｉｎａｔｉｏｎ：ＣＩＥ）が規定するＸＹＺ三刺激値及びＬ＊ａ＊ｂ＊均等色空間が一般的である。
【０００６】
本来、色は、各波長の色の含有割合である分光放射輝度で表現される。ここで、ＸＹＺ三刺激値は、その分光放射輝度を、ＣＩＥが定めた等色関数で重み付けし、その後可視域で積分して得られる３つのスカラー量である。
【０００７】
また、等色関数とは、可視波長の全域にわたり、それぞれ等しい放射パワーを持つ単色光刺激の三刺激値である。
【０００８】
更に、Ｌ＊ａ＊ｂ＊均等色空間は、ＸＹＺ三刺激値をＣＩＥが定めた標準光源で正規化した３次元の色空間である。
【０００９】
色は、本来４００ｎｍから７００ｎｍの可視域を、１０ｎｍ毎にサンプリングしても、３１次元という膨大なデータであるのに対し、これらＸＹＺ三刺激値及びＬ＊ａ＊ｂ＊均等色空間のように、色を３次元で表現する形式は、人間の目が３つの錐体により知覚される理由から、分光色に代わり人間が知覚する色を表現する量としてこれまで広く使われできた。
【００１０】
しかし、近年の情報処理能力の向上により、ＸＹＺ三刺激値及びＬ＊ａ＊ｂ＊均等色空間に代わり、分光的に色を扱う色補正装置が開発されるようになった。
【００１１】
ここで、分光表現には、分光反射率、分光透過率、分光放射輝度及び分光光度がある。分光放射輝度は、被写体から反射又は透過される光の量を表すため、照明の情報を含んでおり、分光放射輝度から照明の分光分布を除いたものを、それぞれ分光反射率又は分光透過率という。
【００１２】
更に、分光表現の中でも、分光反射率又は分光透過率を独立色空間として扱う色補正装置が広まりつつある。すなわち、画像入力手段から被写体の分光反射率又は分光透過率を推定し、これを画像とした分光反射率画像又は分光透過率画像を扱う色補正装置である。以下、分光反射率及び分光透過率に代表される分光表現の色を総称して、分光色と呼ぶ。
【００１３】
分光色を扱った装置には、例えば特開平７−１７４６３１号公報記載の皮膚の測色方法及び分光反射スペクトルの推定方法及び特開平９−２３３４９０号公報記載の色シミュレーション装置がある。これらの装置では、被写体の物理的特性の解析又は様々な照明下で撮影した被写体のシミュレーションを行うために、画像入力手段より入力された画像データから、被写体の分光反射率画像を推定している。
【００１４】
以下に、推定方法を簡単に述べる。まず、分光反射率を主成分分析し、３１次よりも低次元なｍ次の基底で表現する。そして、入力画像データから基底係数ｍ次ベクトルをニューラルネットワーク（神経回路網：Ｎｅｕｒａｌ　Ｎｅｔｗｏｒｋ）で推定する。
【００１５】
ここで、ニューラルネットワークとは、多数の入出力関係の学習からどのような入力があっても、学習が反映された出力がなされるもの、つまり、人間の脳の構造を参考にして構成された情報処理機構である。
【００１６】
更に、推定されたｍ次ベクトルから分光反射率を算出することで、分光反射率の推定を行う。
【００１７】
また、分光色の画像データを色補正する装置には、特開平７−１６２６９４号公報記載の色処理装置がある。図１１は、公報記載の色補正用ディスプレイである。
【００１８】
ここで、図１１（ａ）は、分光色を表示する分光色表示枠を示し、図１１（ｂ）は、ｘｙ色度図表示枠を示し、図１１（ｃ）は、ＨＣＶ表示枠を示す。また、ＨＣＶとは、Ｍｕｎｓｅｌｌ表色系による色相（ｈｕｅ：Ｈ）、彩度（ｃｈｒｏｍａ：Ｃ）及び明度（ｖａｌｕｅ：Ｖ）である。
【００１９】
色補正の手順を説明する。第１に、補正する色見本の分光色が、グラフとして図１１（ａ）の分光色表示枠に表示される。第２に、グラフ上の曲線を、アドビ社製のイラストレーターなどの補正手段により、任意の形状に補正を行う。第３に、補正された分光色は、ｘｙ色度及びＨＣＶ信号に変換され、ｘｙ色度は、図１１（ｂ）に示したｘｙ色度図表示枠に表示され、ＨＣＶ信号は、図１１（ｃ）に示したＨＣＶ表示枠に表示される。第４に、ｘｙ色度図表示枠及びＨＣＶ表示枠の値を参照して、再度、第２の手段から第４の手段を繰り返すことで分光色の補正を行う。
【００２０】
【発明が解決しようとする課題】
しかし、従来は、色編集及び色補正は、画像入力手段及び画像出力手段固有の色空間、例えば、カメラのＲＧＢ色空間、プリンタのＣＭＹＫ色空間及びディスプレイのＲＧＢ色空間上で行われてきた。しかし、分光色を用いたカラーマネジメント装置の普及により、独立色空間による分光色の画像データの色編集及び色補正が要望されている。
【００２１】
更に、分光反射率をある基底で表現すると、基底に波長幅の狭いデルタ関数的なものが含まれていなければ、そのような色は推定することができない。即ち、波長幅の狭い、極めて彩やかな色は推定値として得られないことになる。従って、推定された画像は、比較的彩度が低い印象の画像となってしまうという問題がある。
【００２２】
本発明は、かかる点に鑑みてなされたもので、分光反射率、分光透過率、分光放射輝度及び分光光度を含む分光色の画像データを補正する色補正装置及び色補正方法を、提供することを目的とする。
【００２３】
【課題を解決するための手段】
本発明は、分光色の画像データを色補正する際に、分光色を元の次元よりも低次元の色空間に変換し、低次元色空間内で色補正を行い、変換後の低次元色から適切な次元の分光色を生成するようにしたものである。
【００２４】
また、色補正が鮮やかさを向上する場合には、分光色の波長形状を波長幅の狭い形状に変換するようにしたものであり、分光色の画像データの色補正を簡易に行うことができる。
【００２５】
【発明の実施の形態】
本発明の請求項１に記載の発明は、画像信号が入力され、前記画像信号が分光的に表現した第１の分光色に変換される分光画像推定部と、前記第１の分光色を前記第１の分光色の次元よりも低次元の色空間に変換する空間変換部と、前記低次元色空間から色補正を行う領域が抽出される変更色領域設定部と、前記抽出領域の前記低次元色空間を補正する色補正部と、前記補正された低次元色空間を第２の分光色に変換し、出力する色空間逆変換部とを含む色補正装置としたものであり、分光色を低次元の色空間に変換するため、色補正が容易であり、更に、分光色に逆変換するため、画像出力手段に容易に出力できるという作用を有する。
【００２６】
請求項２に記載の発明は、抽出領域が、予め決められた値以上の色空間の画素である請求項１記載の色補正装置としたものであり、抽出領域を容易に決定し、色補正装置の処理負荷を軽減することができるという作用を有する。
【００２７】
請求項３に記載の発明は、色空間逆変換部が、補正された低次元色空間に、第２の分光色を算出するための変化係数を追加し、前記変化係数を変動させ、第１の分光色に近似した分光色に変換し、出力する請求項１又は２記載の色補正装置としたものであり、低次元の色空間を高次元の分光色に変換することができるという作用を有する。
【００２８】
請求項４に記載の発明は、変化係数が、分割された波長領域毎に設定され、前記分割波長領域の境界では、線形補間される係数である請求項３記載の色補正装置としたものであり、急激な変化係数の変動を抑制し、自然な色補正ができるという作用を有する。
【００２９】
請求項５に記載の発明は、色空間逆変換部が、補正された低次元色空間から三刺激値を算出し、前記三刺激値に変化係数を乗算し、前記変化係数を変動させ、第２の分光色を算出し、出力する請求項１又は２記載の色補正装置としたものであり、低次元の色空間を高次元の分光色に変換することができるという作用を有する。
【００３０】
請求項６に記載の発明は、低次元色空間が、色相、明度及び彩度から構成されるＨＬＳ空間である請求項１ないし５のいずれか記載の色補正装置としたものであり、色補正を容易に行うことができるという作用を有する。
【００３１】
請求項７に記載の発明は、低次元色空間が、明度及び知覚色度から構成されるＬ＊ａ＊ｂ＊空間である請求項１ないし５のいずれか記載の色補正装置としたものであり、色補正を容易に行うことができるという作用を有する。
【００３２】
請求項８に記載の発明は、色補正が、抽出領域の彩度の値と、予め決められた閾値との差分を求め、前記差分に応じて、前記彩度の値に補正値を加える補正である請求項１ないし７のいずれか記載の色補正装置としてものであり、補正値を差分により変化させることで、より自然な彩度変換を行うことができるという作用を有する。
【００３３】
請求項９に記載の発明は、画像信号が入力され、前記画像信号が分光的に表現した分光色に変換される分光画像推定部と、前記分光色から色補正を行う領域が抽出される変更色領域設定部と、前記抽出領域の前記分光色を補正する色補正部とを含む色補正装置としたものであり、分光色の色補正ができるという作用を有する。
【００３４】
請求項１０に記載の発明は、分光色の補正が、抽出領域の波長幅の縮小である請求項９記載の色補正装置としたものであり、分光色の色補正を容易に行うことができるという作用を有する。
【００３５】
請求項１１に記載の発明は、抽出領域が、分光色の２次差分値から極大値及び極小値を算出し、最も大きい前記極大値の両側に存在する２つの前記極小値に囲まれた波長である請求項８又は１０記載の色補正装置としたものであり、色補正する領域が容易に抽出できるという作用を有する。
【００３６】
請求項１２に記載の発明は、画像信号が入力され、前記画像信号が分光的に表現した第１の分光色に変換される第１のステップと、前記第１の分光色を前記第１の分光色の次元よりも低次元の色空間に変換する第２のステップと、予め決められた値以上の色空間の画素が、前記低次元色空間から色補正を行う領域として抽出される第３のステップと、前記抽出領域の前記低次元色空間を補正する第４のステップと前記補正された低次元色空間を第２の分光色に変換し、出力する第５のステップとを含む色補正方法としたものであり、分光色を低次元の色空間に変換するため、色補正が容易であり、更に、分光色に逆変換するため、画像出力手段に容易に出力できるという作用を有する。
【００３７】
請求項１３に記載の発明は、第５のステップが、分割された波長領域毎に設定され、前記分割波長領域の境界では、前記変化係数を線形補間される変化係数を、補正された低次元色空間に加えた第３の分光色に変換し、前記変化係数を変動させ、第１の分光色に近似した第２の分光色に変換し、出力する請求項１２記載の色補正方法としたものであり、急激な変化係数の変動を抑制し、自然な色補正ができ、更に、低次元の色空間を高次元の分光色に変換することができるという作用を有する。
【００３８】
請求項１４に記載の発明は、第５のステップが、補正された低次元色空間を、三刺激値に変換し、前記三刺激値から算出される第２の分光色を出力する請求項１２記載の色補正方法としたものであり、低次元の色空間を高次元の分光色に変換することができるという作用を有する。
【００３９】
請求項１５に記載の発明は、低次元色空間が、色相、明度及び彩度から構成されるＨＬＳ空間である請求項１２ないし１４のいずれか記載の色補正方法としたものであり、色補正を容易に行うことができるという作用を有する。
【００４０】
請求項１６に記載の発明は、低次元色空間が、明度及び知覚色度から構成されるＬ＊ａ＊ｂ＊空間である請求項１２ないし１４のいずれか記載の色補正方法としたものであり、色補正を容易に行うことができるという作用を有する。
【００４１】
請求項１７に記載の発明は、画像信号が入力され、前記画像信号が分光的に表現した分光色に変換される第１のステップと、前記分光色の２次差分値から算出される極大値の両側に存在する２つの極小値に囲まれた波長が、前記分光色から色補正を行う領域として抽出される第２のステップと、前記抽出領域の波長幅の短縮により、前記抽出領域の前記分光色を補正する第３のステップとを含む色補正方法としたものであり、分光色の色補正ができるという作用を有する。
【００４２】
以下、本発明の実施の形態について図面を参照して詳細に説明する。
【００４３】
（実施の形態１）
第１の実施の形態では、入力画像データから推定して得られた分光色に対し、好ましい色補正を施してから出力する色補正装置を開示する。
【００４４】
第１の実施の形態は、入力された画像から分光色の画像データを推定し、得られた分光色の画像データに対し、自動的に彩度を向上させる色補正を施すものである。
【００４５】
また、入力画像データから分光色を完全に推定できた場合においても、人間が好ましいと感じる画像は、彩度が強調されたものである場合が多いため、本実施の形態で開示する色補正を適用することができる。
【００４６】
図１は、本発明第１の実施の形態による分光色の色補正装置の構成を示す図である。
【００４７】
ここで、１０１は、入力画像データから分光色の画像データを推定する分光画像推定部を示し、１０２は、分光色の画像データの各画素値をＨＬＳ空間に変換するＨＬＳ空間変換部を示し、１０３は、ＨＬＳ空間に変換されたＨＬＳ画像において彩度が、予め決められた閾値よりも高い領域を変更すべき色として設定するＨＬＳ空間変更色領域設定部を示し、１０４は、ＨＬＳ画像の変更色領域における彩度の値を高くする色補正を行うＨＬＳ色補正部を示し、１０５は、色補正後のＨＬＳ画像の各画素を元の分光空間に逆変換するＨＬＳ分光空間逆変換部を示す。
【００４８】
本実施の形態では、分光色の画像データとして、各画素値が分光反射率である場合を取り扱う。その他の分光色の画像データとして、各画素値が分光放射輝度、分光光度又は分光透過率であった場合にも、本発明の実施の形態を適用することができる。
【００４９】
以下に、本発明第１の実施の形態による色補正装置について、その動作を説明する。
【００５０】
分光画像推定部１０１では、入力された画像データから分光色の画像データを推定する。ＨＬＳ空間変換部１０２は、分光色の画像データにおける各画素の可視領域のサンプリング次元数を（数１）、（数２）及び（数３）を用いて、より低いサンプリング次元数のＨＬＳ空間に変換する。
【００５１】
【数１】
2004045189.tif 000003
【００５２】
【数２】
2004045189.tif 000004
【００５３】
【数３】
2004045189.tif 000005
【００５４】
ここで、（数１）において、ｒ（λ）は、分光色の画像データによる各画素の値である分光反射率に、指定照明を乗算して得られる分光放射輝度である。また、指定照明とは、撮影対象を照らしている照明の分光放射輝度であり、既知、センサにより計測又は推定方法により推定されているものとする。
【００５５】
更に、Ｓ（λ）は、ＣＩＥにより定められた色温度６５００Ｋの標準光源であるＤ６５照明の分光分布を示し、ｘ（λ）、ｙ（λ）及びｚ（λ）は、ＣＩＥで規定されている視角が２度である２度視野等色関数を示し、Ｘ、Ｙ及びＺは、分光放射輝度から計算される三刺激値である。
【００５６】
次に、（数２）において、Ｌ＊、ａ＊及びｂ＊は、ＣＩＥにより定められた（Ｌ＊ａ＊ｂ＊）表示系であり、Ｌ＊は、明度を示し、ａ＊及びｂ＊は、知覚色度を示し、Ｘｎ、Ｙｎ及びＺｎは、照明に使用する標準光源の三刺激値である。
【００５７】
更に、（数３）において、Ｈ、Ｌ及びＳは、それぞれ求めるべき色相（ｈｕｅ：Ｈ）、明度（ｌｉｇｈｔｎｅｓｓ：Ｌ）及び彩度（ｓａｔｕｒａｔｉｏｎ：Ｓ）を示す。
【００５８】
ＨＬＳ空間変更色領域設定部１０３は、ＨＬＳ画像において、予め決められた閾値Ｃｔｈｒｅ以上の彩度を持つ画素を変更色領域として設定し出力する。
【００５９】
ＨＬＳ色補正部１０４は、変更色領域の画素の彩度に対し、彩度を補正する補正値であるΔＣを加算して、彩度変更を行う。このとき、画素の彩度が、Ｃｔｈｒｅに近いときは、ΔＣの値を０にし、Ｃｔｈｒｅより大きい場合は、ΔＣの値を大きくすることで、より自然な彩度変換を行うことができる。
【００６０】
ＨＬＳ分光空間逆変換部１０５は、ＨＬＳ色補正部１０４から出力された彩度変更が行われたＨＬＳ画像を、分光空間に逆変換し、出力画像データとして出力する。
【００６１】
ＨＬＳ空間から分光空間への逆変換のアルゴリズムは複数存在する。これらすべてをシステムに組み込んでおいて、選択してもらってもよいし、いずれか１種類のアルゴリズムが動作するようにしておいてもよい。
【００６２】
まず、第１の逆変換方法は、色補正後のＨＬＳ画像から得られる分光色が、色補正前の分光色の分布に、最も近いものを選択する方法である。色補正前の分光色と、色補正後のＨＬＳ画像から得られる分光色との距離の近さを拘束条件として、（数４）を最小化する分光色を選択すればよい。
【００６３】
【数４】
2004045189.tif 000006
【００６４】
ここで、ｓｂ（λ）は、色補正前の分光色であり、ｓａ（λ）は、出力する色補正後分光色であり、（Ｈｔａｒｇｅｔ、Ｌｔａｒｇｅｔ、Ｓｔａｒｇｅｔ）は、色補正後のＨＬＳ値であり、（Ｈａ、Ｌａ、Ｓａ）は、初期値設定したｓａ（λ）から算出されるＨＬＳ値であり、αはスカラー係数であり、予め決めておく。
【００６５】
更に、（数４）に求めるべき分光色の形状は滑らかであるという拘束条件を加えてもよい。即ち、ｓａ（λ）の２次差分が０に近いという条件を加味する。スペクトルの２次差分は（数５）となる。
【００６６】
【数５】
2004045189.tif 000007
【００６７】
（数５）のΔｓａ（λｉ）は、波長Ｉによるスペクトルｓａ（λｉ）の２次差分値であり、（数５）の二乗和は（数６）となり、これはスペクトルの滑らかさを表す。
【００６８】
【数６】
2004045189.tif 000008
【００６９】
（数６）を（数４）に加味すると、（数７）が得られ、最小となるｓａ（λ）を求めればよい。
【００７０】
【数７】
2004045189.tif 000009
【００７１】
これらの式を最小化するｓａ（λ）を計算するには、例えば、シンプレックス法など既存の最適化手法で求める。この場合、係数α及びβは経験値から予め決定しておく。又は、係数α及びβを未定係数として、ラグランジェの未定係数法として解くことも可能である。
【００７２】
ここで、計算手順を図２に示す。図２は、（数４）を最適化手法にシンプレックス法を用いて、最小となるｓａ（λ）を求める場合を示す。
【００７３】
まず、色補正前分光色ｓｂ（λ）及び色補正後ＨＬＳ色である（Ｈｔａｒｇｅｔ、Ｌｔａｒｇｅｔ、Ｓｔａｒｇｅｔ）が入力され（Ｓ２０１）、出力すべき色補正後分光色ｓａ（λ）の変化係数の初期値を設定する（Ｓ２０２）。色補正後分光色ｓａ（λ）から（数１）により（Ｈａ、Ｌａ、ｓａ）を求める（Ｓ２０３）。
【００７４】
次に、（Ｈａ、Ｌａ、ｓａ）を（数４）に代入し、分光値が最小となる距離を求める（Ｓ２０４）。最適化手法であるシンプレックス法に基づく収束条件により、ｓａ（λ）の収束状態を判断し（Ｓ２０６）、収束していなければ、ｓａ（λ）の値を、シンプレックス法に基づいて、変化係数を変動させ（Ｓ２０５）、Ｓ２０３に入力する。また、収束していればｓａ（λ）を解として出力する（Ｓ２０７）。
【００７５】
次に、第２の逆変換方法について説明する。第２の逆変換方法は、第１の逆変換方法をより高速化したものである。第１の逆変換方法では、３１次元の分光ベクトルの各値を変数として、数式を最小化する値を求めていた。しかし、最適化手法で求めるには、変数の数が多いほど収束までに時間がかかる。また、各変数を独立に変化させるため、波長形状が不自然なものとなる可能性もある。
【００７６】
ここで、図３に分光反射率と波長の関係を示す。縦軸は分光反射率を示し、横軸は波長を示し、３０１は、色補正前の分光色ｓｂ（λ）を示し、３０２は、可視波長領域を分割した分割波長領域を示す。
【００７７】
また、３０２ａは、３８０ｎｍ以上４８０ｎｍ未満を示し、３０２ｂは、４８０ｎｍ以上５８０ｎｍ未満を示し、３０２ｃは、５８０ｎｍ以上６８０ｎｍ未満を示し、３０２ｄは、６８０ｎｍ以上７８０ｎｍ未満を示す。
【００７８】
また、図４に変化係数と波長の関係を示す。縦軸は変化係数Ｗｉを示し、横軸は波長を示し、４０１は、変化係数Ｗｉを線形補間する線形補間領域を示す。
【００７９】
第２の逆変換方法では、図３に示すように、波長を４つの分割波長領域３０２に分割し、領域毎に変化係数を設定する。図４に分割波長領域３０２毎の変化係数Ｗｉを示す。３０２ａは、変化係数Ｗｉを０．７とし、３０２ｂは、変化係数Ｗｉを１．２とし、３０２ｃは、変化係数Ｗｉを０．６とし、３０２ｄは、変化係数Ｗｉを１．３とする。
【００８０】
図４に示す変化係数Ｗｉと波長の関係のように、変化係数Ｗｉが滑らかにつながるように、線形補間領域４０１で変化係数Ｗｉを線形補間し、変化関数を求める。
【００８１】
次に、図４の得られた変化関数を、図３の色補正前分光色ｓｂ（λ）３０１の各波長値に掛け合わせて、求めるべき色補正後分光色ｓａ（λ）の候補とする。
【００８２】
この第２の逆変換方法によれば、最適化する変数は、変化係数の数、即ち、図４の場合は４つの分割波長領域３０２に分割し、更に、変化係数は連続的になるように変化関数を求めて分光色に乗算するため、色補正前の分光色の形状を可能な限り保ち、不自然さを極力排除した分光色を、（数８）を用いて算出することができる。
【００８３】
【数８】
2004045189.tif 000010
【００８４】
なお、（数８）には、（数６）及び（数７）に示したような拘束条件を加えてもよい。
【００８５】
第２の逆変換手法による計算手順を図５に示す。色補正前分光色ｓｂ（λ）３０１及び色補正後ＨＬＳ色をＨＬＳ分光逆変換部に入力する（Ｓ５０１）。変化係数Ｗｉの初期値を設定する（Ｓ５０２）。変化係数Ｗｉから変化関数Ｗｅｉｇｈｔ（λ）を算出する（Ｓ５０３）。変化関数Ｗｅｉｇｈｔ（λ）と色補正前分光色ｓｂ（λ）３０１から（数９）により色補正後分光色の候補となるＣａｎｄｉｄａｔｅ（λ）を求める（Ｓ５０４）。
【００８６】
【数９】
2004045189.tif 000011
【００８７】
Ｃａｎｄｉｄａｔｅ（λ）をＨＬＳ色空間に変換し、（Ｈａ、Ｌａ、Ｓａ）とする（Ｓ５０５）。（Ｈａ、Ｌａ、Ｓａ）を（数８）に代入し、関数値を求める（Ｓ５０６）。最適化法の収束条件から変化係数Ｗｉの値の状態を判断し（Ｓ５０８）、最適化法の収束条件を満たさない場合は、変化係数Ｗｉの値を、最適化手法に基づいて変動させ（Ｓ５０７）、Ｓ５０３ないしＳ５０６の手順を繰り返す。最適化法の収束条件を満たす場合は、Ｃａｎｄｉｄａｔｅ（λ）を解ｓａ（λ）として出力する（Ｓ５０８）。
【００８８】
更に、第３の逆変換方法について説明する。本方法では、色補正後ＨＬＳ色（Ｈｔａｒｇｅｔ、Ｖｔａｒｇｅｔ、Ｃｔａｒｇｅｔ）からＸＹＺ三刺激値（Ｘｔａｒｇｅｔ、Ｙｔａｒｇｅｔ、Ｚｔａｒｇｅｔ）を求め、これを満足する分光色を決定する。ＸＹＺ三刺激値と分光色は（数１）に示す関係がある。
【００８９】
ここで、（数１）を離散的に表現すれば（数１０）となり、（Ｘｔａｒｇｅｔ、Ｙｔａｒｇｅｔ、Ｚｔａｒｇｅｔ）を満たすｒ（λ）を求めればよい。
【００９０】
【数１０】
2004045189.tif 000012
【００９１】
（数１０）を解く場合、既知の値が３次元（Ｘｔａｒｇｅｔ、Ｙｔａｒｇｅｔ、Ｚｔａｒｇｅｔ）であるのに対し、未知の値が３１次元の分光色であり、既知の値よりも未知の値の方が多い不良設定問題である。このような不良設定問題を解くにはいくつかの方法がある。
【００９２】
第１の不良設定解決方法は、（数１１）を用いて解く方法である。
【００９３】
【数１１】
2004045189.tif 000013
【００９４】
ここで、＋記号はムーアペンローズの一般化逆行列である。一般化逆行列とは、非正方行列に対し、近似的に逆行列を与える。複数存在する解の候補のうち、ノルム最小という拘束条件を満たす解が得られることに相当する。
【００９５】
第２の不良設定解決方法は、未知の値のｒ（λ）を既知の値と同次元、即ち３つの基底関数で表現しておき、その基底関数への係数を求める方法である。この方法であれば、既知の値と未知の値の次元数が同次元のため、解が一意に決まる。
【００９６】
次に、この第２の不良設定解決方法の手順について、図６を用いて説明する。（数１）、（数２）及び（数３）を用いて色補正後ＨＬＳ色（Ｈｔａｒｇｅｔ、Ｖｔａｒｇｅｔ、Ｃｔａｒｇｅｔ）からＸＹＺ三刺激値（Ｘｔａｒｇｅｔ、Ｙｔａｒｇｅｔ、Ｚｔａｒｇｅｔ）を求める（Ｓ６０１）。
【００９７】
分光色を表現するのに適切な基底関数を３つ予め設定する（Ｓ６０２）。基底関数の決定は、例えば、自然界に存在する様々な物体の分光色データベースを主成分分析し、上位３位のベクトルを基底関数とする方法が一般的である。これらの基底関数をｂａｓｅ１（λ）、ｂａｓｅ２（λ）及びｂａｓｅ３（λ）とおく。基底関数により、（数１０）を（数１２）に書き換えることができる。
【００９８】
【数１２】
2004045189.tif 000014
【００９９】
また、（数１２）を解くには、（数１３）に変換し、３次元の基底係数ベクトル（ｂ１、ｂ２、ｂ３）を求める（Ｓ６０３）。
【０１００】
【数１３】
2004045189.tif 000015
【０１０１】
ここで、−１は通常の逆行列を意味する。
【０１０２】
（数１４）により基底係数を分光色ｒ（λ）に変換する（Ｓ６０４）。得られた分光色が（Ｘｔａｒｇｅｔ、Ｙｔａｒｇｅｔ、Ｚｔａｒｇｅｔ）を満たす色補正後分光色である。
【０１０３】
【数１４】
2004045189.tif 000016
【０１０４】
次に、分光画像推定部１０１の詳細な動作を説明する。画像入力手段が分光色の画像データを直接取得する場合には、分光画像推定部１０１は必要なく、入力画像を分光色の画像データとして用いればよい。これは、ハイパースペクトルカメラと呼ばれるカメラで、撮影を行えばよい。しかし、それ以外の画像入力手段であった場合には、入力画像から分光色の画像データを推定する処理が必要となり、この処理を行うのが分光画像推定部１０１である。
【０１０５】
入力画像をＲＧＢ信号とすると、ＲＧＢ信号と分光反射率には（数１５）の関係が成り立つ。
【０１０６】
【数１５】
2004045189.tif 000017
【０１０７】
ここで、（Ｒｄ、Ｇｄ、Ｂｄ）は、実際のカメラ信号（Ｒ、Ｇ、Ｂ）から非線形性を排除したＲＧＢ信号を示し、Ｓ（λ）は、入力画像撮影時の照明の分光分布を示し、ｉがＲ、Ｇ及びＢとなるＣｉ（λ）は、カメラの各フィルタの分光感度を示し、ｒ（λ）は、求めるべき分光反射率を示す。
【０１０８】
分光反射率を予め設定されたｉが、１、２及び３となる３つの基底関数Ｏｉ（λ）で表すと、（数１５）は（数１６）に変換され、分光反射率を求める問題は基底係数（ａ、ｂ、ｃ）を求めることが必要となる。基底関数は通常、自然界に存在するさまざまな物体の分光反射率を主成分分析し、その上位３位の表現ベクトルを使用する。
【０１０９】
【数１６】
2004045189.tif 000018
【０１１０】
（数１６）を離散的な表現に変換することで、（数１７）となるため、（Ｒｄ、Ｇｄ、Ｂｄ）から（ａ、ｂ、ｃ）を行列Ｂの逆行列で求めることがでる。ここで行列Ｂは、分光感度、照明分光色及び基底関数から構成される行列である。
【０１１１】
【数１７】
2004045189.tif 000019
【０１１２】
以上、（数１７）得られた基底変数（ａ、ｂ、ｃ）を（数１８）に代入することで、分光反射率を求めることができる。
【０１１３】
【数１８】
2004045189.tif 000020
【０１１４】
以上の処理を入力画像の各画素に対して行うことにより、分光色の画像データを推定することができる。
【０１１５】
なお、色補正を行うのは、彩度に限らず、色相及び明度を補正しても同様の効果を得ることができる。
【０１１６】
なお、本発明は、図１による色補正装置の動作を、ＣＤ−ＲＯＭなどの記憶媒体に記憶し、記憶媒体のプログラムをＰＣ上のＲＡＭにダウンロードし、ＰＣ上のＣＰＵに処理を行わせるものである。
【０１１７】
また、色補正装置の動作を、スキャナ、ディジタルスチルカメラ、ディジタルムービー、マルチスペクトルカメラ又はハイパースペクトルカメラなどの画像入力手段及びディスプレイ又はプリンタなどの画像出力手段に内蔵して処理を行っても良い。
【０１１８】
なお、記憶媒体はＣＤ−ＲＯＭに限らず、ＩＣカード、メモリーカード、フレキシブルディスク、ＤＶＤ−ＲＯＭ、ハードディスク及び磁気テープなど、プログラムを記録できるものであれば、同様に実施することができる。
【０１１９】
以上のように、本実施の形態によれば、分光色の画像データを鮮やかに色補正する際に、分光色の画像データを低次元の色空間に変換して、低次元色空間上で色補正してから、分光色の画像データに逆変換することにより、簡単に分光色の画像データを好ましい画像に色補正することができる。
【０１２０】
（実施の形態２）
第２の実施の形態は、第１の実施の形態と同様に、画像入力手段から推定して得られた分光色の画像データに対し、好ましい色補正を施してから出力する色補正装置を開示する。
【０１２１】
図７は、第２の実施の形態による分光色の画像データを用いた色補正装置の構成を示す図である。
【０１２２】
図７において、７０１は、分光色の画像データにおいて、顕著なピークを持つ分光色の画素を変更色領域として設定し出力する。分光色変更色領域設定部を示し、７０２は、分光色の画像データの変更色領域において、彩度の値を高くする分光色補正を行う分光色補正部を示す。
【０１２３】
第２の実施の形態による分光色の画像データを用いた色補正装置は、分光画像推定部１０１、分光色変更色領域設定部７０１及び分光色補正部７０２を含む構成である。
【０１２４】
以下に、第２の実施の形態による分光色の画像データを用いた色補正装置の動作について説明する。
【０１２５】
分光画像推定部１０１では、入力画像から分光色の画像データを推定する。その詳細な動作は、第１の実施の形態に記載の分光画像推定部１０１と同様である。
【０１２６】
分光色変更色領域設定部７０１は、分光色の画像データより顕著なピークを持つ分光色の画素を変更色領域として設定し出力する。顕著なピークを持つ分光色は、分光色の２次差分値を算出し、２次差分値から極大、極小点を求める。ここで、分光反射率の大きい極大値を求め、それらを囲む２つの極小点の幅が予め設定した閾値内であれば、顕著なピークを持つと判断する。
【０１２７】
分光色補正部は、変更色領域として設定されている画素の分光色を、ピーク幅が狭く、純度の高い色に彩度変換した画像データを出力する。
【０１２８】
ここで、彩度の高い色は、分光色の形状が比較的波長幅が狭い。これを純度が高いという。図８に分光反射率と波長の関係を示す。縦軸は、分光反射率を示し、横軸は波長を示す。８０１は、色補正後の分光色である色補正後分光反射率を示し、８０３は、色補正前のピークを持つ分光色の波長幅を定義する波長幅定義位置を示し、８０４は、色補正後分光反射率８０１の新波長幅定義位置を示す。
【０１２９】
従って、分光反射率の彩度を上げたい場合に、図８に示すように、分光反射率のピークを検出し、そのピークの波長幅を狭める処理を行えばよい。
【０１３０】
分光色補正部の処理の手順を図９に示す。分光反射率のピークを検出し、その波長幅をλｐｅａｋとする（Ｓ９０１）。検出方法はいくつか存在するが、例えば、分光色の２次差分値を順次求め、０に近い値を持つ極大値のうち、最も大きい分光反射率であるピークの波長幅を定義する波長幅定義位置８０３を求める（Ｓ９０２）。
【０１３１】
波長幅定義位置８０３は、λ１及びλ２からなる。これらはピークを囲む２つの極小点位置を求めればよい。波長幅定義位置８０３を、λ１及びλ２からΔλだけピーク中心に移動させた位置を、新波長幅低地位置とする（Ｓ９０３）。
【０１３２】
Δλの値は、予め設定しておく。鮮やかさを上昇させるには、Δλの値を大きくすればよい。この新波長幅定義位置８０４を、それぞれλ１ｎｅｗ、λ２ｎｅｗとする。λ１による分光反射率ｒ１、λ１ｎｅｗによる分光反射率ｒ２を求め、その比であるｒ１／ｒ２により分光反射率を変化させた割合値ｒａｔｉｏ１とする。λ２及びλ２ｎｅｗについても、同様にして変化割合値ｒａｔｉｏ２を求める（Ｓ９０４）。
【０１３３】
λ１ｎｅｗからλｐｅａｋまでの波長領域の分光反射率に、λ１ｎｅｗでｒａｔｉｏ１、λｐｅａｋで１となるように線形に割り振られた重みを乗算することで、分光反射率を変換する。λｐｅａｋからλ２ｎｅｗまでの波長領域の分光反射率についても同様にｒａｔｉｏ２を用いて分光反射率を変換する。λ１からλ１ｎｅｗまでの波長領域は、λ１による分光反射率に変換する。同様にして、λ２からλ２ｎｅｗまでの波長領域は、λ２による分光反射率に変換する。上記以外の波長領域は、色補正前の分光反射率を色補正後分光反射率８０１とする（Ｓ９０５）。
【０１３４】
以上の手順により、分光色の波長幅を狭めた鮮やかな分光反射率が求められる。
【０１３５】
ただし、本処理では、色の純度は上昇するが、明るさが減少する。従って、明るさは、明度変化で指定、又は、Ｓ９０５において、λｐｅａｋでの重みが１以上となる値を割り振るようにして、明るさを上昇させても良い。
【０１３６】
なお、ピークにおける波長幅の狭い分光色に変換するために、曲線を代表する複数の制御点の座標を元にして、曲線全体はその制御点間を滑らかに接続するスプライン曲線により色補正前の分光色を近似しておき、その制御点を制御する方法を用いても良い。
【０１３７】
なお、分光色を分光反射率としているが、分光放射輝度、分光光度及び分光透過率についても同様に実施することができる。
【０１３８】
なお、色補正を行うのは、波長幅の形状変化に限らず、彩度、色相及び明度を補正しても同様の効果を得ることができる。
【０１３９】
なお、本発明は、図７による色補正装置の動作を、ＣＤ−ＲＯＭなどの記憶媒体に記憶し、記憶媒体のプログラムをＰＣ上のＲＡＭにダウンロードし、ＰＣ上のＣＰＵに処理を行わせるものである。
【０１４０】
また、色補正装置の動作を、画像入力手段及び画像出力手段に内蔵して処理を行っても良い。
【０１４１】
なお、記憶媒体はＣＤ−ＲＯＭに限らず、ＩＣカード、メモリーカード、フレキシブルディスク、ＤＶＤ−ＲＯＭ、ハードディスク及び磁気テープなど、プログラムを記録できるものであれば、同様に実施することができる。
【０１４２】
以上のように、本実施の形態によれば、分光色の画像データをより鮮やかに色補正する際に、分光色の形状を波長幅の狭い形状に変化させることで、簡単に分光色の画像データを好ましい画像に色補正することができる。
【０１４３】
（実施の形態３）
本第３の実施の形態では、ＨＬＳ空間の代わりにＬ＊ａ＊ｂ＊空間により色補正を行う。
【０１４４】
図１０は、第３の実施の形態による分光色の画像データを用いた色補正装置の構成を示す図である。
【０１４５】
図１０において、１００１は、分光色の画像データの各画素値を、Ｌ＊ａ＊ｂ＊空間に変換するＬ＊ａ＊ｂ＊空間変換部を示し、１００２は、Ｌ＊ａ＊ｂ＊空間において、彩度が予め決められた閾値よりも高い領域を変更すべき色として設定するＬ＊ａ＊ｂ＊空間変更色領域設定部を示し、１００３は、変更色領域において彩度の値を高くする色補正を行うＬ＊ａ＊ｂ＊色補正部を示し、１００４は、色補正後のＬ＊ａ＊ｂ＊空間の各画素を元の分光空間に逆変換するＬ＊ａ＊ｂ＊分光空間逆変換部を示す。
【０１４６】
また、分光画像推定部１０１は、他の実施の形態と同様の構成であるので、ここでの説明は省略する。
【０１４７】
なお、本実施の形態では、分光色の画像データとして、各画素値が分光反射率である場合を取り扱う。その他の分光色の画像データとして、各画素値が分光放射輝度、分光光度又は分光透過率であった場合にも、同様に本実施の形態を適用することができる。
【０１４８】
以下に、図１０を用いて動作を説明する。
【０１４９】
分光画像推定部１０１は、他の実施の形態と同様の動作であり、ここでの説明は省略する。
【０１５０】
Ｌ＊ａ＊ｂ＊空間変換部１００１は、分光色の画像データの各画素を（数１）によって元の次元より低次元なＬ＊ａ＊ｂ＊空間に変換して、Ｌ＊ａ＊ｂ＊画像を算出する。
【０１５１】
Ｌ＊ａ＊ｂ＊空間変更色領域設定部１００２は、Ｌ＊ａ＊ｂ＊画像において、予め決められた閾値ａ＊ｔｈｒｅ及びｂ＊ｔｈｒｅ以上のａ＊及びｂ＊の値を持つ画素を変更色領域として設定し出力する。
【０１５２】
Ｌ＊ａ＊ｂ＊色補正部１００３は、変更色領域として設定されている画素の彩度に対し、彩度を予め決められた量Δａ＊及びΔｂ＊を加算する。このとき、画素のａ＊及びｂ＊の値が、それぞれ閾値ａ＊ｔｈｒｅ及びｂ＊ｔｈｒｅに近い場合は、Δａ＊及びΔｂ＊の値を０にし、閾値ａ＊ｔｈｒｅ及びｂ＊ｔｈｒｅより大きい場合は、Δａ＊及びΔｂ＊の値を大きい値にすることで、より自然な色変換を行うことができる。
【０１５３】
Ｌ＊ａ＊ｂ＊分光空間逆変換部１００４は、Ｌ＊ａ＊ｂ＊色補正部１００３から出力された色補正後のＬ＊ａ＊ｂ＊画像をより高次元の分光空間に逆変換し、色補正後の分光色の画像データとして出力する。
【０１５４】
ここで、Ｌ＊ａ＊ｂ＊空間を分光空間に逆変換するアルゴリズムは、第１の実施の形態に記載したＨＬＳ空間を分光空間に逆変換するアルゴリズムを、ＨＬＳをＬ＊ａ＊ｂ＊に置き換えれば同様の動作で実現できる。
【０１５５】
なお、色補正を行うのは、彩度に限らず、色相及び明度を補正しても同様の効果を得ることができる。
【０１５６】
なお、本発明は、図１０による色補正装置の動作を、ＣＤ−ＲＯＭなどの記憶媒体に記憶し、記憶媒体のプログラムをＰＣ上のＲＡＭにダウンロードし、ＰＣ上のＣＰＵに処理を行わせるものである。
【０１５７】
また、色補正装置の動作を、画像入力手段及び画像出力手段に内蔵して処理を行っても良い。
【０１５８】
なお、記憶媒体はＣＤ−ＲＯＭに限らず、ＩＣカード、メモリーカード、フレキシブルディスク、ＤＶＤ−ＲＯＭ、ハードディスク及び磁気テープなど、プログラムを記録できるものであれば、同様に実施することができる。
【０１５９】
以上のように、本実施の形態によれば、分光色の画像データをより鮮やかに色補正する際に、分光色の画像データを低次元の色空間に変換して、低次元色空間上で色補正してから、分光色の画像データに逆変換することにより、簡単に分光色の画像データを好ましい画像に色補正することができる。
【０１６０】
【発明の効果】
以上のように本発明によれば、分光色の補正が簡易な色補正装置及び色補正方法を提供することができる。
【図面の簡単な説明】
【図１】本発明第１の実施の形態による色補正装置を示す図
【図２】同第１の実施の形態による第１の分光空間への変換方法を示す図
【図３】同第１の実施の形態による分光反射輝度と波長の関係を示す図
【図４】同第１の実施の形態による変化関数と波長の関係を示す図
【図５】同第１の実施の形態による第２の分光空間への変換方法を示す図
【図６】同第１の実施の形態による第３の分光空間への変換方法を示す図
【図７】同第２の実施の形態による色補正装置を示す図
【図８】同第２の実施の形態による分光反射率の補正を示す図
【図９】同第２の実施の形態による分光反射率の補正方法を示す図
【図１０】同第３の実施の形態による色補正装置を示す図
【図１１】（ａ）従来の色補正装置による分光色表示枠を示す図
（ｂ）従来の色補正装置によるｘｙ色度図表示枠を示す図
（ｃ）従来の色補正装置によるＨＣＶ表示枠を示す図
【符号の説明】
１０１　分光画像推定部
１０２　ＨＬＳ空間変換部
１０３　ＨＬＳ空間変更色領域設定部
１０４　ＨＬＳ色補正部
１０５　ＨＬＳ分光空間逆変換部
３０１　ｓｂ（λ）
３０２　分割波長領域
４０１　線形補間領域
７０１　分光色変更色領域設定部
７０２　分光色補正部
８０１　色補正後分光反射率
８０３　波長幅定義位置
８０４　新波長幅定義位置
１００１　Ｌ＊ａ＊ｂ＊空間変換部
１００２　Ｌ＊ａ＊ｂ＊空間変更色領域設定部
１００３　Ｌ＊ａ＊ｂ＊色補正部
１００４　Ｌ＊ａ＊ｂ＊分光空間逆変換部
【図１】
2004045189.tif 000021
【図２】
2004045189.tif 000022
【図３】
2004045189.tif 000023
【図４】
2004045189.tif 000024
【図５】
2004045189.tif 000025
【図６】
2004045189.tif 000026
【図７】
2004045189.tif 000027
【図８】
2004045189.tif 000028
【図９】
2004045189.tif 000029
【図１０】
2004045189.tif 000030
【図１１】
2004045189.tif 000031