地形起伏画像生成方法及び地形起伏画像生成装置
【課題】地形の標高を含む三次元地形データを基に、地形の比較的大きな起伏、尾根・谷等の比較的小さな凹凸、及び地形の傾斜度合などを判読し易い地形起伏画像を生成できる地形起伏画像生成方法及び地形起伏画像生成装置を提供する。
【解決手段】対象点の標高値Ｚに応じた色彩を付与して生成した標高段彩データＣ１と（Ｓ７０，Ｓ８０）、対象点の尾根谷度αに応じた色彩を付与して生成した尾根谷度段彩データＣ２（Ｓ９０，Ｓ１００）とを合成して第１合成データＧ１を取得する（Ｓ１１０）。対象点の傾斜度θを計算し、傾斜度θに応じた陰影Ｇθを取得する（Ｓ１２０〜１４０）。第１合成データＧ１に陰影Ｇθを乗算して第２合成データＧ２を生成する。上記各処理を全ての点（地点）に対して繰り返し行って段彩地形データＧを生成する。
【解決手段】対象点の標高値Ｚに応じた色彩を付与して生成した標高段彩データＣ１と（Ｓ７０，Ｓ８０）、対象点の尾根谷度αに応じた色彩を付与して生成した尾根谷度段彩データＣ２（Ｓ９０，Ｓ１００）とを合成して第１合成データＧ１を取得する（Ｓ１１０）。対象点の傾斜度θを計算し、傾斜度θに応じた陰影Ｇθを取得する（Ｓ１２０〜１４０）。第１合成データＧ１に陰影Ｇθを乗算して第２合成データＧ２を生成する。上記各処理を全ての点（地点）に対して繰り返し行って段彩地形データＧを生成する。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、例えば数値標高モデルデータ（ＤＥＭ：Digital Elevation Model）など標高を含む３次元地形データを基に、地形の比較的大きな起伏、尾根・谷などの比較的小さな凹凸（微地形）を判読し易く可視化された地形起伏画像を生成する地形起伏画像生成方法及び地形起伏画像生成装置に関する。
【背景技術】
【０００２】
近年、航空レーザ測量や人工衛星による地形測量技術の進歩により高密度な数値標高データが得られるようになってきた。数値標高データは任意の位置の座標を示し、各地点の座標はＸ、Ｙ、Ｚの座標情報から構成される。特に水平面位置（Ｘ及びＹ）が正規間隔（メッシュ）のデータをメッシュデータと呼ぶ。
【０００３】
従来より、地形の起伏を平面で表現してディスプレイや紙図に出力するための方法として、等高線図、段彩図、陰影図などが知られているが、いずれも地形表現に不得意な範疇がそれぞれあり、結果として正確な地形判読及び理解の妨げとなっていた。
【０００４】
また、例えば段彩陰影図と等高線図とを組み合わせたもの、標高段彩と傾斜区分とを組み合わせたものも知られている。
特許文献１には、ＤＥＭデータを基に、傾斜度が大きいほど濃い赤色となるように最大傾斜度に応じた赤色の階調表示をするとともに、凹凸度を表す浮沈度に応じて赤色の明度を階調表示することで、尾根が白く谷底が黒く表現された傾斜赤色化立体画像を生成する視覚化処理システムが開示されている。また、この特許文献１には、地上開度に応じたグレースケールと地下開度に応じたグレースケールとの合成画像と、傾斜度に応じた赤色画像とを合成することで、赤色化立体画像を生成する構成も開示されている。
【０００５】
また、特許文献２には、凹凸部を検出し、凹部に寒色系の色を付与した凹部の地形データと、凸部に暖色系の色を付与した凸部の地形データと、平坦部の地形データとを１つにまとめた陰影地形データを作成する。この陰影地形データに、同じ地形を表現する別のベース画像（例えば陰影画像、カラー段彩画像、等高線図画像、傾斜勾配画像など）を適当な混合比で混合して立体画像を作成する立体画像作成装置が開示されている。この立体画像では、尾根と谷が強調表示される。
【０００６】
さらに特許文献３には、数値標高データから色相などに割り当てた高度段彩図（カラー標高図）と、傾斜が大きいほど濃い濃度を割り当てたグレースケール傾斜図とを透過合成してカラー標高傾斜図を作成するカラー標高傾斜図作成方法が開示されている。高度段彩図は標高に応じて青〜赤の色相が割り当てられていた。また、特許文献３には、カラー標高図とラプラシアン図との合成も開示されている。なお、カラー標高図における標高に応じた段彩には青（低標高）〜赤（高標高）の色彩範囲が設定され、ラプラシアン図における傾斜変化率に応じた段彩には青（凹部）〜赤（凸部）の色彩範囲が設定されていた。
【０００７】
特許文献４では、標高データに基づいて、色相決定部が色相化標高データを生成し、輝度決定部が太陽光反射強度データを生成し、凹凸補正部が微細な凹凸情報の補正データを生成する。また、彩度決定部は、凹凸補正部の処理結果に基づいて彩度データを決定する。画像合成部は、色相化標高データ、太陽光反射強度データ、彩度データを合成し、対象となっている地表面の画像データを生成する地形図作成方法が開示されている。
【先行技術文献】
【特許文献】
【０００８】
【特許文献１】国際公開第２００４／０４２６７５号パンフレット（特許第３６７０２７４号公報）（例えば明細書段落［００５０］［００５１］、［００９５］−［０１０８］等、図１１、図１９、図２０等）
【特許文献２】特開２００６−７２８５７号公報（特許第４３７９２６４号公報）（例えば明細書段落［００１４］−［００１７］等、図１）
【特許文献３】特開２００７−４８１８５号公報（例えば明細書段落［００３６］［００４２］［００４７］［００５５］［００６４］）
【特許文献４】特開２００８−２４２２９８号公報（例えば明細書段落［００１４］−［００１８］［００２７］［００３３］等、図１、図５等）
【発明の概要】
【発明が解決しようとする課題】
【０００９】
しかし、特許文献１における赤色化立体画像では、谷部に陰影が施されることになり、谷部が暗く判読し難いという問題があった。
また、特許文献２における立体画像では、尾根と谷を判読し易いが、例えばベース画像としてカラー段彩画像と混合した場合、段彩画像の色彩と凹凸部の色彩とが同系色であると、尾根や谷が判読し難くなるという問題があった。また、凹部と凸部と平坦部とが３階調で階調表現（３階調表現）されるだけなので、凹部や凸部の度合までは判読し難かった。また、細い凹部と大きな窪みとを区別した表現ができるように空間周波数の設定を変えて別々に凹部の地形データを生成する必要があった。
【００１０】
さらに特許文献３では、カラー標高傾斜図では、標高と傾斜度を判読し易いものの、尾根や谷を判読し難かった。また、カラー標高図とラプラシアン図との合成図であれば、尾根や谷を判読し易いが、この合成図では傾斜度が判読し難いという問題があった。なお、カラー標高図における標高に応じた段彩と、ラプラシアン図における傾斜変化率に応じた段彩は、共に青（低標高、凹部）〜赤（高標高、凸部）の同系色の色彩範囲に設定されていたので、相対的に標高と凹凸度（尾根谷度）とを判別し難いという問題もある。
【００１１】
また、特許文献４によると、凹凸の標高を大きく変化させる補正を行って凹凸の色彩を強調する方法なので、凹凸を表現した色彩が実際の標高と違った色彩で表示されてしまうという問題があった。
【００１２】
本発明は、上記の課題に着目してなされたものであって、その目的は、地形の標高を含む三次元地形データを基に、地形の比較的大きな起伏、尾根・谷等の比較的小さな凹凸、及び地形の傾斜度合などを判読し易い地形起伏画像を生成できる地形起伏画像生成方法及び地形起伏画像生成装置を提供することにある。
【課題を解決するための手段】
【００１３】
上記目的を達成するために、本発明では、標高を含む三次元地形データを基に地形起伏画像を生成する地形起伏画像生成方法であって、入力した前記三次元地形データに基づき画素となる各地点に当該地点の標高に応じた色彩を付与する段彩処理を行って標高段彩データを生成する標高段彩処理手順と、前記三次元地形データに基づいて前記各地点における地形の凹凸度を演算する凹凸度演算手順と、前記各地点に当該地点毎の前記凹凸度に応じた色彩を付与する段彩処理を行って凹凸度段彩データを生成する凹凸度段彩処理手順と、前記三次元地形データに基づいて前記各地点における傾斜度を演算する傾斜度演算手順と、前記各地点に当該地点毎の前記傾斜度に応じた陰影を付与する陰影処理を行って陰影データを生成する陰影生成手順と、前記標高段彩データと前記凹凸度段彩データと前記陰影データとを合成して段彩地形データを生成する合成処理手順と、前記段彩地形データに基づく地形起伏画像を出力する出力手順と、を備えたことを要旨とする。
【００１４】
この発明によれば、標高に応じた段彩が施されるとともに、凹凸度に応じた段彩が施され、さらに傾斜度に応じた陰影が施される。よって、標高を含む三次元地形データを基に、地形の比較的大きな起伏、尾根・谷等の比較的小さな凹凸（微地形）、及び傾斜の度合などを判読し易い地形起伏画像を生成できる。
【００１５】
請求項２に記載の発明では、請求項１に記載の地形起伏画像生成方法において、地形起伏画像生成方法では、前記標高段彩処理手順では、前記標高の最小値と最大値との間に設定された第１の色彩範囲に従って標高に応じた色彩を付与して前記標高段彩データを生成し、前記凹凸度段彩処理手順では、前記凹凸度の最小値と最大値との間に設定された第２の色彩範囲に従って凹凸度に応じた色彩を付与して前記凹凸度段彩データを生成し、前記第１の色彩範囲と前記第２の色彩範囲は設定された色彩範囲が互いに異なるとともに、前記標高と前記凹凸度との各最小値に対応する色彩が互いに異なる寒色、最大値に対応する色彩が互いに異なる暖色に設定される。
【００１６】
この発明によれば、標高に応じた段彩と、凹凸度に応じた段彩は、それぞれ最小値が互いに異なる寒色、最大値が互いに異なる暖色に設定された互いに色彩範囲の異なる第１の色彩範囲と第２の色彩範囲に従ってそれぞれ施される。よって、地形の比較的大きな起伏と比較的小さな凹凸（微地形）とを、施された色の違いから判別し易い。
【００１７】
請求項３に記載の発明では、請求項１又は２に記載の地形起伏画像生成方法において、前記凹凸度演算手順では、前記三次元地形データにおける前記地点毎に、当該地点から複数の方向に地上角と地下角とをそれぞれ求め、前記地上角の平均である地上開度と、前記地下角の平均である地下開度とに基づき、前記凹凸度としての尾根谷度を演算する。
【００１８】
この発明によれば、着目する地点（対象点）から複数の方向に地上角と地下角とをそれぞれ求め、地上角の平均である地上開度と、地下角の平均である地下開度とに基づき演算した尾根谷度という地形量パラメータを凹凸度として用いて、尾根谷度に応じた段彩が施される。よって、尾根や谷などの微地形が一層適切に判読し易くなる。
【００１９】
請求項４に記載の発明では、請求項１乃至３のいずれか一項に記載の地形起伏画像生成方法において、前記傾斜度演算手順では、前記三次元地形データに基づいて前記各地点におけるｘ方向傾斜度とｙ方向傾斜度とを演算し、当該ｘ方向傾斜度とｙ方向傾斜度とを合成して前記傾斜度を演算する。
【００２０】
この構成によれば、ｘ方向傾斜度とｙ方向傾斜度とを合成して傾斜度を演算するので、傾斜面の向きによらず適切な傾斜度を取得できる。よって、傾斜度に応じた適切な陰影を付与できる。
【００２１】
請求項５に記載の発明では、請求項１乃至４のいずれか一項に記載の地形起伏画像生成方法において、前記合成処理手順では、前記三次元地形データにおける前記各地点に対応する三次元の空間座標点に、前記合成処理手順で取得した色情報を三成分の色座標で付加することにより、三次元の空間座標と三成分の色座標とを含む六次元以上の前記段彩地形データを生成し、前記出力手順では、前記段彩地形データの空間座標点を指定の視野方向にある出力面に写像した点を画素として当該画素に当該空間座標に対応する色座標を画素値として与えることで前記地形起伏画像を生成する。
【００２２】
この構成によれば、段彩地形データの空間座標点を指定の視野方向の出力面に写像した点を画素とし、当該画素に対応する色座標を画素値として与えることで地形起伏画像が生成される。よって、指定の視野方向に応じて写像方向を変えるだけで、指定の視野方向に応じた地形起伏画像を比較的簡単に生成できる。例えばレンダリング処理を不要にできる。
【００２３】
請求項６に記載の発明は、標高を含む三次元地形データを基に地形起伏画像を生成する地形起伏画像生成装置であって、入力した前記三次元地形データに基づき画素となる各地点に当該地点の標高に応じた色彩を付与する段彩処理を行って標高段彩データを生成する標高段彩処理手段と、前記三次元地形データに基づいて前記各地点における地形の凹凸度を演算する凹凸度演算手段と、前記各地点に当該地点毎の前記凹凸度に応じた色彩を付与する段彩処理を行って凹凸度段彩データを生成する凹凸度段彩処理手段と、前記三次元地形データに基づいて前記各地点における傾斜度を演算する傾斜度演算手段と、前記各地点に当該地点毎の前記傾斜度に応じた陰影を付与する陰影処理を行って陰影データを生成する陰影生成手段と、前記標高段彩データと前記凹凸度段彩データと前記陰影データとを合成して段彩地形データを生成する合成処理手段と、前記段彩地形データに基づく地形起伏画像を出力する出力手段と、を備えたことを要旨とする。この発明によれば、上記地形起伏画像生成方法に係る発明と同様の効果を得ることができる。
【発明の効果】
【００２４】
本発明によれば、標高を含む三次元地形データを基に、地形の比較的大きな起伏、尾根・谷等の比較的小さな凹凸、及び傾斜の度合などを判読し易い地形起伏画像を生成できる、という優れた効果が得られる。
【図面の簡単な説明】
【００２５】
【図１】地図生成装置の機能構成を示すブロック図。
【図２】航空レーザ測量の方法を説明する模式図。
【図３】対象点における地上角と地下角とを説明するグラフ。
【図４】（ａ）地上開度ΦLの説明図、（ｂ）は地下開度ΨLの説明図。
【図５】（ａ），（ｂ）凸部と凹部における地上開度ΦLと地下開度ΨLを示す説明図。
【図６】尾根谷度αと地上開度ΦL及び地下開度ΨLとの関係を示すグラフ。
【図７】傾斜度θの算出に用いる対象点の周囲近傍８点の標高を示す模式図。
【図８】（ａ）標高色彩変換テーブルＴ１、（ｂ）尾根谷度色彩変換テーブルＴ２の模式図。
【図９】地図生成装置の処理フロー図。
【図１０】標高段彩画像Ｃ１を示す図。
【図１１】尾根谷度段彩画像Ｃ２を示す図。
【図１２】第１合成画像Ｇ１を示す図。
【図１３】傾斜陰影画像Ｇθを示す図。
【図１４】第２合成画像Ｇ２（段彩地形画像Ｇ）を示す図。
【図１５】段彩地形画像ＧＰを示す図。
【図１６】平面図の地形起伏画像を示す図。
【図１７】鳥瞰図の地形起伏画像を示す図。
【発明を実施するための形態】
【００２６】
以下、本発明を具体化した一実施形態を図１〜図１７に基づいて説明する。
図１は、地形起伏画像生成装置の一例である地図生成装置の機能構成を示すブロック図である。図１に示す地図生成装置１００は、地形起伏画像生成用プログラムがパーソナルコンピュータＰＣの本体１０にインストールされて構築される。地図生成装置１００は、通信インターフェイス（例えばＵＳＢポート、ＬＡＮポート等）や光ディスク読取装置などからなる入力部１１と、キーボードやマウス等からなる入力操作部１２とに接続されている。地図生成装置１００には、入力部１１を通じて三次元地形データ（デジタル標高モデルデータ）ＤＥＭ（Ｘ，Ｙ，Ｚ）が入力される。地図生成装置１００は、三次元地形データＤＥＭ（Ｘ，Ｙ，Ｚ）（ＤＥＭ：Digital Elevation Model（数値標高モデル））を基に、標高Ｚ及び尾根谷度に応じて段彩され、かつ傾斜度に応じて陰影を施した地形起伏画像データＧ（Ｒ，Ｇ，Ｂ）を生成する。地図生成装置１００が生成した地形起伏画像データＧに基づく地形起伏画像（立体地形図）は、本体１０内の表示ドライバ１３を介してモニタ１４に表示されたり、本体１０内の印刷ドライバ１５を介してプリンタ１６により用紙に印刷されたりする。
【００２７】
三次元地形データＤＥＭ（Ｘ，Ｙ，Ｚ）は、例えば図２に示すような航空レーザ測量を利用して取得される。図２に示すように、レーザ測量器２０１を搭載した航空機２０２により地形測量対象エリアの上空を飛行しながらレーザを地上へ発信するレーザスキャンを行い、その反射パルスを計測して航空レーザ測量データを取得する。航空機２０２はＧＰＳ衛星２０３から受信するＧＰＳ信号に基づく現在の測量位置の緯度・経度と、そのときの反射レーザの計測値に基づく標高Ｚとを取得し、緯度・経度に標高が対応付けられた航空レーザ測量データを取得する。この航空レーザ測量データから必要に応じて建物や樹木、車両等の地形形状を表さない地物を分離処理し、さらにモデリングして発生させたＴＩＮ（三角形網：Triangle Irregular Network）等から、所定距離（例えば０．１ｍ〜５０ｍの範囲内の所定値）メッシュの交点の標高Ｚを読み取って、三次元地形データＤＥＭ（Ｘ，Ｙ，Ｚ）は生成される。もちろん、三次元地形データＤＥＭ（Ｘ,Ｙ,Ｚ）は、空中写真測量、衛星画像による測量、レーダによる測量等から取得することもできる。
【００２８】
以下、図１に示す地図生成装置１００の詳細を説明する。地図生成装置１００は、制御部２１、入力バッファ２２、標高段彩処理手段の一例である第１段彩演算部２３、尾根谷度演算手段の一例である尾根谷度演算部２４、尾根谷度段彩処理手段の一例である第２段彩演算部２５、傾斜度演算部２６、陰影演算部２７、合成処理部２８、ポリゴン処理部２９及び画像生成部３０を備えている。
【００２９】
制御部２１は、地図生成装置１００の全体制御を司り、各部２２〜３０に必要な指示を与える。入力バッファ２２には、入力部１１からの三次元地形データＤＥＭ（Ｘ,Ｙ,Ｚ）や、ユーザが入力操作部１２の操作により設定した各種の設定値などが格納される。各部２２〜２８は、三次元地形データＤＥＭ（Ｘ,Ｙ,Ｚ）を構成する１つの座標点（Ｘ，Ｙ，Ｚ）毎に、標高Ｚに応じた色彩を付与する段彩処理（第１の段彩処理）、尾根谷度αに応じた色彩を付与する段彩処理（第２の段彩処理）、傾斜度（傾斜角）θに応じた陰影を付与する陰影処理を施す。そして、合成処理部２８は、第１の段彩処理、第２の段彩処理及び陰影処理の各処理結果を必要に応じて重み付けをして合成することにより、三次元空間座標（Ｘ,Ｙ,Ｚ）と色値（Ｒ，Ｇ，Ｂ）との計６次元で表現される点（Ｘ,Ｙ,Ｚ，Ｒ，Ｇ，Ｂ）を生成する。さらに合成処理部２８は、三次元地形データＤＥＭ（Ｘ，Ｙ，Ｚ）の全ての点に対応する６次元の点（Ｘ,Ｙ,Ｚ，Ｒ，Ｇ，Ｂ）を生成すると、その６次元の点の集まりで表現された６次元のメッシュデータからなる段彩地形データＧ（Ｘ,Ｙ,Ｚ，Ｒ，Ｇ，Ｂ）を生成する。
【００３０】
第１段彩演算部２３は、標高値Ｚを色彩に変換するための標高色彩変換テーブルＴ１（図８（ａ）参照）を作成するための第１変換テーブル作成部３１を備え、第１変換テーブル作成部３１が作成した標高色彩変換テーブルＴ１を参照して、各点の標高値Ｚに応じた色彩を取得する。この色彩を表す点のデータは標高段彩データＣ１（Ｒ，Ｇ，Ｂ）として合成処理部２８に入力される。
【００３１】
尾根谷度演算部２４は、対象点（画素となる地点）の空間座標（Ｘ，Ｙ，Ｚ）から尾根谷度αを演算する。尾根谷度αの詳細な計算方法は後述する。
第２段彩演算部２５は、尾根谷度αを色彩に変換するための尾根谷度色彩変換テーブルＴ２（図８（ｂ）参照）を作成するための第２変換テーブル作成部３２を備え、第２変換テーブル作成部３２が作成した尾根谷度色彩変換テーブルＴ２を参照して、各点の標高値Ｚに応じた色彩を取得する。この色彩を表す点のデータは尾根谷度段彩データＣ２（Ｒ，Ｇ，Ｂ）として合成処理部２８に入力される。
【００３２】
傾斜度演算部２６は、対象点の空間座標（Ｘ，Ｙ，Ｚ）から傾斜度θ（傾斜角）を演算する。傾斜度θの詳細な計算方法は後述する。
陰影演算部２７は、傾斜度θを陰影（明暗）に変換するための傾斜度陰影変換テーブル（図示せず）又は傾斜度陰影変換用の計算式を用いて、傾斜度θに応じた陰影値を取得する。この陰影を表す点のデータは陰影データＧθとして合成処理部２８に入力される。
【００３３】
合成処理部２８は、標高段彩データＣ１（Ｒ，Ｇ，Ｂ）と尾根谷度段彩データＣ２（Ｒ，Ｇ，Ｂ）とを色合成して第１合成データＧ１を生成する色合成処理部３３と、第１合成データＧ１と陰影データＧθとを合成して第２合成データＧ２を生成する陰影合成処理部３４とを備えている。色合成処理部３３による色合成は、例えば重み付け合成（加重平均）により行われる。また、陰影合成処理部３４は、第１合成データＧ１（Ｒ１，Ｇ１，Ｂ１）に、陰影データの陰影値Ｇθ（＝０〜１．０）を乗算することにより、第２合成データＧ２（Ｒ２，Ｇ２，Ｂ２）＝（Ｒ１・Ｇθ，Ｇ１・Ｇθ，Ｂ１・Ｇθ）を求める。
【００３４】
第２合成データＧ２はバッファ３５に蓄積され、三次元地形データＤＥＭ（Ｘ，Ｙ，Ｚ）の全点について処理を終えると、合成処理部２８は、バッファ３５に格納された第２合成データＧ２（Ｒ２，Ｇ２，Ｂ２）の点の集まりとして表現される段彩地形データＧ（Ｘ，Ｙ，Ｚ，Ｒ，Ｇ，Ｂ）をポリゴン処理部２９へ出力する。
【００３５】
ポリゴン処理部２９は、段彩地形データＧ（Ｘ，Ｙ，Ｚ，Ｒ，Ｇ，Ｂ）中において水部ポリゴンデータＰＤに基づき特定される水部ポリゴン内の画素の色を青色に変更し、池、沼、河川、海等の水面が青色で彩色された段彩地形データＧＰ（Ｘ，Ｙ，Ｚ，Ｒ，Ｇ，Ｂ）を画像生成部３０へ出力する。ここで、水部ポリゴンデータＰＤとは、三次元地形データＤＥＭ（Ｘ，Ｙ，Ｚ）又はその他資料から作成され、作業者が水部であると認識した領域をポリゴン（多角形）として指定したデータである。
【００３６】
画像生成部３０は、段彩地形データＧＰ（Ｘ，Ｙ，Ｚ，Ｒ，Ｇ，Ｂ）を、ユーザが指定した視方向から見た地形起伏画像データＧ（Ｒ，Ｇ，Ｂ）を生成する。例えば視方向が真上からであるときの地形起伏画像（平面図）や、斜め上方からであるときの地形起伏画像（鳥瞰図）を生成する。そして、地形起伏画像は、ユーザが選択した出力方法に応じて、モニタ１４に表示されたり、プリンタ１６により用紙に印刷されたりする。
【００３７】
次に、図９に示す処理フローに従って地図生成装置１００による地形起伏画像生成方法の詳細を説明する。以下の説明の中で、尾根谷度α、傾斜度θの具体的な演算方法についても説明する。
【００３８】
本実施形態の処理フローは、テーブル生成処理と起伏画像生成処理とを含む。起伏画像生成処理は、大きく３つの処理に分かれる。つまり、起伏画像生成処理は、対象点の標高値Ｚに応じた色彩を求める第１段彩処理、対象点の尾根谷度αに応じた色彩を求める第２段彩処理、対象点における傾斜度θに応じた陰影を求める陰影処理とを含む。
【００３９】
まず各変換テーブルＴ１，Ｔ２を作成するテーブル生成処理について説明する。
ステップＳ１０では、三次元地形データＤＥＭ（Ｘ，Ｙ，Ｚ）の全点の標高Ｚを取得する。図８（ａ）に示す標高色彩変換テーブルＴ１の段彩設定から分かるように、低標高を緑（厳密には最小値が薄青緑色）、高標高を茶色（厳密には最大値が濃茶色）とすることを基本とする。
【００４０】
ステップＳ２０では、標高の最小値Ｚminと最大値Ｚmaxに基づき、標高色彩変換テーブルＴ１を作成する。緑（薄緑青）〜茶（濃茶）の色彩範囲で徐々に色が変わる複数の階調（本例では２２階調）に色分けされた色彩テーブル（図８（ａ）の右側部分）が予め用意されている。標高の最小値Ｚmin〜最大値Ｚmaxの標高範囲を、色彩テーブルと階調数と同数の複数段階（例えば２２段階）に分け、最小値Ｚmin〜最大値Ｚmaxの各標高に、緑〜茶の２２階調の各色彩を割り当てる。よって、対象エリア（表示エリア）の標高範囲が異なっても、最低標高を緑色、最高標高を茶色とする段彩処理が可能である。このため、比較的大きな起伏は、緑〜茶の段彩によって判読し易くなる。この段彩は、一般的な地図帳の表現と似て把握し易く、見た目も自然な印象を与えながら正確な地形起伏を捉えることが可能となる。
【００４１】
ステップＳ３０では、対象点（着目点）に対する設定距離内の点のデータ（Ｘ，Ｙ，Ｚ）を取得する。
ステップＳ４０では、対象点から設定距離内の点のデータ（Ｘ，Ｙ，Ｚ）を基に、尾根谷度αを計算する。ここで、図３〜図６を用いて、尾根谷度αの算出方法を説明する。尾根谷度αとは、地上開度と地下開度の平均角度より計算される凹凸の度合いを表す地形量である。
【００４２】
図３は、地上開度を計算する際に求める地上角と、地下開度を計算する際に求める地下角とを説明するグラフである。
地上開度とは、各メッシュが周辺に比べて地上に突き出た度合いを表現する地形量パラメータである。地上開度は、着目する地点（対象点）から設定距離Ｌの範囲内で見える空の広さを表し、一般に周囲から突出している地点ほどその値が大きくなる。地上開度は、山頂や尾根で大きな値、窪地や谷底では小さい値を示し、突出した山頂や尾根が強調される。
【００４３】
また、地下開度とは、各メッシュが周辺に比べて地下に突き出た度合いを表現する地形量パラメータである。地下開度は、地表面から地下を見渡すときに設定距離Ｌの範囲における地下の広さを表し、一般に地下にくい込んでいる地点ほどその値が大きくなる。地下開度は、窪地や谷底で大きな値、山頂や尾根では小さい値を示し、窪地や谷地が強調される。
【００４４】
地上開度、地下開度を求めるときには、まず図３に示す地上角Φ1と地下角Φ2をそれぞれ演算する。図３に示すように、対象点から８方向に格子点の標高Ｚを探索し（但し図３では１つの探索方向のみ示す）、各探索方向において設定距離Ｌ内で対象点Ｐjと標高点とを結ぶ線分と、ｚ軸（標高軸）において対象点Ｐjより高標高側に延びる軸部とのなす角度の最小値を、地上角Φ1として求める。この結果、８方向における８つの地上角Φ1が算出される。同様に、図３に示すように、対象点から８方向に格子点の標高Ｚを探索し、各探索方向において設定距離Ｌ内で対象点Ｐjと標高点とを結ぶ線分と、ｚ軸（標高軸）の対象点Ｐjより低標高側に延びる軸部とのなす角度の最小値を、地下角Φ2として求める。この結果、８方向における８つの地下角Φ2が算出される。なお、本実施形態では、０．５ｍ〜５．０ｍ程度の地形メッシュデータの場合は、設定距離Ｌは一例として２５ｍを標準とした。もちろん、設定距離Ｌは、メッシュサイズに応じてその数倍〜数１０倍の値を設定することが望ましく、５ｍ、１０ｍ、５０ｍ等の他の距離を設定してもよい。さらに設定距離Ｌはメッシュ数で設定してもよい。また、対象点からの探索方向は８方向に限らず、４方向又は１６方向や３２方向などさらに細分化してもよい。
【００４５】
次に８方向の各地上角Φ1の平均を演算して図４（ａ）に示す地上開度ΦLを求める。また、８方向の各地下角Φ2の平均を演算して図４（ｂ）に示す地下開度ΨLを求める。なお、図４（ａ）は谷や窪地等の凹部の地上開度ΦLを示した例、図４（ｂ）は尾根等の凸部の地下開度ΨLを示した例である。実際は、図５（ａ），（ｂ）に示すように、凹部と凸部のそれぞれに地上開度ΦLと地下開度ΨLが求められる。
【００４６】
図５（ａ）に示す凸部と図５（ｂ）に示す凹部に対して、地上開度ΦＬで規定される円錐面（地上円錐面）と、地下開度ΨＬで規定される円錐面（地下円錐面）ができる。図５（ａ），（ｂ）において、対象点Ｐjを通る水平面から下側（正）へ地上円錐面までの角度（以下、「地上開角度」と称す）（＝ΦL−９０）と、水平面から地下円錐面までの角度（以下、「地下開角度」と称す）（＝９０−ΨL）を規定する。
【００４７】
図６は、図５（ａ）のような凸部の例で、水平面から地上開角度（＝−９０＋ΦL）と、水平面から地下円錐面までの地下開角度（＝９０−ΨL）とを示した。本実施形態では、地上開角度（＝−９０＋ΦL）と、水平面から地下円錐面までの地下開角度（＝９０−ΨL）との平均（＝（（−９０＋ΦL）＋（９０−ΨL））／２）を、尾根谷度αとして定義する。この尾根谷度αによって、地上開度ΦLと地下開度ΨLの両方の地形量を加味した１つの指標で凹凸度を表現できる。よって、尾根谷度αは、地上開度ΦＬと地下開度ΨLを用いて次式で与えられる。
尾根谷度α＝（ΦL−ΨL）／２ …（１）
図５（ａ）に示す凸部では、相対的に地上開度ΦＬが大きな値、地下開度ΨＬが小さな値をとるので、尾根谷度αは正（α＞０）となる。一方、図５（ｂ）に示す凹部では、相対的に地上開度ΦＬが小さな値、地下開度ΨＬが大きな値をとるので、尾根谷度αは負（α＜０）となる。こうして尾根谷度αは、最小値αmin（＜０）〜最大値αmax（＞０）の範囲の値として求められる。
【００４８】
ステップＳ５０では、尾根谷度の最小値αminと最大値αmaxに基づき、尾根谷度色彩変換テーブルＴ２を作成する。青色〜赤色の範囲で徐々に色が変わる複数階調（本例では２２階調）に色分けされた色彩テーブル（図８（ｂ）の変換テーブルＴ２中の右側部分）が予め用意されている。尾根谷度の最小値αmin〜最大値αmaxの尾根谷度範囲を、色彩テーブルと同数の複数段階（例えば２２段階）に分け、最小値αmin〜最大値αmaxの各尾根谷度に、青〜赤の各色彩を割り当てる。よって、尾根谷度αに応じて、谷部（凹部）が青系、尾根部（凸部）が赤系になるように段彩される。この段彩により、比較的小さな凹凸部（微地形）でも判読し易くなる。こうして図８（ａ），（ｂ）に示す各変換テーブルＴ１，Ｔ２の作成により、テーブル生成処理は終了する。
【００４９】
図８（ａ）に示す各変換テーブルＴ１に設定された段彩のための第１の色彩範囲は「緑系〜茶系」であり、図８（ｂ）に示す各変換テーブルＴ２に設定された段彩のための第２の色彩範囲は「青系〜赤系」である。第１の色彩範囲と第２の色彩範囲は、共に最小値側が「寒色」、最大値側が「暖色」である点で共通しているが、寒色の具体的な色彩（緑系と青系）と、暖色の具体的な色彩（茶と赤）とが互いに異なっている。
【００５０】
本実施形態において、「寒色」とは、色成分の赤（Ｒ）・緑（Ｇ）・青（Ｂ）のうち緑・青のうち少なくとも一方の成分が赤の成分よりも大きな色彩を指す（Ｇ＞ＲとＢ＞Ｒのうち少なくとも一方が成立）。また、「暖色」とは、色成分の赤（Ｒ）・緑（Ｇ）・青（Ｂ）のうち赤の成分が緑・青の各成分よりも大きな色彩を指す（Ｒ＞ＧかつＲ＞Ｂ）。
【００５１】
標高の低い領域に緑系の色彩が付与され、標高の高い領域に茶系の色彩が付与される。また、谷部が青系の色彩に付与され、尾根部が赤系の色彩に付与される。標高Ｚと尾根谷度αで段彩表示のため付与される色彩が異なるので、比較的大きな起伏と、尾根部・谷部などの比較的小さな凹凸（微地形）とを、段彩による色彩の違いから判別し易い。また、標高の低い平地などを緑系、標高の高い山岳の山頂などを茶系に彩色し、一般の地図帳等に使用される馴染みの色彩を使用している。谷部が青系、尾根部が赤系に設定された第２の色彩範囲と、標高用に設定された第１の色彩範囲とを異ならせることで、尾根や谷も判読し易くしている。なお、第１の色彩範囲と前記第２の色彩範囲における最小値側の色彩は、一方が緑系、他方が青系であればよく、例えば標高用と尾根谷度用で逆の色彩を採用してもよい。
【００５２】
次に、三次元地形データＤＥＭ（Ｘ，Ｙ，Ｚ）の全対象点に、各変換テーブルＴ１，Ｔ２及び陰影テーブルを用いて求めた、標高Ｚに応じた緑〜茶の色彩、尾根谷度αに応じた青〜赤の色彩、傾斜度θに応じた黒〜白の陰影を合成して施し、段彩地形データＧ（Ｘ，Ｙ，Ｚ，Ｒ，Ｇ，Ｂ，Ｚ）を生成する起伏画像生成処理（Ｓ６０〜Ｓ１７０）について説明する。
【００５３】
ステップＳ６０では、対象点Ｐj（但しj＝1,2,3,…n（ｎはＤＥＭにおける表示エリアの全点数））を設定する。三次元地形データＤＥＭ（Ｘ，Ｙ，Ｚ）を構成する全ての点（画素とすべき点）を１点ずつ順番に対象点として設定する。以下の対象点の色生成処理（Ｓ７０〜Ｓ１５０）では、対象点に対して第１段彩処理、第２段彩処理、陰影処理、及び合成処理を行って色情報（Ｒ，Ｇ，Ｂ）を付与する。
【００５４】
まずステップＳ７０では、対象点の標高値Ｚを取得する。ステップＳ８０では、標高値Ｚに応じた色彩（標高段彩値）（Ｒz，Ｇz，Ｂz）を、標高色彩変換テーブルＴ１（図８（ａ））を参照して計算する。標高色彩変換テーブルＴ１（図８（ａ））の参照により複数階調（２２階調）のうちどの階調に属するかを判断し、その属する階調範囲における標高及び色彩のそれぞれの下限値（第１境界値）Ｚ１，（Ｒ１，Ｇ１，Ｂ１）とそれぞれの上限値（第２境界値）Ｚ２，（Ｒ２，Ｇ２，Ｂ２）とを取得する。そして、標高に応じた色彩がスムーズになるように、以下の計算式に従ってＲＧＢ成分ごとに加重平均を算出することで、標高Ｚに応じた緑〜茶の色彩を示す標高段彩データＣ１（Ｒz，Ｇz，Ｂz）を計算する。
Ｒz＝Ｒ２＋（Ｒ１−Ｒ２）・（Ｚ−Ｚ２）／（Ｚ１−Ｚ２） …（２）
Ｇz＝Ｇ２＋（Ｇ１−Ｇ２）・（Ｚ−Ｚ２）／（Ｚ１−Ｚ２） …（３）
Ｂz＝Ｂ２＋（Ｂ１−Ｂ２）・（Ｚ−Ｚ２）／（Ｚ１−Ｚ２） …（４）
この標高段彩データＣ１（Ｒz，Ｇz，Ｂz）を空間座標（Ｘ，Ｙ，Ｚ）に付与した全点を画素としてカラー表示したものが、図１０に示す標高段彩画像Ｃ１である。図１０の標高段彩画像Ｃ１に示すように、この表示エリアでは標高Ｚに応じて薄青緑〜茶色に段彩され、段彩された色の違いから比較的大きな起伏の判読が可能である。
【００５５】
ステップＳ９０では、対象点の尾根谷度αを取得する。ステップＳ１００では、尾根谷度αに応じた色彩（尾根谷度段彩値）（Ｒα，Ｇα，Ｂα）を、尾根谷度色彩変換テーブルＴ２（図８（ｂ））を参照して計算する。尾根谷度色彩変換テーブルＴ２（図８（ｂ））の参照により複数階調（２２階調）のうちどの階調に属するかを判断し、その属する階調範囲における尾根谷度α及び色彩のそれぞれの下限値（第１境界値）α１，（Ｒ１，Ｇ１，Ｂ１）とそれぞれの上限値（第２境界値）α２，（Ｒ２，Ｇ２，Ｂ２）とを取得する。そして、尾根谷度αに応じた色彩がスムーズになるように、以下の計算式に従ってＲＧＢ成分ごとに加重平均を算出することで、尾根谷度αに応じた青〜赤の色彩を示す尾根谷度段彩データＣ２（Ｒα，Ｇα，Ｂα）を計算する。
Ｒα＝Ｒ２＋（Ｒ１−Ｒ２）・（α−α２）／（α１−α２） …（５）
Ｇα＝Ｇ２＋（Ｇ１−Ｇ２）・（α−α２）／（α１−α２） …（６）
Ｂα＝Ｂ２＋（Ｂ１−Ｂ２）・（α−α２）／（α１−α２） …（７）
この尾根谷度段彩データＣ２（Ｒα，Ｇα，Ｂα）を空間座標（Ｘ，Ｙ，Ｚ）に付与した全点を画素としてカラー表示したものが、図１１に示す尾根谷度段彩画像Ｃ２である。図１１の尾根谷度段彩画像Ｃ２に示すように、この表示エリアでは尾根谷度αに応じて青〜赤のうちの青緑〜黄色に段彩され、段彩された色の違いから尾根や谷などの比較的小さな凹凸の判読が可能である。
【００５６】
ステップＳ１１０では、標高段彩データＣ１と尾根谷度段彩データＣ２との重み付け平均（加重平均）を計算し、第１合成データＧ１を生成する。第１合成データＧ１（Ｒzα，Ｇzα，Ｂzα）は、以下の計算式に従って算出される。
Ｒzα＝（ｗ１・Ｒz＋ｗ２・Ｒα）／（ｗ１＋ｗ２） …（８）
Ｇzα＝（ｗ１・Ｇz＋ｗ２・Ｇα）／（ｗ１＋ｗ２） …（９）
Ｂzα＝（ｗ１・Ｂz＋ｗ２・Ｂα）／（ｗ１＋ｗ２） …（１０）
ここで、ｗ１は標高段彩の重み係数、ｗ２は尾根谷度段彩の重み係数である。この第１合成データＧ１（Ｒzα，Ｇzα，Ｂzα）を空間座標（Ｘ，Ｙ，Ｚ）に付与した全点を画素としてカラー表示したものが、図１２に示す第１合成画像Ｇ１である。この第１合成画像Ｇ１は、標高段彩の重み係数ｗ１＝０．５、尾根谷度段彩の重み係数ｗ２＝０．５の例を示している。各変換テーブルＴ１，Ｔ２に設定された色彩範囲（緑系〜茶系）（青系〜赤系）の組合せに応じて、標高Ｚによる段彩と尾根谷度αによる段彩とを識別できるように、適宜な重み係数ｗ１，ｗ２を設定することができる。
【００５７】
この第１合成画像Ｇ１の表示エリアでは、標高Ｚに応じた緑〜茶の範囲の段彩と、尾根谷度αに応じた青〜赤の範囲の段彩とが色合成され、平地が薄緑、標高の高い部分が茶鼠でその中に茶色の尾根部が脈状に延び、標高の次に高い部分が黄茶でその中に濃黄色の尾根部が脈状に延び、さらに脈状に延びる谷部が青緑で表現されている。この第１合成画像Ｇ１では、段彩された色の違いから、比較的大きな起伏と比較的小さな凹凸の両方の判読が可能である。
【００５８】
ステップＳ１２０では、対象点の近傍８点のデータ（Ｘ，Ｙ，Ｚ）を取得する。ステップＳ１３０では、対象点の近傍８点のデータ（Ｘ，Ｙ，Ｚ）を用いて傾斜度θを計算する。図７は、対象点の周囲近傍に位置する８点の標高値を示す。まず対象点の周囲近傍８点の標高値Ｚ１〜Ｚ８を用いて、ｘ方向の傾斜度Ｓｘとｙ方向の傾斜度Ｓｙとを次式により計算する。
Ｓｘ＝（（Ｚ１+Ｚ４+Ｚ６）−（Ｚ３+Ｚ５+Ｚ８））／６／Ｍsize …（１１）
Ｓｙ＝（（Ｚ１+Ｚ２+Ｚ３）−（Ｚ６+Ｚ７+Ｚ８））／６／Ｍsize …（１２）
ここで、Ｍsizeはメッシュサイズ（１画素１辺のサイズ（長さ））である。次にｘ方向の傾斜度Ｓｘとｙ方向の傾斜度Ｓｙとを用いて、次式により傾斜度（傾斜角）θを計算する。
傾斜度θ＝arctan（√（Ｓx・Ｓx ＋Ｓy・Ｓy） …（１３）
このように本実施形態では、ｘ方向の傾斜度Ｓｘとｙ方向の傾斜度Ｓｙとを合成して傾斜度θを計算している。この傾斜度θの使用によって、対象点Ｐjの周囲８方位の傾斜角のうち最大傾斜角を使用して陰影を付与する方法に比べ、周囲８方位の傾斜角が適度にスムージングされ、傾斜度θに応じた自然な陰影が付与される。もちろん、「（Ｚ１+Ｚ４+Ｚ６）／３−（Ｚ３+Ｚ５+Ｚ８）／３」のようなｘ方向両側の各３点の標高平均の差をメッシュサイズＭsizeで除した値を傾斜度θに用いてもよい。なお、傾斜度θとして式 θ＝√（Ｓx・Ｓx ＋Ｓy・Ｓy）を使用することもできる。
【００５９】
ステップＳ１４０では、傾斜度θにより陰影Ｇθを計算する。最大傾斜度を「黒」、最小傾斜度を「白」とする予め用意された傾斜度陰影変換テーブル（グレースケールテーブル）を参照して、傾斜度θに応じた陰影Ｇθを求める。本例では、陰影Ｇθは、０〜１．０の範囲内の値として算出される。この陰影Ｇθを空間座標（Ｘ，Ｙ，Ｚ）に付与した全点を画素としてグレースケール表示したものが、図１３に示す傾斜陰影画像Ｇθである。この傾斜陰影画像Ｇθの表示エリアでは、傾斜度θに応じた陰影が表現され、傾斜が大きい領域が濃色（低明度）、傾斜が小さい領域が淡色（高明度）のグレースケールで表現されている。例えば脈状に延びる尾根部や谷部、平地などは淡灰色（明灰色）で表現され、尾根と谷の間の斜面は傾斜に応じた濃さの濃灰色（暗灰色）で表現される。
【００６０】
ステップＳ１５０では、第１合成データＧ１と陰影Ｇθとを乗算し、第２合成データＧ２を生成する。第２合成データＧ２（Ｒzαθ，Ｇzαθ，Ｂzαθ）は、以下の計算式に従って算出される。
Ｒzαθ＝Ｒzα・Ｇθ …（１４）
Ｇzαθ＝Ｇzα・Ｇθ …（１５）
Ｂzαθ＝Ｂzα・Ｇθ …（１６）
この第２合成データＧ２（Ｒzαθ，Ｇzαθ，Ｂzαθ）を空間座標（Ｘ，Ｙ，Ｚ）に付与した全点を画素としてカラー表示したものが、図１４に示す第２合成画像Ｇ２である。この第２合成画像Ｇ２の表示エリアでは、標高Ｚに応じた緑〜茶の範囲の段彩と、尾根谷度αに応じた青〜赤の範囲の段彩と、傾斜度θに応じた陰影とが色合成され、平地が薄緑、標高の高い部分が焦茶でその中に明るい茶色で尾根部が脈状に延び、標高の次に高い部分の中に黄茶の尾根部が脈状に延び、さらに脈状に延びる谷部が青緑で表現されている。そして、尾根部の両側では傾斜度の大きな斜面が濃灰色の混色、傾斜度の小さな斜面が明灰色の混色で表現される。よって、この第２合成画像Ｇ２では、段彩された色の違いから、比較的大きな起伏と比較的小さな凹凸の両方の判読が可能であるうえ、陰影から傾斜の大きさも判読可能である。なお、第２合成データＧ２（Ｒzαθ，Ｇzαθ，Ｂzαθ）は、空間座標（Ｘ，Ｙ，Ｚ）と共に、空間座標と色値とを合わせた６次元のデータ（Ｘ，Ｙ，Ｚ，Ｒzαθ，Ｇzαθ，Ｂzαθ）として、バッファ３５に格納される。
【００６１】
こうして１つの対象点についてＳ７０〜Ｓ１５０の処理を終え、対象点に２種類の段彩及び陰影が合成されて付与されると、ステップＳ１６０において、全点について処理を終了したか否かを判断する。全点の処理終了前であって処理すべき座標点が残っていれば、ステップＳ１７０において次の点の番号を１つ繰り上げて（ｊ＝ｊ＋１）、ステップＳ６０において次の対象点Ｐｊを設定する。以下、同様に、各対象点についてＳ７０〜Ｓ１５０の処理を繰り返し、全ての座標点について処理を終えると、Ｓ１６０において全点の処理を終了したと判断される。このときバッファ３５には、全点の第２合成データＧ２（Ｘ，Ｙ，Ｚ，Ｒzαθ，Ｇzαθ，Ｂzαθ）が格納される。
【００６２】
ステップＳ１８０では、段彩地形データＧ（Ｘ，Ｙ，Ｚ，Ｒ，Ｇ，Ｂ）を生成する。すなわち、バッファ３５に格納された第２合成データＧ２（Ｘ，Ｙ，Ｚ，Ｒzαθ，Ｇzαθ，Ｂzαθ）の全点分の集まりが段彩地形データＧ（Ｘ，Ｙ，Ｚ，Ｒ，Ｇ，Ｂ）として生成される。この段彩地形データＧ（Ｘ，Ｙ，Ｚ，Ｒ，Ｇ，Ｂ）は空間座標値と色値とを成分に含む６次元の数値データであり、各点は空間座標点であると共に画素として機能する。この段彩地形データＧ（Ｘ，Ｙ，Ｚ，Ｒ，Ｇ，Ｂ）の各点を画素としてカラー表現すると、段彩地形データＧが第２合成データＧ２の集まりであることから、図１４と同じ画像となる。
【００６３】
ステップＳ１９０では、水部ポリゴンデータＰＤ（別途Ｘ，Ｙ座標列で表される閉じた多角形面）を取得し、段彩地形データＧ（Ｘ，Ｙ，Ｚ，Ｒ，Ｇ，Ｂ）において水部ポリゴンデータＰＤで規定されるポリゴン内の画素を色変更して段彩地形データＧＰ（Ｘ，Ｙ，Ｚ，Ｒ，Ｇ，Ｂ）を生成する。このポリゴン処理はポリゴン処理部２９が行う。このポリゴン処理の結果、池、沼、河川、海等の水面領域に設定された水部ポリゴン内に存在する全ての点（画素）の色値（Ｒ，Ｇ，Ｂ）が、例えば青色の値に変更される。この段彩地形データＧＰ（Ｘ，Ｙ，Ｚ，Ｒ，Ｇ，Ｂ）の各点を画素としてカラー表現したものが、図１５に示す段彩地形画像ＧＰである。
【００６４】
ステップＳ２００では、段彩地形データＧＰ（Ｘ，Ｙ，Ｚ，Ｒ，Ｇ，Ｂ）を基に、指定の方向視からの地形起伏画像データＧ（Ｒ，Ｇ，Ｂ）を生成する。この地形起伏画像データＧ（Ｒ，Ｇ，Ｂ）の生成は、画像生成部３０が行う。例えば指定方向視が真上からである平面視の場合は、段彩地形データＧＰ（Ｘ，Ｙ，Ｚ，Ｒ，Ｇ，Ｂ）をｚ方向から見たときの出力面（表示面）に指定エリア内の各点（画素）を写像し、出力面上の各写像点に対応する画素値（Ｒ，Ｇ，Ｂ）を与えることで、平面視の地形起伏画像データＧ（Ｒ，Ｇ，Ｂ）を生成する。この結果、生成された地形起伏画像データＧ（Ｒ，Ｇ，Ｂ）に基づき、図１６に示す平面視（ｚ方向視）の地形起伏画像が出力（表示又は印刷）される。また、指定方向視が斜め方向である場合は、段彩地形データＧＰ（Ｘ，Ｙ，Ｚ，Ｒ，Ｇ，Ｂ）をｚ方向に対して斜めから見たときの出力面（表示面）に指定エリア内の各点（画素）を写像し、出力面上の各写像点（画素）に、対応する画素値（Ｒ，Ｇ，Ｂ）を与えることで、斜めから見た地形起伏画像データＧ（Ｒ，Ｇ，Ｂ）を生成する。この結果、生成された地形起伏画像データＧ（Ｒ，Ｇ，Ｂ）に基づき、図１７に示す地形起伏画像（鳥瞰図）が出力（表示又は印刷）される。
【００６５】
従来、三次元地形データＤＥＭ（Ｘ，Ｙ，Ｚ）を基に鳥瞰図を生成する場合、まず地形表面を生成してからレンダリング処理で地形表面に色彩を付与する方法がとられていたので、処理負担の高いレンダリング処理により表示速度が遅くなってしまう。これに対し本実施形態では、三次元座標（Ｘ，Ｙ，Ｚ）に色値（Ｒ，Ｇ，Ｂ）を付加して生成される段彩地形データＧＰには三次元情報が保たれているので、段彩地形データＧＰを出力面（表示面）へ写像する際の向き（視方向）を変えるだけで、鳥瞰図などの画像データを比較的高速に生成できる。この結果、鳥瞰図のような三次元の画像も比較的高速に表示できる。
【００６６】
なお、図９において、Ｓ１０，Ｓ２０，Ｓ７０，Ｓ８０の処理が標高段彩処理手順に相当し、Ｓ３０，Ｓ４０の処理が凹凸度演算手順に相当し、さらにＳ５０，Ｓ９０，Ｓ１００の処理が凹凸度段彩処理手順に相当する。また、Ｓ１２０，Ｓ１３０の処理が傾斜度演算手順に相当し、Ｓ１４０の処理が陰影生成手順に相当し、さらにＳ１１０，Ｓ１５０，Ｓ１８０の処理が合成処理手順に相当する。また、Ｓ２００の処理が出力手順に相当する。
【００６７】
以上詳述したようにこの実施形態によれば、以下の効果が得られる。
（１）標高段彩と尾根谷度段彩（凹凸度段彩）との合成により、標高段彩により比較的大きな起伏を判読し易くなり、尾根谷度段彩により尾根・谷などの比較的小さな凹凸（微地形）を判読し易くなる。
【００６８】
（２）微地形を強調する陰影を傾斜度θから作成するので、斜めの入射光を想定して陰影を施した際に問題となる陰影の光入射方向依存性がないうえ、特許文献１のように凹凸度（浮沈度）に応じた陰影に起因し赤色立体地図において谷部が暗くなる問題も回避できる。尾根谷度段彩画像に傾斜度θによる陰影を重ね合せることで、明暗だけでは区別し難い尾根と谷を一層明確に判読できる。また、傾斜度θに応じて、緩傾斜部分の明度を高く、急傾斜部分の明度を低くすることで立体感を与えることができる。
【００６９】
（３）標高段彩は、標高値Ｚに応じた第１色彩範囲「緑（最低標高Ｚmin）〜茶（最高標高Ｚmax）」を着色することで強調表現し、尾根谷度段彩は、尾根谷度αに応じた第２色彩範囲「青（最小値αmin）〜赤（最大値αmax）」を着色することで強調表現される。第１色彩範囲と第２色彩範囲が異なる色彩範囲（つまりパラメータ（標高Ｚ、尾根谷度α）の最小値の色彩と最大値の色彩が互いに異なるよう）に設定するとともに、共にパラメータの最小値が互いに異なる寒色、最大値が互いに異なる暖色になるように設定した。よって、標高段彩と尾根谷度段彩を合成表示しても、比較的大きな起伏の判読とその中にある比較的小さな凹凸（微地形）の判読が一層し易い。
【００７０】
（４）特許文献１における赤色化立体画像のように浮沈度に応じて赤色の明度で階調表示する方法であると、尾根などの凸部が膨張、谷などの凹部が収縮したイメージ（谷が細く感じる）、谷内が暗く見え難くなる。しかし、本実施形態によれば、尾根や谷などの凹凸が尾根谷度αに応じて青〜赤の色彩範囲の段彩で強調表示されるので、尾根などの凸部が膨張、谷などの凹部が収縮したイメージになることを回避できる。
【００７１】
（５）段彩地形データＧＰは、３次元座標（Ｘ，Ｙ，Ｚ）の付加情報として色情報（Ｒ，Ｇ，Ｂ）を付加することにより生成され、三次元座標（Ｘ，Ｙ，Ｚ）を保持している。このため、空間的に三次元の段彩地形データＧＰを斜めから見ることでより、立体感の効いた地形起伏画像を作成できる。通常、三次元データから面（surface）を発生させてからレンダリング処理を行って地形起伏画像を表現することが一般的であり、表示速度が遅くなってしまう。しかし、本実施形態によれば、三次元地形データにおける着色された点を表示する処理で済むので、地形起伏画像を高速に生成できる。
【００７２】
（６）従来の地形画像は単色や赤〜青の段彩が多いが、本実施形態は、緑〜茶色と一般的な地図帳の表現に近い馴染み易い色調で地形起伏画像を表現できる。また、大きな起伏は緑〜茶色の段彩、凹凸などの微地形は青〜赤の段彩で強調表現されるので、その合成色上に施された注記（赤、青、黒色）が背景に溶け込まず判読し易い効果も得られる。また、別途作成された水部ポリゴン（または水部属性）を使用して水部の色付けを行うことにより、さらに地図帳に近いより馴染み易い地形起伏画像を提供できる。
【００７３】
（７）凹凸度演算手順（Ｓ３０，Ｓ４０）では、三次元地形データＤＥＭにおける対象点毎に、対象点から複数の方向（例えば８方向）に地上角Φ1と地下角Φ2とをそれぞれ求め、地上角Φ1の平均である地上開度ΦLと、地下角Φ2の平均である地下開度ΨLとに基づき、式 α＝（ΦL−ΨL）／２により、凹凸度の一例として尾根谷度αを計算する。尾根などの凸部を表現できる地上開度ΦLと、谷などの凹部を表現できる地下開度ΨLとから、凹凸度を表現できる１つの地形量パラメータとして尾根谷度αを求め、尾根谷度αに応じて変換テーブルＴ２から決まる色彩を施すので、微地形を適切に段彩表現できる。
【００７４】
（８）傾斜度演算手順（Ｓ１２０，Ｓ１３０）では、対象点の周囲近傍８点の標高Ｚを用いて（１１）式及び（１２）式によりそれぞれ計算したｘ方向傾斜度Ｓxとｙ方向傾斜度Ｓyとを（１３）式により合成して傾斜度θを求める。よって、対象点から８方向の傾斜角が適度にスムージングされた傾斜度θを採用できる。この結果、適度にスムージングされた傾斜度θに応じた自然な陰影を付与できる。
【００７５】
（９）標高段彩データＣ１と尾根谷度段彩データＣ２との重み付け平均（加重平均）を、（８）式〜（１０）式により計算し、第１合成データＧ１を生成する。よって、係数ｗ１，ｗ２による重み付けを調整することにより、比較的大きな起伏と比較的小さな凹凸とを判読し易い条件で、段彩地形データＧＰを生成できる。
【００７６】
前記実施形態は上記に限定されず、以下の態様に変更することもできる。
・変換テーブルＴ１，Ｔ２の色彩範囲は前記実施形態に限定されず、適宜な色彩範囲を設定することができる。例えば各変換テーブルＴ１，Ｔ２で使用する色彩範囲を入れ替えてもよい。さらに、変換テーブルＴ１，Ｔ２のうち少なくとも一方について、標高Ｚ又は尾根谷度αの最小値側が暖色、最大値側が寒色となるような色彩範囲を設定してもよい。また、標高用と尾根谷度用で同一の変換テーブルを用いて、使用する色彩範囲を同じとしてもよい。
【００７７】
・標高Ｚに応じた第１段彩処理、尾根谷度αに応じた第２段彩処理、傾斜度に応じた陰影処理及び合成処理を、対象点毎に全点終了まで繰り返す方法を採用したが、これに替えて、全点の処理を終えた後に次の処理へ進む方法も採用できる。
【００７８】
・傾斜度を求める計算式は適宜変更してもよい。最大傾斜角を求める計算式を用いてもよい。
・尾根谷度αを求める計算式は適宜変更してもよい。また、探索方向は８方向に限らず、４方向や１６方向、あるいは３２方向等でもよい。
【００７９】
・凹凸度は尾根谷度αに限定されない。凹凸度として、ラプラシアンを用いてもよい。さらに、曲線上のある点を通る全ての測地線（曲面上で２点間を結ぶ最短距離の曲線）の曲率の最大値と最小値の平均として定義される平均曲率を凹凸度として用いてもよい。
【００８０】
・三次元地形データは、ＴＩＮ（triangulated irregular network，不整三角形網）でもよい。さらには、オルソ画像、ＤＴＭ（Digital Terrain Model：数値地形モデル）、ＤＳＭ（Digital Surface Model：数値表層モデル）でもよい。
【００８１】
・段彩地形データは６次元データに限定されない。他の情報を含めた７次元以上のデータであってもよい。また、三次元空間座標（Ｘ，Ｙ，Ｚ）を無くした三次元の段彩地形データＧＰ（Ｒ，Ｇ，Ｂ）であってもよい。この場合、指定の視野方向毎に段彩地形データの生成処理を行うことになる。
【符号の説明】
【００８２】
１０…本体、１１…入力部、１２…入力操作部、１３…表示ドライバ、１４…モニタ、１５…印刷ドライバ、１６…プリンタ、２１…制御部、２２…入力バッファ、２３…標高段彩処理手段の一例である第１段彩演算部、２４…凹凸度演算手段の一例である尾根谷度演算部、２５…凹凸度段彩処理手段の一例である第２段彩演算部、２６…傾斜度演算手段の一例である傾斜度演算部、２７…陰影生成手段の一例である陰影演算部、２８…合成処理手段の一例である合成処理部、２９…ポリゴン処理部、３０…出力手段の一例である画像生成部、３１…第１変換テーブル作成部、３２…第２変換テーブル作成部、３３…色合成処理部、３４…陰影合成処理部、３５…バッファ、１００…地形起伏画像生成装置の一例である地図生成装置、ＰＣ…パーソナルコンピュータ、ＤＥＭ（Ｘ，Ｙ，Ｚ）…三次元地形データ、Ｚ…標高、Φ1…地上角、Φ2…地下角、ΦL…地上開度、ΨL…地下開度、α…凹凸度の一例である尾根谷度、Ｓx…ｘ方向傾斜度、Ｓy…ｙ方向傾斜度、θ…傾斜度、Ｔ１…高色彩変換テーブル、Ｔ２…尾根谷度色彩変換テーブル、Ｚmin…標高の最小値、Ｚmax…標高の最大値、αmin…尾根谷度の最小値、αmax…尾根谷度の最大値、Ｐj…対象点（着目の地点）Ｃ１…標高段彩データ、Ｃ２…凹凸度段彩データの一例である尾根谷度段彩データ、Ｇθ…陰影データ（陰影）、Ｇ１…第１合成データ、Ｇ２…第２合成データ、Ｇ（Ｘ，Ｙ，Ｚ，Ｒ，Ｇ，Ｂ）…段彩地形データ（水部ポリゴン色付与前）、ＧＰ（Ｘ，Ｙ，Ｚ，Ｒ，Ｇ，Ｂ）…段彩地形データ（水部ポリゴン色付与後）、Ｇ（Ｒ，Ｇ，Ｂ）…地形起伏画像データ。
【技術分野】
【０００１】
本発明は、例えば数値標高モデルデータ（ＤＥＭ：Digital Elevation Model）など標高を含む３次元地形データを基に、地形の比較的大きな起伏、尾根・谷などの比較的小さな凹凸（微地形）を判読し易く可視化された地形起伏画像を生成する地形起伏画像生成方法及び地形起伏画像生成装置に関する。
【背景技術】
【０００２】
近年、航空レーザ測量や人工衛星による地形測量技術の進歩により高密度な数値標高データが得られるようになってきた。数値標高データは任意の位置の座標を示し、各地点の座標はＸ、Ｙ、Ｚの座標情報から構成される。特に水平面位置（Ｘ及びＹ）が正規間隔（メッシュ）のデータをメッシュデータと呼ぶ。
【０００３】
従来より、地形の起伏を平面で表現してディスプレイや紙図に出力するための方法として、等高線図、段彩図、陰影図などが知られているが、いずれも地形表現に不得意な範疇がそれぞれあり、結果として正確な地形判読及び理解の妨げとなっていた。
【０００４】
また、例えば段彩陰影図と等高線図とを組み合わせたもの、標高段彩と傾斜区分とを組み合わせたものも知られている。
特許文献１には、ＤＥＭデータを基に、傾斜度が大きいほど濃い赤色となるように最大傾斜度に応じた赤色の階調表示をするとともに、凹凸度を表す浮沈度に応じて赤色の明度を階調表示することで、尾根が白く谷底が黒く表現された傾斜赤色化立体画像を生成する視覚化処理システムが開示されている。また、この特許文献１には、地上開度に応じたグレースケールと地下開度に応じたグレースケールとの合成画像と、傾斜度に応じた赤色画像とを合成することで、赤色化立体画像を生成する構成も開示されている。
【０００５】
また、特許文献２には、凹凸部を検出し、凹部に寒色系の色を付与した凹部の地形データと、凸部に暖色系の色を付与した凸部の地形データと、平坦部の地形データとを１つにまとめた陰影地形データを作成する。この陰影地形データに、同じ地形を表現する別のベース画像（例えば陰影画像、カラー段彩画像、等高線図画像、傾斜勾配画像など）を適当な混合比で混合して立体画像を作成する立体画像作成装置が開示されている。この立体画像では、尾根と谷が強調表示される。
【０００６】
さらに特許文献３には、数値標高データから色相などに割り当てた高度段彩図（カラー標高図）と、傾斜が大きいほど濃い濃度を割り当てたグレースケール傾斜図とを透過合成してカラー標高傾斜図を作成するカラー標高傾斜図作成方法が開示されている。高度段彩図は標高に応じて青〜赤の色相が割り当てられていた。また、特許文献３には、カラー標高図とラプラシアン図との合成も開示されている。なお、カラー標高図における標高に応じた段彩には青（低標高）〜赤（高標高）の色彩範囲が設定され、ラプラシアン図における傾斜変化率に応じた段彩には青（凹部）〜赤（凸部）の色彩範囲が設定されていた。
【０００７】
特許文献４では、標高データに基づいて、色相決定部が色相化標高データを生成し、輝度決定部が太陽光反射強度データを生成し、凹凸補正部が微細な凹凸情報の補正データを生成する。また、彩度決定部は、凹凸補正部の処理結果に基づいて彩度データを決定する。画像合成部は、色相化標高データ、太陽光反射強度データ、彩度データを合成し、対象となっている地表面の画像データを生成する地形図作成方法が開示されている。
【先行技術文献】
【特許文献】
【０００８】
【特許文献１】国際公開第２００４／０４２６７５号パンフレット（特許第３６７０２７４号公報）（例えば明細書段落［００５０］［００５１］、［００９５］−［０１０８］等、図１１、図１９、図２０等）
【特許文献２】特開２００６−７２８５７号公報（特許第４３７９２６４号公報）（例えば明細書段落［００１４］−［００１７］等、図１）
【特許文献３】特開２００７−４８１８５号公報（例えば明細書段落［００３６］［００４２］［００４７］［００５５］［００６４］）
【特許文献４】特開２００８−２４２２９８号公報（例えば明細書段落［００１４］−［００１８］［００２７］［００３３］等、図１、図５等）
【発明の概要】
【発明が解決しようとする課題】
【０００９】
しかし、特許文献１における赤色化立体画像では、谷部に陰影が施されることになり、谷部が暗く判読し難いという問題があった。
また、特許文献２における立体画像では、尾根と谷を判読し易いが、例えばベース画像としてカラー段彩画像と混合した場合、段彩画像の色彩と凹凸部の色彩とが同系色であると、尾根や谷が判読し難くなるという問題があった。また、凹部と凸部と平坦部とが３階調で階調表現（３階調表現）されるだけなので、凹部や凸部の度合までは判読し難かった。また、細い凹部と大きな窪みとを区別した表現ができるように空間周波数の設定を変えて別々に凹部の地形データを生成する必要があった。
【００１０】
さらに特許文献３では、カラー標高傾斜図では、標高と傾斜度を判読し易いものの、尾根や谷を判読し難かった。また、カラー標高図とラプラシアン図との合成図であれば、尾根や谷を判読し易いが、この合成図では傾斜度が判読し難いという問題があった。なお、カラー標高図における標高に応じた段彩と、ラプラシアン図における傾斜変化率に応じた段彩は、共に青（低標高、凹部）〜赤（高標高、凸部）の同系色の色彩範囲に設定されていたので、相対的に標高と凹凸度（尾根谷度）とを判別し難いという問題もある。
【００１１】
また、特許文献４によると、凹凸の標高を大きく変化させる補正を行って凹凸の色彩を強調する方法なので、凹凸を表現した色彩が実際の標高と違った色彩で表示されてしまうという問題があった。
【００１２】
本発明は、上記の課題に着目してなされたものであって、その目的は、地形の標高を含む三次元地形データを基に、地形の比較的大きな起伏、尾根・谷等の比較的小さな凹凸、及び地形の傾斜度合などを判読し易い地形起伏画像を生成できる地形起伏画像生成方法及び地形起伏画像生成装置を提供することにある。
【課題を解決するための手段】
【００１３】
上記目的を達成するために、本発明では、標高を含む三次元地形データを基に地形起伏画像を生成する地形起伏画像生成方法であって、入力した前記三次元地形データに基づき画素となる各地点に当該地点の標高に応じた色彩を付与する段彩処理を行って標高段彩データを生成する標高段彩処理手順と、前記三次元地形データに基づいて前記各地点における地形の凹凸度を演算する凹凸度演算手順と、前記各地点に当該地点毎の前記凹凸度に応じた色彩を付与する段彩処理を行って凹凸度段彩データを生成する凹凸度段彩処理手順と、前記三次元地形データに基づいて前記各地点における傾斜度を演算する傾斜度演算手順と、前記各地点に当該地点毎の前記傾斜度に応じた陰影を付与する陰影処理を行って陰影データを生成する陰影生成手順と、前記標高段彩データと前記凹凸度段彩データと前記陰影データとを合成して段彩地形データを生成する合成処理手順と、前記段彩地形データに基づく地形起伏画像を出力する出力手順と、を備えたことを要旨とする。
【００１４】
この発明によれば、標高に応じた段彩が施されるとともに、凹凸度に応じた段彩が施され、さらに傾斜度に応じた陰影が施される。よって、標高を含む三次元地形データを基に、地形の比較的大きな起伏、尾根・谷等の比較的小さな凹凸（微地形）、及び傾斜の度合などを判読し易い地形起伏画像を生成できる。
【００１５】
請求項２に記載の発明では、請求項１に記載の地形起伏画像生成方法において、地形起伏画像生成方法では、前記標高段彩処理手順では、前記標高の最小値と最大値との間に設定された第１の色彩範囲に従って標高に応じた色彩を付与して前記標高段彩データを生成し、前記凹凸度段彩処理手順では、前記凹凸度の最小値と最大値との間に設定された第２の色彩範囲に従って凹凸度に応じた色彩を付与して前記凹凸度段彩データを生成し、前記第１の色彩範囲と前記第２の色彩範囲は設定された色彩範囲が互いに異なるとともに、前記標高と前記凹凸度との各最小値に対応する色彩が互いに異なる寒色、最大値に対応する色彩が互いに異なる暖色に設定される。
【００１６】
この発明によれば、標高に応じた段彩と、凹凸度に応じた段彩は、それぞれ最小値が互いに異なる寒色、最大値が互いに異なる暖色に設定された互いに色彩範囲の異なる第１の色彩範囲と第２の色彩範囲に従ってそれぞれ施される。よって、地形の比較的大きな起伏と比較的小さな凹凸（微地形）とを、施された色の違いから判別し易い。
【００１７】
請求項３に記載の発明では、請求項１又は２に記載の地形起伏画像生成方法において、前記凹凸度演算手順では、前記三次元地形データにおける前記地点毎に、当該地点から複数の方向に地上角と地下角とをそれぞれ求め、前記地上角の平均である地上開度と、前記地下角の平均である地下開度とに基づき、前記凹凸度としての尾根谷度を演算する。
【００１８】
この発明によれば、着目する地点（対象点）から複数の方向に地上角と地下角とをそれぞれ求め、地上角の平均である地上開度と、地下角の平均である地下開度とに基づき演算した尾根谷度という地形量パラメータを凹凸度として用いて、尾根谷度に応じた段彩が施される。よって、尾根や谷などの微地形が一層適切に判読し易くなる。
【００１９】
請求項４に記載の発明では、請求項１乃至３のいずれか一項に記載の地形起伏画像生成方法において、前記傾斜度演算手順では、前記三次元地形データに基づいて前記各地点におけるｘ方向傾斜度とｙ方向傾斜度とを演算し、当該ｘ方向傾斜度とｙ方向傾斜度とを合成して前記傾斜度を演算する。
【００２０】
この構成によれば、ｘ方向傾斜度とｙ方向傾斜度とを合成して傾斜度を演算するので、傾斜面の向きによらず適切な傾斜度を取得できる。よって、傾斜度に応じた適切な陰影を付与できる。
【００２１】
請求項５に記載の発明では、請求項１乃至４のいずれか一項に記載の地形起伏画像生成方法において、前記合成処理手順では、前記三次元地形データにおける前記各地点に対応する三次元の空間座標点に、前記合成処理手順で取得した色情報を三成分の色座標で付加することにより、三次元の空間座標と三成分の色座標とを含む六次元以上の前記段彩地形データを生成し、前記出力手順では、前記段彩地形データの空間座標点を指定の視野方向にある出力面に写像した点を画素として当該画素に当該空間座標に対応する色座標を画素値として与えることで前記地形起伏画像を生成する。
【００２２】
この構成によれば、段彩地形データの空間座標点を指定の視野方向の出力面に写像した点を画素とし、当該画素に対応する色座標を画素値として与えることで地形起伏画像が生成される。よって、指定の視野方向に応じて写像方向を変えるだけで、指定の視野方向に応じた地形起伏画像を比較的簡単に生成できる。例えばレンダリング処理を不要にできる。
【００２３】
請求項６に記載の発明は、標高を含む三次元地形データを基に地形起伏画像を生成する地形起伏画像生成装置であって、入力した前記三次元地形データに基づき画素となる各地点に当該地点の標高に応じた色彩を付与する段彩処理を行って標高段彩データを生成する標高段彩処理手段と、前記三次元地形データに基づいて前記各地点における地形の凹凸度を演算する凹凸度演算手段と、前記各地点に当該地点毎の前記凹凸度に応じた色彩を付与する段彩処理を行って凹凸度段彩データを生成する凹凸度段彩処理手段と、前記三次元地形データに基づいて前記各地点における傾斜度を演算する傾斜度演算手段と、前記各地点に当該地点毎の前記傾斜度に応じた陰影を付与する陰影処理を行って陰影データを生成する陰影生成手段と、前記標高段彩データと前記凹凸度段彩データと前記陰影データとを合成して段彩地形データを生成する合成処理手段と、前記段彩地形データに基づく地形起伏画像を出力する出力手段と、を備えたことを要旨とする。この発明によれば、上記地形起伏画像生成方法に係る発明と同様の効果を得ることができる。
【発明の効果】
【００２４】
本発明によれば、標高を含む三次元地形データを基に、地形の比較的大きな起伏、尾根・谷等の比較的小さな凹凸、及び傾斜の度合などを判読し易い地形起伏画像を生成できる、という優れた効果が得られる。
【図面の簡単な説明】
【００２５】
【図１】地図生成装置の機能構成を示すブロック図。
【図２】航空レーザ測量の方法を説明する模式図。
【図３】対象点における地上角と地下角とを説明するグラフ。
【図４】（ａ）地上開度ΦLの説明図、（ｂ）は地下開度ΨLの説明図。
【図５】（ａ），（ｂ）凸部と凹部における地上開度ΦLと地下開度ΨLを示す説明図。
【図６】尾根谷度αと地上開度ΦL及び地下開度ΨLとの関係を示すグラフ。
【図７】傾斜度θの算出に用いる対象点の周囲近傍８点の標高を示す模式図。
【図８】（ａ）標高色彩変換テーブルＴ１、（ｂ）尾根谷度色彩変換テーブルＴ２の模式図。
【図９】地図生成装置の処理フロー図。
【図１０】標高段彩画像Ｃ１を示す図。
【図１１】尾根谷度段彩画像Ｃ２を示す図。
【図１２】第１合成画像Ｇ１を示す図。
【図１３】傾斜陰影画像Ｇθを示す図。
【図１４】第２合成画像Ｇ２（段彩地形画像Ｇ）を示す図。
【図１５】段彩地形画像ＧＰを示す図。
【図１６】平面図の地形起伏画像を示す図。
【図１７】鳥瞰図の地形起伏画像を示す図。
【発明を実施するための形態】
【００２６】
以下、本発明を具体化した一実施形態を図１〜図１７に基づいて説明する。
図１は、地形起伏画像生成装置の一例である地図生成装置の機能構成を示すブロック図である。図１に示す地図生成装置１００は、地形起伏画像生成用プログラムがパーソナルコンピュータＰＣの本体１０にインストールされて構築される。地図生成装置１００は、通信インターフェイス（例えばＵＳＢポート、ＬＡＮポート等）や光ディスク読取装置などからなる入力部１１と、キーボードやマウス等からなる入力操作部１２とに接続されている。地図生成装置１００には、入力部１１を通じて三次元地形データ（デジタル標高モデルデータ）ＤＥＭ（Ｘ，Ｙ，Ｚ）が入力される。地図生成装置１００は、三次元地形データＤＥＭ（Ｘ，Ｙ，Ｚ）（ＤＥＭ：Digital Elevation Model（数値標高モデル））を基に、標高Ｚ及び尾根谷度に応じて段彩され、かつ傾斜度に応じて陰影を施した地形起伏画像データＧ（Ｒ，Ｇ，Ｂ）を生成する。地図生成装置１００が生成した地形起伏画像データＧに基づく地形起伏画像（立体地形図）は、本体１０内の表示ドライバ１３を介してモニタ１４に表示されたり、本体１０内の印刷ドライバ１５を介してプリンタ１６により用紙に印刷されたりする。
【００２７】
三次元地形データＤＥＭ（Ｘ，Ｙ，Ｚ）は、例えば図２に示すような航空レーザ測量を利用して取得される。図２に示すように、レーザ測量器２０１を搭載した航空機２０２により地形測量対象エリアの上空を飛行しながらレーザを地上へ発信するレーザスキャンを行い、その反射パルスを計測して航空レーザ測量データを取得する。航空機２０２はＧＰＳ衛星２０３から受信するＧＰＳ信号に基づく現在の測量位置の緯度・経度と、そのときの反射レーザの計測値に基づく標高Ｚとを取得し、緯度・経度に標高が対応付けられた航空レーザ測量データを取得する。この航空レーザ測量データから必要に応じて建物や樹木、車両等の地形形状を表さない地物を分離処理し、さらにモデリングして発生させたＴＩＮ（三角形網：Triangle Irregular Network）等から、所定距離（例えば０．１ｍ〜５０ｍの範囲内の所定値）メッシュの交点の標高Ｚを読み取って、三次元地形データＤＥＭ（Ｘ，Ｙ，Ｚ）は生成される。もちろん、三次元地形データＤＥＭ（Ｘ,Ｙ,Ｚ）は、空中写真測量、衛星画像による測量、レーダによる測量等から取得することもできる。
【００２８】
以下、図１に示す地図生成装置１００の詳細を説明する。地図生成装置１００は、制御部２１、入力バッファ２２、標高段彩処理手段の一例である第１段彩演算部２３、尾根谷度演算手段の一例である尾根谷度演算部２４、尾根谷度段彩処理手段の一例である第２段彩演算部２５、傾斜度演算部２６、陰影演算部２７、合成処理部２８、ポリゴン処理部２９及び画像生成部３０を備えている。
【００２９】
制御部２１は、地図生成装置１００の全体制御を司り、各部２２〜３０に必要な指示を与える。入力バッファ２２には、入力部１１からの三次元地形データＤＥＭ（Ｘ,Ｙ,Ｚ）や、ユーザが入力操作部１２の操作により設定した各種の設定値などが格納される。各部２２〜２８は、三次元地形データＤＥＭ（Ｘ,Ｙ,Ｚ）を構成する１つの座標点（Ｘ，Ｙ，Ｚ）毎に、標高Ｚに応じた色彩を付与する段彩処理（第１の段彩処理）、尾根谷度αに応じた色彩を付与する段彩処理（第２の段彩処理）、傾斜度（傾斜角）θに応じた陰影を付与する陰影処理を施す。そして、合成処理部２８は、第１の段彩処理、第２の段彩処理及び陰影処理の各処理結果を必要に応じて重み付けをして合成することにより、三次元空間座標（Ｘ,Ｙ,Ｚ）と色値（Ｒ，Ｇ，Ｂ）との計６次元で表現される点（Ｘ,Ｙ,Ｚ，Ｒ，Ｇ，Ｂ）を生成する。さらに合成処理部２８は、三次元地形データＤＥＭ（Ｘ，Ｙ，Ｚ）の全ての点に対応する６次元の点（Ｘ,Ｙ,Ｚ，Ｒ，Ｇ，Ｂ）を生成すると、その６次元の点の集まりで表現された６次元のメッシュデータからなる段彩地形データＧ（Ｘ,Ｙ,Ｚ，Ｒ，Ｇ，Ｂ）を生成する。
【００３０】
第１段彩演算部２３は、標高値Ｚを色彩に変換するための標高色彩変換テーブルＴ１（図８（ａ）参照）を作成するための第１変換テーブル作成部３１を備え、第１変換テーブル作成部３１が作成した標高色彩変換テーブルＴ１を参照して、各点の標高値Ｚに応じた色彩を取得する。この色彩を表す点のデータは標高段彩データＣ１（Ｒ，Ｇ，Ｂ）として合成処理部２８に入力される。
【００３１】
尾根谷度演算部２４は、対象点（画素となる地点）の空間座標（Ｘ，Ｙ，Ｚ）から尾根谷度αを演算する。尾根谷度αの詳細な計算方法は後述する。
第２段彩演算部２５は、尾根谷度αを色彩に変換するための尾根谷度色彩変換テーブルＴ２（図８（ｂ）参照）を作成するための第２変換テーブル作成部３２を備え、第２変換テーブル作成部３２が作成した尾根谷度色彩変換テーブルＴ２を参照して、各点の標高値Ｚに応じた色彩を取得する。この色彩を表す点のデータは尾根谷度段彩データＣ２（Ｒ，Ｇ，Ｂ）として合成処理部２８に入力される。
【００３２】
傾斜度演算部２６は、対象点の空間座標（Ｘ，Ｙ，Ｚ）から傾斜度θ（傾斜角）を演算する。傾斜度θの詳細な計算方法は後述する。
陰影演算部２７は、傾斜度θを陰影（明暗）に変換するための傾斜度陰影変換テーブル（図示せず）又は傾斜度陰影変換用の計算式を用いて、傾斜度θに応じた陰影値を取得する。この陰影を表す点のデータは陰影データＧθとして合成処理部２８に入力される。
【００３３】
合成処理部２８は、標高段彩データＣ１（Ｒ，Ｇ，Ｂ）と尾根谷度段彩データＣ２（Ｒ，Ｇ，Ｂ）とを色合成して第１合成データＧ１を生成する色合成処理部３３と、第１合成データＧ１と陰影データＧθとを合成して第２合成データＧ２を生成する陰影合成処理部３４とを備えている。色合成処理部３３による色合成は、例えば重み付け合成（加重平均）により行われる。また、陰影合成処理部３４は、第１合成データＧ１（Ｒ１，Ｇ１，Ｂ１）に、陰影データの陰影値Ｇθ（＝０〜１．０）を乗算することにより、第２合成データＧ２（Ｒ２，Ｇ２，Ｂ２）＝（Ｒ１・Ｇθ，Ｇ１・Ｇθ，Ｂ１・Ｇθ）を求める。
【００３４】
第２合成データＧ２はバッファ３５に蓄積され、三次元地形データＤＥＭ（Ｘ，Ｙ，Ｚ）の全点について処理を終えると、合成処理部２８は、バッファ３５に格納された第２合成データＧ２（Ｒ２，Ｇ２，Ｂ２）の点の集まりとして表現される段彩地形データＧ（Ｘ，Ｙ，Ｚ，Ｒ，Ｇ，Ｂ）をポリゴン処理部２９へ出力する。
【００３５】
ポリゴン処理部２９は、段彩地形データＧ（Ｘ，Ｙ，Ｚ，Ｒ，Ｇ，Ｂ）中において水部ポリゴンデータＰＤに基づき特定される水部ポリゴン内の画素の色を青色に変更し、池、沼、河川、海等の水面が青色で彩色された段彩地形データＧＰ（Ｘ，Ｙ，Ｚ，Ｒ，Ｇ，Ｂ）を画像生成部３０へ出力する。ここで、水部ポリゴンデータＰＤとは、三次元地形データＤＥＭ（Ｘ，Ｙ，Ｚ）又はその他資料から作成され、作業者が水部であると認識した領域をポリゴン（多角形）として指定したデータである。
【００３６】
画像生成部３０は、段彩地形データＧＰ（Ｘ，Ｙ，Ｚ，Ｒ，Ｇ，Ｂ）を、ユーザが指定した視方向から見た地形起伏画像データＧ（Ｒ，Ｇ，Ｂ）を生成する。例えば視方向が真上からであるときの地形起伏画像（平面図）や、斜め上方からであるときの地形起伏画像（鳥瞰図）を生成する。そして、地形起伏画像は、ユーザが選択した出力方法に応じて、モニタ１４に表示されたり、プリンタ１６により用紙に印刷されたりする。
【００３７】
次に、図９に示す処理フローに従って地図生成装置１００による地形起伏画像生成方法の詳細を説明する。以下の説明の中で、尾根谷度α、傾斜度θの具体的な演算方法についても説明する。
【００３８】
本実施形態の処理フローは、テーブル生成処理と起伏画像生成処理とを含む。起伏画像生成処理は、大きく３つの処理に分かれる。つまり、起伏画像生成処理は、対象点の標高値Ｚに応じた色彩を求める第１段彩処理、対象点の尾根谷度αに応じた色彩を求める第２段彩処理、対象点における傾斜度θに応じた陰影を求める陰影処理とを含む。
【００３９】
まず各変換テーブルＴ１，Ｔ２を作成するテーブル生成処理について説明する。
ステップＳ１０では、三次元地形データＤＥＭ（Ｘ，Ｙ，Ｚ）の全点の標高Ｚを取得する。図８（ａ）に示す標高色彩変換テーブルＴ１の段彩設定から分かるように、低標高を緑（厳密には最小値が薄青緑色）、高標高を茶色（厳密には最大値が濃茶色）とすることを基本とする。
【００４０】
ステップＳ２０では、標高の最小値Ｚminと最大値Ｚmaxに基づき、標高色彩変換テーブルＴ１を作成する。緑（薄緑青）〜茶（濃茶）の色彩範囲で徐々に色が変わる複数の階調（本例では２２階調）に色分けされた色彩テーブル（図８（ａ）の右側部分）が予め用意されている。標高の最小値Ｚmin〜最大値Ｚmaxの標高範囲を、色彩テーブルと階調数と同数の複数段階（例えば２２段階）に分け、最小値Ｚmin〜最大値Ｚmaxの各標高に、緑〜茶の２２階調の各色彩を割り当てる。よって、対象エリア（表示エリア）の標高範囲が異なっても、最低標高を緑色、最高標高を茶色とする段彩処理が可能である。このため、比較的大きな起伏は、緑〜茶の段彩によって判読し易くなる。この段彩は、一般的な地図帳の表現と似て把握し易く、見た目も自然な印象を与えながら正確な地形起伏を捉えることが可能となる。
【００４１】
ステップＳ３０では、対象点（着目点）に対する設定距離内の点のデータ（Ｘ，Ｙ，Ｚ）を取得する。
ステップＳ４０では、対象点から設定距離内の点のデータ（Ｘ，Ｙ，Ｚ）を基に、尾根谷度αを計算する。ここで、図３〜図６を用いて、尾根谷度αの算出方法を説明する。尾根谷度αとは、地上開度と地下開度の平均角度より計算される凹凸の度合いを表す地形量である。
【００４２】
図３は、地上開度を計算する際に求める地上角と、地下開度を計算する際に求める地下角とを説明するグラフである。
地上開度とは、各メッシュが周辺に比べて地上に突き出た度合いを表現する地形量パラメータである。地上開度は、着目する地点（対象点）から設定距離Ｌの範囲内で見える空の広さを表し、一般に周囲から突出している地点ほどその値が大きくなる。地上開度は、山頂や尾根で大きな値、窪地や谷底では小さい値を示し、突出した山頂や尾根が強調される。
【００４３】
また、地下開度とは、各メッシュが周辺に比べて地下に突き出た度合いを表現する地形量パラメータである。地下開度は、地表面から地下を見渡すときに設定距離Ｌの範囲における地下の広さを表し、一般に地下にくい込んでいる地点ほどその値が大きくなる。地下開度は、窪地や谷底で大きな値、山頂や尾根では小さい値を示し、窪地や谷地が強調される。
【００４４】
地上開度、地下開度を求めるときには、まず図３に示す地上角Φ1と地下角Φ2をそれぞれ演算する。図３に示すように、対象点から８方向に格子点の標高Ｚを探索し（但し図３では１つの探索方向のみ示す）、各探索方向において設定距離Ｌ内で対象点Ｐjと標高点とを結ぶ線分と、ｚ軸（標高軸）において対象点Ｐjより高標高側に延びる軸部とのなす角度の最小値を、地上角Φ1として求める。この結果、８方向における８つの地上角Φ1が算出される。同様に、図３に示すように、対象点から８方向に格子点の標高Ｚを探索し、各探索方向において設定距離Ｌ内で対象点Ｐjと標高点とを結ぶ線分と、ｚ軸（標高軸）の対象点Ｐjより低標高側に延びる軸部とのなす角度の最小値を、地下角Φ2として求める。この結果、８方向における８つの地下角Φ2が算出される。なお、本実施形態では、０．５ｍ〜５．０ｍ程度の地形メッシュデータの場合は、設定距離Ｌは一例として２５ｍを標準とした。もちろん、設定距離Ｌは、メッシュサイズに応じてその数倍〜数１０倍の値を設定することが望ましく、５ｍ、１０ｍ、５０ｍ等の他の距離を設定してもよい。さらに設定距離Ｌはメッシュ数で設定してもよい。また、対象点からの探索方向は８方向に限らず、４方向又は１６方向や３２方向などさらに細分化してもよい。
【００４５】
次に８方向の各地上角Φ1の平均を演算して図４（ａ）に示す地上開度ΦLを求める。また、８方向の各地下角Φ2の平均を演算して図４（ｂ）に示す地下開度ΨLを求める。なお、図４（ａ）は谷や窪地等の凹部の地上開度ΦLを示した例、図４（ｂ）は尾根等の凸部の地下開度ΨLを示した例である。実際は、図５（ａ），（ｂ）に示すように、凹部と凸部のそれぞれに地上開度ΦLと地下開度ΨLが求められる。
【００４６】
図５（ａ）に示す凸部と図５（ｂ）に示す凹部に対して、地上開度ΦＬで規定される円錐面（地上円錐面）と、地下開度ΨＬで規定される円錐面（地下円錐面）ができる。図５（ａ），（ｂ）において、対象点Ｐjを通る水平面から下側（正）へ地上円錐面までの角度（以下、「地上開角度」と称す）（＝ΦL−９０）と、水平面から地下円錐面までの角度（以下、「地下開角度」と称す）（＝９０−ΨL）を規定する。
【００４７】
図６は、図５（ａ）のような凸部の例で、水平面から地上開角度（＝−９０＋ΦL）と、水平面から地下円錐面までの地下開角度（＝９０−ΨL）とを示した。本実施形態では、地上開角度（＝−９０＋ΦL）と、水平面から地下円錐面までの地下開角度（＝９０−ΨL）との平均（＝（（−９０＋ΦL）＋（９０−ΨL））／２）を、尾根谷度αとして定義する。この尾根谷度αによって、地上開度ΦLと地下開度ΨLの両方の地形量を加味した１つの指標で凹凸度を表現できる。よって、尾根谷度αは、地上開度ΦＬと地下開度ΨLを用いて次式で与えられる。
尾根谷度α＝（ΦL−ΨL）／２ …（１）
図５（ａ）に示す凸部では、相対的に地上開度ΦＬが大きな値、地下開度ΨＬが小さな値をとるので、尾根谷度αは正（α＞０）となる。一方、図５（ｂ）に示す凹部では、相対的に地上開度ΦＬが小さな値、地下開度ΨＬが大きな値をとるので、尾根谷度αは負（α＜０）となる。こうして尾根谷度αは、最小値αmin（＜０）〜最大値αmax（＞０）の範囲の値として求められる。
【００４８】
ステップＳ５０では、尾根谷度の最小値αminと最大値αmaxに基づき、尾根谷度色彩変換テーブルＴ２を作成する。青色〜赤色の範囲で徐々に色が変わる複数階調（本例では２２階調）に色分けされた色彩テーブル（図８（ｂ）の変換テーブルＴ２中の右側部分）が予め用意されている。尾根谷度の最小値αmin〜最大値αmaxの尾根谷度範囲を、色彩テーブルと同数の複数段階（例えば２２段階）に分け、最小値αmin〜最大値αmaxの各尾根谷度に、青〜赤の各色彩を割り当てる。よって、尾根谷度αに応じて、谷部（凹部）が青系、尾根部（凸部）が赤系になるように段彩される。この段彩により、比較的小さな凹凸部（微地形）でも判読し易くなる。こうして図８（ａ），（ｂ）に示す各変換テーブルＴ１，Ｔ２の作成により、テーブル生成処理は終了する。
【００４９】
図８（ａ）に示す各変換テーブルＴ１に設定された段彩のための第１の色彩範囲は「緑系〜茶系」であり、図８（ｂ）に示す各変換テーブルＴ２に設定された段彩のための第２の色彩範囲は「青系〜赤系」である。第１の色彩範囲と第２の色彩範囲は、共に最小値側が「寒色」、最大値側が「暖色」である点で共通しているが、寒色の具体的な色彩（緑系と青系）と、暖色の具体的な色彩（茶と赤）とが互いに異なっている。
【００５０】
本実施形態において、「寒色」とは、色成分の赤（Ｒ）・緑（Ｇ）・青（Ｂ）のうち緑・青のうち少なくとも一方の成分が赤の成分よりも大きな色彩を指す（Ｇ＞ＲとＢ＞Ｒのうち少なくとも一方が成立）。また、「暖色」とは、色成分の赤（Ｒ）・緑（Ｇ）・青（Ｂ）のうち赤の成分が緑・青の各成分よりも大きな色彩を指す（Ｒ＞ＧかつＲ＞Ｂ）。
【００５１】
標高の低い領域に緑系の色彩が付与され、標高の高い領域に茶系の色彩が付与される。また、谷部が青系の色彩に付与され、尾根部が赤系の色彩に付与される。標高Ｚと尾根谷度αで段彩表示のため付与される色彩が異なるので、比較的大きな起伏と、尾根部・谷部などの比較的小さな凹凸（微地形）とを、段彩による色彩の違いから判別し易い。また、標高の低い平地などを緑系、標高の高い山岳の山頂などを茶系に彩色し、一般の地図帳等に使用される馴染みの色彩を使用している。谷部が青系、尾根部が赤系に設定された第２の色彩範囲と、標高用に設定された第１の色彩範囲とを異ならせることで、尾根や谷も判読し易くしている。なお、第１の色彩範囲と前記第２の色彩範囲における最小値側の色彩は、一方が緑系、他方が青系であればよく、例えば標高用と尾根谷度用で逆の色彩を採用してもよい。
【００５２】
次に、三次元地形データＤＥＭ（Ｘ，Ｙ，Ｚ）の全対象点に、各変換テーブルＴ１，Ｔ２及び陰影テーブルを用いて求めた、標高Ｚに応じた緑〜茶の色彩、尾根谷度αに応じた青〜赤の色彩、傾斜度θに応じた黒〜白の陰影を合成して施し、段彩地形データＧ（Ｘ，Ｙ，Ｚ，Ｒ，Ｇ，Ｂ，Ｚ）を生成する起伏画像生成処理（Ｓ６０〜Ｓ１７０）について説明する。
【００５３】
ステップＳ６０では、対象点Ｐj（但しj＝1,2,3,…n（ｎはＤＥＭにおける表示エリアの全点数））を設定する。三次元地形データＤＥＭ（Ｘ，Ｙ，Ｚ）を構成する全ての点（画素とすべき点）を１点ずつ順番に対象点として設定する。以下の対象点の色生成処理（Ｓ７０〜Ｓ１５０）では、対象点に対して第１段彩処理、第２段彩処理、陰影処理、及び合成処理を行って色情報（Ｒ，Ｇ，Ｂ）を付与する。
【００５４】
まずステップＳ７０では、対象点の標高値Ｚを取得する。ステップＳ８０では、標高値Ｚに応じた色彩（標高段彩値）（Ｒz，Ｇz，Ｂz）を、標高色彩変換テーブルＴ１（図８（ａ））を参照して計算する。標高色彩変換テーブルＴ１（図８（ａ））の参照により複数階調（２２階調）のうちどの階調に属するかを判断し、その属する階調範囲における標高及び色彩のそれぞれの下限値（第１境界値）Ｚ１，（Ｒ１，Ｇ１，Ｂ１）とそれぞれの上限値（第２境界値）Ｚ２，（Ｒ２，Ｇ２，Ｂ２）とを取得する。そして、標高に応じた色彩がスムーズになるように、以下の計算式に従ってＲＧＢ成分ごとに加重平均を算出することで、標高Ｚに応じた緑〜茶の色彩を示す標高段彩データＣ１（Ｒz，Ｇz，Ｂz）を計算する。
Ｒz＝Ｒ２＋（Ｒ１−Ｒ２）・（Ｚ−Ｚ２）／（Ｚ１−Ｚ２） …（２）
Ｇz＝Ｇ２＋（Ｇ１−Ｇ２）・（Ｚ−Ｚ２）／（Ｚ１−Ｚ２） …（３）
Ｂz＝Ｂ２＋（Ｂ１−Ｂ２）・（Ｚ−Ｚ２）／（Ｚ１−Ｚ２） …（４）
この標高段彩データＣ１（Ｒz，Ｇz，Ｂz）を空間座標（Ｘ，Ｙ，Ｚ）に付与した全点を画素としてカラー表示したものが、図１０に示す標高段彩画像Ｃ１である。図１０の標高段彩画像Ｃ１に示すように、この表示エリアでは標高Ｚに応じて薄青緑〜茶色に段彩され、段彩された色の違いから比較的大きな起伏の判読が可能である。
【００５５】
ステップＳ９０では、対象点の尾根谷度αを取得する。ステップＳ１００では、尾根谷度αに応じた色彩（尾根谷度段彩値）（Ｒα，Ｇα，Ｂα）を、尾根谷度色彩変換テーブルＴ２（図８（ｂ））を参照して計算する。尾根谷度色彩変換テーブルＴ２（図８（ｂ））の参照により複数階調（２２階調）のうちどの階調に属するかを判断し、その属する階調範囲における尾根谷度α及び色彩のそれぞれの下限値（第１境界値）α１，（Ｒ１，Ｇ１，Ｂ１）とそれぞれの上限値（第２境界値）α２，（Ｒ２，Ｇ２，Ｂ２）とを取得する。そして、尾根谷度αに応じた色彩がスムーズになるように、以下の計算式に従ってＲＧＢ成分ごとに加重平均を算出することで、尾根谷度αに応じた青〜赤の色彩を示す尾根谷度段彩データＣ２（Ｒα，Ｇα，Ｂα）を計算する。
Ｒα＝Ｒ２＋（Ｒ１−Ｒ２）・（α−α２）／（α１−α２） …（５）
Ｇα＝Ｇ２＋（Ｇ１−Ｇ２）・（α−α２）／（α１−α２） …（６）
Ｂα＝Ｂ２＋（Ｂ１−Ｂ２）・（α−α２）／（α１−α２） …（７）
この尾根谷度段彩データＣ２（Ｒα，Ｇα，Ｂα）を空間座標（Ｘ，Ｙ，Ｚ）に付与した全点を画素としてカラー表示したものが、図１１に示す尾根谷度段彩画像Ｃ２である。図１１の尾根谷度段彩画像Ｃ２に示すように、この表示エリアでは尾根谷度αに応じて青〜赤のうちの青緑〜黄色に段彩され、段彩された色の違いから尾根や谷などの比較的小さな凹凸の判読が可能である。
【００５６】
ステップＳ１１０では、標高段彩データＣ１と尾根谷度段彩データＣ２との重み付け平均（加重平均）を計算し、第１合成データＧ１を生成する。第１合成データＧ１（Ｒzα，Ｇzα，Ｂzα）は、以下の計算式に従って算出される。
Ｒzα＝（ｗ１・Ｒz＋ｗ２・Ｒα）／（ｗ１＋ｗ２） …（８）
Ｇzα＝（ｗ１・Ｇz＋ｗ２・Ｇα）／（ｗ１＋ｗ２） …（９）
Ｂzα＝（ｗ１・Ｂz＋ｗ２・Ｂα）／（ｗ１＋ｗ２） …（１０）
ここで、ｗ１は標高段彩の重み係数、ｗ２は尾根谷度段彩の重み係数である。この第１合成データＧ１（Ｒzα，Ｇzα，Ｂzα）を空間座標（Ｘ，Ｙ，Ｚ）に付与した全点を画素としてカラー表示したものが、図１２に示す第１合成画像Ｇ１である。この第１合成画像Ｇ１は、標高段彩の重み係数ｗ１＝０．５、尾根谷度段彩の重み係数ｗ２＝０．５の例を示している。各変換テーブルＴ１，Ｔ２に設定された色彩範囲（緑系〜茶系）（青系〜赤系）の組合せに応じて、標高Ｚによる段彩と尾根谷度αによる段彩とを識別できるように、適宜な重み係数ｗ１，ｗ２を設定することができる。
【００５７】
この第１合成画像Ｇ１の表示エリアでは、標高Ｚに応じた緑〜茶の範囲の段彩と、尾根谷度αに応じた青〜赤の範囲の段彩とが色合成され、平地が薄緑、標高の高い部分が茶鼠でその中に茶色の尾根部が脈状に延び、標高の次に高い部分が黄茶でその中に濃黄色の尾根部が脈状に延び、さらに脈状に延びる谷部が青緑で表現されている。この第１合成画像Ｇ１では、段彩された色の違いから、比較的大きな起伏と比較的小さな凹凸の両方の判読が可能である。
【００５８】
ステップＳ１２０では、対象点の近傍８点のデータ（Ｘ，Ｙ，Ｚ）を取得する。ステップＳ１３０では、対象点の近傍８点のデータ（Ｘ，Ｙ，Ｚ）を用いて傾斜度θを計算する。図７は、対象点の周囲近傍に位置する８点の標高値を示す。まず対象点の周囲近傍８点の標高値Ｚ１〜Ｚ８を用いて、ｘ方向の傾斜度Ｓｘとｙ方向の傾斜度Ｓｙとを次式により計算する。
Ｓｘ＝（（Ｚ１+Ｚ４+Ｚ６）−（Ｚ３+Ｚ５+Ｚ８））／６／Ｍsize …（１１）
Ｓｙ＝（（Ｚ１+Ｚ２+Ｚ３）−（Ｚ６+Ｚ７+Ｚ８））／６／Ｍsize …（１２）
ここで、Ｍsizeはメッシュサイズ（１画素１辺のサイズ（長さ））である。次にｘ方向の傾斜度Ｓｘとｙ方向の傾斜度Ｓｙとを用いて、次式により傾斜度（傾斜角）θを計算する。
傾斜度θ＝arctan（√（Ｓx・Ｓx ＋Ｓy・Ｓy） …（１３）
このように本実施形態では、ｘ方向の傾斜度Ｓｘとｙ方向の傾斜度Ｓｙとを合成して傾斜度θを計算している。この傾斜度θの使用によって、対象点Ｐjの周囲８方位の傾斜角のうち最大傾斜角を使用して陰影を付与する方法に比べ、周囲８方位の傾斜角が適度にスムージングされ、傾斜度θに応じた自然な陰影が付与される。もちろん、「（Ｚ１+Ｚ４+Ｚ６）／３−（Ｚ３+Ｚ５+Ｚ８）／３」のようなｘ方向両側の各３点の標高平均の差をメッシュサイズＭsizeで除した値を傾斜度θに用いてもよい。なお、傾斜度θとして式 θ＝√（Ｓx・Ｓx ＋Ｓy・Ｓy）を使用することもできる。
【００５９】
ステップＳ１４０では、傾斜度θにより陰影Ｇθを計算する。最大傾斜度を「黒」、最小傾斜度を「白」とする予め用意された傾斜度陰影変換テーブル（グレースケールテーブル）を参照して、傾斜度θに応じた陰影Ｇθを求める。本例では、陰影Ｇθは、０〜１．０の範囲内の値として算出される。この陰影Ｇθを空間座標（Ｘ，Ｙ，Ｚ）に付与した全点を画素としてグレースケール表示したものが、図１３に示す傾斜陰影画像Ｇθである。この傾斜陰影画像Ｇθの表示エリアでは、傾斜度θに応じた陰影が表現され、傾斜が大きい領域が濃色（低明度）、傾斜が小さい領域が淡色（高明度）のグレースケールで表現されている。例えば脈状に延びる尾根部や谷部、平地などは淡灰色（明灰色）で表現され、尾根と谷の間の斜面は傾斜に応じた濃さの濃灰色（暗灰色）で表現される。
【００６０】
ステップＳ１５０では、第１合成データＧ１と陰影Ｇθとを乗算し、第２合成データＧ２を生成する。第２合成データＧ２（Ｒzαθ，Ｇzαθ，Ｂzαθ）は、以下の計算式に従って算出される。
Ｒzαθ＝Ｒzα・Ｇθ …（１４）
Ｇzαθ＝Ｇzα・Ｇθ …（１５）
Ｂzαθ＝Ｂzα・Ｇθ …（１６）
この第２合成データＧ２（Ｒzαθ，Ｇzαθ，Ｂzαθ）を空間座標（Ｘ，Ｙ，Ｚ）に付与した全点を画素としてカラー表示したものが、図１４に示す第２合成画像Ｇ２である。この第２合成画像Ｇ２の表示エリアでは、標高Ｚに応じた緑〜茶の範囲の段彩と、尾根谷度αに応じた青〜赤の範囲の段彩と、傾斜度θに応じた陰影とが色合成され、平地が薄緑、標高の高い部分が焦茶でその中に明るい茶色で尾根部が脈状に延び、標高の次に高い部分の中に黄茶の尾根部が脈状に延び、さらに脈状に延びる谷部が青緑で表現されている。そして、尾根部の両側では傾斜度の大きな斜面が濃灰色の混色、傾斜度の小さな斜面が明灰色の混色で表現される。よって、この第２合成画像Ｇ２では、段彩された色の違いから、比較的大きな起伏と比較的小さな凹凸の両方の判読が可能であるうえ、陰影から傾斜の大きさも判読可能である。なお、第２合成データＧ２（Ｒzαθ，Ｇzαθ，Ｂzαθ）は、空間座標（Ｘ，Ｙ，Ｚ）と共に、空間座標と色値とを合わせた６次元のデータ（Ｘ，Ｙ，Ｚ，Ｒzαθ，Ｇzαθ，Ｂzαθ）として、バッファ３５に格納される。
【００６１】
こうして１つの対象点についてＳ７０〜Ｓ１５０の処理を終え、対象点に２種類の段彩及び陰影が合成されて付与されると、ステップＳ１６０において、全点について処理を終了したか否かを判断する。全点の処理終了前であって処理すべき座標点が残っていれば、ステップＳ１７０において次の点の番号を１つ繰り上げて（ｊ＝ｊ＋１）、ステップＳ６０において次の対象点Ｐｊを設定する。以下、同様に、各対象点についてＳ７０〜Ｓ１５０の処理を繰り返し、全ての座標点について処理を終えると、Ｓ１６０において全点の処理を終了したと判断される。このときバッファ３５には、全点の第２合成データＧ２（Ｘ，Ｙ，Ｚ，Ｒzαθ，Ｇzαθ，Ｂzαθ）が格納される。
【００６２】
ステップＳ１８０では、段彩地形データＧ（Ｘ，Ｙ，Ｚ，Ｒ，Ｇ，Ｂ）を生成する。すなわち、バッファ３５に格納された第２合成データＧ２（Ｘ，Ｙ，Ｚ，Ｒzαθ，Ｇzαθ，Ｂzαθ）の全点分の集まりが段彩地形データＧ（Ｘ，Ｙ，Ｚ，Ｒ，Ｇ，Ｂ）として生成される。この段彩地形データＧ（Ｘ，Ｙ，Ｚ，Ｒ，Ｇ，Ｂ）は空間座標値と色値とを成分に含む６次元の数値データであり、各点は空間座標点であると共に画素として機能する。この段彩地形データＧ（Ｘ，Ｙ，Ｚ，Ｒ，Ｇ，Ｂ）の各点を画素としてカラー表現すると、段彩地形データＧが第２合成データＧ２の集まりであることから、図１４と同じ画像となる。
【００６３】
ステップＳ１９０では、水部ポリゴンデータＰＤ（別途Ｘ，Ｙ座標列で表される閉じた多角形面）を取得し、段彩地形データＧ（Ｘ，Ｙ，Ｚ，Ｒ，Ｇ，Ｂ）において水部ポリゴンデータＰＤで規定されるポリゴン内の画素を色変更して段彩地形データＧＰ（Ｘ，Ｙ，Ｚ，Ｒ，Ｇ，Ｂ）を生成する。このポリゴン処理はポリゴン処理部２９が行う。このポリゴン処理の結果、池、沼、河川、海等の水面領域に設定された水部ポリゴン内に存在する全ての点（画素）の色値（Ｒ，Ｇ，Ｂ）が、例えば青色の値に変更される。この段彩地形データＧＰ（Ｘ，Ｙ，Ｚ，Ｒ，Ｇ，Ｂ）の各点を画素としてカラー表現したものが、図１５に示す段彩地形画像ＧＰである。
【００６４】
ステップＳ２００では、段彩地形データＧＰ（Ｘ，Ｙ，Ｚ，Ｒ，Ｇ，Ｂ）を基に、指定の方向視からの地形起伏画像データＧ（Ｒ，Ｇ，Ｂ）を生成する。この地形起伏画像データＧ（Ｒ，Ｇ，Ｂ）の生成は、画像生成部３０が行う。例えば指定方向視が真上からである平面視の場合は、段彩地形データＧＰ（Ｘ，Ｙ，Ｚ，Ｒ，Ｇ，Ｂ）をｚ方向から見たときの出力面（表示面）に指定エリア内の各点（画素）を写像し、出力面上の各写像点に対応する画素値（Ｒ，Ｇ，Ｂ）を与えることで、平面視の地形起伏画像データＧ（Ｒ，Ｇ，Ｂ）を生成する。この結果、生成された地形起伏画像データＧ（Ｒ，Ｇ，Ｂ）に基づき、図１６に示す平面視（ｚ方向視）の地形起伏画像が出力（表示又は印刷）される。また、指定方向視が斜め方向である場合は、段彩地形データＧＰ（Ｘ，Ｙ，Ｚ，Ｒ，Ｇ，Ｂ）をｚ方向に対して斜めから見たときの出力面（表示面）に指定エリア内の各点（画素）を写像し、出力面上の各写像点（画素）に、対応する画素値（Ｒ，Ｇ，Ｂ）を与えることで、斜めから見た地形起伏画像データＧ（Ｒ，Ｇ，Ｂ）を生成する。この結果、生成された地形起伏画像データＧ（Ｒ，Ｇ，Ｂ）に基づき、図１７に示す地形起伏画像（鳥瞰図）が出力（表示又は印刷）される。
【００６５】
従来、三次元地形データＤＥＭ（Ｘ，Ｙ，Ｚ）を基に鳥瞰図を生成する場合、まず地形表面を生成してからレンダリング処理で地形表面に色彩を付与する方法がとられていたので、処理負担の高いレンダリング処理により表示速度が遅くなってしまう。これに対し本実施形態では、三次元座標（Ｘ，Ｙ，Ｚ）に色値（Ｒ，Ｇ，Ｂ）を付加して生成される段彩地形データＧＰには三次元情報が保たれているので、段彩地形データＧＰを出力面（表示面）へ写像する際の向き（視方向）を変えるだけで、鳥瞰図などの画像データを比較的高速に生成できる。この結果、鳥瞰図のような三次元の画像も比較的高速に表示できる。
【００６６】
なお、図９において、Ｓ１０，Ｓ２０，Ｓ７０，Ｓ８０の処理が標高段彩処理手順に相当し、Ｓ３０，Ｓ４０の処理が凹凸度演算手順に相当し、さらにＳ５０，Ｓ９０，Ｓ１００の処理が凹凸度段彩処理手順に相当する。また、Ｓ１２０，Ｓ１３０の処理が傾斜度演算手順に相当し、Ｓ１４０の処理が陰影生成手順に相当し、さらにＳ１１０，Ｓ１５０，Ｓ１８０の処理が合成処理手順に相当する。また、Ｓ２００の処理が出力手順に相当する。
【００６７】
以上詳述したようにこの実施形態によれば、以下の効果が得られる。
（１）標高段彩と尾根谷度段彩（凹凸度段彩）との合成により、標高段彩により比較的大きな起伏を判読し易くなり、尾根谷度段彩により尾根・谷などの比較的小さな凹凸（微地形）を判読し易くなる。
【００６８】
（２）微地形を強調する陰影を傾斜度θから作成するので、斜めの入射光を想定して陰影を施した際に問題となる陰影の光入射方向依存性がないうえ、特許文献１のように凹凸度（浮沈度）に応じた陰影に起因し赤色立体地図において谷部が暗くなる問題も回避できる。尾根谷度段彩画像に傾斜度θによる陰影を重ね合せることで、明暗だけでは区別し難い尾根と谷を一層明確に判読できる。また、傾斜度θに応じて、緩傾斜部分の明度を高く、急傾斜部分の明度を低くすることで立体感を与えることができる。
【００６９】
（３）標高段彩は、標高値Ｚに応じた第１色彩範囲「緑（最低標高Ｚmin）〜茶（最高標高Ｚmax）」を着色することで強調表現し、尾根谷度段彩は、尾根谷度αに応じた第２色彩範囲「青（最小値αmin）〜赤（最大値αmax）」を着色することで強調表現される。第１色彩範囲と第２色彩範囲が異なる色彩範囲（つまりパラメータ（標高Ｚ、尾根谷度α）の最小値の色彩と最大値の色彩が互いに異なるよう）に設定するとともに、共にパラメータの最小値が互いに異なる寒色、最大値が互いに異なる暖色になるように設定した。よって、標高段彩と尾根谷度段彩を合成表示しても、比較的大きな起伏の判読とその中にある比較的小さな凹凸（微地形）の判読が一層し易い。
【００７０】
（４）特許文献１における赤色化立体画像のように浮沈度に応じて赤色の明度で階調表示する方法であると、尾根などの凸部が膨張、谷などの凹部が収縮したイメージ（谷が細く感じる）、谷内が暗く見え難くなる。しかし、本実施形態によれば、尾根や谷などの凹凸が尾根谷度αに応じて青〜赤の色彩範囲の段彩で強調表示されるので、尾根などの凸部が膨張、谷などの凹部が収縮したイメージになることを回避できる。
【００７１】
（５）段彩地形データＧＰは、３次元座標（Ｘ，Ｙ，Ｚ）の付加情報として色情報（Ｒ，Ｇ，Ｂ）を付加することにより生成され、三次元座標（Ｘ，Ｙ，Ｚ）を保持している。このため、空間的に三次元の段彩地形データＧＰを斜めから見ることでより、立体感の効いた地形起伏画像を作成できる。通常、三次元データから面（surface）を発生させてからレンダリング処理を行って地形起伏画像を表現することが一般的であり、表示速度が遅くなってしまう。しかし、本実施形態によれば、三次元地形データにおける着色された点を表示する処理で済むので、地形起伏画像を高速に生成できる。
【００７２】
（６）従来の地形画像は単色や赤〜青の段彩が多いが、本実施形態は、緑〜茶色と一般的な地図帳の表現に近い馴染み易い色調で地形起伏画像を表現できる。また、大きな起伏は緑〜茶色の段彩、凹凸などの微地形は青〜赤の段彩で強調表現されるので、その合成色上に施された注記（赤、青、黒色）が背景に溶け込まず判読し易い効果も得られる。また、別途作成された水部ポリゴン（または水部属性）を使用して水部の色付けを行うことにより、さらに地図帳に近いより馴染み易い地形起伏画像を提供できる。
【００７３】
（７）凹凸度演算手順（Ｓ３０，Ｓ４０）では、三次元地形データＤＥＭにおける対象点毎に、対象点から複数の方向（例えば８方向）に地上角Φ1と地下角Φ2とをそれぞれ求め、地上角Φ1の平均である地上開度ΦLと、地下角Φ2の平均である地下開度ΨLとに基づき、式 α＝（ΦL−ΨL）／２により、凹凸度の一例として尾根谷度αを計算する。尾根などの凸部を表現できる地上開度ΦLと、谷などの凹部を表現できる地下開度ΨLとから、凹凸度を表現できる１つの地形量パラメータとして尾根谷度αを求め、尾根谷度αに応じて変換テーブルＴ２から決まる色彩を施すので、微地形を適切に段彩表現できる。
【００７４】
（８）傾斜度演算手順（Ｓ１２０，Ｓ１３０）では、対象点の周囲近傍８点の標高Ｚを用いて（１１）式及び（１２）式によりそれぞれ計算したｘ方向傾斜度Ｓxとｙ方向傾斜度Ｓyとを（１３）式により合成して傾斜度θを求める。よって、対象点から８方向の傾斜角が適度にスムージングされた傾斜度θを採用できる。この結果、適度にスムージングされた傾斜度θに応じた自然な陰影を付与できる。
【００７５】
（９）標高段彩データＣ１と尾根谷度段彩データＣ２との重み付け平均（加重平均）を、（８）式〜（１０）式により計算し、第１合成データＧ１を生成する。よって、係数ｗ１，ｗ２による重み付けを調整することにより、比較的大きな起伏と比較的小さな凹凸とを判読し易い条件で、段彩地形データＧＰを生成できる。
【００７６】
前記実施形態は上記に限定されず、以下の態様に変更することもできる。
・変換テーブルＴ１，Ｔ２の色彩範囲は前記実施形態に限定されず、適宜な色彩範囲を設定することができる。例えば各変換テーブルＴ１，Ｔ２で使用する色彩範囲を入れ替えてもよい。さらに、変換テーブルＴ１，Ｔ２のうち少なくとも一方について、標高Ｚ又は尾根谷度αの最小値側が暖色、最大値側が寒色となるような色彩範囲を設定してもよい。また、標高用と尾根谷度用で同一の変換テーブルを用いて、使用する色彩範囲を同じとしてもよい。
【００７７】
・標高Ｚに応じた第１段彩処理、尾根谷度αに応じた第２段彩処理、傾斜度に応じた陰影処理及び合成処理を、対象点毎に全点終了まで繰り返す方法を採用したが、これに替えて、全点の処理を終えた後に次の処理へ進む方法も採用できる。
【００７８】
・傾斜度を求める計算式は適宜変更してもよい。最大傾斜角を求める計算式を用いてもよい。
・尾根谷度αを求める計算式は適宜変更してもよい。また、探索方向は８方向に限らず、４方向や１６方向、あるいは３２方向等でもよい。
【００７９】
・凹凸度は尾根谷度αに限定されない。凹凸度として、ラプラシアンを用いてもよい。さらに、曲線上のある点を通る全ての測地線（曲面上で２点間を結ぶ最短距離の曲線）の曲率の最大値と最小値の平均として定義される平均曲率を凹凸度として用いてもよい。
【００８０】
・三次元地形データは、ＴＩＮ（triangulated irregular network，不整三角形網）でもよい。さらには、オルソ画像、ＤＴＭ（Digital Terrain Model：数値地形モデル）、ＤＳＭ（Digital Surface Model：数値表層モデル）でもよい。
【００８１】
・段彩地形データは６次元データに限定されない。他の情報を含めた７次元以上のデータであってもよい。また、三次元空間座標（Ｘ，Ｙ，Ｚ）を無くした三次元の段彩地形データＧＰ（Ｒ，Ｇ，Ｂ）であってもよい。この場合、指定の視野方向毎に段彩地形データの生成処理を行うことになる。
【符号の説明】
【００８２】
１０…本体、１１…入力部、１２…入力操作部、１３…表示ドライバ、１４…モニタ、１５…印刷ドライバ、１６…プリンタ、２１…制御部、２２…入力バッファ、２３…標高段彩処理手段の一例である第１段彩演算部、２４…凹凸度演算手段の一例である尾根谷度演算部、２５…凹凸度段彩処理手段の一例である第２段彩演算部、２６…傾斜度演算手段の一例である傾斜度演算部、２７…陰影生成手段の一例である陰影演算部、２８…合成処理手段の一例である合成処理部、２９…ポリゴン処理部、３０…出力手段の一例である画像生成部、３１…第１変換テーブル作成部、３２…第２変換テーブル作成部、３３…色合成処理部、３４…陰影合成処理部、３５…バッファ、１００…地形起伏画像生成装置の一例である地図生成装置、ＰＣ…パーソナルコンピュータ、ＤＥＭ（Ｘ，Ｙ，Ｚ）…三次元地形データ、Ｚ…標高、Φ1…地上角、Φ2…地下角、ΦL…地上開度、ΨL…地下開度、α…凹凸度の一例である尾根谷度、Ｓx…ｘ方向傾斜度、Ｓy…ｙ方向傾斜度、θ…傾斜度、Ｔ１…高色彩変換テーブル、Ｔ２…尾根谷度色彩変換テーブル、Ｚmin…標高の最小値、Ｚmax…標高の最大値、αmin…尾根谷度の最小値、αmax…尾根谷度の最大値、Ｐj…対象点（着目の地点）Ｃ１…標高段彩データ、Ｃ２…凹凸度段彩データの一例である尾根谷度段彩データ、Ｇθ…陰影データ（陰影）、Ｇ１…第１合成データ、Ｇ２…第２合成データ、Ｇ（Ｘ，Ｙ，Ｚ，Ｒ，Ｇ，Ｂ）…段彩地形データ（水部ポリゴン色付与前）、ＧＰ（Ｘ，Ｙ，Ｚ，Ｒ，Ｇ，Ｂ）…段彩地形データ（水部ポリゴン色付与後）、Ｇ（Ｒ，Ｇ，Ｂ）…地形起伏画像データ。
【特許請求の範囲】
【請求項１】
標高を含む三次元地形データを基に地形起伏画像を生成する地形起伏画像生成方法であって、
入力した前記三次元地形データに基づき画素となる各地点に当該地点の標高に応じた色彩を付与する段彩処理を行って標高段彩データを生成する標高段彩処理手順と、
前記三次元地形データに基づいて前記各地点における地形の凹凸度を演算する凹凸度演算手順と、
前記各地点に当該地点毎の前記凹凸度に応じた色彩を付与する段彩処理を行って凹凸度段彩データを生成する凹凸度段彩処理手順と、
前記三次元地形データに基づいて前記各地点における傾斜度を演算する傾斜度演算手順と、
前記各地点に当該地点毎の前記傾斜度に応じた陰影を付与する陰影処理を行って陰影データを生成する陰影生成手順と、
前記標高段彩データと前記凹凸度段彩データと前記陰影データとを合成して段彩地形データを生成する合成処理手順と、
前記段彩地形データに基づく地形起伏画像を出力する出力手順と、
を備えたことを特徴とする地形起伏画像生成方法。
【請求項２】
前記標高段彩処理手順では、前記標高の最小値と最大値との間に設定された第１の色彩範囲に従って標高に応じた色彩を付与して前記標高段彩データを生成し、
前記凹凸度段彩処理手順では、前記凹凸度の最小値と最大値との間に設定された第２の色彩範囲に従って凹凸度に応じた色彩を付与して前記凹凸度段彩データを生成し、
前記第１の色彩範囲と前記第２の色彩範囲は設定された色彩範囲が互いに異なるとともに、前記標高と前記凹凸度との各最小値に対応する色彩が互いに異なる寒色、最大値に対応する色彩が互いに異なる暖色に設定されることを特徴とする請求項１に記載の地形起伏画像生成方法。
【請求項３】
前記凹凸度演算手順では、前記三次元地形データにおける前記地点毎に、当該地点から複数の方向に地上角と地下角とをそれぞれ求め、前記地上角の平均である地上開度と、前記地下角の平均である地下開度とに基づき、前記凹凸度としての尾根谷度を演算することを特徴とする請求項１又は２に記載の地形起伏画像生成方法。
【請求項４】
前記傾斜度演算手順では、前記三次元地形データに基づいて前記各地点におけるｘ方向傾斜度とｙ方向傾斜度とを演算し、当該ｘ方向傾斜度とｙ方向傾斜度とを合成して前記傾斜度を演算することを特徴とする請求項１乃至３のいずれか一項に記載の地形起伏画像生成方法。
【請求項５】
前記合成処理手順では、前記三次元地形データにおける前記各地点に対応する三次元の空間座標点に、前記合成処理手順で取得した色情報を三成分の色座標で付加することにより、三次元の空間座標と三成分の色座標とを含む六次元以上の前記段彩地形データを生成し、
前記出力手順では、前記段彩地形データの空間座標点を指定の視野方向にある出力面に写像した点を画素として当該画素に当該空間座標に対応する色座標を画素値として与えることで前記地形起伏画像を生成することを特徴とする請求項１乃至４のいずれか一項に記載の地形起伏画像生成方法。
【請求項６】
標高を含む三次元地形データを基に地形起伏画像を生成する地形起伏画像生成装置であって、
入力した前記三次元地形データに基づき画素となる各地点に当該地点の標高に応じた色彩を付与する段彩処理を行って標高段彩データを生成する標高段彩処理手段と、
前記三次元地形データに基づいて前記各地点における地形の凹凸度を演算する凹凸度演算手段と、
前記各地点に当該地点毎の前記凹凸度に応じた色彩を付与する段彩処理を行って凹凸度段彩データを生成する凹凸度段彩処理手段と、
前記三次元地形データに基づいて前記各地点における傾斜度を演算する傾斜度演算手段と、
前記各地点に当該地点毎の前記傾斜度に応じた陰影を付与する陰影処理を行って陰影データを生成する陰影生成手段と、
前記標高段彩データと前記凹凸度段彩データと前記陰影データとを合成して段彩地形データを生成する合成処理手段と、
前記段彩地形データに基づく地形起伏画像を出力する出力手段と、
を備えたことを特徴とする地形起伏画像生成装置。
【請求項１】
標高を含む三次元地形データを基に地形起伏画像を生成する地形起伏画像生成方法であって、
入力した前記三次元地形データに基づき画素となる各地点に当該地点の標高に応じた色彩を付与する段彩処理を行って標高段彩データを生成する標高段彩処理手順と、
前記三次元地形データに基づいて前記各地点における地形の凹凸度を演算する凹凸度演算手順と、
前記各地点に当該地点毎の前記凹凸度に応じた色彩を付与する段彩処理を行って凹凸度段彩データを生成する凹凸度段彩処理手順と、
前記三次元地形データに基づいて前記各地点における傾斜度を演算する傾斜度演算手順と、
前記各地点に当該地点毎の前記傾斜度に応じた陰影を付与する陰影処理を行って陰影データを生成する陰影生成手順と、
前記標高段彩データと前記凹凸度段彩データと前記陰影データとを合成して段彩地形データを生成する合成処理手順と、
前記段彩地形データに基づく地形起伏画像を出力する出力手順と、
を備えたことを特徴とする地形起伏画像生成方法。
【請求項２】
前記標高段彩処理手順では、前記標高の最小値と最大値との間に設定された第１の色彩範囲に従って標高に応じた色彩を付与して前記標高段彩データを生成し、
前記凹凸度段彩処理手順では、前記凹凸度の最小値と最大値との間に設定された第２の色彩範囲に従って凹凸度に応じた色彩を付与して前記凹凸度段彩データを生成し、
前記第１の色彩範囲と前記第２の色彩範囲は設定された色彩範囲が互いに異なるとともに、前記標高と前記凹凸度との各最小値に対応する色彩が互いに異なる寒色、最大値に対応する色彩が互いに異なる暖色に設定されることを特徴とする請求項１に記載の地形起伏画像生成方法。
【請求項３】
前記凹凸度演算手順では、前記三次元地形データにおける前記地点毎に、当該地点から複数の方向に地上角と地下角とをそれぞれ求め、前記地上角の平均である地上開度と、前記地下角の平均である地下開度とに基づき、前記凹凸度としての尾根谷度を演算することを特徴とする請求項１又は２に記載の地形起伏画像生成方法。
【請求項４】
前記傾斜度演算手順では、前記三次元地形データに基づいて前記各地点におけるｘ方向傾斜度とｙ方向傾斜度とを演算し、当該ｘ方向傾斜度とｙ方向傾斜度とを合成して前記傾斜度を演算することを特徴とする請求項１乃至３のいずれか一項に記載の地形起伏画像生成方法。
【請求項５】
前記合成処理手順では、前記三次元地形データにおける前記各地点に対応する三次元の空間座標点に、前記合成処理手順で取得した色情報を三成分の色座標で付加することにより、三次元の空間座標と三成分の色座標とを含む六次元以上の前記段彩地形データを生成し、
前記出力手順では、前記段彩地形データの空間座標点を指定の視野方向にある出力面に写像した点を画素として当該画素に当該空間座標に対応する色座標を画素値として与えることで前記地形起伏画像を生成することを特徴とする請求項１乃至４のいずれか一項に記載の地形起伏画像生成方法。
【請求項６】
標高を含む三次元地形データを基に地形起伏画像を生成する地形起伏画像生成装置であって、
入力した前記三次元地形データに基づき画素となる各地点に当該地点の標高に応じた色彩を付与する段彩処理を行って標高段彩データを生成する標高段彩処理手段と、
前記三次元地形データに基づいて前記各地点における地形の凹凸度を演算する凹凸度演算手段と、
前記各地点に当該地点毎の前記凹凸度に応じた色彩を付与する段彩処理を行って凹凸度段彩データを生成する凹凸度段彩処理手段と、
前記三次元地形データに基づいて前記各地点における傾斜度を演算する傾斜度演算手段と、
前記各地点に当該地点毎の前記傾斜度に応じた陰影を付与する陰影処理を行って陰影データを生成する陰影生成手段と、
前記標高段彩データと前記凹凸度段彩データと前記陰影データとを合成して段彩地形データを生成する合成処理手段と、
前記段彩地形データに基づく地形起伏画像を出力する出力手段と、
を備えたことを特徴とする地形起伏画像生成装置。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図９】
【図８】
【図１０】
【図１１】
【図１２】
【図１３】
【図１４】
【図１５】
【図１６】
【図１７】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図９】
【図８】
【図１０】
【図１１】
【図１２】
【図１３】
【図１４】
【図１５】
【図１６】
【図１７】
【公開番号】特開２０１２−３４００（Ｐ２０１２−３４００Ａ）
【公開日】平成２４年１月５日（２０１２．１．５）
【国際特許分類】
【出願番号】特願２０１０−１３６１９４（Ｐ２０１０−１３６１９４）
【出願日】平成２２年６月１５日（２０１０．６．１５）
【出願人】（５９３０２６９５９）中日本航空株式会社 (8)
【Ｆターム（参考）】
【公開日】平成２４年１月５日（２０１２．１．５）
【国際特許分類】
【出願日】平成２２年６月１５日（２０１０．６．１５）
【出願人】（５９３０２６９５９）中日本航空株式会社 (8)
【Ｆターム（参考）】
[ Back to top ]