(57)【要約】【課題】谷線・尾根線を含む線形構造が明瞭で、山地、平野、段丘崖の区分も可能であり、地形の違いが明瞭に識別できる画像を作成する。【解決手段】傾斜方向と傾斜量の空間的分布を同時に表現するために、傾斜面の法線ベクトルをＨＳＩ円錐カラーモデルの底面に描いたウルフ網上の点に対応させることにより、斜面ベクトル図を作成した。この図は、単独で傾斜方向、傾斜量、地形の変化を表現することが可能で、線形構造などの特徴的地形の抽出に適している。さらに、斜面ベクトル図の立体視対を作成し立体視することで、断層や環状構造などの特徴的地形が容易に抽出できる。【選択図】        図８
2004362087.tif 000002
【特許請求の範囲】
【請求項１】
主制御部は、地図上の格子点の位置を識別する位置識別子に対応して格子点における標高値を記憶した標高モデルファイルを参照して、各格子点の位置識別子に基づき、該当する格子点の位置識別子における標高値と、該格子点の周りの複数の格子点の位置識別子における標高値を読み取り、これら標高値に基づき該格子点における斜面の傾斜方向と傾斜角を表す斜面ベクトルデータを求めるステップと、
主制御部は、求められた斜面ベクトルデータに基づき、色相Ｈ、彩度Ｓ及び明度Ｉを求めるステップと、
主制御部は、求められた色相Ｈ、彩度Ｓ及び明度Ｉに基づき、ＨＳＩ円錐カラーモデルからＲＧＢカラーモデルに変換するステップと、
主制御部は、変換されたＲＧＢカラーモデルの各原色を表す色画像データを、位置識別子に対応して色画像ファイルに記憶するステップと、
主制御部は、色画像ファイルから各位置識別子に基づいて色画像データを読み取り、該色画像データを地図上の位置識別子が示す格子点に対応させることにより斜面ベクトル図を作成し、表示部に表示又は出力部に出力するステップと、
を含む画像作成方法。
【請求項２】
主制御部は、傾斜角の最大値及び／又は最小値を入力するステップと、
主制御部は、求められた彩度Ｓを入力された最大値及び／又は最小値により値域幅を伸張して変換するステップと
をさらに含む請求項１に記載の画像作成方法。
【請求項３】
主制御部は、さらに、色画像ファイルから読み取った色画像データに基づき、立体視対の画像データを求め、位置識別子に対応して該画像データを立体視画像ファイルに記憶し、該画像データにより立体視対の斜面ベクトル図を作成し、表示部に表示又は出力部に出力する請求項１に記載の画像作成方法。
【請求項４】
主制御部は、各位置識別子に基づき斜面ベクトルデータ中の傾斜方向を読み取り、該傾斜方向に基づきＨＳＩ円錐カラーモデルの断面図又はＲＧＢカラーモデルの断面図の外周近傍の色に対応する色画像データを求め、位置識別子に対応して色画像ファイルに記憶するステップをさらに含み、
主制御部は、該色画像データを用いることにより傾斜方位図を作成し、表示部に表示又は出力部に出力する請求項１に記載の画像作成方法。
【請求項５】
主制御部は、各位置識別子に基づき斜面ベクトルデータ中の傾斜角を読み取り、該傾斜角に基づき濃淡に対応する色画像データを求め、位置識別子に対応して色画像ファイルに記憶するステップをさらに含み、
主制御部は該色画像データを用いることにより傾斜量図を表示部に表示又は出力部に出力する請求項１に記載の画像作成方法。
【請求項６】
前記斜面ベクトルデータを求めるステップは、周りの格子点の標高値から近似平面を求め該近似平面の法線ベクトルの傾斜方向θ及び傾斜角φを求める請求項１に記載の画像作成方法。
【請求項７】
前記色相Ｈ、彩度Ｓ及び明度Ｉを求めるステップは、以下の式により各値を求める請求項１に記載の画像作成方法。
Ｈ＝θ
Ｓ＝ｔａｎ（（９０−φ）／２）
Ｉ＝Ｉ０＋（１−Ｉ０）Ｓ  （０≦Ｉ０≦１．０）
【請求項８】
主制御部は、傾斜方向θ及び傾斜角φに基づき斜面ベクトルの彩度色相円への投影点を求め、その投影点の色に従い格子点の色画像データを求める請求項１に記載の画像作成方法。
【請求項９】
主制御部は、各格子点の斜面ベクトルデータに基づき、傾斜角の最大値及び／又は最小値を算出するステップと、
主制御部は、求められた彩度Ｓを算出された最大値及び／又は最小値により値域幅を伸張して変換するステップと
をさらに含む請求項１に記載の画像作成方法。
【請求項１０】
主制御部は、使用する色相及び／又は色相の変化幅を調整するステップをさらに含む請求項１に記載の画像作成方法。
【請求項１１】
主制御部は、地図上の格子点の位置を識別する位置識別子に対応して格子点における標高値を記憶した標高モデルファイルを参照して、各格子点の位置識別子に基づき、該当する格子点の位置識別子における標高値と、該格子点の周りの複数の格子点の位置識別子における標高値を読み取り、これら標高値に基づき該格子点における斜面の傾斜方向と傾斜角を表す斜面ベクトルデータを求めるステップと、
主制御部は、求められた斜面ベクトルデータに基づき、色相Ｈ、彩度Ｓ及び明度Ｉを求めるステップと、
主制御部は、求められた色相Ｈ、彩度Ｓ及び明度Ｉに基づき、ＨＳＩ円錐カラーモデルからＲＧＢカラーモデルに変換するステップと、
主制御部は、変換されたＲＧＢカラーモデルの各原色を表す色画像データを、位置識別子に対応して色画像ファイルに記憶するステップと、
主制御部は、色画像ファイルから各位置識別子に基づいて色画像データを読み取り、該色画像データを地図上の位置識別子が示す格子点に対応させることにより斜面ベクトル図を作成し、表示部に表示又は出力部に出力するステップと、
をコンピュータに実行させるため画像作成プログラム。
【請求項１２】
主制御部は、地図上の格子点の位置を識別する位置識別子に対応して格子点における標高値を記憶した標高モデルファイルを参照して、各格子点の位置識別子に基づき、該当する格子点の位置識別子における標高値と、該格子点の周りの複数の格子点の位置識別子における標高値を読み取り、これら標高値に基づき該格子点における斜面の傾斜方向と傾斜角を表す斜面ベクトルデータを求めるステップと、
主制御部は、求められた斜面ベクトルデータに基づき、色相Ｈ、彩度Ｓ及び明度Ｉを求めるステップと、
主制御部は、求められた色相Ｈ、彩度Ｓ及び明度Ｉに基づき、ＨＳＩ円錐カラーモデルからＲＧＢカラーモデルに変換するステップと、
主制御部は、変換されたＲＧＢカラーモデルの各原色を表す色画像データを、位置識別子に対応して色画像ファイルに記憶するステップと、
主制御部は、色画像ファイルから各位置識別子に基づいて色画像データを読み取り、該色画像データを地図上の位置識別子が示す格子点に対応させることにより斜面ベクトル図を作成し、表示部に表示又は出力部に出力するステップと、
をコンピュータに実行させるため画像作成プログラムを記録したコンピュータ読み取り可能な記録媒体。
【請求項１３】
地図上の格子点の位置を識別する位置識別子に対応して格子点における標高値を記憶した標高モデルファイルと、
格子点の位置識別子に対応して格子点におけるＲＧＢカラーモデルの各原色の値を表す色画像データを記憶した色画像ファイルと、
各々の前記ファイルをアクセスし、斜面ベクトル図を作成し、表示部に表示又は出力部に出力するための主制御部と、
を備え、
主制御部は、前記標高モデルファイルを参照して、各格子点の位置識別子に基づき、該当する格子点の位置識別子における標高値と、該格子点の周りの複数の格子点の位置識別子における標高値を読み取り、これら標高値に基づき該格子点における斜面の傾斜方向と傾斜角を表す斜面ベクトルデータを求め、
主制御部は、求められた斜面ベクトルデータに基づき、色相Ｈ、彩度Ｓ及び明度Ｉを求め、
主制御部は、求められた色相Ｈ、彩度Ｓ及び明度Ｉに基づき、ＨＳＩ円錐カラーモデルからＲＧＢカラーモデルに変換し、
主制御部は、変換されたＲＧＢカラーモデルの各原色を表す色画像データを、位置識別子に対応して色画像ファイルに記憶し、
主制御部は、前記色画像ファイルから各位置識別子に基づいて色画像データを読み取り、該色画像データを地図上の位置識別子が示す格子点に対応させることにより斜面ベクトル図を作成し、表示部に表示又は出力部に出力する
画像作成装置。
【請求項１４】
予め定められた明度ＩにおけるＨＳＩ円錐カラーモデルの断面図における点により、色相Ｈで斜面の傾斜方向を表し、彩度Ｓで斜面の傾斜角を表した斜面ベクトル図であって、前記斜面ベクトル図は、
主制御部により、地図上の格子点の位置を識別する位置識別子に対応して格子点における標高値を記憶した標高モデルファイルを参照して、各格子点の位置識別子に基づき、該当する格子点の位置識別子における標高値と、該格子点の周りの複数の格子点の位置識別子における標高値を読み取り、これら標高値に基づき該格子点における斜面の傾斜方向と傾斜角を表す斜面ベクトルデータを求め、
主制御部により、求められた斜面ベクトルデータに基づき、色相Ｈ、彩度Ｓ及び明度Ｉを求め、
主制御部により、求められた色相Ｈ、彩度Ｓ及び明度Ｉに基づき、ＨＳＩ円錐カラーモデルからＲＧＢカラーモデルに変換し、
主制御部により、変換されたＲＧＢカラーモデルの各原色を表す色画像データを、位置識別子に対応して色画像ファイルに記憶し、
主制御部により、色画像ファイルから各位置識別子に基づいて色画像データを読み取り、該色画像データを地図上の位置識別子が示す格子点に対応させることにより斜面ベクトル図を作成し、表示部に表示又は出力部に出力する
ことにより作成された前記斜面ベクトル図。
【請求項１５】
立体視対の画像を含む請求項１４に記載の斜面ベクトル図。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、画像作成装置、画像作成方法及びプログラム及び記録媒体、斜面ベクトル図に係り、特に、特徴的地形抽出のための斜面ベクトルのカラー表示を可能とし、また、その立体視を可能とする画像作成装置、画像作成方法及びプログラム及び記録媒体、斜面ベクトル図に関する。
【０００２】
【従来の技術】
地形図から断裂帯や環状構造などを抽出する（例えば、非特許文献１及び２）とき，地表面の傾斜は重要な情報を提供する。線形構造をはじめとする特徴的地形を把握するために標高図を眺めることや空中写真を立体視することは、以前から知られている有効な手段である。最近では、コンピュータの発達に伴ってデジタル標高モデル（ＤＥＭ：Ｄｉｇｉｔａｌ  Ｅｌｅｖａｔｉｏｎ  Ｍｏｄｅｌ）の整備が進み、それらを用いた地形解析（例えば、非特許文献３、４参照）も活発に行われている。それは、地形情報がアナログからデジタルに変化したことにより、地形の変化量の数値解析や、空間情報の表示システムであるＧＩＳ（Ｇｅｏｇｒａｐｈｉｃ  Ｉｎｆｏｒｍａｔｉｏｎ  Ｓｙｓｔｅｍ）を利用した解析結果の管理・表示が容易に行えるようになったからである。
【０００３】
ところで、一般的な色は、光の３原色である赤（Ｒ）、緑（Ｇ）、青（Ｂ）と、それらの補色であるシアン（Ｃ）、マゼンダ（Ｍ）、黄（Ｙ）、および、黒（Ｋ）、白（Ｗ）の６点を頂点とするＲＧＢ立方体内の一点として表される。今、黒Ｋと白Ｗを結ぶ直線をＩ軸として、この立方体のＲ、Ｙ、Ｇ、Ｃ、Ｂ、Ｍを、白Ｗを含みＩ軸に直交する平面に投影すると、正六角形ができる。この正六角形を円に変換し、その円を底面とし、黒Ｋを頂点とした円錐をＨＳＩ円錐カラーモデルという。ＨＳＩ円錐カラーモデルは、色相Ｈ、彩度Ｓ、明度Ｉで色を表現する。他の表現によると円錐側面の表面の色を底面に投影することでＨＳＩ円錐カラーモデルの各明度Ｉにおける色相Ｈ及び彩度Ｓの色の断面図が得られる。
【０００４】
【非特許文献１】
羽坂・渡辺・池田（１９９１）地形数値データを基にしたリニアメントの抽出−２万５千分の１地形図「余市岳」での応用例−，地質調査所月報，ｖｏｌ．４２，ｎｏ．１，ｐｐ．１１−１７．
【非特許文献２】
田中・幾志・熊谷（１９８５）三宅島の傾斜分布図．国立防災科学技術センター研究速報，ｖｏｌ．６９，ｐｐ．１７−２７．
【非特許文献３】
岩橋（１９９２）ディジタルフィルタを用いた数値地形モデルの地形解析，情報地質，ｖｏｌ．３，ｎｏ．１，ｐｐ．３１−３７．
【非特許文献４】
神谷・黒木・田中（２０００）傾斜量図を用いた地形・地質の判読，情報地質、ｖｏｌ．１１，ｎｏ．１，ｐｐ．１１−２４．
【０００５】
【発明が解決しようとする課題】
一般に、地形を表現する図面には、標高図又は陰影図の他に、傾斜方位図や傾斜量図がある。
傾斜方位図には、傾斜方向のみが表示されている。傾斜量図には、地形の傾斜量のみが表示されている。標高図や陰影図を参考に傾斜方位図を見ると、斜面の連続性の抽出に有効であり、谷線・尾根線を含む線形構造がよく分かる。しかし、平野と山地の区分が不明瞭で、また斜面の凹凸も不明瞭な場合がある。傾斜量図は、平野と山地が明瞭に区分されており、また、谷線・尾根線を含む線形構造や段丘崖も明瞭に認識できる。ただ、山地では、地形の起伏が不明瞭な場合がある。
【０００６】
本発明は、以上の点に鑑み、谷線・尾根線を含む線形構造が明瞭で、山地、平野、段丘崖の区分も可能であり、地形の違いが明瞭に識別できる画像を作成することを目的とする。また、本発明は、作成された斜面ベクトル図により、傾斜方位図と傾斜量図の長所を合わせもち、地形の起伏も表現することを目的とする。
【０００７】
また、一般に、地図上のある点付近の地表面の傾きは、その面に垂直で上向きの３次元単位ベクトル、すなわち法線ベクトルで表現できる。このベクトルをここでは斜面ベクトルと呼ぶ。斜面ベクトルはウルフ網上の１点で表わされる。そこで、この円錐の底面をウルフ網と一致させることで、ウルフ網上の１点として与えられる斜面ベクトルを色で表わすことができる。
【０００８】
なお、図１７に、ウルフ網の説明図を示す。一般に、ウルフ網は、地球の南極に視点をおき、北半球の経線と緯線を投影すると図（ａ）が得られる。また、赤道上に視点をおき、それを極とする子午面に視点と反対側の経線と緯線を投影すると図（ｂ）が得られる。これらをそれぞれ赤道面投影と子午面投影のステレオ網（ｓｔｅｒｅｏｎｅｔ）あるいはウルフ網（Ｗｕｌｆｆ  ｎｅｔ）という。ウルフ網は等角網（ｅｑｕａｌ  ａｎｇｌｅ  ｎｅｔ）ともいわれる。これは、ウルフ網上では、球面上の角度が保たれていることによる。本実施の形態では、特に、赤道面投影のウルフ網が用いられる。
【０００９】
ウルフ網とＨＳＩ円錐カラーモデルの底面を対応させて斜面ベクトルを色で表し、これを地図として表現すると、傾斜方向、傾斜量、地形の変化を同時に表現できる。この図を斜面ベクトル図と名づけた。
本発明は、この斜面ベクトルの色を地図上に表示することにより、傾斜方向と傾斜量の空間分布を同時に表現することを目的とする。さらに、本発明は、このようにして作成した斜面ベクトル図を、地形が立体視できる立体視対として表現することにより、断層帯などの特徴的地形をより容易に抽出できる斜面ベクトル図を形成することを目的とする。
【００１０】
【課題を解決するための手段】
本発明の第１の解決手段によると、
主制御部は、地図上の格子点の位置を識別する位置識別子に対応して格子点における標高値を記憶した標高モデルファイルを参照して、各格子点の位置識別子に基づき、該当する格子点の位置識別子における標高値と、該格子点の周りの複数の格子点の位置識別子における標高値を読み取り、これら標高値に基づき該格子点における斜面の傾斜方向と傾斜角を表す斜面ベクトルデータを求めるステップと、
主制御部は、求められた斜面ベクトルデータに基づき、色相Ｈ、彩度Ｓ及び明度Ｉを求めるステップと、
主制御部は、求められた色相Ｈ、彩度Ｓ及び明度Ｉに基づき、ＨＳＩ円錐カラーモデルからＲＧＢカラーモデルに変換するステップと、
主制御部は、変換されたＲＧＢカラーモデルの各原色を表す色画像データを、位置識別子に対応して色画像ファイルに記憶するステップと、
主制御部は、色画像ファイルから各位置識別子に基づいて色画像データを読み取り、該色画像データを地図上の位置識別子が示す格子点に対応させることにより斜面ベクトル図を作成し、表示部に表示又は出力部に出力するステップと、
を含む画像作成方法、これら各ステップをコンピュータに実行させるため画像作成プログラム及びそのプログラムを記録したコンピュータ読み取り可能な記録媒体が提供される。
【０００１１】
本発明の第２の解決手段によると、
地図上の格子点の位置を識別する位置識別子に対応して格子点における標高値を記憶した標高モデルファイルと、
格子点の位置識別子に対応して格子点におけるＲＧＢカラーモデルの各原色の値を表す色画像データを記憶した色画像ファイルと、
各々の前記ファイルをアクセスし、斜面ベクトル図を作成し、表示部に表示又は出力部に出力するための主制御部と、
を備え、
主制御部は、前記標高モデルファイルを参照して、各格子点の位置識別子に基づき、該当する格子点の位置識別子における標高値と、該格子点の周りの複数の格子点の位置識別子における標高値を読み取り、これら標高値に基づき該格子点における斜面の傾斜方向と傾斜角を表す斜面ベクトルデータを求め、
主制御部は、求められた斜面ベクトルデータに基づき、色相Ｈ、彩度Ｓ及び明度Ｉを求め、
主制御部は、求められた色相Ｈ、彩度Ｓ及び明度Ｉに基づき、ＨＳＩ円錐カラーモデルからＲＧＢカラーモデルに変換し、
主制御部は、変換されたＲＧＢカラーモデルの各原色を表す色画像データを、位置識別子に対応して色画像ファイルに記憶し、
主制御部は、前記色画像ファイルから各位置識別子に基づいて色画像データを読み取り、該色画像データを地図上の位置識別子が示す格子点に対応させることにより斜面ベクトル図を作成し、表示部に表示又は出力部に出力する
画像作成装置が提供される。
【００１２】
本発明の第３の解決手段によると、
予め定められた明度ＩにおけるＨＳＩ円錐カラーモデルの断面図における点により、色相Ｈで斜面の傾斜方向を表し、彩度Ｓで斜面の傾斜角を表した斜面ベクトル図であって、前記斜面ベクトル図は、
主制御部により、地図上の格子点の位置を識別する位置識別子に対応して格子点における標高値を記憶した標高モデルファイルを参照して、各格子点の位置識別子に基づき、該当する格子点の位置識別子における標高値と、該格子点の周りの複数の格子点の位置識別子における標高値を読み取り、これら標高値に基づき該格子点における斜面の傾斜方向と傾斜角を表す斜面ベクトルデータを求め、
主制御部により、求められた斜面ベクトルデータに基づき、色相Ｈ、彩度Ｓ及び明度Ｉを求め、
主制御部により、求められた色相Ｈ、彩度Ｓ及び明度Ｉに基づき、ＨＳＩ円錐カラーモデルからＲＧＢカラーモデルに変換し、
主制御部により、変換されたＲＧＢカラーモデルの各原色を表す色画像データを、位置識別子に対応して色画像ファイルに記憶し、
主制御部により、色画像ファイルから各位置識別子に基づいて色画像データを読み取り、該色画像データを地図上の位置識別子が示す格子点に対応させることにより斜面ベクトル図を作成し、表示部に表示又は出力部に出力する
ことにより作成された前記斜面ベクトル図が提供される。
【００１３】
【発明の実施の形態】
１．序論
（１）使用データ
斜面ベクトル図がさまざまな特徴的地形の抽出に有効であることを示すために、以下では、一例として、１９９７年に国土地理院から刊行された茨城県筑波山周辺の５０ｍＤＥＭを用いて、標高図、陰影図、傾斜方位図、傾斜量図、斜面ベクトル図、斜面ベクトル図の立体視対を作成した。また、特に断裂系と環状構造に対する斜面ベクトル図と斜面ベクトル図の立体視対の効果を検証するために、大分県豊肥地域の万年山周辺と、青森県と秋田県の県境に位置する十和田湖周辺の傾斜量図と斜面ベクトル図の立体視対も作成した。筑波山周辺は、平野と山地が比較的明瞭に分かれた地形を示す。豊肥地域は日本最大の活火山地帯であり、火山地形と東西に走る大小規模の正断層が発達している。十和田湖はカルデラであり、環状構造の例として選んだ。
これらの使用データは一例を示すものであり、これに限られず、様々な地図上の位置とその標記を示すデータを用いることができる。
【００１４】
（２）斜面ベクトル図の作成方法
（斜面ベクトル）
５０ｍＤＥＭから斜面ベクトルを算出した。筑波山周辺のＤＥＭの格子点間隔は、東西（緯線）方向が約５６．０ｍ、南北（経線）方向が約４６．３ｍである。斜面ベクトルの傾斜方向と傾斜角は、対象とする格子点を中心とする９個の格子点に近似平面を当てはめ、その法線から求めた。
直交座標系で求めた法線ベクトルの傾斜方向（θ）と傾斜角（φ）は次式で与えられる。
【００１５】
【数１】
2004362087.tif 000003
ここで、ａとｂは近似平面ｚ＝ａｘ＋ｂｙ＋ｃの係数である。
【００１６】
（斜面ベクトルのステレオ投影とＨＳＩ色空間）
斜面ベクトルの空間分布を地形図上に表現するために、色相（Ｈ：Ｈｕｅ）、彩度（Ｓ：Ｓａｔｕｒａｔｉｏｎ）、明度（Ｉ：Ｉｎｔｅｎｓｉｔｙ）の３成分で構成されるＨＳＩ円錐カラーモデルを利用した。
図１に、斜面ベクトルのステレオ投影点とＨＳＩ円錐カラーモデルの対応関係の説明図を示す。この図では、傾斜量と傾斜方向を示すベクトルのステレオ投影点Ｐ’とＨＳＩ双円錐カラーモデルの関係（Ｒ：赤、Ｙ：黄、Ｇ：緑、Ｃ：シアン、Ｂ：青、Ｍ：マゼンダ、Ｗ：白、Ｂｌ：黒、Ｐ：ベクトルと半球面の交点、Ｈ：色相（０≦Ｈ≦２π）、Ｓ：彩度（０≦Ｓ≦１）、Ｉ：明度（０≦Ｉ≦１）、φ：傾斜角、θ：傾斜方向、Ｒの方向が東（Ｅ））が示される。ウルフ網とＨＳＩ円錐カラーモデルの底面を対応させたとき、斜面ベクトル（θ，φ）のステレオ投影点とＨＳ面の座標（Ｈ，Ｓ）は次のように対応する。
【００１７】
【数２】
2004362087.tif 000004
この場合、ＨＳＩ空間ではＩ＝１．０である。この明度Ｉを、
【００１８】
【数３】
2004362087.tif 000005
とすることもできる。ここで、Ｉ０は円錐の頂点の明度である。適当なＩ０を与えると、傾斜角度が明度でも表現される。Ｉ０は図が見やすくなるように任意に決めればよい。
【００１９】
図２に、ＨＳＩ円錐カラーモデルで着色されたウルフ網の図を示す。左から、円錐カラーモデルの明度がＩ０＝１．０、０．５、０．０であり、経線と緯線の間隔はそれぞれ３０°である。
ところで、実際に斜面ベクトルの傾き角が０°に近くなることはほとんどないので、（２）式だとＨＳＩ円錐カラーモデルの底面を有効に使えない場合がある。そこで、φに最大傾斜を与えるφｍｉｎを代入し、Ｓの最大輝度Ｓｓｔｅｅｐｅｓｔを、
【００２０】
【数４】
2004362087.tif 000006
とすると、ＨＳＩ円錐カラーモデルの底面で高輝度の部分を一層有効に使うことができる。また、φに最小傾斜を与えるφｍａｘを代入し、Ｓの最小輝度Ｓｇｅｎｔｌｅを
【００２１】
【数５】
2004362087.tif 000007
として、これを用いてもよい。さらに、これら両方の値を用いて傾斜角の値域幅を求めることができる。例えば、通常は０≦Ｓ≦１であるが、このスケール幅を０≦Ｓ≦Ｓｓｔｅｅｐｅｓｔ、Ｓｇｅｎｔｌｅ≦Ｓ≦Ｓｓｔｅｅｐｅｓｔ又はＳｇｅｎｔｌｅ≦Ｓ≦１で示す範囲を０≦Ｓ≦１で示す広範囲に変換するようにすればよい。
【００２２】
以上のように、斜面ベクトル図では、斜面ベクトルの色相Ｈで斜面の傾斜方向を表し、彩度Ｓで斜面の傾斜角を表すことができる。図２に示すように、予め定められた明度ＩにおけるＨＳＩ円錐カラーモデルの断面図（ウルフ網と重なる彩度最大色相円）を使用して、図１に示すように、斜面ベクトルと彩度最大色相円を含む球面との交点Ｐを、彩度最大色相円内にステレオ投影した投影点Ｐ’の位置の色により、斜面の傾斜方向及び傾斜角が表現される。すなわち、この例では、地図上の格子点における斜面が急勾配になるにつれて、彩度Ｓが大きくなり（１に近づき）濃い色で表現され、一方、斜面が平坦になるにつれて、彩度Ｓが小さくなり（０に近づき）薄い色で表現される（なお、濃淡が逆でもよい）。地図上の格子点における斜面の方向は、色相Ｈを使用して、図２に示すような円のように、Ｒ：赤、Ｙ：黄、Ｇ：緑、Ｃ：シアン、Ｂ：青、Ｍ：マゼンダ、にそれぞれ対応して東西南北が示される。本実施の形態では、一例として、東（Ｅ）方向に傾いている斜面を赤（Ｒ）で表し、西（Ｗ）方向に傾いている斜面をシアン（Ｃ）で表し、北（Ｎ）方向に傾いている斜面を緑（Ｇ）と黄（Ｙ）の中間の色で表し、南（Ｓ）方向に傾いている斜面を青（Ｂ）とマゼンダ（Ｍ）の中間の色で表している。なお、後述の図１４及びその説明箇所を参照すると、傾斜方向及び傾斜角の色画像との対応の一例が示されている。また、色の傾斜角との対応及び傾斜方向との対応は、この例に限らず適宜定めることができる。
【００２３】
（３）立体視対の作成
斜面ベクトルを表現する図で、さらに、同時に地形の起伏を表現したい場合、斜面ベクトル図と地形の起伏情報との重ね合わせが必要である。地形の起伏は、しばしば陰影図として表示される。そこで、斜面ベクトルはＨＳ面で表し、Ｉを陰影図に利用することが可能である。しかし、この場合、当然のことながら、明度の変化として斜面ベクトル図を表現することはできない。この欠点を除くために、斜面ベクトルを示す色で地表面を着色し、この地形の立体視対を作成した。これを斜面ベクトル図の立体視対と呼ぶ。立体視対は、左目で見た像と、右目で見た像の２枚の画像からなる。この対画像をそれぞれの目で見ることにより、脳内で３次元像が認識される。
【００２４】
２．ハードウェア
図３は、本実施の形態に関するハードウェアの構成図である。
このハードウェアは、中央処理装置（ＣＰＵ）である主制御部１、入力部２、出力部３、表示部４、記憶部５及び入出力制御部６を備える。
【００２５】
主制御部１は、解析プログラム５５及び画像表示用プログラム５６等の各種プログラムを記憶部５からロードして、適宜の処理を実行する。入力部２は、例えば、設定、指示、条件等の画像形成に必要な各種データ等を入力するための入力手段であり、キーボード、マウス、ポインティングディバイス等の適宜の装置が用いられる。出力部３は、カラープリンタ、又は、各種記憶媒体、他の装置、ネットワーク等にデータを出力する出力装置である。表示部４は、例えば、ＲＧＢ画像や立体画像等を画面に表示するためのカラーディスプレイ装置である。記憶部５は、ハードディスク等の適宜のメモリデバイスであり、標高モデルファイル５１、斜面ベクトルファイル５２、ＲＧＢ画像ファイル５３、立体視画像ファイル５４、解析プログラム５５、画像表示用プログラム５６を含む。記憶部５には、初期設定、パラメータ等の各種データや、画像形成の最終結果や中間結果等のデータが記憶される。入出力制御部６は、主制御部１と、入力部２、出力部３、表示部４及び記憶部との間のデータの送受信を制御する。入出力制御部６は、バス、スター等の適宜のネットワークで構成することもできる。
【００２６】
図４に、標高モデルファイル５１の説明図を示す。標高モデルファイル５１は、地図上の格子点の位置を識別する位置識別子（位置ＩＤ）に対応して格子点における標高値を記憶する。
図５に、斜面ベクトルファイル５２の説明図を示す。斜面ベクトルファイル５２は、各格子点の位置識別子に対応して格子点における斜面の傾斜方向θと傾斜角φを表す斜面ベクトルデータを記憶する。
図６に、ＲＧＢ色画像ファイル５３の説明図を示す。ＲＧＢ色画像ファイル５３は、各格子点の位置識別子に対応して格子点におけるＲＧＢカラーモデルの各原色の値を表す色画像データを記憶する。
図７に、立体視画像ファイル５４の説明図を示す。立体視画像ファイル５４は、各格子点の位置識別子に対応して、左画像及び右画像を含む立体視対の画像データを記憶する。
なお、各データファイルでは、位置識別子は、シーケンス番号を用いているが、これに限らず、２次元マトリクス等の適宜位置を識別する番号・記号を用いることができる。
【００２７】
３．ソフトウェア
図８に、画像形成のためのフローチャートを示す。
まず、主制御部１は、入力ファイル及び出力ファイルを特定し、標高モデルファイル５１を参照して、各格子点の位置識別子における標高値と、その格子点の周りの複数の格子点の位置識別子における標高値を読み取り、入力する（Ｓ１０１）。主制御部は、これら標高値等のデータに基づき各位置識別子における格子点の斜面の傾斜方向θ及び傾斜角φを含む斜面ベクトルデータを求め、そのデータを位置識別子に対応して斜面ベクトルファイル５２に記憶する（Ｓ１０３）。斜面ベクトルデータを求めるには、例えば、前述のように、周りの格子点の標高値から近似平面を求め法線ベクトルの傾斜方向θ及び傾斜角φを求めればよい。
【００２８】
まず、主制御部１は、斜面ベクトルファイル５２を参照して、傾斜角の最大値及び／又は最小値を算出し、表示部４にそれを表示する（Ｓ１０５）。これにより、操作者はこれら最大値及び／又は最小値が適当なＳ（彩度）で斜面ベクトル図が表示できる値かどうかを確認することができる。ここで、操作者は、入力部２により傾斜角の値域幅を入力する（Ｓ１０７）。ここで、値域幅として、算出された最大値及び／又は最小値が適当であれば、エンター等によりそれを使用する指示を入力し、変更する場合は最大値及び／又は最小値を入力する。つぎに、主制御部１は、入力部２又は記憶部５から、明度値Ｉ０等の予め定められた値を入力する（Ｓ１０８）。
【００２９】
つぎに、主制御部１は、求められた斜面ベクトルデータに基づき、色相Ｈ、彩度Ｓ及び明度Ｉを求める（Ｓ１０９）。ここで、色相Ｈ、彩度Ｓ及び明度Ｉを求めるためには、例えば、前述「１．序論」で説明したように、上述のＨ、Ｓ、Ｉ式により各値を求めることができる。さらに、主制御部１は、前述のように、求められた彩度Ｓを入力された最大値及び／又は最小値により値域幅を伸張して変換することもできる。
【００３０】
なお、ステップＳ１０５、又は、ステップＳ１０５とＳ１０７の両方を省略して、傾斜角の値域幅として予め記憶部５に定められたデフォルト値を用いることで、ステップＳ１０９における傾斜角の値域幅による修正を省略しても良い。
【００３１】
さらに、色相Ｈも、使用する色相とその変化幅を調整できる。図１８に、色相の変化幅の調整についての説明図を示す。図（ａ）のように、通常、Ｒ、Ｇ、Ｂがそれぞれ１２０度間隔であるが、例えば、図（ｂ）のように、その間隔を９０度・９０度・１８０度としたり、図（ｃ）のように、９０度・１８０度・９０度とすることができる。さらに、Ｒ、Ｇ、Ｂのみを用いて、Ｙ，Ｃ，Ｍを用いないようにしても良い。主制御部１は、入力部２又は記憶部５からこれらの設定を行うことができ、その設定に応じて変化幅を調整する計算を実行する。
【００３２】
つぎに、入力部２により画像形成処理を選択する（Ｓ１１１）。画像形成処理には、例えば、（ａ）傾斜ベクトル図、（ｂ）傾斜方位図、（ｃ）傾斜量図の形成処理が含まれる。
ステップＳ１１１で「（ａ）傾斜ベクトル図作成」が選択された場合、主制御部１は、求められた色相、彩度及び明度に基づき、ＨＳＩ円錐カラーモデルからＲＧＢカラーモデルに変換することにより斜面ベクトル図を作成する（Ｓ１１３）。主制御部１は、変換されたＲＧＢカラーモデルの各色に対応する色画像データを位置識別子に対応して色画像ファイル５３に記憶する（Ｓ１１５）。
【００３３】
さらに、主制御部１は、色画像ファイル５２から読み取った色画像データに基づき、立体視対の画像データを求め位置識別子に対応して立体視画像ファイル５４に記憶し、該画像データを表示部４に表示又は出力部３に出力するようにしてもよい（Ｓ１１７）。なお、ステップＳ１１７の処理は、省略することもできる。
【００３４】
主制御部１は、色画像ファイル５３から各位置識別子に対応した色画像データを読み取り、該色画像データを地図上の位置識別子が示す格子点に貼り付ける等により対応させることで斜面ベクトル図を作成し、表示部４に表示及び／又は出力部３に出力する（Ｓ１１９）。ここで、ステップＳ１１７を実行した際は、立体視画像ファイル５４から各位置識別子に対応した主体部画像データを読み取り、同様に立体視対の斜面ベクトル図を作成し、表示部４に表示及び／又は出力部３に出力する。
【００３５】
一方、ステップＳ１１１で「（ｂ）傾斜方位図作成」が選択された場合、主制御部１は、位置識別子に対応して斜面ベクトルデータ中の傾斜方向θを読み取り、傾斜方向θに基づきＨＳＩ円錐カラーモデルの断面図又はＲＧＢカラーモデルの断面図の外周近傍の色に対応する色画像データを求める（Ｓ１２１）。その後ステップＳ１１５以降の処理により、同様に、色画像データを位置識別子に対応して色画像ファイル５３に記憶し、必要に応じて立体視画像を作成し、得られた傾斜方位図を表示部４に表示又は出力部３に出力する。
【００３６】
また、ステップＳ１１１で「（ｃ）傾斜量図作成」が選択された場合、主制御部１は、位置識別子に対応して斜面ベクトルデータ中の傾斜角φを読み取り、傾斜角φに基づき濃淡に対応する色画像データを求める（Ｓ１２３）。その後ステップＳ１１５以降の処理により、同様に、色画像データを位置識別子に対応して色画像ファイル５３に記憶し、必要に応じて立体視画像を作成し、得られた傾斜量図を表示部４に表示又は出力部３に出力する。このとき、同時に使用した円形カラーチャートも出力できる。
なお、他の画像形成処理をさらに選択する場合、ステップＳ１１１に戻りその後処理を繰り返しても良い。また、（ａ）傾斜ベクトル図、（ｂ）傾斜方位図、（ｃ）傾斜量図の各色画像データをそれぞれ別個のファイルに記憶するようにしても良い。さらに、各ステップにおける所定ファイルへのデータの記憶を適宜行ったり、逆に、それを省略するようにしてもよい。
【００３７】
（ＨＳＩ−ＲＧＢ変換）
以下に、ステップＳ１１３の斜面ベクトル図作成におけるＨＳＩ−ＲＧＢ変換について詳述する。
実際に画像に着色する際には、最小二乗近似平面のベクトルから求めたＨＳＩ（Ｈ＝φ、Ｓ＝θ）の値をＲＧＢ値に変換しなければならない。ＨＳＩ−ＲＧＢ変換式を以下に示す。図１に示すとおり、Ｈ（φ）＝０は、直交座標系のｘ軸方向と、ＨＩＳ双円錐カラーモデルのＲ軸方向と一致する。また、使用するカラーモデルはＩの値によって決まる。
【００３８】
【数６】
2004362087.tif 000008
【００３９】
【数７】
2004362087.tif 000009
（０．０≦ｒ、ｇ、ｂ≦１．０、０≦Ｈ≦２π、０．０≦Ｓ≦１．０、０．０≦Ｉ≦１．０とする）
【００４０】
（立体視画像の作成）
以下に、ステップＳ１１７における立体視画像の作成処理について説明する。傾斜量と傾斜方位を同時に表現する画像に正確な地形の起伏を合わせて表現したい場合、さらにそれらの情報との重ねあわせが必要である。地形の起伏に関しては陰影図が有効であるが、色で情報を表現した合成画像に、さらに影を重ねると合成画像の情報が隠されてしまい全体として意味をなさない画像になってしまう。そこで、作成した合成画像を立体視することで、ここに挙げる地形特性のすべてを表現する。
【００４１】
図９に、立体視対の作成方法の説明図を示す。立体視対の作成とは、両眼視差を考慮した画像を２枚作成することである。つまり、右目で見た像と、左目で見た像を一枚の画像から幾何学的に作成し、それぞれの画像をそれぞれの側の目で見ることにより、脳内で立体画像が合成される。
対象とする画像上の１点、Ｉｍ（ｘ，ｙ，ｚ）を、視点Ｅｒ（右眼）、Ｅｌ（左眼）の位置から見たときの投影面（Ｐ）上の各点をＰｒ、Ｐｌとし、そのときのｘ座標をＸｒ、Ｘｌ、ｙ座標をＹとする、両眼の位置はｚ軸をはさんで左右各３０ｍｍの等間隔、ｘ軸に水平とする。また明視距離（ｘｙｚ座標の原点０から視点の中心Ｅまでの垂直距離０Ｅ）は３００ｍｍとする。Ｐｒ、Ｐｌ、Ｙの座標をそれぞれ以下のように求める。
【００４２】
【数８】
2004362087.tif 000010
【００４３】
【数９】
2004362087.tif 000011
【００４４】
【数１０】
2004362087.tif 000012
【００４５】
４．結果と考察
図１０に、筑波山周辺の地形解析図（東西約２２．４ｋｍ、南北約１８．５ｋｍ）を示す。図中、ａ：標高図、ｂ：陰影図、ｃ：傾斜方位図、ｄ：傾斜量図、ｅ：斜面ベクトル図、が示される。
また、図１１に、筑波山周辺の立体視図（東西約２２．４ｋｍ、南北約１８．５ｋｍ）を示す。
傾斜方位図は、傾斜方向（θ）に円形カラーチャートのＳ＝１．０のときの色を割り当てて表示してある。傾斜量図は、傾斜角（φ）が０°を白、９０°を黒とする２５５〜０の濃淡値に変換して表示してある。これらの図の斜面ベクトル図は、図１に示すＩ０＝１．０のときの円形カラーチャートで表現してある。斜面ベクトルの立体視対（図１１）を立体視すると、傾斜方向、傾斜量、標高の空間的変化を１枚の地図で容易に認識することができる。
【００４６】
まず、筑波山周辺の傾斜方位図および傾斜量図と斜面ベクトル図とを比較する。傾斜方位図には、傾斜方向のみが表示されている。傾斜量図には、地形の傾斜量のみが表示されている。標高図や陰影図を参考に傾斜方位図を見ると、斜面の連続性の抽出に有効であり、谷線・尾根線を含む線形構造がよく分かる。しかし、平野と山地の区分が不明瞭で、また斜面の凹凸も不明瞭である。傾斜量図は平野と山地が明瞭に区分されている。また、谷線・尾根線を含む線形構造や段丘崖も明瞭に認識できる。ただ、山地では、地形の起伏が不明瞭である。対して、斜面ベクトル図を見ると、谷線・尾根線を含む線形構造が明瞭で、山地、平野、段丘崖の区分も可能であり、地形の違いが明瞭に識別できる。これらのことから、斜面ベクトル図は、傾斜方位図と傾斜量図の長所を合わせもつ上に、地形の起伏も表現していることがわかる。
【００４７】
次に、標高図および陰影図を斜面ベクトル図と比較する。標高図は、地形の外観を見るためには一般的なものであるが、標高図から特徴的な地形を判読するためにはある程度の熟練を要する。また、標高図から特徴的地形を定量的に抽出することはさらに難しい。陰影図は、影を付けることで地形の凹凸が立体的に見えるという点で、線形構造の抽出や地形の外観把握に優れている。しかし、光の照射方向に平行する線形構造が捉えられないことや、陰になる部分の地形が不明瞭であるという欠点がある。このため線形構造や地形の変化を方向に関係なく抽出するためには、複数の方向からの光源を想定した多くの画像を見なければならない。その点、斜面ベクトル図は１枚の図で線形構造などの特徴的地形や地形の変化を表現している上に、定量的に扱うことも可能である。また、斜面ベクトル図はカラー図であるため、人間の直感的な認識を助け、見落としが少なくなることが期待できる。
【００４８】
斜面ベクトル図は、単独で傾斜方向、傾斜量、地形の変化などを表現することが可能で、線形構造などの特徴的地形の抽出に適していることがわかった。しかし、正確に地形の変化を認識するためには、円形カラーチャートの色相と傾斜方向の対応付けを参照する必要がある。この煩雑さを防ぐために、斜面ベクトルに対応する色がつけられた地形を想定して、これを立体視することにした。これにより、地形の変化がより容易に認識できる（図１１）。
【００４９】
斜面ベクトル図が線形構造の抽出に適していることを確かめるために、以下に例を示す。
図１２に、大分県万年山周辺の地形解析図（東西約２３．３ｋｍ、南北約１８．５ｋｍ）を示す。図中、ａ：陰影図（光源：北）ｂ：陰影図（光源：南東）ｃ：傾斜方位図ｄ：傾斜量図ｅ：斜面ベクトル図、がそれぞれ示される。
図１３に、大分県万年山周辺の斜面ベクトルの地形立体視図（東西約２３．３ｋｍ、南北約１８．５ｋｍ）を示す。
【００５０】
まず、図１２の斜面ベクトル図を見ると、万年山の西〜南西部（区域Ａ）と北東部（区域Ｂ）において、ほぼ東西に走る線構造（青）が見られる。これら東西に走る線構造は、南〜南西向きの幅の狭い急峻な斜面の連続であり、本地域に多数発達している東西に走る正断層（町田・太田・河名・森脇・長岡（２００１）日本の地形７九州・南西諸島）に対比される。また、南北に走る線構造（赤）が青で示す線構造と交差している。これらの赤で示された線構造は東向きの幅の狭い急斜面の連続である。このことは、図１３の立体視対を用いて立体視して見るとさらに明瞭に認識できる。これらの線構造は、立体視からわかる地形の特徴からも、断裂構造を捉えている可能性が高い。また、区域Ａでは、東西方向の線構造である青い線が南北性の赤の線構造を連続的に分断しているようにも見える。
【００５１】
次に、図１２の陰影図と斜面ベクトル図を比較する。陰影図（図１２のａとｂ）を見ると、光源が北側にある場合は、区域Ａ、区域Ｂに見られる東西方向の線形構造が鮮明に認識できるが、光源が南東側にある場合はかなり不鮮明である。対して、斜面ベクトル図ではさまざまな方向の線形構造が色相の違う線で全て同時に表現されている。さらに、斜面ベクトル図を、傾斜方位図および傾斜量図と比較すると、傾斜方位図は、傾斜方向の連続性の抽出には適しているが、傾斜角が急峻な部分のみを抽出することは不可能であり、抽出される断層の幅が実際より広くなる。一方、傾斜量図では、さまざまな方向性を持つ斜面が混在している場合、傾斜方向の識別ができず線形構造の連続性を追いにくい。
【００５２】
図１４に、傾斜方位図・傾斜量図・斜面ベクトル図における地形断面と表示色の関係図を示す。断層の抽出で斜面ベクトル図が傾斜方位図や傾斜量図より優れているのは、この図を見れば明らかである。例えば、両側が同程度の傾きをもつ谷の断面の場合、傾斜方位図（ａ）では、傾斜方向のみによって色分けされ、同じ傾斜方向を持つ斜面は同じ色相と彩度で表現されるため（この例では、図中左側が赤、右側が青）、抽出された線構造の幅が広くなる。傾斜量図（ｂ）でも谷線の両側の斜面が同じ色相と彩度で表現されるため（この例では、図中上側が淡く、下側が濃い。）、抽出された線構造の幅が広くなる。これに対し、斜面ベクトル図（傾斜ベクトル図）（ｃ）では、傾斜角が大きい部分は彩度が高く（この例では、図中上側が淡く、下側が濃い。）、傾斜方向が色相で区別されるため（この例では、図中左側が赤、右側が青）、推定される断層の幅がシャープになる。また、環状構造の抽出にも斜面ベクトル図は有効である。
【００５３】
図１５に、十和田湖周辺の傾斜量図（左）と傾斜ベクトル図（右）（東西約２１．２ｋｍ、南北約１８．５ｋｍ）を示す。
図１６に、十和田湖周辺の立体視対の図（東西約２１．２ｋｍ、南北約２７．８ｋｍ）を示す。
図１５の傾斜量図と斜面ベクトル図を比較すると、傾斜量図は環状構造を形成している急峻な傾斜面の連続をよく表現している。しかし、環状構造の内側が陥没しているか隆起しているかの判断は難しい。斜面ベクトル図でも、傾斜の急峻な部分は彩度が大きいため環状構造は抽出される。さらに、斜面ベクトル図では、傾斜の方向がわかるので、陥没構造と隆起構造の識別が可能となる。また、斜面ベクトル図を立体視することで陥没構造と隆起構造をより容易に認識できる。
【００５４】
５．付記
本発明の画像作成方法又は画像作成装置・システムは、その各手順をコンピュータに実行させるための画像作成プログラム、画像作成プログラムを記録したコンピュータ読み取り可能な記録媒体、画像作成プログラムを含みコンピュータの内部メモリにロード可能なプログラム製品、そのプログラムを含むサーバ等のコンピュータ、等により提供されることができる。
【００５５】
【発明の効果】
本発明によると、以上のように、谷線・尾根線を含む線形構造が明瞭で、山地、平野、段丘崖の区分も可能であり、地形の違いが明瞭に識別できる画像を作成することができる。また、本発明によると、作成された斜面ベクトル図により、傾斜方位図と傾斜量図の長所を合わせもち、地形の起伏も表現することができる。
【００５６】
また、本発明によると、斜面ベクトルの色を地図上に表示することにより、傾斜方向と傾斜量の空間分布を同時に表現することができる。さらに、このようにして作成した斜面ベクトル図を、地形が立体視できる立体視対として表現することにより、断層帯などの特徴的地形をより容易に抽出できる斜面ベクトル図を形成することができる。
【図面の簡単な説明】
【図１】斜面ベクトルのステレオ投影点とＨＳＩ円錐カラーモデルの対応関係の説明図。
【図２】Ｉ０が１．０、０．５、０．０のときのウルフ網の図。
【図３】本実施の形態に関するハードウェアの構成図。
【図４】標高モデルファイル５１の説明図。
【図５】斜面ベクトルファイル５２の説明図。
【図６】ＲＧＢ色画像ファイル５３の説明図。
【図７】立体視画像ファイル５４の説明図。
【図８】画像形成のためのフローチャート。
【図９】立体視対の作成方法の説明図。
【図１０】筑波山周辺の標高図、陰影図、傾斜方位図、傾斜量図、斜面ベクトル図。
【図１１】筑波山周辺の立体視対の図。
【図１２】大分県万年山周辺の陰影図、傾斜方位図、傾斜量図、斜面ベクトル図。
【図１３】大分県万年山周辺の立体視対の図。
【図１４】傾斜方位図・傾斜量図・斜面ベクトル図における地形断面と表示色の関係図。
【図１５】十和田湖周辺の傾斜量図と斜面ベクトル図。
【図１６】和田湖周辺の立体視対の図。
【図１７】ウルフ網の説明図。
【図１８】色相の変化幅の調整についての説明図。
【符号の説明】
１    主制御部（ＣＰＵ）
２    入力部
３    表示部
４    出力部
５    記憶部
６    入出力制御部
５１  標高モデルファイル
５２  斜面ベクトルファイル
５３  ＲＧＢ色画像ファイル
５４  立体視画像ファイル
５５  解析プログラム
５６  画像表示用プログラム
【図１】
2004362087.tif 000013
【図２】
2004362087.tif 000014
【図３】
2004362087.tif 000015
【図４】
2004362087.tif 000016
【図５】
2004362087.tif 000017
【図６】
2004362087.tif 000018
【図７】
2004362087.tif 000019
【図８】
2004362087.tif 000020
【図９】
2004362087.tif 000021
【図１０】
2004362087.tif 000022
【図１１】
2004362087.tif 000023
【図１２】
2004362087.tif 000024
【図１３】
2004362087.tif 000025
【図１４】
2004362087.tif 000026
【図１５】
2004362087.tif 000027
【図１６】
2004362087.tif 000028
【図１７】
2004362087.tif 000029
【図１８】
2004362087.tif 000030