(57)【特許請求の範囲】
【請求項１】
緯度経度で定義された数値標高モデルを記憶した数値標高モデル用メモリと、
前記数値標高モデルに緯度経度で定義されている一定の大きさのメッシュを順次、平面直角座標系で定義された平面直角メッシュ用メモリに投影変換した平面直角メッシュ群を生成する平面直角座標変換部と、
平面直角座標系で定義されたラスタ用メモリと、
前記平面直角メッシュ用メモリの前記平面直角メッシュ毎に、この平面直角メッシュの
Ｘ軸方向となる横辺の長さを奇数分割し、これらの分割辺の長さを平均化した平均長で前
記ラスタ用メモリのＸ軸、Ｙ軸を分割して前記ラスタ用メモリのＸ−Ｙ平面に縦横が前記
平均長の微細格子の微細格子群を定義する微細格子生成処理部と、
前記ラスタ用メモリの微細格子毎の標高値を内挿補間して割付ける標高値補間部と、
前記ラスタ用メモリの全ての微細格子の内挿補間後の標高値の内で、起点となる前記微
細格子を指定し、この指定した微細格子と同じ標高値を有して閉曲する微細格子群を通る
直線を求め、これらの直線をベクター化し、これを等高線ベクターとして生成する等高線
生成部と、
前記等高線ベクターを画像にして表示用メモリに書き込んで画面に表示する表示処理部
と、
を有することを特徴とする等高線生成システム。
【請求項２】
縦列及び横列ともに前記奇数分割の個数にした移動平均化メッシュを生成し、この移動
平均化メッシュを前記ラスタ用メモリの前記微細格子に順次、かけて前記内挿補間後の標
高値を順次、平滑し、この平滑した滑らか標高値に基づいた滑らか等高線画像を生成し、
この滑らか等高線画像を前記表示処理部によって表示させる滑か処理部と、
を有することを特徴とする請求項１記載の等高線生成システム。
【請求項３】
前記ラスタ用メモリにおける平滑後の微細格子を着目点とし、この着目点毎に、この着
目点からの考慮距離を、この考慮距離に対応する微細格子数で定義し、この微細格子数内
に存在する滑らか標高値に基づいて尾根谷度を求め、この尾根谷度に応じた諧調色の画像
を生成し、画像重ね合わせ部によって前記滑らか等高線画像にこの諧調色の画像を重ねて
、前記表示処理部によって表示させる画像処理部と、
を有することを特徴とする請求項２記載の等高線生成システム。
【請求項４】
前記画像処理部は、
前記微細格子数内に存在する滑らか標高値に基づいて斜度と地上開度と地下開度とを求
め、この地下開度と地上開度とより前記尾根谷度を求め、斜度を赤系の彩度に割り当て、
尾根谷度を明度に割り当てることで前記諧調色の画像を生成することを特徴とする請求項
３記載の等高線生成システム。
【請求項５】
道路、建物、河川、沼がベクター情報で定義された２万５千分の１の標準地図情報が記
憶された標準地図用メモリを備え、
前記画像重ね合わせ部は、
前記標準地図情報を画像化し、この画像を前記滑らか等高線画像又は諧調色の画像に重
ね合わせて前記表示処理部によって表示させることを特徴とする請求項３記載の等高線生
成システム。
【請求項６】
前記平均長のサイズを前記表示用メモリのセルの解像度として定義する解像度定義部と
、
を有することを特徴とする請求項１記載の等高線生成システム。
【請求項７】
前記数値標高モデルは、５ｍＤＥＭであり、前記奇数分割の数は、「９」であることを
特徴とする請求項１記載の等高線生成システム。
【請求項８】
前記内挿補間は、ＴＩＮバイリニア補間であることを特徴とする請求項１記載の等高線
生成システム。
【請求項９】
前記数値標高モデル用メモリから所定のエリアに存在する前記一定の大きさのメッシュ
を構成する点群を読み出して前記平面直角座標変換部に出力するエリア読込部と、
を有することを特徴とする請求項１記載の等高線生成システム。
【請求項１０】
コンピュータに、
緯度経度で定義された数値標高モデルを数値標高モデル用メモリに記憶する手段、
前記数値標高モデルに緯度経度で定義されている一定の大きさのメッシュを順次、平面
直角座標系で定義された平面直角メッシュ用メモリに投影変換した平面直角メッシュ群を
生成する手段、
ラスタ用メモリを平面直角座標系で定義する手段、
前記平面直角メッシュ用メモリの前記平面直角メッシュ毎に、この平面直角メッシュの
Ｘ軸方向となる横辺の長さを奇数分割し、これらの分割辺の長さを平均化した平均長で前
記ラスタ用メモリのＸ軸、Ｙ軸を分割して前記ラスタ用メモリのＸ−Ｙ平面に縦横が前記
平均長の微細格子の微細格子群を定義する手段、
前記ラスタ用メモリの微細格子毎の標高値を内挿補間して割付ける手段、
前記ラスタ用メモリの全ての微細格子の内挿補間後の標高値の内で、起点となる前記微
細格子を指定し、この指定した微細格子と同じ標高値を有して閉曲する微細格子群を通る
直線を求め、これらの直線をベクター化し、これを等高線ベクターとして生成する手段、
前記等高線ベクターを画像にして表示用メモリに書き込んで表示する表示手段、
としての機能を実行させる等高線生成プログラム。
【請求項１１】
コンピュータに、
縦列及び横列ともに前記奇数分割の個数にした移動平均化メッシュを生成し、この移動
平均化メッシュを前記ラスタ用メモリの前記微細格子に順次、かけて前記内挿補間後の標
高値を順次、平滑し、この平滑した滑らか標高値に基づいた滑らか等高線画像を生成し、
この滑らか等高線画像を前記表示手段によって表示させる滑らか等高線性生成手段、
としての機能を実行させる請求項１０記載の等高線生成プログラム。
【請求項１２】
コンピュータに、
前記ラスタ用メモリにおける平滑後の微細格子を着目点とし、この着目点毎に、この着
目点からの考慮距離を、この考慮距離に対応する微細格子数で定義し、この微細格子数内
に存在する滑らか標高値に基づいて尾根谷度を求め、この尾根谷度に応じた諧調色の画像
を生成し、画像重ね合わせ手段によって前記滑らか等高線画像にこの諧調色の画像を重ね
、これを前記表示手段によって表示させる画像処理手段、
としての機能を実行させる請求項１１記載の等高線生成プログラム。
【請求項１３】
コンピュータに、
前記画像処理手段は、
前記微細格子数内に存在する滑らか標高値に基づいて斜度と地上開度と地下開度とを求
め、この地下開度と地上開度とより前記尾根谷度を求め、斜度を赤系の彩度に割り当て、
尾根谷度を明度に割り当てることで前記諧調色の画像を生成させることを特徴とする請求
項１２記載の等高線生成プログラム。
【請求項１４】
コンピュータに、
道路、建物、河川、沼がベクター情報で定義された２万５千分の１の標準地図情報を標
準地図用メモリに記憶する手段、
前記画像重ね合わせ手段として、
前記標準地図情報を画像化し、この画像を前記滑らか等高線画像又は諧調色の画像に重
ね合わせて前記表示手段によって表示させることを特徴とする請求項１２記載の等高線生成プログラム。
【請求項１５】
コンピュータに、
前記平均長のサイズを前記表示用メモリのセルの解像度として定義する手段、
としての機能を実行させるための請求項１０記載の等高線生成プログラム。
【請求項１６】
前記数値標高モデルは、５ｍＤＥＭであり、前記奇数分割の数は、「９」であることを
特徴とする請求項１０記載の等高線生成プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、等高線生成システムに関する。
【背景技術】
【０００２】
近年は国土地理院（以下、地理院という）が、インターネット網でデジタル標高モデル（DEM： Digital Elevation Model）を公開している。ＤＥＭは、５ｍＤＥＭ、１０ｍＭＥＭが主流である。
【０００３】
このＤＥＭは、航空機からレーザ光を地上に照射して得た各々のレーザ計測点をＴＩＮ（triangulated irregular network）で結び、これに経度差、緯度差0.2秒（約5m）間隔若しくは経度差、緯度差0.4秒（約10ｍ）のメッシュで区切った枠を被せ、各々のメッシュの中心点の高さを、その中心点を含むＴＩＮの頂点の高さから内挿補間により求めている。なお、１０ｍＤＥＭは、等高線より手作業で作成するのが一般的である。
【０００４】
このようなＤＥＭ（５ｍＤＥＭ、１０ｍＤＥＭ）を用いて地形の標高を表現する等高線図が作成される。等高線は標高が等しい点の集まりを結んだ線であり、依然として様々な分野で用いられている。
【０００５】
例えば、地形の等高線は、国土地理院でネット上で公開されている２万５０００分の１の等高線図がある。この等高線は手作業によるのがほとんどであり、曲率最大化処理（スプライン曲線、ベジェ曲線、３次エルミート関数等）で等高線の曲線、直線を現わしている。
【０００６】
一方、不動産の価値判断、人、車両等の移動の安全性等の観点から崖部、平野部、山岳等の地形の起伏をさらに詳細に分かるような等高線図が求められる場合もある。
しかし、５ｍＤＥＭのメッシュを用いた場合は、等高線間隔は、５ｍ以下の詳細な等高線を得ることができない。
【０００７】
この詳細な等高線を得るため等高線生成方法として特許文献１が公開されている。
例えば、特許文献１の等高線図作成方法は、１０ｍ間隔の基本等高線図を画面に表示し、この基本等高線図上に縦横１０ｍ間隔で区切った区分線を定義する。
そして、この１０ｍメッシュを複数の格子で分割している（図５においては、サイズが２０ｃｍの格子を縦５、横５に分割している）。
【０００８】
そして、等高線と等高線との間の複数の格子の中から一つを選出して、着目点とし、この着目点を中心にして８方向に延長線を定義して、最も初めに交わる等高線を検索する。そして、一定範囲の等高線同士に交わる基準直線との標高値を得て、これらの標高値を結ぶ曲線を求めて、着目点の標高値を補間していく。そして、これらの補間標高値をを適用して結んで微細間隔の等高線を生成している。
【先行技術文献】
【特許文献】
【０００９】
【特許文献１】特開２００３−１８６３９３号公報
【発明の概要】
【発明が解決しようとする課題】
【００１０】
しかしながら、従来の地形の等高線図は、５ｍＤＥＭ以上のＤＥＭ（例えば、５ｍ、１０ｍ、２０ｍ、３０ｍ、・・）を用いて作成しているので、等高線と等高線との間が広くなり（ＤＥＭ間隔の２倍以上）、かつ曲率最大化処理（スプライン曲線、ベジェ曲線等）を用いて等高線を自動的に描いているので、隣接する等高線同士の平行感が損なわれる。このため、隣接する等高線同士の曲率感が損なわれる。
【００１１】
一方、近年は１ｍＤＥＭも存在するが、この１ｍＤＥＭは地理院では原則非公開であり、また取得されているエリアも限られている。
【００１２】
このため、１／２５００の等高線図から手作業で生成しないといけないので、非常に手間がかかる。
【００１３】
さらに、特許文献１の等高線生成方法は、既に生成されている基準等高線（ベクタ）上に、複数の格子を定義して、任意の格子の標高値を基準等高線の標高値と基準等高線の標高値とで内挿補間して、新たに等高線を生成する方法であるので、新たな等高線もこの基準等高線の影響を受けた曲線となる。
【００１４】
つまり、曲率最大化処理に基づく等高線であるので、なめらか処理を施したとしてもカーブの部分が急激に変化したり、あるいは丸みを帯びすぎることになる。
【００１５】
本発明は以上の課題を鑑みてなされたもので、ＤＥＭを用いて曲率最大化処理を行わなくとも、実地形を詳細に表現でき、かつ等高線同士の曲率感を損なわない高精度な等高線を生成する等高線生成システムを得ることを目的とする。
【課題を解決するための手段】
【００１６】
本発明に係る等高線生成システムは、緯度経度で定義された数値標高モデルを記憶した数値標高モデル用メモリと、
前記数値標高モデルに緯度経度で定義されている一定の大きさのメッシュを順次、平面直角座標系で定義された平面直角メッシュ用メモリに投影変換した平面直角メッシュ群を生成する平面直角座標変換部と、
平面直角座標系で定義されたラスタ用メモリと、
前記平面直角メッシュ用メモリの前記平面直角メッシュ毎に、この平面直角メッシュのＸ軸方向となる横辺の長さを奇数分割し、これらの分割辺の長さを平均化した平均長で前記ラスタ用メモリのＸ軸、Ｙ軸を分割して前記ラスタ用メモリのＸ−Ｙ平面に縦横が前記平均長の微細格子の微細格子群を定義する微細格子生成処理部と、
前記ラスタ用メモリの微細格子毎の標高値を内挿補間して割付ける標高値補間部と、
前記ラスタ用メモリの全ての微細格子の内挿補間後の標高値の内で、起点となる前記微細格子を指定し、この指定した微細格子と同じ標高値を有して閉曲する微細格子群を通る直線を求め、これらの直線をベクター化し、これを等高線ベクターとして生成する等高線生成部と、
前記等高線ベクターを画像にして表示用メモリに書き込んで画面に表示する表示処理部と、を備えたことを要旨とする。
【発明の効果】
【００１７】
以上のように本発明によれば、緯度経度で定義された一定の大きさのメッシュ（５ｍＤＥＭ）を平面直角座標系に順次、投影変換し、これらの平面直角メッシュのＸ方向を奇数分割して、これらの分割長の平均長（例えば、０．５５９ｍ）でラスター用メモリのＸ軸、Ｙ軸を分割して微細格子群（ｍｉ：０．５５９ｍ）を定義する。
【００１８】
そして、これらの微細格子（ｍｉ）毎の標高値を内挿補間して平滑化し、微細格子（ｍｉ）と同じ平滑後の標高値を有して閉曲する微細格子（ｍｉ）群を通る直線群の塊をベクター化した等高線情報を生成し、これを表示する。
【００１９】
すなわち、表示された等高線は、微細格子（例えば、０．５５９ｍ）と通る閉曲する直線の塊をベクター化した等高線（等高線）であるので、曲率最大化処理を行わなくとも、この等高線（等高線）は、実地形の形状に非常に近い曲線となり、かつ滑らかな曲線となる。このため、等高線（等高線）同士の平行感が得られるので、曲率感が損なわれない。
【００２０】
また、一般に公開されている５ｍＤＥＭのメッシュを約５０ｃｍ〜６０ｃｍの微細分割した格子の滑らか標高値をベクタ化しているので、１万分の１の等高線地図として利用できる。
【００２１】
また、等高線（等高線）と等高線（等高線）との間の高度差に応じた色で表示しても、表示用メモリのセルは、微細格子（０．５５９ｍ）と同じサイズであるので、ある程度、拡大したとしてもジャギーは抑えられる。
【００２２】
また、エリアをこれらの平面直角メッシュのＸ方向を奇数分割（例えば、９分割）して、これらの分割長の平均長（例えば、０．５５９ｍ）の微細格子にしているので、等高線（等高線）と等高線（等高線）との間の高度差に応じた色で表示してもアーチファクトが抑えられる。
【００２３】
さらに、標準地図（レベル１６）と合成しても違和感もない。
【図面の簡単な説明】
【００２４】
【図１】本実施の形態の等高線生成システムの概念を説明する説明図である。
【図２】５ｍＤＥＭで生成した等高線図画像ＧＯｉの説明図である。
【図３】本実施の形態の乗算合成画像ＧＳＪｉ（標準地図画像ＧＯｉ＋滑か等高線画像ＧＪｉ）の説明図である。
【図４】本実施の形態の等高線生成システムの概略構成図である。
【図５】地理院の標準地図画像ＧＯｉ（レベル１６）の拡大図である。
【図６】緯度経度５ｍメッシュＭａｉと等高線の関係の説明図である。
【図７】５ｍＤＥＭ緯度経度ポイントデータＰａｉの具体例の説明図である。
【図８】平面直角５ｍメッシュＭｂｉと等高線の関係の説明図である。
【図９】緯度経度変換と平面直角座標変換との関係を説明する説明図である。
【図１０】Ｘ方向分割平均距離Ｂｄｉ（０．５５９ｍ）で分割されたＸ−Ｙ平面の説明図である    。
【図１１】分割後の微細格子ｍｉの三角形の補間後の座標例の説明図である。
【図１２】分割後の点群と平面直角メッシュと等高線との関係を説明する説明図である。
【図１３】平面直角５ｍメッシュＭｂｉと微細格子ｍｉ（０．５５９ｍ）との関係を説明する説明図である。
【図１４】ＴＩＮバイナリ補間後の平面直角５ｍメッシュと等高線ｃｊとの関係を説明する説明図である。
【図１５】従来の５ｍＤＥＭによる標高値と本実施の形態の微細格子による標高値の説明図である。
【図１６】本実施の形態の移動平均化用メッシュＨＭｉの説明図である。
【図１７】本実施の形態の移動平均化用メッシュＨＭｉを１回かけた場合の等高線段彩図の説明図である。
【図１８】本実施の形態の移動平均化用メッシュＨＭｉを２回かけた場合の等高線段彩図の説明図である。
【図１９】移動平均化用メッシュＨＭｉを２回かけた場合の等高線段彩図の拡大図である。
【図２０】滑か微細標高値ラスタ画像データＲＧｉの説明図である。
【図２１】等高線の生成を説明する説明図である。
【図２２】実施の形態の滑らか滑か等高線画像ＧＪｉとスムージング処理を行わなかった場合の滑か等高線画像ＧＪｉの違いの説明図（１）である。
【図２３】実施の形態の滑らか滑か等高線画像ＧＪｉとスムージング処理を行わなかった場合の滑か等高線画像ＧＪｉの違いの説明図（２）である。
【図２４】実施の形態の滑らか滑か等高線画像ＧＪｉとスムージング処理を行わなかった場合の滑か等高線画像ＧＪｉの違いの説明図（３）である。
【図２５】実施の形態の滑らか滑か等高線画像ＧＪｉとスムージング処理を行わなかった場合の滑か等高線画像ＧＪｉの違いの説明図（４）である。
【図２６】スムージング処理の前の補間標高値ｈｚｂｉに対して超解像度赤色立体画像作成処理を行った場合の説明図である。
【図２７】スムージング処理後の値に対して超解像度赤色立体画像作成処理を行った場合の説明図である。
【図２８】本実施の形態の乗算合成画像（赤色重ね）の説明図である。
【発明を実施するための形態】
【００２５】
本実施の形態では、地理院の５ｍＤＥＭ（Ａ：Ａはレーザを意味する）の基盤地図（以下、５ｍＤＥＭ基盤地図Ｆａという）を一例として等高線図を生成するとして説明する。
【００２６】
＜実施の形態＞
図１は本実施の形態の等高線生成システムの概念を説明する説明図であり、コンピュータが行う処理である。本実施の形態の等高線生成システムは、一般に広く提供されている国土地理院（以下、地理院という）の５ｍＤＥＭ（５ｍメッシュ）を微細分割し、この微細メッシュ（以下、微細格子ｍｉという）に割り付けられている標高値と同じ微細格子ｍｉを繋げたベクタ情報の等高線を生成する。この等高線は本実施の形態では滑か処理が施される（以下、滑か等高線ｃｊｉと称している）。
【００２７】
以下に図１を用いてこの生成方法の概略を説明する。
図１に示すように、任意のエリアＥｉ（例えば、石巻市）の緯度及び経度等で定義された５ｍＤＥＭのデータ（以下、緯度経度５ｍＤＥＭデータｐａｉという）を読込む（Ｓ１０）。
【００２８】
そして、これを等緯度経度ＸＹＺポイントファイルに順次、エクスポートする（Ｓ２０）。この緯度経度５ｍＤＥＭデータＰａｉ（Ｐａ１、Ｐａ２、・・）を繋げたメッシュ（以下、緯度経度５ｍＤＥＭメッシュＭａｉという）の例を図６に示している。
【００２９】
そして、これらの緯度経度５ｍＤＥＭデータＰａｉ（ｘａｉ、ｙａｉ、ｚａｉ）を平面直角座標系の座標に変換し（以下、平面直角５ｍＤＥＭデータＰｂｉ（ｘｂｉ、ｙｂｉ、ｚｂｉ））（Ｓ２０）、これを平面直角ＸＹＺポイントファイルに順次、エクスポートする（Ｓ３０）。この平面直角５ｍＤＥＭデータＰｂｉを繋げたメッシュ（以下、平面直角５ｍメッシュＭｂｉという）の例を図８に示している。
【００３０】
そして、ラスター化処理を行う（Ｓ５０）。このラスター化処理は、平面直角５ｍメッシュＭｂｉ（Ｍｂ１、Ｍｂ２・・、Ｍｂ００、・・）を順に指定し、この指定毎に、この平面直角５ｍメッシュＭｂｉの横辺（以下、単にＸ方向とも称する）を９分割するＸ方向分割距離Ｂｉ（例えば、０．５８４４ｍ又は０．５５９ｍ：小数点４以下は切り捨て）を求める。
【００３１】
そして、これらのＸ方向分割距離Ｂｉを平均化し（以下、Ｘ方向分割平均距離Ｂｄｉという）、後述するラスタ用メモリにエリアＥｉのＸ−Ｙ平面を定義し、このＸ方向分割平均距離Ｂｄｉ（例えば、０．５５９ｍ）でエリアＥｉのＸ軸、Ｙ軸を分割する。この分割格子を微細格子ｍｉと称している。
【００３２】
なお、図示しないが表示用メモリの解像度をＸ方向分割平均距離Ｂｄｉ（例えば、０．５５９ｍ）にする（セルサイズは、０．５５９ｍに対応するサイズ）。
そして、ラスタ用メモリに平面直角５ｍメッシュＭｂｉを順次定義して、各々の平面直角５ｍメッシュＭｂｉに対してＴＩＮ（triangulated irregular network）バイナリ補間を行って、ラスタ用メモリの微細格子ｍｉの標高値を補間する。
【００３３】
そして、Ｘ方向分割平均距離Ｂｄｉ（例えば、０．５５９ｍ）と同じサイズの格子を横９列、縦９行にした移動平均化用メッシュＨＭｉ（相加平均ともい）を生成して、順次、ラスタ用メモリの微細格子ｍｉにかけて滑か処理を行う（Ｓ６０）。
【００３４】
これらは、オペレータが滑らかと判断されるまで、移動平均化用メッシュＨＭｉ（９×９：スムージングフィルタともいう）をかけて行うのが好ましい。滑らかにされたラスタ用メモリの微細格子ｍｉに割り付けられている標高値を滑か処理後標高値ｚｒｉと称する。
【００３５】
そして、ラスタ用メモリの微細格子ｍｉに割り付けられている滑か処理後標高値ｚｒｉと同じ滑か処理後標高値ｚｒｉを有する微細格子ｍｉを直線で繋げ、これらをベクターに変換し、これを滑か等高線情報Ｊｉとして生成する（Ｓ７０）。この滑か等高線情報Ｊｉを画像化したものを滑か等高線画像ＧＪｉと称する（図３参照）。
【００３６】
次に、地理院の標準地図情報ＳＤｉ（レベル１６：ベクター）を読み込み、この標準地図情報ＳＤｉ（２万５０００分の１）を画像化し（以下、標準地図画像ＧＯｉという）、滑か等高線画像ＧＪｉと乗算合成して標準地図画像ＧＯｉ＋滑か等高線画像ＧＪｉ（総称して乗算合成画像ＧＳＪｉともいう）を生成（表示）する（Ｓ９０）。
なお、標準地図画像ＧＯｉ（図２及び図３参照）は、等高線、建物、道路等のベクターデータを画像化したものである。図３は図２の範囲ｕｃｉ付近を拡大した図である。
【００３７】
上記の処理によって得られた本実施の形態の図３の標準地図画像ＧＯｉ＋滑か等高線画像ＧＪｉ（乗算合成画像ＧＳＪｉ）と、図２に示す５ｍＤＥＭで生成した等高線図画像Ｇｏｉと比較すると、本実施の滑か等高線は各段に実際の地形の形状にあった等高線になっている（細かくなっている）。なお、図２には、後述する滑か等高線ｃｊｉを説明する範囲Ｕａｉを示し、図３にも図２の範囲Ｕａｉを示している。
【００３８】
図２と図３とを比較すると、従来の５ｍＤＥＭに基づいて生成した従来の等高線図画像Ｇｏｉは間隔が広い。これに対して図３の滑か等高線ｃｊｉは、間隔が狭く、かつ滑か等高線ｃｊｉ同士の曲率感が合っている。
【００３９】
（具体的な構成）
図４は本実施の形態の等高線生成システムの概略構成図である。図４に示すように、実施の形態１の等高線生成システム３００は、コンピュータ本体部１００と、表示部２００等で構成されている。
【００４０】
コンピュータ本体部１００は、地理座標（緯度経度）で定義された５ｍＤＥＭ基盤地図Ｆａを記憶した基盤地図用データベース１１０と、５ｍＤＥＭメッシュ読込部１１２と、地理座標用ＸＹＺポイントファイル１１４（メモリ）と、平面直角用ＸＹＺポイントファイル１１８（メモリ）と、ラスター化処理部１３０と、滑か処理部１４０と、等高線生成部１４５と、重ね合せ部１４３と、解像度定義部１５２と、表示処理部１５０等を備えて、図３に示す滑か等高線画像ＧＪｉを表示部２００の画面に得る。
【００４１】
また、記憶手段として基盤地図用データベース１１０と、地理院の標準地図（レベル１６）を記憶した標準地図用データベース１４９と、第１のラスタ用メモリ１４２と、第２のラスタ用メモリ１４７等を備えている。第１のラスタ用メモリ１４２と第２のラスタ用メモリ１４７とを総称して単にラスタ用メモリという。
【００４２】
（各部の説明）
標準地図用データベース１４９は、２万５０００分の１の建物、河川、道路等（総称して都市図ともいう）のベクター情報（２万５０００分の１）を記憶している。
なお、図５は、ベクター情報をラスター化した画像を標準地図画像ＧＯｉとして示している。但し、図５の標準地図画像ＧＯｉは等高線も示している。本実施の形態ではこの等高線（ベクター）は用いない。但し、図５は図２のＵａｉの付近である。
【００４３】
基盤地図用データベース１１０は、緯度経度５ｍＤＥＭ（地形）を記憶している。この緯度経度５ｍＤＥＭは、航空レーザ測量で計測して得た点群（例えば、５０ｃｍ間隔：ｘａｉ、ｙａｉ、ｚａｉ）に対して経度0.2秒×緯度0.2秒で区切ったメッシュ（以下、緯度経度５ｍメッシュＭａｉという）で定義している。この、緯度経度５ｍメッシュＭａｉの角の点を５ｍＤＥＭ緯度経度ポイントデータＰａｉ（緯度ｘａｉ、緯度ｙａｉ及び標高値ｚａｉ）と称している（図６、図７参照）。
【００４４】
図６においては、緯度軸Ｎ、経度軸Ｅを０．１秒間隔で区切った枠を示し、ＴＩＮ補間用の斜線を示す。また、サイズスケールＱｉ（０．０ｍ〜６．０ｍ）と、等高線の標高を知らせる色バーＶｉを示す。さらに、図６には理解を容易にするために、背景として後述する等高線の段彩画像データを緯度、経度座標に変換して示している。
【００４５】
図７は、５ｍＤＥＭ緯度経度ポイントデータＰａｉの具体例である。但し、標高値は示さない。
【００４６】
５ｍＤＥＭメッシュ読込部１１２は、オペレータにより入力（指定）されたエリアＥｉ（例えば、石巻市）の範囲に対応する全ての緯度経度５ｍメッシュＭａｉを構成する５ｍＤＥＭ緯度経度ポイントデータＰａｉを５ｍＤＥＭの基盤地図用データベース１１０より地理座標用ＸＹＺポイントファイル１１４（メモリ）にエクスポートする（図６及び図７参照）。図６においては、理解を容易にするために、背景として後述する等高線ｃｊｉの段彩図（ラスタ画像）を経度座標に変換して示している。
【００４７】
平面直角座標変換部１１６は、地理座標用ＸＹＺポイントファイル１１４（メモリ）の５ｍＤＥＭ緯度経度ポイントデータＰａｉを平面直角座標に投影変換し（高さはそのまま用いる）、平面直角用ＸＹＺポイントファイル１１８（メモリ）にエクスポートする（図８参照）。
【００４８】
これを５ｍＤＥＭ平面直角ポイントデータＰｂｉ（ｘｂｉ、ｙｂｉ、ｚｂｉ（ｚａｉ））と称する。
【００４９】
図８においては、理解を容易にするために、背景として後述する等高線ｃｊｉの段彩図）を平面直角座標に変換して示している。
図８に示すように、緯度経度５ｍメッシュＭａｉの５ｍＤＥＭ緯度経度ポイントデータＰａｉ（緯度ｘａｉ、緯度ｙａｉ及び標高値ｚａｉ）に対応する４点の５ｍＤＥＭ平面直角ポイントデータＰｂｉ（ｘｂｉ、ｙｂｉ、ｚｂｉ（ｚａｉ））を繋ぐメッシュ（以下、平面直角５ｍメッシュＭｂｉという）は、Ｘ方向が狭く、かつＹ方向が長い長方形になっている。
【００５０】
また、図８に示すように、緯度経度５ｍメッシュＭａｉは平面直角座標に変換されるので、等高線の段彩図の画像もＸ方向が短くなり、Ｙ方向が長くなる。
この理由を説明する。
【００５１】
５ｍＤＥＭ緯度経度ポイントデータＰａｉは、図９（ａ）に示すように、地球中心「Ｏ」から緯度経度座標に定義した緯度経度５ｍメッシュＭａｉ（Ｍａ１、Ｍａ２、・・）ものである。
【００５２】
そして、平面直角座標は、経緯線を円筒に投影してから円筒を開いて生成した「正角円筒図法」であり、極に近づくほど緯線の間隔が広くなる。具体的には、エリアＥｉによって、Ｘ方向、Ｙ方向が歪んでおり、例えば図９（ｂ）の平面直角５ｍメッシュＭｂｉは略平行四辺形状になるためである。
【００５３】
ラスター化処理部１３０は、微細格子生成処理部１３２と、ＴＩＮバイナリ補間部１３９等を備えている。
微細格子生成処理部１３２は、Ｘ方向９分割距離算出部１３２と、メッシュ分割部１３４等を備えて、第１のラスタ用メモリ１４２のＸ軸（Ｘ方向ともいう）、Ｙ軸（Ｙ方向ともいう）を例えば０．５５９ｍの間隔で奇数分割してＸ−Ｙ平面に０．５５９ｍのサイズの微細格子ｍｉ（ｍ１、ｍ２、・・・）を生成する（図１０参照）。
【００５４】
Ｘ方向９分割距離算出部１３２は、平面直角用ＸＹＺポイントファイル１１８（メモリ）の４点の５ｍＤＥＭ平面直角ポイントデータＰｂｉ（ｘｂｉ、ｙｂｉ、ｚｂｉ（ｚａｉ））を繋ぐ平面直角５ｍメッシュＭｂｉ（Ｐｂ１、Ｐｂ２、Ｐｂ３、Ｐｂ４、・・・）を順次指定する。
【００５５】
そして、この指定毎に平面直角５ｍメッシュＭｂｉ（Ｐｂ１、Ｐｂ２、Ｐｂ３、Ｐｂ４、・・・）のＸ方向の２点の５ｍＤＥＭ平面直角ポイントデータＰｂｉ（例えば、Ｐｂｉ、Ｐｂｉ＋１）との間の距離Ｄｘｉ（Ｄｘ１、Ｄｘ２、・・・）を順次求め、この距離Ｄｘｉ毎に９分割（メッシュの大きさによっては、３分割、５分割、７分割、１１分割、１３分割、・・・奇数分割）する。これらをＸ方向分割距離Ｂｉと称する。本実施の形態では９分割とする。
【００５６】
そして、これらのＸ方向分割距離Ｂｉを平均化する（以下、Ｘ方向分割平均距離Ｂｄｉという）。本実施の形態では、Ｘ方向分割平均距離Ｂｄｉを０．５５９ｍとする。
一般に格子間隔を適切にとらないと、平面直角座標系に対する緯度、経度の投影変換と、内挿補間の干渉による、格子状のアーチファクトが生じする。
【００５７】
このため、平面直角５ｍメッシュＭｂｉのＸ方向を９分割している。この分割は、できるだけ正方形又は長方形のメッシュになるように分割するのが好ましい。なお、平面直角５ｍメッシュＭｂｉは、９分割にした場合は、画像を画面に表示させたときに、格子状のアーチファクトが最も生じないことが実験によって判明したことによる。
【００５８】
そして、Ｘ方向９分割距離算出部１３２は、このＸ方向分割平均距離Ｂｄｉをメッシュ分割部１３４及び滑か処理部１４０に出力すると共に、ＴＩＮバイナリ補間部１３９に、指定された平面直角５ｍメッシュＭｂｉ（Ｐｂ１、Ｐｂ２、Ｐｂ３、Ｐｂ４、・・・）の番号（単に、平面直角５ｍメッシュＭｂｉと記載する）と、Ｘ方向分割平均距離Ｂｄｉと、平面直角５ｍメッシュＭｂｉを構成する４点の５ｍＤＥＭ平面直角ポイントデータＰｂｉ（ｘｂｉ、ｙｂｉ、ｚｂｉ（ｚａｉ））等を出力する。これらを総称してＴＩＮ補間用データＦｉと称する。
【００５９】
メッシュ分割部１３４は、ラスタ用メモリ１４２にエリアＥｉの範囲のＸ−Ｙ平面を定義する。そして、このＸ−Ｙ平面の基準位置（例えば、Ｘ軸）からＸ方向に向かってＸ方向分割平均距離Ｂｄｉ（０．５５９ｍ）で分割すると共に、Ｙ方向をＸ方向分割平均距離Ｂｄｉ（０．５５９ｍ）で分割する微細格子生成処理を行う（図１０参照）。このＸ−Ｙ平面に定義された格子を微細格子ｍｉ（ｍ１、ｍ２、ｍ３、・・・）と称している。
【００６０】
図１１は分割後の微細格子ｍｉ（ｍ１、ｍ２、ｍ３、・・・）の三角形の補間後の座標例である。図１１は、。おおむね9×11で分割した場合の例である。
【００６１】
Idx,X,Y,Elevation (m),Length,Total Length,Heading
1,-10835.893,-32871.056,41.274,0.559 m,---,269° 55' 48.4"
2,-10836.452,-32871.056,41.412,0.79 m,0.559 m,134° 52' 44.3"
3,-10835.893,-32871.614,41.214,---,1.349 m,---
となる。
【００６２】
ＴＩＮバイナリ補間部１３９（単に標高値補間部ともいう）は、ラスタ用メモリの微細格子毎の標高値を内挿補間して割付ける。
すなわち、ＴＩＮ補間用データＦｉが入力する毎に、これに含まれている平面直角５ｍメッシュＭｂｉを構成する４点の５ｍＤＥＭ平面直角ポイントデータＰｂｉ（ｘｂｉ、ｙｂｉ、ｚｂｉ（ｚａｉ））を第１のラスタ用メモリ１４２に順次定義してこれらを結ぶ（図１２参照）。
【００６３】
つまり、Ｘ方向、Ｙ方向が０．５５９ｍ間隔で微細に分割された第１のラスタ用メモリ１４２のＸ−Ｙ平面に、平面直角５ｍメッシュＭｂｉ（Ｍｂ１、Ｍｂ２、・・・）を定義する。
【００６４】
そして、この平面直角５ｍメッシュＭｂｉ毎にＴＩＮを定義して、このＴＩＮに含まれる微細格子ｍｉ（ｍ１、ｍ２、ｍ３、・・・）の標高値ｚｂｉを補間する（図１２参照）。これを補間標高値ｈｚｂｉと称する。
【００６５】
また、図１２は、Ｘ方向を９等分に分割し、Ｙ方向をおおむね１１分割して、背景を等高線のラスタ画像にしている。補間標高値ｈｚｂｉは黒点で示している。
そして、ＴＩＮバイナリ補間部１３９は、さらにＴＩＮバイナリ補間処理により順次補間する。この補間されたデータをＴＩＮバイナリ補間標高値ｚｒｉ（ｚｒ１、ｚｒ２、・・）と称する。そして、滑か処理部１４０（スムージング処理ともいう）を起動する。
【００６６】
平面直角５ｍメッシュＭｂｉは、歪んでいるので、図１３に示すように、Ｘ、Ｙ方向がずれる。ＴＩＮバイナリ補間処理による補間値の決定は、平面直角５ｍメッシュＭｂｉに隣接する平面直角５ｍメッシュＭｂｉの内で、微細格子ｍｉ（ｍ１、ｍ２、ｍ３、・・・）の面積が大きい方の平面直角５ｍメッシュＭｂｉの補間値を採用するのが好ましい。
【００６７】
図１４はＴＩＮバイナリ補間後の等高線のラスタ画像の例であり、第１のラスタ用メモリ１４２のＴＩＮバイナリ補間標高値ｚｒｉを色別で示した例である。
すなわち、ラスター化処理部１３０は、図１５（ａ）に示す平面直角５ｍメッシュＭｂｉで代表される標高値を０．５５９ｍ（小数点４位以下省略）単位の微細間隔の標高値にして、これをバイナリ補間している。
【００６８】
しかし、バイナリ補間は、図１５（ｂ）に示すように、縁が急激に吐出（ｈｉ）したり、あるは谷が急激に下がったり（ｈｉ）する。このため、本実施の形態では、滑か処理部１４０を備えている。
【００６９】
滑か処理部１４０は、移動平均化処理を行う。この移動平均化処理は、オペレータによって滑か処理指示が入力する毎に、Ｘ方向分割平均距離Ｂｄｉ（例えば、０．５５９ｍ）を読み込む。そして、縦横のサイズをこのＸ方向分割平均距離Ｂｄｉ（例えば、０．５５９ｍ）にした格子を、横列に９個、縦に９個にした図１６に示す移動平均化用メッシュＨＭｉを生成する。
【００７０】
そして、第２のラスタ用メモリ１４７（等高線用メモリともいう）に、第１のラスタ用メモリ１４２と同じエリアＥｉのＸ−Ｙ平面を定義し、このＸ−Ｙ平面のＸ方向、Ｙ方向をＸ方向分割平均距離Ｂｄｉ（例えば、０．５５９ｍ）で分割する上記の微細格子生成処理を行う。
【００７１】
微細格子生成処理によって生成された第２のラスタ用メモリ１４７（等高線用メモリともいう）のレイヤーをて生成された等高線用レイヤと称する。
一回目の滑か処理指示に伴って生成された等高線用レイヤを第１の等高線用レイヤーと称し、二回目の滑か処理指示に伴って生成された等高線用レイヤーを第１の等高線用レイヤーと称する。これらを総称して等高線用レイヤーと称する。
【００７２】
そして、滑か処理部１４０は、この移動平均化用メッシュＨＭｉの中心格子を第１のラスタ用メモリ１４２の微細メッシュｍｉ（ｍ１、ｍ２、ｍ３、・・・）に順次定義し、この中心格子における移動平均値（以下、滑か処理後標高値ｚｒｉ）等を第２のラスタ用メモリ１４７の等高線用レイヤーに順次記憶する（以下、移動平均化処理（加重平均）という）。
【００７３】
滑か処理後標高値ｚｒｉ（ｚｒ１、ｚｒ２、・・・）には各種データが関連付けられる（以下、滑か微細標高値ラスタ画像データＲＧｉという（図２０参照））。
図１７は一回目の移動平均化処理による例であり、図１８は二回目の移動平均化処理である。図１９は二回目の移動平均化処理の拡大表示例である。図１７〜図１９に示すように、等高線の鋭い箇所が滑らかになっている。
【００７４】
つまり、図１５（ｂ）のｈｉが平滑される。
オペレータは、画面に表示された滑か微細標高値画像ＧＲｉの等高線の滑らか具合を確認して、必要により移動平均化処理（スムージングともいう）を指示している。この指示は、表示処理部１５０から滑か処理部１４０に出力される。
【００７５】
滑か微細標高値ラスタ画像データＲＧｉは、図２０に示すように、エリアＥｉと、平面直角５ｍメッシュＭｂｉと、平面直角５ｍメッシュＭｂｉと、平面直角５ｍメッシュＭｂｉの四角の座標（例えば、Ｐｂ１、Ｐｂ２、Ｐｂ３、Ｐｂ４）と、微細格子ｍｉ（番号）と、Ｘ方向分割平均距離Ｂｄｉ（例えば、０．５５９ｍ）と、微細格子ｍｉの四角の座標（例えば、ＰＬ１、ＰＬ２、ＰＬ３、ＰＬ４：図２１に示している）と、補間標高値ｈｚｂｉと、ＴＩＮバイナリ補間標高値ｚｒｉと、一回目の滑か微細標高値Ｒｉと、２回目の滑か微細標高値Ｒｉ´等よりなる。例えば、ＰＬ１とＰＬ２とが入り口の線、ＰＬ３とＰＬ４とが出口の線である。
【００７６】
本実施の形態では一回目の滑か微細標高値Ｒｉとして説明する（以下、単に滑か微細標高値Ｒｉと称する）。
等高線生成部１４５について図２１を用いて説明する。図２１においては、同じ滑か微細標高値Ｒｉの微細格子ｍｉは１個として説明する。
【００７７】
等高線生成部１４５は、等高線用メモリ１８０にエリアＥｉの全ての滑か微細標高値ラスタ画像データＲＧｉが記憶さる毎に、このエリアＥｉの基準点（既知となる点）に該当する微細格子ｍｉ（滑か微細標高値ラスタ画像データＲＧｉ）を指定する。
【００７８】
そして、指定した微細格子ｍｉ（０．５５９ｍサイズ）毎に、一定範囲（例えば、５ｍ、１０ｍ、２０ｍ、・・）を定義し、指定した微細格子ｍｉ（滑か微細標高値ラスタ画像データＲＧｉ）の滑か微細標高値Ｒｉと同じ標高値を有する微細格子ｍｉ（微細標高値ラスタ画像データＲＧｉ）を検索する。
【００７９】
そして、これらの微細格子ｍｉ（微細標高値ラスタ画像データＲＧｉ）に対して標準偏差算出処理等により、繋ぐ微細格子ｍｉ（微細標高値ラスタ画像データＲＧｉ）を決定して行って、閉曲させる。
【００８０】
そして、この閉曲となる微細格子ｍｉ（微細標高値ラスタ画像データＲＧｉ）の集合をベクター化（関数）し、これを滑か等高線情報Ｊｉとする。
【００８１】
このベクター化は、繋ぐべき、隣の微細メッシュｍｉがＸ方向又はＹ方向の場合は、中心座標同士（ｘ、ｙ）を直線で結ぶ、また繋ぐ隣の微細メッシュｍｉが斜め方向の場合は、繋ぐ方向の微細メッシュｍｉ側の角２点を座標の中心と、繋ぐ斜め方向の微細メッシュｍｉの２点間の中心座標とを結んで直線にする。
【００８２】
そして、これらの直線の集合を関数（近似関数にしてもよい）にする。
すなわち、滑か等高線情報Ｊｉは、従来のように、スプライン曲線、ベジェ曲線等の曲率最大化処理を行わない、微細格子ｍｉ（０．５５９ｍ）を通る直線を繋げた等高線となっている。
【００８３】
このとき、色値を割り付ける（例えば、橙色）。すなわち、滑か等高線情報Ｊｉは、エリアＥｉと、微細格子ｍｉと、サイズ（０．５５９ｍ）と、滑か微細標高値Ｒｉと、接続方向（Ｘ方向上（又は下）、Ｙ方向上（又は下）又は右斜め或いは左斜め）等よりなる。
【００８４】
等高線生成部１４５は、このようにして生成した滑か等高線情報Ｊｉを等高線用メモリ１４８に記憶する。なお、滑か等高線ｃｊｉの間隔は、１ｍ、２ｍ、３ｍ、・・でも構わない。
【００８５】
重ね合せ部１４３は、メモリ１４２の微細格子ｍｉの格子サイズ（例えば、０．５５９ｍ）及びＸ−Ｙ平面の範囲（エリアＥｉ）を読み込み、メモリ１８０のＸ−Ｙ平面にエリアＥｉを定義して、このエリアＥｉのＸ方向及びＹ方向を０．５５９ｍで分割した、標準地図用レイヤと、等高線用レイヤと、標準地図等高線レイヤとを生成する。
そして、等高線用メモリ１４８から滑か等高線情報Ｊｉを読み込む。
【００８６】
そして、これをラスタ変換して等高線用レイヤに書き込む（滑か等高線画像ＣＪｉ）。また、標準地図用データベース１４９の標準地図情報ＳＤｉをラスタ変換して標準地図用レイヤに書き込む（標準地図画像ＧＯｉ）。
【００８７】
そして、標準地図画像ＧＯｉと滑か等高線画像ＣＪｉとを重ねた乗算合成画像ＧＳＪｉ（標準地図画像ＧＯｉ＋滑か等高線画像ＧＪｉ）を標準地図等高線レイヤに書き込む。
【００８８】
そして、この乗算合成画像ＧＳＪｉ（標準地図画像ＧＯｉ＋滑か等高線画像ＧＪｉ）を表示処理部１５０に出力する。
表示処理部１５０は、表示用メモリ（図示せず）を解像度定義部１５２によって設定された解像度（１セル：０．５５９ｍサイズ）にする。
【００８９】
また、等高線と等高線との間を、予め設定入力された色値の色で表示する（例えば、緑、黄色、赤、白、橙等）。色の表示は、高度差に応じて、色を濃くしたり、又は薄くしたりするのが好ましい。
解像度定義部１５２は、Ｘ方向９分割距離算出部１３２からのＸ方向分割平均距離Ｂｄｉ（例えば、０．５５９ｍ）を読込み、表示用メモリのセルサイズの縦横を０．５５９ｍに設定する。
【００９０】
そして、表示処理部１５０は、乗算合成画像ＧＳＪｉ（標準地図画像ＧＯｉ＋滑か等高線画像ＧＪｉ）を表示用メモリに書き込んで表示部２００の画面に乗算合成画像ＧＳＪｉ（標準地図画像ＧＯｉ＋滑か等高線画像ＧＪｉ）を表示する（図２８参照）。
【００９１】
図２８に示すように、等高線の解像度（線の太さ）と、標準地図の解像度（線の太さ）とが同程度となり、かつ隣接する等高線同士の曲率程度は違和感のない曲率となる。
また、５ｍＤＥＭを用いた０．５５９ｍサイズで示される細かな等高線（図３参照）であるので、５ｍＤＥＭ（１／２５０００対応）であっても１万分の１地図として使用できる。
【００９２】
なお、上記の滑か等高線ｃｊｉと滑か等高線ｃｊｉとの間の高低差に応じてその間を色表示（青、緑、黄色、赤、白等）しても構わない。
【００９３】
また、等高線は、陸地の他に、海底に適用させても構わない。また、上記実施の形態では、５ｍＤＥＭのメッシュとして説明したが１０ｍＤＥＭの場合は、Ｘ方向を９分割するのがよい。
【００９４】
次に、図２２〜図２５を用いて本実施の形態の滑らか滑か等高線画像ＧＪｉと、スムージング処理を行わなかった場合の滑か等高線画像ＧＪｉの違いを説明する。
図２２はＴＩＮバイナリ補間後の滑か等高線画像ＧＪｉを表示した例であり、図２３は図２２の範囲Ｕｂｉの拡大図である。図２４は図２２と同様なエリアであり、スムージング処理後の滑か等高線画像ＧＪｉを表示した例である。図２５は図３のＵａｉ付近の拡大図である。但し、図２５は、図３の範囲Ｕａｉに対応する領域である。
【００９５】
図２２に示すように、ＴＩＮバイナリ補間後の滑か等高線画像ＧＪｉは等高線の湾曲箇所は鋭く湾曲しているが、スムージング処理を施した場合は、図２５に湾曲部分は滑らかになっている。図２５は図２のＵｃｉのあたりである。図２５と図２とを比較すると、ＵｃｉのＵｃａｉに対応するエリアの図２５に示す等高線は、計曲線も共に滑らかになっており、その間隔も微細な間隔である。
【００９６】
すなわち、図２５に示すように、本実施の形態の等高線は１万分の１の等高線図として利用できる。
【００９７】
さらに、上記実施の形態では、滑か処理部１４０を備えて、第１のラスタ用メモリ１４２のＴＩＮバイナリ補間標高値ｚｒｉ（ｚｒ１、ｚｒ２、・・）を移動平均化処理したが、移動平均化処理をしない等高線を表示させても構わない。
すなわち、等高線生成部１４５が第１のラスタ用メモリ１４２のＴＩＮバイナリ補間標高値ｚｒｉ（ｚｒ１、ｚｒ２、・・）を等高線用メモリ１８０にエリアＥｉの全ての滑か微細標高値ラスタ画像データＲＧｉが記憶さる毎に、このエリアＥｉの基準点に該当する微細格子ｍｉ（滑か微細標高値ラスタ画像データＲＧｉ）を指定する。
【００９８】
そして、指定した微細格子ｍｉ（０．５５９ｍサイズ）毎に、一定範囲（例えば、５ｍ、１０ｍ、２０ｍ、・・）を定義し、指定した微細格子ｍｉ（滑か微細標高値ラスタ画像データＲＧｉ）の滑か微細標高値Ｒｉと同じ標高値を有する微細格子ｍｉ（微細標高値ラスタ画像データＲＧｉ）を検索する。
【００９９】
そして、これらの微細格子ｍｉ（微細標高値ラスタ画像データＲＧｉ）に対して標準偏差算出処理等により、繋ぐ微細格子ｍｉ（微細標高値ラスタ画像データＲＧｉ）を決定して行って、閉曲させる。
【０１００】
そして、この閉曲となる微細格子ｍｉ（微細標高値ラスタ画像データＲＧｉ）の集合をベクター化（関数）し、これを滑か等高線情報Ｊｉとし、これを表示処理部１５０に出力して画面に表示させる（図２２、図２３参照）。
【０１０１】
つまり、繋ぐべき、隣の微細メッシュｍｉがＸ方向又はＹ方向の場合は、中心座標同士（ｘ、ｙ）を直線で結ぶ、また繋ぐ隣の微細メッシュｍｉが斜め方向の場合は、繋ぐ方向の微細メッシュｍｉ側の角２点を座標の中心と、繋ぐ斜め方向の微細メッシュｍｉの２点間の中心座標とを結んで直線にしている。そして、これらの直線の集合を関数（近似関数にしてもよい）にする。
【０１０２】
前述の図２８は、特許第３６７０２７４号公報に記載の赤色立体地図（画像）の技術を用いた場合の赤色立体地図（本実施の形態では、超解像度赤色画像と称する）を重ねた乗算合成画像ＧＳＪｉとして表示させた例である。
【０１０３】
この超解像度赤色画像の生成の概略を下記に記載する。超解像度赤色画像の生成は、超解像度赤色立体視画像生成部（図示せず）が行う。超解像度赤色立体視画像生成部（以下、超解像度赤色立体画像作成処理という）。
超解像度赤色立体画像作成処理は、ラスタ用メモリ（第１のラスタ用メモリ１４２又は第２のラスタ用メモリ１４７）の微細格子ｍｉを指定し、この指定毎に、この指定した微細格子ｍｉを着目点とする。
【０１０４】
そして、この着目点から考慮距離Ｌ内に存在する滑か微細標高値を用いて斜度と地上開度と地下開度とを求め、地上開度と地下開度とから尾根谷度を求め、斜度を赤の彩度に割りあて、尾根谷度を明度に割り付けた平滑超解像度赤色立体視化画像を生成し、図示しない赤色立体画像記憶用メモリに記憶する。前述の考慮距離Ｌは、５ｍＤＥＭを用いた超解像度赤色立体画像作成処理が考慮距離Ｌを５０ｍとしている場合は、着目点から（５０ｍ／９）×９個＝Ｌとしている。なお、平滑超解像度赤色立体視化画像の色値は、５０％程度低下させている。
【０１０５】
そして、重ね合せ部１４３がこの平滑超解像度赤色立体視化画像と、標準地図用レイヤの標準地図画像ＧＯｉと、等高線用メモリ１４８の滑か等高線情報Ｊｉ（ベクタ）に基づく滑か等高線画像ＧＪｉとを乗算合成して図２８の乗算合成画像を得ている。
【０１０６】
図２６はスムージング処理の前の第１のラスタ用メモリ１４２における補間標高値ｈｚｂｉ（ｈｚｂ１、ｈｚｂ２、・・・）に対して超解像度赤色立体画像作成処理を行った例である。
【０１０７】
図２７はスムージング処理後の滑か処理後標高値ｚｒｉ（ｚｒ１、ｚｒ２、・・）に対して超解像度赤色立体画像作成処理を行った例である。但し、図２６のｍｉにおける黒点は、表示上のずれであり、実際は図２７に示すように、ｍｉの中央に位置している。
図２６と図２７とを比較すると、スムージング処理を施した方が全体的に色（橙）が均一化している。
【符号の説明】
【０１０８】
１１０  基盤地図用データベース
１１２  ５ｍＤＥＭ読込部
１１４  地理座標用ＸＹＺポイントファイル１１４
１１８  平面直角用ＸＹＺポイントファイル１１８
１３０  微細間隔標高値ラスタ化処理部
１４０  滑か処理部
１４５  等高線生成部
１４３  重ね合せ部
１５２  解像度定義部
１５０  表示処理部
【図１】
0006949926.tif 000002
【図２】
0006949926.tif 000003
【図３】
0006949926.tif 000004
【図４】
0006949926.tif 000005
【図５】
0006949926.tif 000006
【図６】
0006949926.tif 000007
【図７】
0006949926.tif 000008
【図８】
0006949926.tif 000009
【図９】
0006949926.tif 000010
【図１０】
0006949926.tif 000011
【図１１】
0006949926.tif 000012
【図１２】
0006949926.tif 000013
【図１３】
0006949926.tif 000014
【図１４】
0006949926.tif 000015
【図１５】
0006949926.tif 000016
【図１６】
0006949926.tif 000017
【図１７】
0006949926.tif 000018
【図１８】
0006949926.tif 000019
【図１９】
0006949926.tif 000020
【図２０】
0006949926.tif 000021
【図２１】
0006949926.tif 000022
【図２２】
0006949926.tif 000023
【図２３】
0006949926.tif 000024
【図２４】
0006949926.tif 000025
【図２５】
0006949926.tif 000026
【図２６】
0006949926.tif 000027
【図２７】
0006949926.tif 000028
【図２８】
0006949926.tif 000029