JP3535598B2 - 図形データ整合方法 - Google Patents
図形データ整合方法Info
- Publication number
- JP3535598B2 JP3535598B2 JP33927594A JP33927594A JP3535598B2 JP 3535598 B2 JP3535598 B2 JP 3535598B2 JP 33927594 A JP33927594 A JP 33927594A JP 33927594 A JP33927594 A JP 33927594A JP 3535598 B2 JP3535598 B2 JP 3535598B2
- Authority
- JP
- Japan
- Prior art keywords
- map
- data
- mesh
- affine transformation
- map data
- 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.)
- Expired - Lifetime
Links
Landscapes
- Processing Or Creating Images (AREA)
Description
【発明の詳細な説明】
【０００１】
【産業上の利用分野】本発明は、図形データ整合方法に
関するものである。
関するものである。
【０００２】
【従来の技術】一般に、マッピングシステムでは、たと
えば、道路管理システムの地図データを都市計画図に整
合させるような場合が生じる。
えば、道路管理システムの地図データを都市計画図に整
合させるような場合が生じる。
【０００３】
【発明が解決しようとする課題】しかしながら、道路管
理システムの地図はメッシュ化されており、この地図デ
ータを都市計画図に整合させる場合、通常各メッシュ毎
にアフィン変換パラメータを求めてアフィン変換を行う
ので、メッシュの境界においては１本の道路が一致しな
いような事態が生じることがある。
理システムの地図はメッシュ化されており、この地図デ
ータを都市計画図に整合させる場合、通常各メッシュ毎
にアフィン変換パラメータを求めてアフィン変換を行う
ので、メッシュの境界においては１本の道路が一致しな
いような事態が生じることがある。
【０００４】本発明は、このような問題に鑑みてなされ
たもので、その目的とするところは、メッシュ化された
地図データをメッシュ毎にアフィン変換を行っても、各
メッシュの境界で地図データが整合される図形データ整
合方法を提供することにある
たもので、その目的とするところは、メッシュ化された
地図データをメッシュ毎にアフィン変換を行っても、各
メッシュの境界で地図データが整合される図形データ整
合方法を提供することにある
【０００５】。
【課題を解決するための手段】前述した目的を達成する
ために、本発明は、メッシュ化された第１の地図データ
を変換して第２の地図データに整合させる図形データ整
合方法であって、（ａ）第１の地図データと第２の地図
データの対応する点を複数個指定する工程と、（ｂ）前
記対応する点を基にして、アフィン変換用パラメータを
算出する工程と、（ｃ）前記第１の地図データの各メッ
シュの４隅の図郭を前記アフィン変換用パラメータを用
いてアフィン変換する工程と、（ｄ）あるメッシュの１
つの図郭は、隣接する他の３つのメッシュの図郭でもあ
り、隣接する４つのメッシュの同一点を表す４つの図郭
が、前記工程（ｃ）により夫々アフィン変換され、変換
後の４つの図郭が一致しない場合には代表点を算出する
工程と、（ｅ）前記第１の地図データの図郭が前記代表
点と一致するように、前記第１の地図データを逆正規化
変換する工程と、を具備することを特徴とする図形デー
タ整合方法である。
ために、本発明は、メッシュ化された第１の地図データ
を変換して第２の地図データに整合させる図形データ整
合方法であって、（ａ）第１の地図データと第２の地図
データの対応する点を複数個指定する工程と、（ｂ）前
記対応する点を基にして、アフィン変換用パラメータを
算出する工程と、（ｃ）前記第１の地図データの各メッ
シュの４隅の図郭を前記アフィン変換用パラメータを用
いてアフィン変換する工程と、（ｄ）あるメッシュの１
つの図郭は、隣接する他の３つのメッシュの図郭でもあ
り、隣接する４つのメッシュの同一点を表す４つの図郭
が、前記工程（ｃ）により夫々アフィン変換され、変換
後の４つの図郭が一致しない場合には代表点を算出する
工程と、（ｅ）前記第１の地図データの図郭が前記代表
点と一致するように、前記第１の地図データを逆正規化
変換する工程と、を具備することを特徴とする図形デー
タ整合方法である。
【０００６】
【作用】本発明では、（ａ）第１の地図データと第２の
地図データの対応する点を複数個指定し、（ｂ）対応す
る点を基にして、アフィン変換用パラメータを算出し、
（ｃ）第１の地図データの各メッシュの４隅の図郭を前
記アフィン変換用パラメータを用いてアフィン変換し、
（ｄ）あるメッシュの１つの図郭は、隣接する他の３つ
のメッシュの図郭でもあり、隣接する４つのメッシュの
同一点を表す４つの図郭が、工程（ｃ）により夫々アフ
ィン変換され、変換後の４つの図郭が一致しない場合に
は代表点を算出し、（ｅ）第１の地図データの図郭が代
表点と一致するように、第１の地図データを逆正規化変
換する。
地図データの対応する点を複数個指定し、（ｂ）対応す
る点を基にして、アフィン変換用パラメータを算出し、
（ｃ）第１の地図データの各メッシュの４隅の図郭を前
記アフィン変換用パラメータを用いてアフィン変換し、
（ｄ）あるメッシュの１つの図郭は、隣接する他の３つ
のメッシュの図郭でもあり、隣接する４つのメッシュの
同一点を表す４つの図郭が、工程（ｃ）により夫々アフ
ィン変換され、変換後の４つの図郭が一致しない場合に
は代表点を算出し、（ｅ）第１の地図データの図郭が代
表点と一致するように、第１の地図データを逆正規化変
換する。
【０００７】
【実施例】図１は、本発明の１実施例に係る図形データ
整合方法に用いられる装置の構成例を示す図である。１
はホストコンピュータ、２は内部記憶装置、３は地図デ
ータベース、５はアプリケーションプログラム、９はデ
ィスプレイ、１１はデジタイザ（またはタブレット
等）、１３ａはスタイラスペン、カーソル等の指示具、
１３ｂはマウス等の指示具、１５はデータ入力、コマン
ド入力、ディスプレイ９の表示カーソル移動等に用いる
キーボード、１７はプリンタ（または静電プロッタ等）
である。
整合方法に用いられる装置の構成例を示す図である。１
はホストコンピュータ、２は内部記憶装置、３は地図デ
ータベース、５はアプリケーションプログラム、９はデ
ィスプレイ、１１はデジタイザ（またはタブレット
等）、１３ａはスタイラスペン、カーソル等の指示具、
１３ｂはマウス等の指示具、１５はデータ入力、コマン
ド入力、ディスプレイ９の表示カーソル移動等に用いる
キーボード、１７はプリンタ（または静電プロッタ等）
である。
【０００８】地図データベース３には、道路台帳図から
入力された道路管理システムの道路データおよび都市計
画図の街区線データが保存される。道路管理システムの
道路データは、図８に示すようにメッシュ化されてい
る。３３ａ、３３ｂ……、は道路データの１メッシュを
示す。また、３１は都市計画図の地図データの仮想メッ
シュを示す。
入力された道路管理システムの道路データおよび都市計
画図の街区線データが保存される。道路管理システムの
道路データは、図８に示すようにメッシュ化されてい
る。３３ａ、３３ｂ……、は道路データの１メッシュを
示す。また、３１は都市計画図の地図データの仮想メッ
シュを示す。
【０００９】図２は、本実施例の動作を示すフローチャ
ートである。まず、道路管理システムの各メッシュ３３
毎に道路データを抽出し、ワーク領域７内の第１のワー
クファイルに展開する。また、メッシュ３３をカバーす
る範囲に含まれる都市計画図３１の街区線データを抽出
し、ワーク領域７内の第２のワークファイルに展開する
（ステップ３０１）。
ートである。まず、道路管理システムの各メッシュ３３
毎に道路データを抽出し、ワーク領域７内の第１のワー
クファイルに展開する。また、メッシュ３３をカバーす
る範囲に含まれる都市計画図３１の街区線データを抽出
し、ワーク領域７内の第２のワークファイルに展開する
（ステップ３０１）。
【００１０】次に、グラフィックディスプレイ９に道路
データ４１と街区線データ４３とを色を区別して重ねて
表示する（ステップ３０２）。図４は、グラフィックデ
ィスプレイ９上に表示された道路データ４１と街区線４
３とを示すものである。オペレータはグラフィックディ
スプレイ９上で道路データ４１上の交点等の２つの点Ｇ
１、Ｇ２を指定し、点Ｇ１、Ｇ２に対応する街区線デー
タ４３上の点Ｆ１、Ｆ２を指定する（ステップ３０
３）。そして、点Ｇ１、Ｇ２および点Ｆ１、Ｆ２の座標
値を用いて道路管理システムの道路データ４１およびメ
ッシュ３３の図郭の４点の座標を平行移動又は回転によ
り、座標変換する（ステップ３０４）。
データ４１と街区線データ４３とを色を区別して重ねて
表示する（ステップ３０２）。図４は、グラフィックデ
ィスプレイ９上に表示された道路データ４１と街区線４
３とを示すものである。オペレータはグラフィックディ
スプレイ９上で道路データ４１上の交点等の２つの点Ｇ
１、Ｇ２を指定し、点Ｇ１、Ｇ２に対応する街区線デー
タ４３上の点Ｆ１、Ｆ２を指定する（ステップ３０
３）。そして、点Ｇ１、Ｇ２および点Ｆ１、Ｆ２の座標
値を用いて道路管理システムの道路データ４１およびメ
ッシュ３３の図郭の４点の座標を平行移動又は回転によ
り、座標変換する（ステップ３０４）。
【００１１】すなわち、点Ｇ１、Ｇ２の座標を（Ｘ１、
Ｙ１）、（Ｘ２、Ｙ２）、点Ｆ１、Ｆ２の座標を（Ｘ１
´、Ｙ１´）、（Ｘ２´、Ｙ２´）とすると、道路デー
タ（Ｘ´、Ｙ´）は、次式のように変換される。
Ｙ１）、（Ｘ２、Ｙ２）、点Ｆ１、Ｆ２の座標を（Ｘ１
´、Ｙ１´）、（Ｘ２´、Ｙ２´）とすると、道路デー
タ（Ｘ´、Ｙ´）は、次式のように変換される。
【００１２】ＸＸ´＝（Ｘ´−Ｘ１´）cos θ＋（Ｙ´
−Ｙ１´）sin θ＋Ｘ１ ＹＹ´＝（Ｙ´−Ｙ１´）sin θ＋（Ｘ´−Ｘ１´）co
n θ＋Ｙ１ ただし、 θ＝tan -1［（Ｙ１−Ｙ２）／（Ｘ１−Ｘ２）］−tan
-1［（Ｙ１´−Ｙ２´）／（Ｘ１´−Ｘ２´）］ 次に、オペレータはアフィン変換を行うため、ディスプ
レイ９上で交差点を構成する２つのベクトルの組を指定
し（ステップ３０５）、２つのベクトルから交点を計算
する（ステップ３０６）。すなわち、図５に示すように
道路データｇai、ｇbiを指定することにより交点Ｇi を
計算し、同様に街区線データｆai、ｆbiを指定すること
により交点Ｆi を計算する。そして、ベクトルの指定が
終了するまで、ステップ３０５、３０６の処理を繰り返
す（ステップ３０７）。
−Ｙ１´）sin θ＋Ｘ１ ＹＹ´＝（Ｙ´−Ｙ１´）sin θ＋（Ｘ´−Ｘ１´）co
n θ＋Ｙ１ ただし、 θ＝tan -1［（Ｙ１−Ｙ２）／（Ｘ１−Ｘ２）］−tan
-1［（Ｙ１´−Ｙ２´）／（Ｘ１´−Ｘ２´）］ 次に、オペレータはアフィン変換を行うため、ディスプ
レイ９上で交差点を構成する２つのベクトルの組を指定
し（ステップ３０５）、２つのベクトルから交点を計算
する（ステップ３０６）。すなわち、図５に示すように
道路データｇai、ｇbiを指定することにより交点Ｇi を
計算し、同様に街区線データｆai、ｆbiを指定すること
により交点Ｆi を計算する。そして、ベクトルの指定が
終了するまで、ステップ３０５、３０６の処理を繰り返
す（ステップ３０７）。
【００１３】このようにして、道路データ４１と街区線
データ４３とで対応する点を多数指定して、アフィン変
換用パラメータおよび変位の標準偏差を求めて計算結果
を保存する（ステップ３０８）。
データ４３とで対応する点を多数指定して、アフィン変
換用パラメータおよび変位の標準偏差を求めて計算結果
を保存する（ステップ３０８）。
【００１４】ここで、アフィン変換用パラメータの算出
について説明する。まず、点Ｇi （Ｘi ，Ｙi ）を基準
とした場合の、２つの交点Ｇi 、Ｆi の座標のＸ方向の
差分δＸi およびＹ方向の差分δＹi を算出する。図６
は、交点Ｇi 、Ｆi の位置関係およびアフィン変換を説
明する図であり、実際にはディスプレイ９には表示され
ない第２の地図の線分も合わせて示されている。図６
（ａ）は、交点Ｇi と、交点Ｆi との位置関係を示し、
図６（ｂ）は、Ｘ方向の差分δＸi およびＹ方向の差分
δＹi 、ずれδＬi を示している。つまり、ディスプレ
イ９上のＦi の座標は（Ｘi ＋δＸi ，Ｙi ＋δＹi ）
となる。
について説明する。まず、点Ｇi （Ｘi ，Ｙi ）を基準
とした場合の、２つの交点Ｇi 、Ｆi の座標のＸ方向の
差分δＸi およびＹ方向の差分δＹi を算出する。図６
は、交点Ｇi 、Ｆi の位置関係およびアフィン変換を説
明する図であり、実際にはディスプレイ９には表示され
ない第２の地図の線分も合わせて示されている。図６
（ａ）は、交点Ｇi と、交点Ｆi との位置関係を示し、
図６（ｂ）は、Ｘ方向の差分δＸi およびＹ方向の差分
δＹi 、ずれδＬi を示している。つまり、ディスプレ
イ９上のＦi の座標は（Ｘi ＋δＸi ，Ｙi ＋δＹi ）
となる。
【００１５】なお、これまでは、算定された交点が地図
上の交点と一致する場合についてのみ説明したが、図６
（ｃ）に示すように交差点の隅切り等がなされている場
合には、２本の線分の延長線の交点を計算して求める。
このように、２つの地図において対応する点を正確に認
識することが困難な場合にも、より認識の容易な対応す
る線分を指定することにより対応する点を算定すること
ができるためオペレータの指定作業が容易になり、また
対応する点の組を多数求めることができるため的確なア
フィン変換用パラメータを算出することができる。
上の交点と一致する場合についてのみ説明したが、図６
（ｃ）に示すように交差点の隅切り等がなされている場
合には、２本の線分の延長線の交点を計算して求める。
このように、２つの地図において対応する点を正確に認
識することが困難な場合にも、より認識の容易な対応す
る線分を指定することにより対応する点を算定すること
ができるためオペレータの指定作業が容易になり、また
対応する点の組を多数求めることができるため的確なア
フィン変換用パラメータを算出することができる。
【００１６】また、図６（ｂ）、（ｃ）に示したずれδ
Ｌi の大きさは、Ｘ方向の差分δＸi およびＹ方向の差
分δＹi のデータから次式により求められる。
Ｌi の大きさは、Ｘ方向の差分δＸi およびＹ方向の差
分δＹi のデータから次式により求められる。
【００１７】
【数１】
次に、ホストコンピュータ１は、Ｎ組のＸ方向の差分δ
Ｘi およびＹ方向の差分δＹi のデータを用いて、Ｘ方
向とＹ方向につき別々に最小二乗法によりアフィン変換
用パラメータを算出する。
Ｘi およびＹ方向の差分δＹi のデータを用いて、Ｘ方
向とＹ方向につき別々に最小二乗法によりアフィン変換
用パラメータを算出する。
【００１８】つまり、一般に第１の地図の座標データ
（Ｘ，Ｙ）を、第２の地図の座標に適合する座標
（ｘ* ，ｙ* ）に補正するアフィン変換の式を、 ｘ* ＝−（ａＸ＋ｂＹ＋ｃ）＋Ｘ ｙ* ＝−（ｄＸ＋ｅＹ＋ｆ）＋Ｙ とするための、第１の地図全体に対して統一して用いる
アフィン変換用パラメータａ、ｂ、ｃ、ｄ、ｅ、ｆを決
定する必要がある。
（Ｘ，Ｙ）を、第２の地図の座標に適合する座標
（ｘ* ，ｙ* ）に補正するアフィン変換の式を、 ｘ* ＝−（ａＸ＋ｂＹ＋ｃ）＋Ｘ ｙ* ＝−（ｄＸ＋ｅＹ＋ｆ）＋Ｙ とするための、第１の地図全体に対して統一して用いる
アフィン変換用パラメータａ、ｂ、ｃ、ｄ、ｅ、ｆを決
定する必要がある。
【００１９】個々のデータの組については、第１の地図
のｉ番目の点Ｇi の座標（Ｘi 、Ｙi ）を第２の地図の
座標に適合する座標（ｘ* i ，ｙ* i ）に補正するため
のアフィン変換を、 ｘ* i ＝−（ａＸi ＋ｂＹi ＋ｃ）＋Ｘi ｙ* i ＝−（ｄＸi ＋ｅＹi ＋ｆ）＋Ｙi とすると、第２の地図から求めた、第１の地図の点Ｇi
に対応する点Ｆi の座標（ｘi ，ｙi ）と、アフィン変
換後の点Ｇi の座標（ｘ* i ，ｙ* i ）との、Ｘ方向の
誤差εxi、Ｙ方向の誤差εyiは、各々、次式で表され
る。
のｉ番目の点Ｇi の座標（Ｘi 、Ｙi ）を第２の地図の
座標に適合する座標（ｘ* i ，ｙ* i ）に補正するため
のアフィン変換を、 ｘ* i ＝−（ａＸi ＋ｂＹi ＋ｃ）＋Ｘi ｙ* i ＝−（ｄＸi ＋ｅＹi ＋ｆ）＋Ｙi とすると、第２の地図から求めた、第１の地図の点Ｇi
に対応する点Ｆi の座標（ｘi ，ｙi ）と、アフィン変
換後の点Ｇi の座標（ｘ* i ，ｙ* i ）との、Ｘ方向の
誤差εxi、Ｙ方向の誤差εyiは、各々、次式で表され
る。
【００２０】εxi＝ｘi −ｘ* i
＝δｘi −（ａＸi ＋ｂＹi ＋ｃ）
εyi＝ｙi −ｙ* i
＝δｙi −（ｄＸi ＋ｅＹi ＋ｆ）
なお、このとき、Ｘ方向の誤差εxi、Ｙ方向の誤差εyi
から、アフィン変換後の誤差の大きさδｌi は、
から、アフィン変換後の誤差の大きさδｌi は、
【００２１】
【数２】
により求めることができる。図７は、アフィン変換前後
の各点の位置関係等を説明する図であり、点Ｇi の座標
（Ｘi 、Ｙi ）、点Ｇi に対応する点Ｆi の座標
（ｘi ，ｙi ）、点Ｇi と点Ｆi との距離δＬi 、Ｘ方
向の差分δＸi およびＹ方向の差分δＹi 、アフィン変
換により点Ｇi を第２の地図の座標に適合させた後の座
標（ｘ* i ，ｙ* i ）、アフィン変換後の誤差の大きさ
δｌi 、Ｘ方向の誤差εxi、Ｙ方向の誤差εyiの関係を
示す。
の各点の位置関係等を説明する図であり、点Ｇi の座標
（Ｘi 、Ｙi ）、点Ｇi に対応する点Ｆi の座標
（ｘi ，ｙi ）、点Ｇi と点Ｆi との距離δＬi 、Ｘ方
向の差分δＸi およびＹ方向の差分δＹi 、アフィン変
換により点Ｇi を第２の地図の座標に適合させた後の座
標（ｘ* i ，ｙ* i ）、アフィン変換後の誤差の大きさ
δｌi 、Ｘ方向の誤差εxi、Ｙ方向の誤差εyiの関係を
示す。
【００２２】このアフィン変換後のＸ方向、Ｙ方向の誤
差εxiおよびεyiの二乗和をそれぞれ最小とするために
満足すべき条件は、
差εxiおよびεyiの二乗和をそれぞれ最小とするために
満足すべき条件は、
【００２３】
【数３】
である。したがって、誤差εxiおよびεyiの二乗和をそ
れぞれ最小とするアフィン変換用パラメータａ、ｂ、
ｃ、ｄ、ｅ、ｆを求めるには、
れぞれ最小とするアフィン変換用パラメータａ、ｂ、
ｃ、ｄ、ｅ、ｆを求めるには、
【００２４】
【数４】
の連立方程式を解けばよい。この連立方程式を解くに
は、たとえば「ＦＯＲＴＲＡＮ７７による科学計算サブ
ルーチンライブラリ」（啓学出版Ｐ８４〜８６）をプロ
グラムに準備しておけば、ホストコンピュータ１により
計算することができる。なお、アフィン変換用パラメー
タの内、ａはＸ方向の伸縮率、ｂ＋ｄは偏平率、ｃはＸ
方向の移動量、ｅはＹ方向の伸縮率、ｆはＹ方向の移動
量に関するものである。
は、たとえば「ＦＯＲＴＲＡＮ７７による科学計算サブ
ルーチンライブラリ」（啓学出版Ｐ８４〜８６）をプロ
グラムに準備しておけば、ホストコンピュータ１により
計算することができる。なお、アフィン変換用パラメー
タの内、ａはＸ方向の伸縮率、ｂ＋ｄは偏平率、ｃはＸ
方向の移動量、ｅはＹ方向の伸縮率、ｆはＹ方向の移動
量に関するものである。
【００２５】このようにして、アフィン変換用パラメー
タを算出したのち、道路管理システムのメッシュ３３の
４隅の各点（図郭）の座標についてアフィン変換を行い
保存する（ステップ３０９）。そして、ステップ３０１
からステップ３０９までの処理を道路管理システムの対
象とする全てのメッシュ３３について繰り返して行う
（ステップ３１０）。すなわち、図８に示すような複数
のメッシュ３３ａ、３３ｂ、３３ｃ、３３ｄ毎にメッシ
ュの図郭のアフィン変換を行う。
タを算出したのち、道路管理システムのメッシュ３３の
４隅の各点（図郭）の座標についてアフィン変換を行い
保存する（ステップ３０９）。そして、ステップ３０１
からステップ３０９までの処理を道路管理システムの対
象とする全てのメッシュ３３について繰り返して行う
（ステップ３１０）。すなわち、図８に示すような複数
のメッシュ３３ａ、３３ｂ、３３ｃ、３３ｄ毎にメッシ
ュの図郭のアフィン変換を行う。
【００２６】次に全てのメッシュ３３についてステップ
３０９の処理を行ったため、メッシュの４隅の各点（図
郭）についての座標値が１点につき４つ求められている
ので、単純平均又は、標準偏差による加重平均等の方法
により座標値を計算する（ステップ３１１）。
３０９の処理を行ったため、メッシュの４隅の各点（図
郭）についての座標値が１点につき４つ求められている
ので、単純平均又は、標準偏差による加重平均等の方法
により座標値を計算する（ステップ３１１）。
【００２７】図９は、ステップ３１１の説明図であり、
点Ａはメッシュ３３ａ、３３ｂ、３３ｃ、３３ｄの図郭
である。メッシュ３３ａをアフィン変換すると、点Ａは
点５１ａに変換され、メッシュ３３ｂをアフィン変換す
ると、点Ａは点５１ｂに変換され、メッシュ３３ｃをア
フィン変換すると、点Ａは点５１ｃに変換され、メッシ
ュ３３ｄをアフィン変換すると、点Ａは点５１ｄに変換
される。このように、通常図郭は４つの異なる点に変換
されるので、点５１ａ、５１ｂ、５１ｃ、５１ｄの座標
値の平均値を求める。図１０はこのようにして求められ
た代表点Ａ１、Ｂ１、Ｃ１、Ｄ１を示すものである。
点Ａはメッシュ３３ａ、３３ｂ、３３ｃ、３３ｄの図郭
である。メッシュ３３ａをアフィン変換すると、点Ａは
点５１ａに変換され、メッシュ３３ｂをアフィン変換す
ると、点Ａは点５１ｂに変換され、メッシュ３３ｃをア
フィン変換すると、点Ａは点５１ｃに変換され、メッシ
ュ３３ｄをアフィン変換すると、点Ａは点５１ｄに変換
される。このように、通常図郭は４つの異なる点に変換
されるので、点５１ａ、５１ｂ、５１ｃ、５１ｄの座標
値の平均値を求める。図１０はこのようにして求められ
た代表点Ａ１、Ｂ１、Ｃ１、Ｄ１を示すものである。
【００２８】次に、ステップ３１１までの処理で得られ
た道路データの座標変換結果、メッシュ３３の位置、標
準偏差等のデータを分かりやすく出図し、目視等により
適否を判断する（ステップ３１２）。そして、道路デー
タのメッシュ３３の図郭がステップ３１１の代表点に一
致するようにメッシュ内の道路データを逆正規化変換す
る。
た道路データの座標変換結果、メッシュ３３の位置、標
準偏差等のデータを分かりやすく出図し、目視等により
適否を判断する（ステップ３１２）。そして、道路デー
タのメッシュ３３の図郭がステップ３１１の代表点に一
致するようにメッシュ内の道路データを逆正規化変換す
る。
【００２９】図１０、図１１、図１２を用いて、逆正規
化変換について説明する。道路データのメッシュ３３毎
にメッシュの四隅の図郭Ａ、Ｂ、Ｃ、Ｄと変換すべき点
ＡＡの座標データを用いて、比例配分式により中間デー
タρ、μを算出する。
化変換について説明する。道路データのメッシュ３３毎
にメッシュの四隅の図郭Ａ、Ｂ、Ｃ、Ｄと変換すべき点
ＡＡの座標データを用いて、比例配分式により中間デー
タρ、μを算出する。
【００３０】次に、図１２に示すように中間データρ、
μを用いて、点Ｅ、Ｆ、Ｇ、Ｈを求め、更に、直線ＥＦ
と直線ＧＨの交点ＡＡ´を求める。この交点ＡＡ´が逆
正規化座標である。このように、メッシュ３３内の全て
のベクトルについて逆正規化変換を行う。この逆正規化
変換は、ｂｉｔ別冊「計算幾何学と地理情報処理」１９
８６年、９月号別冊、共立出版株式会社、３１頁に記載
された正規化座標変換の逆変換を行うものである。
μを用いて、点Ｅ、Ｆ、Ｇ、Ｈを求め、更に、直線ＥＦ
と直線ＧＨの交点ＡＡ´を求める。この交点ＡＡ´が逆
正規化座標である。このように、メッシュ３３内の全て
のベクトルについて逆正規化変換を行う。この逆正規化
変換は、ｂｉｔ別冊「計算幾何学と地理情報処理」１９
８６年、９月号別冊、共立出版株式会社、３１頁に記載
された正規化座標変換の逆変換を行うものである。
【００３１】図１３はこのようにして逆正規化変換され
た道路データ６１を示すもので、６１ａはメッシュ３３
ａ内で逆正規化変換された道路データを示し、６１ｂは
メッシュ３３ｅ内で逆正規化変換された道路データを示
す。このように、逆正規化変換することにより、道路デ
ータ６１は２つのメッシュの境界でも連続している。な
お、６１ｃはメッシュ３３ａを単純にアフィン変換した
場合の道路データであり、６１ｄはメッシュ３３ｅを単
純にアフィン変換した場合の道路データであり、メッシ
ュの境界で連続していない。
た道路データ６１を示すもので、６１ａはメッシュ３３
ａ内で逆正規化変換された道路データを示し、６１ｂは
メッシュ３３ｅ内で逆正規化変換された道路データを示
す。このように、逆正規化変換することにより、道路デ
ータ６１は２つのメッシュの境界でも連続している。な
お、６１ｃはメッシュ３３ａを単純にアフィン変換した
場合の道路データであり、６１ｄはメッシュ３３ｅを単
純にアフィン変換した場合の道路データであり、メッシ
ュの境界で連続していない。
【００３２】このように、本実施例では複数のメッシュ
の境界においても道路や設備等のデータを連続させるこ
とができる。
の境界においても道路や設備等のデータを連続させるこ
とができる。
【００３３】なお、代表点等を算出せず、目分量で二つ
の地図の四隅の図郭を合致させ、一報の地図の逆正規化
変換を行ってもよい。
の地図の四隅の図郭を合致させ、一報の地図の逆正規化
変換を行ってもよい。
【００３４】
【発明の効果】以上、詳細に説明したように本発明によ
れば、メッシュ化された地図データの各メッシュの境界
で地図データが整合される。
れば、メッシュ化された地図データの各メッシュの境界
で地図データが整合される。
【図１】 本発明の１実施例に係る図形データ整合方法
に用いられる装置の構成図
に用いられる装置の構成図
【図２】 実施例の動作を示すフローチャート
【図３】 実施例の動作を示すフローチャート
【図４】 ディスプレイ９上の表示を示す図
【図５】 対応する点を指定する場合の説明図
【図６】 交点Ｇｉ、Ｆ２ｉとの位置およびアフィン変
換の説明図
換の説明図
【図７】 アフィン変換前後の各点の位置関係等の説明
図
図
【図８】 道路管理システムのメッシュ３３と都市計画
図３１を示す図
図３１を示す図
【図９】 メッシュの図郭を変換した場合の説明図
【図１０】 代表点を示す図
【図１１】 逆正規化変換の説明図
【図１２】 逆正規化変換の説明図
【図１３】 逆正規化変換された道路データを示す図
１………ホストコンピュータ
２………内部記憶装置
３………地図データベース
５………アプリケーションプログラム
７………ワーク領域
９………ディスプレイ
１１………デジダイザ
１３ａ、１３ｂ………指示具
１５………キーボード
１７………プリンタ
─────────────────────────────────────────────────────
フロントページの続き
(56)参考文献 特開 平４−158388（ＪＰ，Ａ)
大西啓介・加藤誠巳，区分入力された
ディジタル地図データの相互接続，電子
情報通信学会論文誌Ｄ−ＩＩ，電子情報
通信学会，1992年 ７月25日，第Ｊ75−
Ｄ−ＩＩ，第７号，ｐ．1202−1212
(58)調査した分野(Int.Cl.7，ＤＢ名)
G06T 11/60 - 11/80
G06F 17/50
Claims (4)
- 【請求項１】 メッシュ化された第１の地図データを変
換して第２の地図データに整合させる図形データ整合方
法であって、 （ａ）第１の地図データと第２の地図データの対応する
点を複数個指定する工程と、 （ｂ）前記対応する点を基にして、アフィン変換用パラ
メータを算出する工程と、 （ｃ）前記第１の地図データの各メッシュの４隅の図郭
を前記アフィン変換用パラメータを用いてアフィン変換
する工程と、 （ｄ）あるメッシュの１つの図郭は、隣接する他の３つ
のメッシュの図郭でもあり、隣接する４つのメッシュの
同一点を表す４つの図郭が、前記工程（ｃ）により夫々
アフィン変換され、変換後の４つの図郭が一致しない場
合には代表点を算出する工程と、 （ｅ）前記第１の地図データの図郭が前記代表点と一致
するように、前記第１の地図データを逆正規化変換する
工程と、 を具備することを特徴とする図形データ整合方法。 - 【請求項２】 前記工程（ｄ）は、４個の図郭の座標値
の平均値を代表点の座標として算出することを特徴とす
る請求項１記載の図形データ整合方法。 - 【請求項３】 前記第１の地図データは、道路台帳図デ
ータであることを特徴とする請求項１又は２記載の図形
データ整合方法。 - 【請求項４】 前記第２の地図データは、国土基本図デ
ータであることを特徴とする請求項１又は２記載の図形
データ整合方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP33927594A JP3535598B2 (ja) | 1994-12-31 | 1994-12-31 | 図形データ整合方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP33927594A JP3535598B2 (ja) | 1994-12-31 | 1994-12-31 | 図形データ整合方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH08190617A JPH08190617A (ja) | 1996-07-23 |
| JP3535598B2 true JP3535598B2 (ja) | 2004-06-07 |
Family
ID=18325916
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP33927594A Expired - Lifetime JP3535598B2 (ja) | 1994-12-31 | 1994-12-31 | 図形データ整合方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3535598B2 (ja) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN100433584C (zh) * | 2005-02-25 | 2008-11-12 | 北京市东城区信息中心 | 城市管理信息采集方法及其所用的装置 |
-
1994
- 1994-12-31 JP JP33927594A patent/JP3535598B2/ja not_active Expired - Lifetime
Non-Patent Citations (1)
| Title |
|---|
| 大西啓介・加藤誠巳，区分入力されたディジタル地図データの相互接続，電子情報通信学会論文誌Ｄ−ＩＩ，電子情報通信学会，１９９２年 ７月２５日，第Ｊ７５−Ｄ−ＩＩ，第７号，ｐ．１２０２−１２１２ |
Also Published As
| Publication number | Publication date |
|---|---|
| JPH08190617A (ja) | 1996-07-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| KR0156270B1 (ko) | 지도데이터 작성방법 및 작성장치 | |
| CN113297650B (zh) | 一种基于bim技术的单元式玻璃幕墙施工方法及系统 | |
| Doytsher | A rubber sheeting algorithm for non-rectangular maps | |
| US6266676B1 (en) | Link information management method | |
| CN108711356B (zh) | 矢量地理pdf地图制作中地理目标与符号图形配准方法 | |
| JPH04232579A (ja) | 近傍データに基づく画像形状の比較方法 | |
| WO2009104218A1 (ja) | 地図表示装置 | |
| JP3535598B2 (ja) | 図形データ整合方法 | |
| JP2920195B2 (ja) | 形状変換方法及び装置 | |
| CN111898183B (zh) | 基于三维场景的架空电力线路辅助设计方法及装置 | |
| JP2662856B2 (ja) | 形状特徴計測装置および方法 | |
| JP3535597B2 (ja) | 図形処理装置 | |
| JP3535594B2 (ja) | 図形処理方法 | |
| CA2012406C (en) | Graphics system | |
| WO2024042702A1 (ja) | コンピュータシステム、寸法計測方法、および半導体装置製造システム | |
| JP3259981B2 (ja) | 図形データの一部修正システムおよび図形データの一部修正方法 | |
| JP3259998B2 (ja) | 設備断面図の自動作成システムおよび設備断面図の自動作成方法 | |
| JPH05324921A (ja) | 曲線形成方法および装置 | |
| CN110569620B (zh) | 一种混合模型的动态可视化系统及方法 | |
| JP3168117B2 (ja) | ３次元形状モデル処理装置 | |
| JP3402651B2 (ja) | 地図データベース検定方法 | |
| JP2571178B2 (ja) | ３次元図形処理装置 | |
| Qian et al. | Method for continuous multiscale representation of terrain contours using T‐spline technique | |
| JPH07141508A (ja) | 形状記述装置 | |
| CN117115371A (zh) | 一种高精地图的绘制方法、系统、电子设备及存储介质 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| 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: 20040309 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20040312 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080319 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080319 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090319 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090319 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100319 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100319 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110319 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110319 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120319 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120319 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130319 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130319 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140319 Year of fee payment: 10 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| EXPY | Cancellation because of completion of term |