JPH06251120A - 図形処理方法 - Google Patents
図形処理方法Info
- Publication number
- JPH06251120A JPH06251120A JP5062926A JP6292693A JPH06251120A JP H06251120 A JPH06251120 A JP H06251120A JP 5062926 A JP5062926 A JP 5062926A JP 6292693 A JP6292693 A JP 6292693A JP H06251120 A JPH06251120 A JP H06251120A
- Authority
- JP
- Japan
- Prior art keywords
- map
- point
- display
- affine transformation
- line segment
- 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
Landscapes
- Instructional Devices (AREA)
- Processing Or Creating Images (AREA)
Abstract
(57)【要約】 （修正有）
【構成】 ワーク領域７にある大縮尺の地図をディスプ
レイ９上に表示させ（ステップ２１）、たとえば点Ｇに
カーソルを合わせ、点Ｇを通る線分ｇを指定する（ステ
ップ２２）。つぎにデジタイザ１１の地図を見て、ディ
スプレイ９上の線分ｇに対応する線分ｇ１上の点Ｆを指
定する（ステップ２３）と、ディスプレイ９上では、点
Ｆ´となって表示される。アプリケーションプログラム
５は点Ｆ´と線分ｇとの距離と角度を算出する（ステッ
プ２４）。ステップ２２からステップ２４までの動作を
繰り返して点Ｇ、点Ｆを異ならせて、線分ｇと点Ｆの組
を複数組指定する（ステップ２５）と、プログラム５は
ステップ２４で算出された複数組の距離と角度を用いて
アフィン変換用パラメータを算出する（ステップ２
６）。ワーク領域７に格納されている地図データをアフ
ィン変換し（ステップ２７）、地図データベース３に格
納する。 【効果】 一方の地図を高精度でかつ効率よく他の地図
に整合できる。
レイ９上に表示させ（ステップ２１）、たとえば点Ｇに
カーソルを合わせ、点Ｇを通る線分ｇを指定する（ステ
ップ２２）。つぎにデジタイザ１１の地図を見て、ディ
スプレイ９上の線分ｇに対応する線分ｇ１上の点Ｆを指
定する（ステップ２３）と、ディスプレイ９上では、点
Ｆ´となって表示される。アプリケーションプログラム
５は点Ｆ´と線分ｇとの距離と角度を算出する（ステッ
プ２４）。ステップ２２からステップ２４までの動作を
繰り返して点Ｇ、点Ｆを異ならせて、線分ｇと点Ｆの組
を複数組指定する（ステップ２５）と、プログラム５は
ステップ２４で算出された複数組の距離と角度を用いて
アフィン変換用パラメータを算出する（ステップ２
６）。ワーク領域７に格納されている地図データをアフ
ィン変換し（ステップ２７）、地図データベース３に格
納する。 【効果】 一方の地図を高精度でかつ効率よく他の地図
に整合できる。
Description
【０００１】
【産業上の利用分野】本発明は、地図データをデータベ
ース化して利用するマッピングシステムに用いられる図
形処理方法に関するものである。
ース化して利用するマッピングシステムに用いられる図
形処理方法に関するものである。
【０００２】
【従来の技術】一般にマッピングシステムでは、地図デ
ータベースとして国家基準点を利用して作成した地図を
用いることが望まれる。
ータベースとして国家基準点を利用して作成した地図を
用いることが望まれる。
【０００３】しかしながら、国家基準点以外の基準点を
基に作成した近代の地図は、精度の高い測量方法に基づ
いており、部分的には国家基準点を利用して作成した地
図よりも高い精度を有することがある。たとえば他に先
駆けて地図を整備した川崎市は国家基準点を利用せず、
独自に縮尺１／５００の大縮尺の地図を作成しており、
この地図には道路の幅員、境界杭等が国家基準点に基づ
く地図（例えば縮尺１／２５００の国土基本図）よりも
正確に記されている。
基に作成した近代の地図は、精度の高い測量方法に基づ
いており、部分的には国家基準点を利用して作成した地
図よりも高い精度を有することがある。たとえば他に先
駆けて地図を整備した川崎市は国家基準点を利用せず、
独自に縮尺１／５００の大縮尺の地図を作成しており、
この地図には道路の幅員、境界杭等が国家基準点に基づ
く地図（例えば縮尺１／２５００の国土基本図）よりも
正確に記されている。
【０００４】
【発明が解決しようとする課題】このため、このような
部分的に高い精度を有する地図をマッピングシステムに
取り入れることが考えられるが、国家基準点以外の基準
点を基に作成した地図を地図データベース化する場合、
行政界で他の地図と整合しないことが多い。
部分的に高い精度を有する地図をマッピングシステムに
取り入れることが考えられるが、国家基準点以外の基準
点を基に作成した地図を地図データベース化する場合、
行政界で他の地図と整合しないことが多い。
【０００５】本発明は、このような問題に鑑みてなされ
たもので、その目的とするところは、一方の地図を高精
度でかつ効率よく他の地図に整合させることができる図
形処理方法を提供することにある。
たもので、その目的とするところは、一方の地図を高精
度でかつ効率よく他の地図に整合させることができる図
形処理方法を提供することにある。
【０００６】
【課題を解決するための手段】図１は、本発明の原理説
明図である。図１において１はホストコンピュータ、３
は地図データベース、５はアプリケーションプログラ
ム、７はワーク領域、９はディスプレイ、１１はデジタ
イザ、１３は指示具である。
明図である。図１において１はホストコンピュータ、３
は地図データベース、５はアプリケーションプログラ
ム、７はワーク領域、９はディスプレイ、１１はデジタ
イザ、１３は指示具である。
【０００７】ホストコンピュータ１はアプリケーション
プログラム５にしたがって各種処理を行う。地図データ
ベース３には例えば縮尺１／５００の大縮尺の地図デー
タが格納される。ディスプレイ９にはワーク領域７に格
納された、たとえば縮尺１／５００の大縮尺の地図デー
タが表示される。デジタイザ１１には、たとえば縮尺１
／２５００の中縮尺の地図が貼られる。この地図はディ
スプレイ９に表示された地図に記載された領域を含む。
デジタイザ１１において指示具１３により地図データを
入力することができる。指示具１３を動かすと、この動
きにつれて、ディスプレイ９上でカーソルが動き、この
指示具１３をクリックすることによりディスプレイ９上
の点、あるいは線分を指定することができる。
プログラム５にしたがって各種処理を行う。地図データ
ベース３には例えば縮尺１／５００の大縮尺の地図デー
タが格納される。ディスプレイ９にはワーク領域７に格
納された、たとえば縮尺１／５００の大縮尺の地図デー
タが表示される。デジタイザ１１には、たとえば縮尺１
／２５００の中縮尺の地図が貼られる。この地図はディ
スプレイ９に表示された地図に記載された領域を含む。
デジタイザ１１において指示具１３により地図データを
入力することができる。指示具１３を動かすと、この動
きにつれて、ディスプレイ９上でカーソルが動き、この
指示具１３をクリックすることによりディスプレイ９上
の点、あるいは線分を指定することができる。
【０００８】Ａはアプリケーションプログラム５の内容
およびオペレータの操作を示すフローチャートである。
まずワーク領域７にある大縮尺の地図をディスプレイ９
上に表示させる（ステップ２１）。オペレータはディス
プレイ９を見ながら指示具１３を動かしディスプレイ９
上の、たとえば点Ｇにカーソルを合わせ、指示具１３を
用いてＧを通る線分ｇを指定する（ステップ２２）。つ
ぎにオペレータはデジタイザ１１に貼られた中縮尺の地
図に目を向け、ディスプレイ９上の線分ｇに対応する線
分ｇ１上の点Ｆを指示具１３で指定する（ステップ２
３）。デジタイザ１１により点Ｆを指定すると、ディス
プレイ９上では、それが点Ｆ´となって表示される。点
Ｆ´が線分ｇ上に存在しないのは、中縮尺の地図と大縮
尺の地図の対応する部分で微妙なずれがあり、両者が整
合していないからである。
およびオペレータの操作を示すフローチャートである。
まずワーク領域７にある大縮尺の地図をディスプレイ９
上に表示させる（ステップ２１）。オペレータはディス
プレイ９を見ながら指示具１３を動かしディスプレイ９
上の、たとえば点Ｇにカーソルを合わせ、指示具１３を
用いてＧを通る線分ｇを指定する（ステップ２２）。つ
ぎにオペレータはデジタイザ１１に貼られた中縮尺の地
図に目を向け、ディスプレイ９上の線分ｇに対応する線
分ｇ１上の点Ｆを指示具１３で指定する（ステップ２
３）。デジタイザ１１により点Ｆを指定すると、ディス
プレイ９上では、それが点Ｆ´となって表示される。点
Ｆ´が線分ｇ上に存在しないのは、中縮尺の地図と大縮
尺の地図の対応する部分で微妙なずれがあり、両者が整
合していないからである。
【０００９】アプリケーションプログラム５はディスプ
レイ９上における点Ｆに対応する点Ｆ´と線分ｇとの距
離と角度を算出する（ステップ２４）。オペレータはス
テップ２２からステップ２４までの動作を繰り返して点
Ｇ、点Ｆを異ならしめて、線分ｇと点Ｆの組を複数組指
定する（ステップ２５）。このように線分ｇと点Ｆが複
数組指定されると、アプリケーションプログラム５はス
テップ２４で算出された複数組の距離と角度を用いてア
フィン変換用パラメータを算出する（ステップ２６）。
そしてワーク領域７に格納されている大縮尺の地図デー
タをアフィン変換用パラメータを用いてアフィン変換す
る（ステップ２７）。アフィン変換された地図データは
地図データベース３に格納される。
レイ９上における点Ｆに対応する点Ｆ´と線分ｇとの距
離と角度を算出する（ステップ２４）。オペレータはス
テップ２２からステップ２４までの動作を繰り返して点
Ｇ、点Ｆを異ならしめて、線分ｇと点Ｆの組を複数組指
定する（ステップ２５）。このように線分ｇと点Ｆが複
数組指定されると、アプリケーションプログラム５はス
テップ２４で算出された複数組の距離と角度を用いてア
フィン変換用パラメータを算出する（ステップ２６）。
そしてワーク領域７に格納されている大縮尺の地図デー
タをアフィン変換用パラメータを用いてアフィン変換す
る（ステップ２７）。アフィン変換された地図データは
地図データベース３に格納される。
【００１０】
【作用】本発明では、ディスプレイ９に表示されている
大縮尺の地図とデジタイザ１１上に貼り付けられた中縮
尺の地図の対応する部分が整合していない場合でも、複
数組の線分と点を指定することにより、ワーク領域７に
ある大縮尺の地図データをアフィン変換し、大縮尺の地
図を中縮尺の地図に整合するように変換することができ
る。
大縮尺の地図とデジタイザ１１上に貼り付けられた中縮
尺の地図の対応する部分が整合していない場合でも、複
数組の線分と点を指定することにより、ワーク領域７に
ある大縮尺の地図データをアフィン変換し、大縮尺の地
図を中縮尺の地図に整合するように変換することができ
る。
【００１１】
【実施例】以下、図面に基づいて本発明の実施例を詳細
に説明する。図２は、本発明の１実施例に係る図形処理
方法に用いられる装置の構成を示す図である。１はホス
トコンピュータ、３は地図データベース、５はアプリケ
ーションプログラム、９はディスプレイ、１１はデジタ
イザ、１３は指示具、１５はキーボード、１７はプリン
タである。
に説明する。図２は、本発明の１実施例に係る図形処理
方法に用いられる装置の構成を示す図である。１はホス
トコンピュータ、３は地図データベース、５はアプリケ
ーションプログラム、９はディスプレイ、１１はデジタ
イザ、１３は指示具、１５はキーボード、１７はプリン
タである。
【００１２】図３、図４は本実施例の処理を示すフロー
チャートである。まずオペレータは大縮尺の地図（たと
えば縮尺１／５００の川崎市の地図）をデジタイザ１１
に貼り付け、指示具１３を用いて、この地図上の街区線
等の地図データを入力してワーク領域７に保管する。
チャートである。まずオペレータは大縮尺の地図（たと
えば縮尺１／５００の川崎市の地図）をデジタイザ１１
に貼り付け、指示具１３を用いて、この地図上の街区線
等の地図データを入力してワーク領域７に保管する。
【００１３】図５は、デジタイザ１１に貼り付けられた
大縮尺の地図を示す。図５においてＥ１部は現状の道路
幅員で記載されており、国土基本図よりもＥ１部におい
ては正確に記載されている。
大縮尺の地図を示す。図５においてＥ１部は現状の道路
幅員で記載されており、国土基本図よりもＥ１部におい
ては正確に記載されている。
【００１４】この地図をデジタイザ１１により入力する
場合、まずＡ１を基準点として入力し、４隅Ａ１、Ｂ
１、Ｃ１、Ｄ１を指示具１３で指定して、その座標値を
入力し、たとえば特開平１−３２０５７９号公報に記載
された方法により図郭の補正を行うためのパラメータを
算出する。そして指示具１３を用いて街区線等のデータ
を入力する。さらに入力されたデータに縮尺の逆数（こ
の場合５００）を掛けて、現地寸法に換算し、前述した
パラメータにより補正し、ワーク領域７に保管する。そ
して、その地図をディスプレイ９上に表示させる（ステ
ップ３０１）。
場合、まずＡ１を基準点として入力し、４隅Ａ１、Ｂ
１、Ｃ１、Ｄ１を指示具１３で指定して、その座標値を
入力し、たとえば特開平１−３２０５７９号公報に記載
された方法により図郭の補正を行うためのパラメータを
算出する。そして指示具１３を用いて街区線等のデータ
を入力する。さらに入力されたデータに縮尺の逆数（こ
の場合５００）を掛けて、現地寸法に換算し、前述した
パラメータにより補正し、ワーク領域７に保管する。そ
して、その地図をディスプレイ９上に表示させる（ステ
ップ３０１）。
【００１５】つぎにオペレータは、中縮尺の地図（縮尺
１／２５００の国土基本図）をデジタイザ１１に貼り付
け、図郭の４隅Ａ２、Ｂ２、Ｃ２、Ｄ２を入力し、ディ
スプレイ９に表示された領域に該当する領域の左下隅の
点Ａの座標値を入力する（ステップ３０２）。
１／２５００の国土基本図）をデジタイザ１１に貼り付
け、図郭の４隅Ａ２、Ｂ２、Ｃ２、Ｄ２を入力し、ディ
スプレイ９に表示された領域に該当する領域の左下隅の
点Ａの座標値を入力する（ステップ３０２）。
【００１６】このステップ３０２の処理を詳細に説明す
る。図６は、デジタイザ１１に貼り付けられた中縮尺の
地図を示す。この中縮尺の地図においてＡＢＣＤがディ
スプレイ９に表示されている大縮尺の地図に予定される
部分である。
る。図６は、デジタイザ１１に貼り付けられた中縮尺の
地図を示す。この中縮尺の地図においてＡＢＣＤがディ
スプレイ９に表示されている大縮尺の地図に予定される
部分である。
【００１７】図７は、図６の領域ＡＢＣＤを拡大したも
のである。この図で示されるように、Ｅ部では道路幅員
が狭いので所定の値、例えば１mm幅で描かれている。
のである。この図で示されるように、Ｅ部では道路幅員
が狭いので所定の値、例えば１mm幅で描かれている。
【００１８】図８は、ステップ３０２の入力の説明図で
ある。まず指示具１３を用いて４隅Ａ２、Ｂ２、Ｃ２、
Ｄ２を指定して、この座標値を入力する。つぎにＡ２の
座標（ＸA ，ＹA ）を入力する。つぎに指示具１３によ
り点Ａを指定して、点Ａの座標（ｌx ，ｌy ）を入力す
る。ここでＸA 、ＹA は絶対座標である。またｌx 、ｌ
y は図６の地図上において、Ａ２を原点とした場合の点
Ａの座標である。このようにして、点Ａの座標値を入力
する。
ある。まず指示具１３を用いて４隅Ａ２、Ｂ２、Ｃ２、
Ｄ２を指定して、この座標値を入力する。つぎにＡ２の
座標（ＸA ，ＹA ）を入力する。つぎに指示具１３によ
り点Ａを指定して、点Ａの座標（ｌx ，ｌy ）を入力す
る。ここでＸA 、ＹA は絶対座標である。またｌx 、ｌ
y は図６の地図上において、Ａ２を原点とした場合の点
Ａの座標である。このようにして、点Ａの座標値を入力
する。
【００１９】つぎに組数を表わすパラメータｉをクリア
し（ステップ３０３）、ディスプレイ９上の点Ｇを指示
具１３により入力する（ステップ３０４）。
し（ステップ３０３）、ディスプレイ９上の点Ｇを指示
具１３により入力する（ステップ３０４）。
【００２０】ステップ３０４の処理について説明する。
図９はディスプレイ９の現状の表示内容と、デジタイザ
１１上の領域ＡＢＣＤをそのままディスプレイ９上に表
示したと仮定した表示内容を重畳させて示したものであ
る。図９において実線はディスプレイ９の実際の表示内
容を示し、破線はデジタイザ１１上の図面をそのまま表
示したと仮定した場合の表示内容を示す。
図９はディスプレイ９の現状の表示内容と、デジタイザ
１１上の領域ＡＢＣＤをそのままディスプレイ９上に表
示したと仮定した表示内容を重畳させて示したものであ
る。図９において実線はディスプレイ９の実際の表示内
容を示し、破線はデジタイザ１１上の図面をそのまま表
示したと仮定した場合の表示内容を示す。
【００２１】オペレータはディスプレイ９を見ながら指
示具１３を操作する。指示具１３の動きに連動してディ
スプレイ９上でカーソルが動く。オペレータはカーソル
が図９の点Ｇ付近まで移動すると、指示具１３により入
力を行う。
示具１３を操作する。指示具１３の動きに連動してディ
スプレイ９上でカーソルが動く。オペレータはカーソル
が図９の点Ｇ付近まで移動すると、指示具１３により入
力を行う。
【００２２】このように点Ｇが指定されると、ホストコ
ンピュータ１は指定された点の最も近傍にある線分ｇ
（本実施例の場合点Ｇは線分ｇ上に存在する）をたとえ
ば赤で表示する（ステップ３０５）。図１０は、このと
きのディスプレイ９の表示内容を示すもので、線分ｇが
赤で表示されている。
ンピュータ１は指定された点の最も近傍にある線分ｇ
（本実施例の場合点Ｇは線分ｇ上に存在する）をたとえ
ば赤で表示する（ステップ３０５）。図１０は、このと
きのディスプレイ９の表示内容を示すもので、線分ｇが
赤で表示されている。
【００２３】つぎにオペレータはデジタイザ１１に目を
向け、デジタイザ１１に貼られた中縮尺の地図上で線分
ｇ１に対応する線分上の点Ｆを指示具１３を用いて入力
する（ステップ３０６）。
向け、デジタイザ１１に貼られた中縮尺の地図上で線分
ｇ１に対応する線分上の点Ｆを指示具１３を用いて入力
する（ステップ３０６）。
【００２４】Ｆ点を入力した後、以下の変換を行う。
【００２５】Ｘ＝ｂｘ＋ＸA Ｙ＝ｂｙ＋ＹA ここでｘ、ｙは点Ｆのデジタイザ１１上での座標、ｂは
縮尺の逆数を表わす。
縮尺の逆数を表わす。
【００２６】図９に示すように大縮尺の地図と中縮尺の
地図とは、その対応する部分においても微妙なずれが存
在する。デジタイザ１１に貼られた中縮尺の地図上で点
Ｆを入力すると、図１０に示すようにディスプレイ９上
では点Ｆ´が入力されたことになる。
地図とは、その対応する部分においても微妙なずれが存
在する。デジタイザ１１に貼られた中縮尺の地図上で点
Ｆを入力すると、図１０に示すようにディスプレイ９上
では点Ｆ´が入力されたことになる。
【００２７】つぎにホストコンピュータ１は、線分ｇと
点Ｆ´との距離Ｌi と角度θi を計算し、点Ｆ´の座標
（Ｘi ，Ｙi ）とともに格納する（ステップ３０７）。
点Ｆ´との距離Ｌi と角度θi を計算し、点Ｆ´の座標
（Ｘi ，Ｙi ）とともに格納する（ステップ３０７）。
【００２８】距離Ｌi と角度θi を計算はつぎのように
して行う。図１１は線分ｇと点Ｆ´を示す図である。点
Ｆ´の座標を今（Ｘ1 ，Ｙ1 ）とし、線分ｇ上の任意の
２点の座標（Ｘ2 ，Ｙ2 ，）、（Ｘ3 ，Ｙ3 ）とする。
点Ｆ´から直線ｇに垂線をおろし、その垂線の長さを距
離Ｌi 、その垂線の水平軸に対する傾きを角度θi とす
る。
して行う。図１１は線分ｇと点Ｆ´を示す図である。点
Ｆ´の座標を今（Ｘ1 ，Ｙ1 ）とし、線分ｇ上の任意の
２点の座標（Ｘ2 ，Ｙ2 ，）、（Ｘ3 ，Ｙ3 ）とする。
点Ｆ´から直線ｇに垂線をおろし、その垂線の長さを距
離Ｌi 、その垂線の水平軸に対する傾きを角度θi とす
る。
【００２９】ａ＝Ｙ2 −Ｙ3 、ｂ＝Ｘ3 −Ｘ2 とすると θi ＝ｓｇｎ（ｂ／ａ）・ｔａｎ-1｜ｂ／ａ｜ （ただし−π／２＜θi ＜π／２、ａ＝０のときθi ＝
π／２） となる。ここでｓｇｎ（ｂ／ａ）は線分ｇのベクトル方
向が第２象限、第４象限にあるときは正、第１象限、第
３象限にあるときは負となることを示す。
π／２） となる。ここでｓｇｎ（ｂ／ａ）は線分ｇのベクトル方
向が第２象限、第４象限にあるときは正、第１象限、第
３象限にあるときは負となることを示す。
【００３０】一方、距離Ｌi を求めるには ｃ＝（Ｘ3 −Ｘ2 ）・Ｙ2 ＋（Ｙ2 −Ｙ3 ）・Ｘ2 とすると Ｌi ＝ｓｇｎ（ａ）・（ａＸ1 ＋ｂＹ1 −ｃ）／（ａ2
＋ｂ2 ）1/2 となる。
＋ｂ2 ）1/2 となる。
【００３１】つぎにパラメータｉを「１」増加させ（ス
テップ３０８）、オペレータが終了の指示を行うまでス
テップ３０４からステップ３０８までの処理が繰り返さ
れる（ステップ３０９）。すなわち異なる点Ｇおよび点
Ｆの組が複数組指定され、各組ごとに、その距離Ｌ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 の２乗和が最小となるための条件は、
【００３５】
【数１】
【００３６】である。計算結果は未定数パラメータを
Ａ、Ｂ、Ｃ、Ｄ、Ｅ、Ｆとした連立方程式となり、行列
で示すと次式の通りである。
Ａ、Ｂ、Ｃ、Ｄ、Ｅ、Ｆとした連立方程式となり、行列
で示すと次式の通りである。
【００３７】 ［Ｋhj］［Ａ，Ｂ，Ｃ，Ｄ，Ｅ，Ｆ，ｌ］＝０ ［Ｋhj］はＫ11、Ｋ12、…Ｋ17、…Ｋ61、…Ｋ67という
ように４２個の係数の行列になるが、たとえばこの中の
Ｋ11からＫ17を示すと、つぎのようになる。
ように４２個の係数の行列になるが、たとえばこの中の
Ｋ11からＫ17を示すと、つぎのようになる。
【００３８】
【数２】
【００３９】この連立方程式を解くには、たとえば「Ｆ
ＯＲＴＲＡＮ７７による科学計算サブルーチンライブラ
リ」（啓学出版Ｐ８４〜８６）をプログラムに準備して
おけば、ホストコンピュータ１により計算することがで
きる。なお、アフィン変換用パラメータの内、ＡはＸ方
向の伸縮率、（Ｂ＋Ｄ）／２は偏平率、ＣはＸ方向の移
動量、ＥはＹ方向の伸縮率、（Ｂ−Ｄ）／２は回転量、
ＦはＹ方向の移動量に関するものである。
ＯＲＴＲＡＮ７７による科学計算サブルーチンライブラ
リ」（啓学出版Ｐ８４〜８６）をプログラムに準備して
おけば、ホストコンピュータ１により計算することがで
きる。なお、アフィン変換用パラメータの内、ＡはＸ方
向の伸縮率、（Ｂ＋Ｄ）／２は偏平率、ＣはＸ方向の移
動量、ＥはＹ方向の伸縮率、（Ｂ−Ｄ）／２は回転量、
ＦはＹ方向の移動量に関するものである。
【００４０】つぎにワーク領域７内のデータをステップ
３１１で求められたアフィン変換用パラメータを用いて
変換する（ステップ３１２）。
３１１で求められたアフィン変換用パラメータを用いて
変換する（ステップ３１２）。
【００４１】すなわちワーク領域７内のデータの座標を
Ｘ、Ｙとし、アフィン変換用パラメータを用いて変換し
た座標をＸ＊、Ｙ＊とすると、 Ｘ＊＝Ｘ＋（ＡＸ＋ＢＹ＋Ｃ） Ｙ＊＝Ｙ＋（ＤＸ＋ＥＹ＋Ｆ）となる。この座標変換を
ワーク領域７内にある地図データについて行う。
Ｘ、Ｙとし、アフィン変換用パラメータを用いて変換し
た座標をＸ＊、Ｙ＊とすると、 Ｘ＊＝Ｘ＋（ＡＸ＋ＢＹ＋Ｃ） Ｙ＊＝Ｙ＋（ＤＸ＋ＥＹ＋Ｆ）となる。この座標変換を
ワーク領域７内にある地図データについて行う。
【００４２】そして変換されたデータを地図データベー
ス３内に格納する（ステップ３１３）。なお、本実施例
では、大縮尺地図を入力しながら修正したが、この大縮
尺地図が既に入力されていて、データベース化されてい
るような場合にも適応される。 このように本実施例で
は、部分的に高い精度を有する大縮尺の地図（たとえば
縮尺１／５００の川崎市の地図）をアフィン変換を用い
て中縮尺の地図（たとえば国土基本図）に整合するよう
に変換することができる。アフィン変換を用いて変換さ
れた地図は部分的に高い精度を有するので、地図データ
ベースとしては好適である。
ス３内に格納する（ステップ３１３）。なお、本実施例
では、大縮尺地図を入力しながら修正したが、この大縮
尺地図が既に入力されていて、データベース化されてい
るような場合にも適応される。 このように本実施例で
は、部分的に高い精度を有する大縮尺の地図（たとえば
縮尺１／５００の川崎市の地図）をアフィン変換を用い
て中縮尺の地図（たとえば国土基本図）に整合するよう
に変換することができる。アフィン変換を用いて変換さ
れた地図は部分的に高い精度を有するので、地図データ
ベースとしては好適である。
【００４３】［第２実施例］つぎに本発明の第２実施例
を図１３、図１４に示すフローチャートに従って説明す
る。この第２実施例は、交差点の隅切り等が行われて、
街区が経年変化しているにも拘らずオペレータが、その
事実に気付かずに経年変化に係る街区線を入力指定する
と、アフィン変換用パラメータの精度が悪くなる恐れが
あるので、その経年変化したと思われるデータを除外し
てアフィン変換用パラメータを求めるようにしたもので
ある。
を図１３、図１４に示すフローチャートに従って説明す
る。この第２実施例は、交差点の隅切り等が行われて、
街区が経年変化しているにも拘らずオペレータが、その
事実に気付かずに経年変化に係る街区線を入力指定する
と、アフィン変換用パラメータの精度が悪くなる恐れが
あるので、その経年変化したと思われるデータを除外し
てアフィン変換用パラメータを求めるようにしたもので
ある。
【００４４】すなわちホストコンピュータ１はステップ
３１１において、つぎのような処理を行う。
３１１において、つぎのような処理を行う。
【００４５】すなわち格納されている座標Ｘｉ、Ｙ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を読み出して、その区域に
ついてのアフィン変換用パラメータを最小二乗法により
決定する（ステップ１６０５）。そして決定したアフィ
ン変換用パラメータにより、その区域のすべての座標デ
ータをアフィン変換する（ステップ１６０６）。 つぎ
に残りの区域についてステップ１６０５、１６０６と全
く同様の処理を行う（ステップ１６０７）。そして最後
にオペレータは各区域間の境界部分の座標データについ
て接合編集を行い（ステップ１６０８）、処理を終了す
る。
【００５７】［第４実施例］つぎに第４実施例について
説明する。第４実施例は街区線の種別、たとえば国道、
都道、区道、市道、私道、歩道等に応じて重み付けをし
てアフィン変換用パラメータを決定するようにする。た
とえば国道、都道、区道、市道、私道、歩道についてそ
れぞれ１つずづ街区線が指定されたとすると、それらが
各々６個、５個、４個、３個、２個、１個ずつ指定され
たものとしてアフィン変換用パラメータを決定する。こ
のようにすれば重要度の高い街区線に沿った座標データ
ほど高精度にアフィン変換される。
説明する。第４実施例は街区線の種別、たとえば国道、
都道、区道、市道、私道、歩道等に応じて重み付けをし
てアフィン変換用パラメータを決定するようにする。た
とえば国道、都道、区道、市道、私道、歩道についてそ
れぞれ１つずづ街区線が指定されたとすると、それらが
各々６個、５個、４個、３個、２個、１個ずつ指定され
たものとしてアフィン変換用パラメータを決定する。こ
のようにすれば重要度の高い街区線に沿った座標データ
ほど高精度にアフィン変換される。
【００５８】［第５実施例］つぎに第５実施例について
説明する。第５実施例では、線分ｇと点Ｆ´との距離と
角度を求める代わりに、点Ｇと点Ｆ´との間で、座標位
置データをＸ成分と、Ｙ成分とに分割することにより線
データ間の距離と角度に基づいて、アフィン変換用のパ
ラメータを決定する図３、図４の処理を全く同様に実行
するようにしたものである。
説明する。第５実施例では、線分ｇと点Ｆ´との距離と
角度を求める代わりに、点Ｇと点Ｆ´との間で、座標位
置データをＸ成分と、Ｙ成分とに分割することにより線
データ間の距離と角度に基づいて、アフィン変換用のパ
ラメータを決定する図３、図４の処理を全く同様に実行
するようにしたものである。
【００５９】すなわち指定された大縮尺地図の点Ｇと中
縮尺地図の点Ｆ´の座標位置データが図１６に示したよ
うに（Ｘ1 ，Ｙ1 ）、（Ｘ2 ，Ｙ2 ）であったとする
と、図３のステップ３０４では、それぞれＸ軸、Ｙ軸が
指定されたものとみなして、（ｉ＋１）番目におけるス
テップ３０７では距離Ｌi+1 を（Ｘ1 −Ｘ2 ）、角度θ
i+1 を０度とし、（ｉ＋２）番目におけるステップ３０
７では距離Ｌi+2 を（Ｙ1 −Ｙ2 ）、角度θi+2 を９０
度として処理するものである。
縮尺地図の点Ｆ´の座標位置データが図１６に示したよ
うに（Ｘ1 ，Ｙ1 ）、（Ｘ2 ，Ｙ2 ）であったとする
と、図３のステップ３０４では、それぞれＸ軸、Ｙ軸が
指定されたものとみなして、（ｉ＋１）番目におけるス
テップ３０７では距離Ｌi+1 を（Ｘ1 −Ｘ2 ）、角度θ
i+1 を０度とし、（ｉ＋２）番目におけるステップ３０
７では距離Ｌi+2 を（Ｙ1 −Ｙ2 ）、角度θi+2 を９０
度として処理するものである。
【００６０】［第６実施例］つぎに第６実施例について
説明する。第６実施例は図３のステップ３０２で図郭の
４隅を入力するさいに、その座標値をアフィン変換用パ
ラメータを用いて変換し、そのデータを他の中縮尺の地
図データ（国土基本図）と識別できるようにして保管
し、出図のさいにはオペレータのコマンド指示により、
その座標およびその座標で結ばれる直線を表示するもの
である。
説明する。第６実施例は図３のステップ３０２で図郭の
４隅を入力するさいに、その座標値をアフィン変換用パ
ラメータを用いて変換し、そのデータを他の中縮尺の地
図データ（国土基本図）と識別できるようにして保管
し、出図のさいにはオペレータのコマンド指示により、
その座標およびその座標で結ばれる直線を表示するもの
である。
【００６１】図１７は、このときのディスプレイ９の表
示内容を示す図である。図１７において、２点鎖線が中
縮尺地図の図郭の４隅を結ぶ直線である。なおディスプ
レイ９からはみだした２点鎖線は他のメッシュに対応し
て記憶されている。
示内容を示す図である。図１７において、２点鎖線が中
縮尺地図の図郭の４隅を結ぶ直線である。なおディスプ
レイ９からはみだした２点鎖線は他のメッシュに対応し
て記憶されている。
【００６２】すなわち、Ａ１´、Ｂ１´、Ｃ１´、Ｄ１
´の座標を（ｘA 、ｙA ）、（ｘB、ｙB ）、（ｘC 、
ｙC ）、（ｘD 、ｙD ）とすると、各座標値を次式のよ
うにアフィン変換して、変換後の座標データ等と、とも
に保管する。
´の座標を（ｘA 、ｙA ）、（ｘB、ｙB ）、（ｘC 、
ｙC ）、（ｘD 、ｙD ）とすると、各座標値を次式のよ
うにアフィン変換して、変換後の座標データ等と、とも
に保管する。
【００６３】Ｘ＝ｘA ＋（ＡｘA ＋ＢｙA ＋Ｃ） Ｙ＝ｙA ＋（ＤｘA ＋ＥＹA ＋Ｆ） Ｘ＝ｘB ＋（ＡｘB ＋ＢｙB ＋Ｃ） Ｙ＝ｙB ＋（ＤｘB ＋ＥＹB ＋Ｆ） Ｘ＝ｘC ＋（ＡｘC ＋ＢｙC ＋Ｃ） Ｙ＝ｙC ＋（ＤｘC ＋ＥＹC ＋Ｆ） Ｘ＝ｘD ＋（ＡｘD ＋ＢｙD ＋Ｃ） Ｙ＝ｙD ＋（ＤｘD ＋ＥＹD ＋Ｆ）
【００６４】
【発明の効果】以上、詳細に説明したように本発明によ
れば、一方の地図を高精度でかつ効率よく他の地図に整
合させることができる。
れば、一方の地図を高精度でかつ効率よく他の地図に整
合させることができる。
【図１】 本発明の原理説明図
【図２】 本発明の１実施例に係る図形処理方法を実行
するための装置構成図
するための装置構成図
【図３】 実施例の動作を示すフローチャート
【図４】 実施例の動作を示すフローチャート
【図５】 ディスプレイ９に表示された大縮尺地図を示
す図
す図
【図６】 デジタイザ１１上に貼られた中縮尺地図を示
す図
す図
【図７】 図６の矩形Ａ、Ｂ、Ｃ、Ｄの拡大図
【図８】 図６の矩形Ａ、Ｂ、Ｃ、Ｄの入力を説明する
図
図
【図９】 ディスプレイ９の表示内容とデジタイザ１１
上の図面の該当部分を、そのままディスプレイ９に表示
したと仮定した場合の説明図
上の図面の該当部分を、そのままディスプレイ９に表示
したと仮定した場合の説明図
【図１０】 ステップ３０５において線分ｇが強調され
て表示されているディスプレイ９の表示内容を示す図
て表示されているディスプレイ９の表示内容を示す図
【図１１】 線分ｇと点Ｆとの位置関係を示す図
【図１２】 アフィン変換を説明する図
【図１３】 アフィン変換用パラメータの決定要素を自
動的に除外する場合の動作を示すフローチャート
動的に除外する場合の動作を示すフローチャート
【図１４】 アフィン変換用パラメータの決定要素を半
自動的に除外する場合の動作を示すフローチャート
自動的に除外する場合の動作を示すフローチャート
【図１５】 区域別にアフィン変換用パラメータを決定
する場合の動作を示すフローチャート
する場合の動作を示すフローチャート
【図１６】 各地図上の基準点が、それぞれ指定された
場合の処理を説明する図
場合の処理を説明する図
【図１７】 中縮尺の地図上に大縮尺の該当する部分の
図郭が示された図
図郭が示された図
１………ホストコンピュータ ３………地図データベース ５………アプリケーションプログラム ７………ワーク領域 ９………ディスプレイ １１………デジタイザ １５………キーボード １７………プリンタ
Claims (8)
- 【請求項１】 デジタイザと、ディスプレイと、記憶装
置とを備えたシステムにおいて、 （ａ）第１の地図を前記デジタイザから入力して、前記
記憶装置に記憶させる工程と、 （ｂ）前記記憶装置に記憶された前記第１の地図を前記
ディスプレイ上に表示させる工程と、 （ｃ）前記ディスプレイ上の第１の線分を指定する工程
と、 （ｄ）前記デジタイザ上に貼りつけられており、前記第
１の地図に対応する部分を有する第２の地図上の前記第
１の線分に対応する第２の線分上の点を指定する工程
と、 （ｅ）前記点に対応する前記第１の地図上の参照点と、
前記第１の線分との距離と角度を算出する工程と、 （ｆ）工程（ｃ）〜（ｅ）を繰り返して、複数組の第１
の線分と参照点との距離と角度を用いて、アフィン変換
用パラメータを算出する工程と、 （ｇ）前記記憶装置に記憶された前記第１の地図のデー
タを前記アフィン変換用パラメータを用いてアフィン変
換する工程と、 を具備する図形処理方法。 - 【請求項２】 デジタイザと、ディスプレイと、記憶装
置とを備えたシステムにおいて、 （ａ）第１の地図を前記デジタイザから入力して、前記
記憶装置に記憶させる工程と、 （ｂ）前記記憶装置に記憶された前記第１の地図を前記
ディスプレイ上に表示させる工程と、 （ｃ）前記ディスプレイ上の第１の点を指定する工程
と、 （ｄ）前記デジタイザ上に貼りつけられており、前記第
１の地図に対応する部分を有する第２の地図上の前記第
１の点に対応する第２の点を指定する工程と、 （ｅ）前記第１の点と前記第２の点を対角線とする矩形
の長辺および短辺のそれぞれの長さおよび傾きを距離と
角度として算出する工程と、 （ｆ）工程（ｃ）〜（ｅ）を繰り返して、複数組の第１
の点と第２の点との距離と角度を用いて、アフィン変換
用パラメータを算出する工程と、 （ｇ）前記記憶装置に記憶された前記第１の地図のデー
タを前記アフィン変換用パラメータを用いてアフィン変
換する工程と、 を具備する図形処理方法。 - 【請求項３】 デジタイザと、ディスプレイと、記憶装
置とを備えたシステムにおいて、 （ａ）第１の地図を前記デジタイザから入力して、前記
記憶装置に記憶させる工程と、 （ｂ）前記記憶装置に記憶された前記第１の地図を前記
ディスプレイ上に表示させる工程と、 （ｃ）前記ディスプレイ上の第１の線分または第３の点
を指定する工程と、 （ｄ）前記デジタイザ上に貼りつけられており、前記第
１の地図に対応する部分を有する第２の地図上の前記第
１の線分に対応する第２の線分上の点または、前記第３
の点に対応する第４の点を指定する工程と、 （ｅ）前記点に対応する前記第１の地図上の参照点と、
前記第１の線分との距離と角度を算出し、または前記第
３の点と前記第４の点を対角線とする矩形の長辺および
短辺のそれぞれの長さおよび傾きを距離と角度として算
出する工程と、 （ｆ）工程（ｃ）〜（ｅ）を繰り返して、複数組の第１
の線分と参照点との距離と角度または、第３の点と第４
の点との距離と角度を用いて、アフィン変換用パラメー
タを算出する工程と、 （ｇ）前記記憶装置に記憶された前記第１の地図のデー
タを前記アフィン変換用パラメータを用いてアフィン変
換する工程と、 を具備する図形処理方法。 - 【請求項４】 前記工程（ｅ）は、前記参照点から前記
第１の線分におろされた垂線の長さを距離とし、前記垂
線の水平軸に対する傾きを角度とすることを特徴とする
請求項１に記載された図形処理方法。 - 【請求項５】 前記工程（ｆ）は、距離が標準偏差から
大きくずれた組の距離と角度を除外して、アフィン変換
用パラメータを算出することを特徴とする請求項１から
請求項３までのいずれかに記載された図形処理方法。 - 【請求項６】 前記工程（ｆ）は、指定された線の種類
に対して重みづけを行ってアフィン変換用パラメータを
算出することを特徴とする請求項１から請求項３までの
いずれかに記載された図形処理方法。 - 【請求項７】 前記第１の地図の図郭の４隅を入力する
際に、その座標値を工程（ｆ）で算出されたアフィン変
換用パラメータを用いてアフィン変換する工程と、 アフィン変換された前記座標値を保持する工程と、 出図を行う場合、前記座標値で結ばれる直線を他の地図
データとともに表示させる工程と、 をさらに具備する請求項１から請求項３までのいずれか
に記載された図形処理方法。 - 【請求項８】 前記第１の地図は、国家基準点以外の基
準点を基に作成した大縮尺の地図であり、前記第２の地
図は、国家基準点を基準として作成した中縮尺の地図で
ある請求項１から請求項７のいずれかに記載された図形
処理方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP06292693A JP3298972B2 (ja) | 1993-02-28 | 1993-02-28 | 図形処理方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP06292693A JP3298972B2 (ja) | 1993-02-28 | 1993-02-28 | 図形処理方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH06251120A true JPH06251120A (ja) | 1994-09-09 |
| JP3298972B2 JP3298972B2 (ja) | 2002-07-08 |
Family
ID=13214376
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP06292693A Expired - Lifetime JP3298972B2 (ja) | 1993-02-28 | 1993-02-28 | 図形処理方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3298972B2 (ja) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH10307907A (ja) * | 1997-05-08 | 1998-11-17 | Hitachi Ltd | リンク情報保守管理方法 |
| JP2004101780A (ja) * | 2002-09-09 | 2004-04-02 | Hitachi Ltd | 位置座標補正処理装置及び処理方法 |
| JP2004355505A (ja) * | 2003-05-30 | 2004-12-16 | Hitachi Ltd | 時空間データ管理装置及び方法 |
| JP2006053179A (ja) * | 2004-08-09 | 2006-02-23 | Hitachi Ltd | 地理情報処理装置及び方法 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0561407A (ja) * | 1991-08-30 | 1993-03-12 | Tokyo Gas Co Ltd | 図形処理装置 |
-
1993
- 1993-02-28 JP JP06292693A patent/JP3298972B2/ja not_active Expired - Lifetime
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0561407A (ja) * | 1991-08-30 | 1993-03-12 | Tokyo Gas Co Ltd | 図形処理装置 |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH10307907A (ja) * | 1997-05-08 | 1998-11-17 | Hitachi Ltd | リンク情報保守管理方法 |
| JP2004101780A (ja) * | 2002-09-09 | 2004-04-02 | Hitachi Ltd | 位置座標補正処理装置及び処理方法 |
| JP2004355505A (ja) * | 2003-05-30 | 2004-12-16 | Hitachi Ltd | 時空間データ管理装置及び方法 |
| JP4630528B2 (ja) * | 2003-05-30 | 2011-02-09 | 株式会社日立製作所 | 時空間データ管理装置 |
| JP2006053179A (ja) * | 2004-08-09 | 2006-02-23 | Hitachi Ltd | 地理情報処理装置及び方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP3298972B2 (ja) | 2002-07-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US8072450B2 (en) | System and method for measuring a three-dimensional object | |
| JP5324240B2 (ja) | 路面標示地図生成方法及び路面標示地図生成装置 | |
| JP6591131B1 (ja) | 構造物計測装置および構造物計測方法 | |
| JP5363752B2 (ja) | 路面標示地図生成方法 | |
| JPH0769958B2 (ja) | 画像変換処理方法 | |
| JP4535337B2 (ja) | 登記異動情報反映装置および方法 | |
| JP2859004B2 (ja) | 図形処理装置 | |
| JP3298972B2 (ja) | 図形処理方法 | |
| WO2009104218A1 (ja) | 地図表示装置 | |
| JP3298973B2 (ja) | 図形処理方法 | |
| JP3535595B2 (ja) | 図形処理方法 | |
| JP2913814B2 (ja) | 地図データベース作成・修正システム | |
| JP3417998B2 (ja) | トラバース変換機能を有する登記測量システム | |
| JP3402651B2 (ja) | 地図データベース検定方法 | |
| JP2005070373A (ja) | 地図座標変換プログラムおよび地図座標変換装置 | |
| JP3535594B2 (ja) | 図形処理方法 | |
| JP2006293809A (ja) | 図形データ処理装置 | |
| CN108614552B (zh) | 一种基于离散位姿的路径规划方法 | |
| JPH09138870A (ja) | 画像処理装置及び方法 | |
| JP4153323B2 (ja) | 写真測量における計測点の対応付け方法及び装置 | |
| JPH09179967A (ja) | 地図データ修正方法 | |
| JPH1166130A (ja) | 二次元データの三次元化処理システム及び二次元データの三次元化処理方法 | |
| JPH08190577A (ja) | 図形処理装置 | |
| JP3087992B2 (ja) | 建物図作成方法及び装置 | |
| Holmes et al. | The MAPEDIT system for automatic map digitization |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090419 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090419 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100419 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110419 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120419 Year of fee payment: 10 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130419 Year of fee payment: 11 |
|
| EXPY | Cancellation because of completion of term |