JP2021082112A - 等高線生成システム及びそのプログラム - Google Patents
等高線生成システム及びそのプログラム Download PDFInfo
- Publication number
- JP2021082112A JP2021082112A JP2019210350A JP2019210350A JP2021082112A JP 2021082112 A JP2021082112 A JP 2021082112A JP 2019210350 A JP2019210350 A JP 2019210350A JP 2019210350 A JP2019210350 A JP 2019210350A JP 2021082112 A JP2021082112 A JP 2021082112A
- Authority
- JP
- Japan
- Prior art keywords
- contour line
- mesh
- memory
- fine
- contour
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Instructional Devices (AREA)
- Processing Or Creating Images (AREA)
Abstract
Description
しかし、5mDEMのメッシュを用いた場合は、等高線間隔は、5m以下の詳細な等高線を得ることができない。
例えば、特許文献1の等高線図作成方法は、10m間隔の基本等高線図を画面に表示し、この基本等高線図上に縦横10m間隔で区切った区分線を定義する。
そして、この10mメッシュを複数の格子で分割している(図5においては、サイズが20cmの格子を縦5、横5に分割している)。
前記数値標高モデルに緯度経度で定義されている一定の大きさのメッシュを順次、平面直角座標系で定義された平面直角メッシュ用メモリに投影変換した平面直角メッシュ群を生成する平面直角座標変換部と、
平面直角座標系で定義されたラスター用メモリと、
前記平面直角メッシュ用メモリの前記平面直角メッシュ毎に、この平面直角メッシュのX軸方向となる横辺長を奇数分割し、これらの分割長を平均化した平均長で前記ラスター用メモリのX軸、Y軸を分割して前記ラスタ用メモリのX−Y平面に縦横が前記平均長の微細格子群を定義する微細格子生成処理部と、
前記ラスタ用メモリの微細格子毎の標高値を内挿補間して割付ける標高値補間部と、
前記ラスタ用メモリの全ての微細格子の内挿補間後の標高値の内で、起点となる前記微細格子を指定し、この指定した微細格子と同じ標高値を有して閉曲する微細格子群を通る直線を求め、これらの直線をベクター化し、これを等高線ベクターとして生成する等高線生成部と、
前記等高線ベクターを画像にして表示用メモリに書き込んで表示する表示処理部とを備えたことを要旨とする。
図1は本実施の形態の等高線生成システムの概念を説明する説明図であり、コンピュータが行う処理である。本実施の形態の等高線生成システムは、一般に広く提供されている国土地理院(以下、地理院という)の5mDEM(5mメッシュ)を微細分割し、この微細メッシュ(以下、微細格子miという)に割り付けられている標高値と同じ微細格子miを繋げたベクタ情報の等高線を生成する。この等高線は本実施の形態では滑か処理が施される(以下、滑か等高線cjiと称している)。
図1に示すように、任意のエリアEi(例えば、石巻市)の緯度及び経度等で定義された5mDEMのデータ(以下、緯度経度5mDEMデータpaiという)を読込む(S10)。
そして、ラスタ用メモリに平面直角5mメッシュMbiを順次定義して、各々の平面直角5mメッシュMbiに対してTIN(triangulated irregular network)バイナリ補間を行って、ラスタ用メモリの微細格子miの標高値を補間する。
なお、標準地図画像GOi(図2及び図3参照)は、等高線、建物、道路等のベクターデータを画像化したものである。図3は図2の範囲uci付近を拡大した図である。
図4は本実施の形態の等高線生成システムの概略構成図である。図4に示すように、実施の形態1の等高線生成システム300は、コンピュータ本体部100と、表示部200等で構成されている。
標準地図用データベース149は、2万5000分の1の建物、河川、道路等(総称して都市図ともいう)のベクター情報(2万5000分の1)を記憶している。
なお、図5は、ベクター情報をラスター化した画像を標準地図画像GOiとして示している。但し、図5の標準地図画像GOiは等高線も示している。本実施の形態ではこの等高線(ベクター)は用いない。但し、図5は図2のUaiの付近である。
図8に示すように、緯度経度5mメッシュMaiの5mDEM緯度経度ポイントデータPai(緯度xai、緯度yai及び標高値zai)に対応する4点の5mDEM平面直角ポイントデータPbi(xbi、ybi、zbi(zai))を繋ぐメッシュ(以下、平面直角5mメッシュMbiという)は、X方向が狭く、かつY方向が長い長方形になっている。
この理由を説明する。
微細格子生成処理部132は、X方向9分割距離算出部132と、メッシュ分割部134等を備えて、第1のラスタ用メモリ142のX軸(X方向ともいう)、Y軸(Y方向ともいう)を例えば0.559mの間隔で奇数分割してX−Y平面に0.559mのサイズの微細格子mi(m1、m2、・・・)を生成する(図10参照)。
一般に格子間隔を適切にとらないと、平面直角座標系に対する緯度、経度の投影変換と、内挿補間の干渉による、格子状のアーチファクトが生じする。
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,---
となる。
すなわち、TIN補間用データFiが入力する毎に、これに含まれている平面直角5mメッシュMbiを構成する4点の5mDEM平面直角ポイントデータPbi(xbi、ybi、zbi(zai))を第1のラスタ用メモリ142に順次定義してこれらを結ぶ(図12参照)。
そして、TINバイナリ補間部139は、さらにTINバイナリ補間処理により順次補間する。この補間されたデータをTINバイナリ補間標高値zri(zr1、zr2、・・)と称する。そして、滑か処理部140(スムージング処理ともいう)を起動する。
すなわち、ラスター化処理部130は、図15(a)に示す平面直角5mメッシュMbiで代表される標高値を0.559m(小数点4位以下省略)単位の微細間隔の標高値にして、これをバイナリ補間している。
一回目の滑か処理指示に伴って生成された等高線用レイヤを第1の等高線用レイヤーと称し、二回目の滑か処理指示に伴って生成された等高線用レイヤーを第1の等高線用レイヤーと称する。これらを総称して等高線用レイヤーと称する。
図17は一回目の移動平均化処理による例であり、図18は二回目の移動平均化処理である。図19は二回目の移動平均化処理の拡大表示例である。図17〜図19に示すように、等高線の鋭い箇所が滑らかになっている。
オペレータは、画面に表示された滑か微細標高値画像GRiの等高線の滑らか具合を確認して、必要により移動平均化処理(スムージングともいう)を指示している。この指示は、表示処理部150から滑か処理部140に出力される。
等高線生成部145について図21を用いて説明する。図21においては、同じ滑か微細標高値Riの微細格子miは1個として説明する。
すなわち、滑か等高線情報Jiは、従来のように、スプライン曲線、ベジェ曲線等の曲率最大化処理を行わない、微細格子mi(0.559m)を通る直線を繋げた等高線となっている。
そして、等高線用メモリ148から滑か等高線情報Jiを読み込む。
表示処理部150は、表示用メモリ(図示せず)を解像度定義部152によって設定された解像度(1セル:0.559mサイズ)にする。
解像度定義部152は、X方向9分割距離算出部132からのX方向分割平均距離Bdi(例えば、0.559m)を読込み、表示用メモリのセルサイズの縦横を0.559mに設定する。
また、5mDEMを用いた0.559mサイズで示される細かな等高線(図3参照)であるので、5mDEM(1/25000対応)であっても1万分の1地図として使用できる。
図22はTINバイナリ補間後の滑か等高線画像GJiを表示した例であり、図23は図22の範囲Ubiの拡大図である。図24は図22と同様なエリアであり、スムージング処理後の滑か等高線画像GJiを表示した例である。図25は図3のUai付近の拡大図である。但し、図25は、図3の範囲Uaiに対応する領域である。
すなわち、等高線生成部145が第1のラスタ用メモリ142のTINバイナリ補間標高値zri(zr1、zr2、・・)を等高線用メモリ180にエリアEiの全ての滑か微細標高値ラスタ画像データRGiが記憶さる毎に、このエリアEiの基準点に該当する微細格子mi(滑か微細標高値ラスタ画像データRGi)を指定する。
超解像度赤色立体画像作成処理は、ラスタ用メモリ(第1のラスタ用メモリ142又は第2のラスタ用メモリ147)の微細格子miを指定し、この指定毎に、この指定した微細格子miを着目点とする。
図26と図27とを比較すると、スムージング処理を施した方が全体的に色(橙)が均一化している。
112 5mDEM読込部
114 地理座標用XYZポイントファイル114
118 平面直角用XYZポイントファイル118
130 微細間隔標高値ラスタ化処理部
140 滑か処理部
145 等高線生成部
143 重ね合せ部
152 解像度定義部
150 表示処理部
前記数値標高モデルに緯度経度で定義されている一定の大きさのメッシュを順次、平面直角座標系で定義された平面直角メッシュ用メモリに投影変換した平面直角メッシュ群を生成する平面直角座標変換部と、
平面直角座標系で定義されたラスタ用メモリと、
前記平面直角メッシュ用メモリの前記平面直角メッシュ毎に、この平面直角メッシュのX軸方向となる横辺長を奇数分割し、これらの分割長を平均化した平均長で前記ラスタ用メモリのX軸、Y軸を分割して前記ラスタ用メモリのX−Y平面に縦横が前記平均長の微細格子群を定義する微細格子生成処理部と、
前記ラスタ用メモリの微細格子毎の標高値を内挿補間して割付ける標高値補間部と、
前記ラスタ用メモリの全ての微細格子の内挿補間後の標高値の内で、起点となる前記微細格子を指定し、この指定した微細格子と同じ標高値を有して閉曲する微細格子群を通る直線を求め、これらの直線をベクター化し、これを等高線ベクターとして生成する等高線生成部と、
前記等高線ベクターを画像にして表示用メモリに書き込んで画面に表示する表示処理部とを備えることを要旨とする。
Claims (12)
- 緯度経度で定義された数値標高モデルを記憶した数値標高モデル用メモリと、
前記数値標高モデルに緯度経度で定義されている一定の大きさのメッシュを順次、平面直角座標系で定義された平面直角メッシュ用メモリに投影変換した平面直角メッシュ群を生成する平面直角座標変換部と、
平面直角座標系で定義されたラスター用メモリと、
前記平面直角メッシュ用メモリの前記平面直角メッシュ毎に、この平面直角メッシュのX軸方向となる横辺長を奇数分割し、これらの分割長を平均化した平均長で前記ラスター用メモリのX軸、Y軸を分割して前記ラスタ用メモリのX−Y平面に縦横が前記平均長の微細格子群を定義する微細格子生成処理部と、
前記ラスタ用メモリの微細格子毎の標高値を内挿補間して割付ける標高値補間部と、
前記ラスタ用メモリの全ての微細格子の内挿補間後の標高値の内で、起点となる前記微細格子を指定し、この指定した微細格子と同じ標高値を有して閉曲する微細格子群を通る直線を求め、これらの直線をベクター化し、これを等高線ベクターとして生成する等高線生成部と、
前記等高線ベクターを画像にして表示用メモリに書き込んで表示する表示処理部と
を有することを特徴とする等高線生成システム。 - 前記微細格子を縦列及び横列ともに前記奇数分割の個数にした移動平均化メッシュを生成し、この移動平均化メッシュを前記ラスター用メモリの前記微細格子に順次、かけて内挿補間後の標高値を順次、平滑する滑か処理部と
を有することを特徴とする請求項1記載の等高線生成システム。 - 道路、建物、河川、沼がベクター情報で定義された2万5千分の1の標準地図情報が記憶された標準地図用メモリを備え、
前記標準地図情報を画像化し、前記等高線情報を画像に重ねた合成画像を生成して表示部に表示させる重ね合せ部と
を有することを特徴とする請求項1又は2記載の等高線生成システム。 - 前記平均長のサイズを前記表示用メモリのセルの解像度として定義する解像度定義部と
を有することを特徴とする請求項1記載の等高線生成システム。 - 前記数値標高モデルは、5mDEMであり、前記奇数分割の数は、「9」であることを特徴とする請求項1記載の等高線生成システム。
- 前記内挿補間は、TINバイリニア補間であることを特徴とする請求項1記載の等高線生成システム。
- 前記数値標高モデル用メモリから所定のエリアに存在する前記一定の大きさのメッシュを構成する点群を読み出して前記平面直角座標変換部に出力するエリア読込部と
を有することを特徴とする請求項1記載の等高線生成システム。 - コンピュータを、
緯度経度で定義された数値標高モデルを数値標高モデル用メモリに記憶する手段、
前記数値標高モデルに緯度経度で定義されている一定の大きさのメッシュを順次、平面直角座標系で定義された平面直角メッシュ用メモリに投影変換した平面直角メッシュ群を生成する手段、
ラスター用メモリを平面直角座標系で定義する手段、
前記平面直角メッシュ用メモリの前記平面直角メッシュ毎に、この平面直角メッシュのX軸方向となる横辺長を奇数分割し、これらの分割長を平均化した平均長で前記ラスター用メモリのX軸、Y軸を分割して前記ラスタ用メモリのX−Y平面に縦横が前記平均長の微細格子群を定義する手段、
前記ラスタ用メモリの微細格子毎の標高値を内挿補間して割付ける手段、
前記ラスタ用メモリの全ての微細格子の内挿補間後の標高値の内で、起点となる前記微細格子を指定し、この指定した微細格子と同じ標高値を有して閉曲する微細格子群を通る直線を求め、これらの直線をベクター化し、これを等高線ベクターとして生成する手段、
前記等高線ベクターを画像にして表示用メモリに書き込んで表示する手段
として機能させるための等高線生成プログラム。 - コンピュータを、
前記微細格子を縦列及び横列ともに前記奇数分割の個数にした移動平均化メッシュを生成し、この移動平均化メッシュを前記ラスター用メモリの前記微細格子に順次、かけて内挿補間後の標高値を順次、平滑する手段
としての機能させるための請求項8記載の等高線生成プログラム。 - コンピュータに、
道路、建物、河川、沼がベクター情報で定義された2万5千分の1の標準地図情報を標準地図用メモリに記憶する手段、
前記標準地図情報を画像化し、前記等高線情報を画像に重ねた合成画像を生成して表示部に表示させる手段
としての機能を実行させる請求項8記載の等高線生成プログラム。 - コンピュータを、
前記平均長のサイズを前記表示用メモリのセルの解像度として定義する手段
としての機能を実行させるための請求項8記載の等高線生成プログラム。 - 前記数値標高モデルは、5mDEMであり、前記奇数分割の数は、「9」であることを特徴とする請求項8記載の等高線生成プログラム。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2019210350A JP6949926B2 (ja) | 2019-11-21 | 2019-11-21 | 等高線生成システム及びそのプログラム |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2019210350A JP6949926B2 (ja) | 2019-11-21 | 2019-11-21 | 等高線生成システム及びそのプログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2021082112A true JP2021082112A (ja) | 2021-05-27 |
JP6949926B2 JP6949926B2 (ja) | 2021-10-13 |
Family
ID=75965324
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2019210350A Active JP6949926B2 (ja) | 2019-11-21 | 2019-11-21 | 等高線生成システム及びそのプログラム |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP6949926B2 (ja) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2021112678A (ja) * | 2016-09-12 | 2021-08-05 | 株式会社三洋物産 | 遊技機 |
CN117876621A (zh) * | 2024-03-07 | 2024-04-12 | 贵州省第一测绘院(贵州省北斗导航位置服务中心) | 一种基于高分辨率遥感图像和地形数据的国土测绘方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2003090184A1 (fr) * | 2002-04-22 | 2003-10-30 | Dgs Computer | Procede et dispositif permettant le dessin d'une carte altimetrique numerique |
WO2004042675A1 (ja) * | 2002-11-05 | 2004-05-21 | Asia Air Survey Co.,Ltd. | 視覚化処理システム、視覚化処理方法、及び視覚化処理プログラム |
-
2019
- 2019-11-21 JP JP2019210350A patent/JP6949926B2/ja active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2003090184A1 (fr) * | 2002-04-22 | 2003-10-30 | Dgs Computer | Procede et dispositif permettant le dessin d'une carte altimetrique numerique |
WO2004042675A1 (ja) * | 2002-11-05 | 2004-05-21 | Asia Air Survey Co.,Ltd. | 視覚化処理システム、視覚化処理方法、及び視覚化処理プログラム |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2021112678A (ja) * | 2016-09-12 | 2021-08-05 | 株式会社三洋物産 | 遊技機 |
CN117876621A (zh) * | 2024-03-07 | 2024-04-12 | 贵州省第一测绘院(贵州省北斗导航位置服务中心) | 一种基于高分辨率遥感图像和地形数据的国土测绘方法 |
Also Published As
Publication number | Publication date |
---|---|
JP6949926B2 (ja) | 2021-10-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105336003B (zh) | 结合gpu技术实时流畅绘制出三维地形模型的方法 | |
US7352369B2 (en) | System and method for approximating an editable surface | |
TWI780945B (zh) | 使用有超解析度影像之等高線產生系統及使用有超解析度影像之等高線產生程式 | |
JP6949926B2 (ja) | 等高線生成システム及びそのプログラム | |
JPS58198712A (ja) | 産業プラント構造物に関する情報を引き出し、蓄積する方法 | |
CN112465976A (zh) | 堆场三维地图建立方法、库存管理方法、设备及介质 | |
AU2018381378A1 (en) | Feature/Ground Height-Based Colored Image Generating Apparatus and Feature Height-Based Colored Image Generating Program | |
KR20200096950A (ko) | 지물 높이별 컬러 화상 생성 장치 및 지물 높이별 컬러 화상 생성 프로그램 | |
CN112102489B (zh) | 导航界面显示方法、装置、计算设备和存储介质 | |
Yoo et al. | Image‐Based Modeling of Urban Buildings Using Aerial Photographs and Digital Maps | |
US6856312B1 (en) | Method for generating three-dimensional form data and apparatus therefor | |
JP2002236019A (ja) | 地表面抽出処理システム | |
JP5241296B2 (ja) | 数値地図データ処理プログラム及び数値地図データ処理装置 | |
CN117274527B (zh) | 一种发电机设备三维可视化模型数据集构建方法 | |
JP4592510B2 (ja) | 立体地図画像生成装置および方法 | |
JP2012073520A (ja) | 立体視画像表示処理装置、立体視画像表示処理方法及びプログラム | |
Frommholz et al. | Inlining 3d reconstruction, multi-source texture mapping and semantic analysis using oblique aerial imagery | |
JPH0636013A (ja) | 地形データの作成方法および装置 | |
JP4543820B2 (ja) | 3次元データ処理装置及びプログラム | |
US20090146998A1 (en) | Computer graphics using coarse level meshes | |
CN113282695B (zh) | 一种基于遥感影像的矢量地理信息采集方法和装置 | |
JP5946369B2 (ja) | 3次元地図画像データ生成システム | |
EP0473152B1 (en) | Topographical data construction system | |
Böhm | From point samples to surfaces-on meshing and alternatives | |
CN106600694A (zh) | 一种地形数据的平滑处理方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A821 Effective date: 20191122 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20210128 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20210128 |
|
A871 | Explanation of circumstances concerning accelerated examination |
Free format text: JAPANESE INTERMEDIATE CODE: A871 Effective date: 20210128 |
|
A975 | Report on accelerated examination |
Free format text: JAPANESE INTERMEDIATE CODE: A971005 Effective date: 20210419 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20210422 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20210427 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20210701 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20210706 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20210916 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20210922 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6949926 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |