JP3298973B2 - 図形処理方法 - Google Patents
図形処理方法Info
- Publication number
- JP3298973B2 JP3298973B2 JP06292793A JP6292793A JP3298973B2 JP 3298973 B2 JP3298973 B2 JP 3298973B2 JP 06292793 A JP06292793 A JP 06292793A JP 6292793 A JP6292793 A JP 6292793A JP 3298973 B2 JP3298973 B2 JP 3298973B2
- Authority
- JP
- Japan
- Prior art keywords
- map
- storage device
- point
- display
- 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.)
- Expired - Lifetime
Links
Landscapes
- Instructional Devices (AREA)
- Processing Or Creating Images (AREA)
Description
【０００１】
【産業上の利用分野】本発明は、地図データをデータベ
ース化して利用するマッピングシステムに用いられる図
形処理方法に関するものである。
ース化して利用するマッピングシステムに用いられる図
形処理方法に関するものである。
【０００２】
【従来の技術】一般にマッピングシステムでは、地図デ
ータベースとして国家基準点を利用して作成した地図を
用いることが望まれる。
ータベースとして国家基準点を利用して作成した地図を
用いることが望まれる。
【０００３】しかしながら、国道台帳図等においては精
度の高い測量方法に基づいて作成されており、部分的に
は国家基準点を利用して作成した地図よりも高い精度を
有している。
度の高い測量方法に基づいて作成されており、部分的に
は国家基準点を利用して作成した地図よりも高い精度を
有している。
【０００４】
【発明が解決しようとする課題】このため、このような
部分的に高い精度を有する地図をマッピングシステムに
取り入れることが考えられる。しかしながら、国道台帳
図は道路に沿って細長い形態であり、長手方向が必ずし
も北を指しておらず、かかる国道台帳図を地図データベ
ース化する場合、国家基準点を利用して作成した地図と
整合させるためには、磁北を合わせた後、道路管理境界
で地図同志の整合を図る必要がある。
部分的に高い精度を有する地図をマッピングシステムに
取り入れることが考えられる。しかしながら、国道台帳
図は道路に沿って細長い形態であり、長手方向が必ずし
も北を指しておらず、かかる国道台帳図を地図データベ
ース化する場合、国家基準点を利用して作成した地図と
整合させるためには、磁北を合わせた後、道路管理境界
で地図同志の整合を図る必要がある。
【０００５】本発明は、このような問題に鑑みてなされ
たもので、その目的とするところは、一方の地図を高精
度でかつ効率よく他の地図に整合させることができる図
形処理方法を提供することにある。
たもので、その目的とするところは、一方の地図を高精
度でかつ効率よく他の地図に整合させることができる図
形処理方法を提供することにある。
【０００６】
【課題を解決するための手段】図１は、本発明の原理説
明図である。図１において１はホストコンピュータ、３
は地図データベース、５はアプリケーションプログラ
ム、７はワークファイル、９はディスプレイ、１１はデ
ジタイザ、１３は指示具、１４は国道台帳図である。
ホストコンピュータ１はアプリケーションプログラム５
にしたがって各種処理を行う。地図データベース３に
は、たとえば国家基準点を利用して作成された市・区道
台帳図データが格納される。ディスプレイ９には地図デ
ータベース３内の地図と、デジタイザ１１に貼られ、指
示具１３で入力された国道台帳図１４とが重畳されて表
示される。この時、両方の地図が識別できるように色が
変えられて表示される。図１においては地図データベー
ス３の地図データが実線として表示され、国道台帳図１
４のほうが破線で表示される。
明図である。図１において１はホストコンピュータ、３
は地図データベース、５はアプリケーションプログラ
ム、７はワークファイル、９はディスプレイ、１１はデ
ジタイザ、１３は指示具、１４は国道台帳図である。
ホストコンピュータ１はアプリケーションプログラム５
にしたがって各種処理を行う。地図データベース３に
は、たとえば国家基準点を利用して作成された市・区道
台帳図データが格納される。ディスプレイ９には地図デ
ータベース３内の地図と、デジタイザ１１に貼られ、指
示具１３で入力された国道台帳図１４とが重畳されて表
示される。この時、両方の地図が識別できるように色が
変えられて表示される。図１においては地図データベー
ス３の地図データが実線として表示され、国道台帳図１
４のほうが破線で表示される。
【０００７】デジタイザ１１には、たとえば国道台帳図
１４が貼られる。この地図はディスプレイ９に表示され
た地図に記載された領域に含まれる。デジタイザ１１に
おいて指示具１３により地図データを入力することがで
きる。指示具１３を動かすと、この動きにつれて、ディ
スプレイ９上でカーソルが動き、この指示具１３をクリ
ックすることによりディスプレイ９上の点、あるいは線
分を指定することができる。
１４が貼られる。この地図はディスプレイ９に表示され
た地図に記載された領域に含まれる。デジタイザ１１に
おいて指示具１３により地図データを入力することがで
きる。指示具１３を動かすと、この動きにつれて、ディ
スプレイ９上でカーソルが動き、この指示具１３をクリ
ックすることによりディスプレイ９上の点、あるいは線
分を指定することができる。
【０００８】Ａはアプリケーションプログラム５の内容
およびオペレータの操作を示すフローチャートである。
まずオペレータは国道台帳図１４をデジタイザ１１から
入力し、ワークファイル７に保管する（ステップ２
１）。そして国道台帳図１４の方位を入力する（ステッ
プ２２）。オペレータは国道台帳図１４に対応する地図
データベース３内の国土基本図にかかる地図をディスプ
レイ９に表示させる旨の指示を行う（ステップ２３）。
ディスプレイ９には実線で示されるような地図が表示さ
れる。
およびオペレータの操作を示すフローチャートである。
まずオペレータは国道台帳図１４をデジタイザ１１から
入力し、ワークファイル７に保管する（ステップ２
１）。そして国道台帳図１４の方位を入力する（ステッ
プ２２）。オペレータは国道台帳図１４に対応する地図
データベース３内の国土基本図にかかる地図をディスプ
レイ９に表示させる旨の指示を行う（ステップ２３）。
ディスプレイ９には実線で示されるような地図が表示さ
れる。
【０００９】ホストコンピュータ１はステップ２２で入
力された方位に基づいて、ワークファイル７内のデータ
の方位変換を行う（ステップ２４）。そして変換された
データをディスプレイ９上に色を変えて表示する（ステ
ップ２５）。図１の場合、このデータは破線で示してあ
る。すなわちデジタイザ１１に貼られた国道台帳図１４
の磁北が上を指すように変換されてディスプレイ９上に
色を変えて表示される。 オペレータはデジタイザ１１
の指示具１３を用いて対応する線分ｇ１、ｇ２を指定す
る（ステップ２６）。そしてホストコンピュータ１は２
つの線分ｇ１、ｇ２間の距離と角度を算出する（ステッ
プ２７）。オペレータはステップ２６、ステップ２７の
動作を繰り返して線分ｇ１、ｇ２を異ならしめて、線分
ｇ１、ｇ２の組をＮ組指定する（ステップ２８）。この
ように線分ｇ１、ｇ２が複数組指定されるとホストコン
ピュータ１はステップ２７で算出された複数組の距離と
角度を用いてアフィン変換用パラメータを算出し（ステ
ップ２９）、ワークファイル７内の国道台帳図１４に係
る地図のデータをアフィン変換用パラメータを用いてア
フィン変換する（ステップ３０）。アフィン変換された
地図データは地図データベース３に格納される。
力された方位に基づいて、ワークファイル７内のデータ
の方位変換を行う（ステップ２４）。そして変換された
データをディスプレイ９上に色を変えて表示する（ステ
ップ２５）。図１の場合、このデータは破線で示してあ
る。すなわちデジタイザ１１に貼られた国道台帳図１４
の磁北が上を指すように変換されてディスプレイ９上に
色を変えて表示される。 オペレータはデジタイザ１１
の指示具１３を用いて対応する線分ｇ１、ｇ２を指定す
る（ステップ２６）。そしてホストコンピュータ１は２
つの線分ｇ１、ｇ２間の距離と角度を算出する（ステッ
プ２７）。オペレータはステップ２６、ステップ２７の
動作を繰り返して線分ｇ１、ｇ２を異ならしめて、線分
ｇ１、ｇ２の組をＮ組指定する（ステップ２８）。この
ように線分ｇ１、ｇ２が複数組指定されるとホストコン
ピュータ１はステップ２７で算出された複数組の距離と
角度を用いてアフィン変換用パラメータを算出し（ステ
ップ２９）、ワークファイル７内の国道台帳図１４に係
る地図のデータをアフィン変換用パラメータを用いてア
フィン変換する（ステップ３０）。アフィン変換された
地図データは地図データベース３に格納される。
【００１０】
【作用】本発明では、デジタイザ１１に貼られた国道台
帳図１４は横長であり、磁北を上にしてデジタイザ１１
に貼ることが困難である。このため国道台帳図１４をデ
ジタイザ１１で入力した後、方位ベクトルを入力し、こ
の方位ベクトルに従って国道台帳図１４のデータを変換
し、磁北が上に来るようにディスプレイ９上に表示させ
る。ディスプレイ９上では、地図データベース３内の市
・区道台帳図に係る地図データと国道台帳図１４にかか
る地図データが重畳されて色違いで表示される。そし
て、この２つの地図データの整合を図るためにアフィン
変換用パラメータが算出され、ワークファイル７内の国
道台帳図１４にかかる地図データがアフィン変換され
る。
帳図１４は横長であり、磁北を上にしてデジタイザ１１
に貼ることが困難である。このため国道台帳図１４をデ
ジタイザ１１で入力した後、方位ベクトルを入力し、こ
の方位ベクトルに従って国道台帳図１４のデータを変換
し、磁北が上に来るようにディスプレイ９上に表示させ
る。ディスプレイ９上では、地図データベース３内の市
・区道台帳図に係る地図データと国道台帳図１４にかか
る地図データが重畳されて色違いで表示される。そし
て、この２つの地図データの整合を図るためにアフィン
変換用パラメータが算出され、ワークファイル７内の国
道台帳図１４にかかる地図データがアフィン変換され
る。
【００１１】
【実施例】以下、図面に基づいて本発明の実施例を詳細
に説明する。図２は、本発明の１実施例に係る図形処理
方法に用いられる装置の構成を示す図である。１はホス
トコンピュータ、３は地図データベース、５はアプリケ
ーションプログラム、９はディスプレイ、１１はデジタ
イザ、１３は指示具、１５はキーボード、１７はプリン
タである。
に説明する。図２は、本発明の１実施例に係る図形処理
方法に用いられる装置の構成を示す図である。１はホス
トコンピュータ、３は地図データベース、５はアプリケ
ーションプログラム、９はディスプレイ、１１はデジタ
イザ、１３は指示具、１５はキーボード、１７はプリン
タである。
【００１２】図３、図４は本実施例の処理を示すフロー
チャートである。まずオペレータは国道台帳図１４をデ
ジタイザ１１に貼りつけ国道部分の道路、歩道、道路施
設を入力してワークファイル７内の所定の領域（以下、
ワークファイルＷ１）に保管する（ステップ３０１）。
チャートである。まずオペレータは国道台帳図１４をデ
ジタイザ１１に貼りつけ国道部分の道路、歩道、道路施
設を入力してワークファイル７内の所定の領域（以下、
ワークファイルＷ１）に保管する（ステップ３０１）。
【００１３】図６は、このときのディスプレイ９および
デジタイザ１１を示す。デジタイザ１１には横長の国道
台帳図１４が貼り付けられている。国道台帳図１４は横
長であり、その磁北を上にしてデジタイザ１１上に貼り
付けることは難しいので、このように横長の状態でデジ
タイザ１１に貼り付けられる。
デジタイザ１１を示す。デジタイザ１１には横長の国道
台帳図１４が貼り付けられている。国道台帳図１４は横
長であり、その磁北を上にしてデジタイザ１１上に貼り
付けることは難しいので、このように横長の状態でデジ
タイザ１１に貼り付けられる。
【００１４】さらにオペレータは仮基準位置Ａおよび方
位ベクトルＢ、Ｃを入力し、ワークファイルＷ２に保管
する（ステップ３０２）。すなわちオペレータは指示具
１３を用いて、国道台帳図１４上の点Ａ、Ｂ、Ｃを指定
する。そしてこれらの点の座標（ｘA ，ｙA ）、
（ｘB ，ｙB ）、（ｘC ，ｙC ）がワークファイルＷ２
に保管される。
位ベクトルＢ、Ｃを入力し、ワークファイルＷ２に保管
する（ステップ３０２）。すなわちオペレータは指示具
１３を用いて、国道台帳図１４上の点Ａ、Ｂ、Ｃを指定
する。そしてこれらの点の座標（ｘA ，ｙA ）、
（ｘB ，ｙB ）、（ｘC ，ｙC ）がワークファイルＷ２
に保管される。
【００１５】つぎにオペレータはこの国道台帳図１４に
相当する区道台帳図のメッシュ番号をキーボード１５か
ら入力して、地図データベース３内の地図データをディ
スプレイ９上に表示させる（ステップ３０３）。図７
は、このとき指定される地図データベース３内の地図デ
ータ４１を示す。そして図６に示されるように、この地
図データ４１がディスプレイ９上に表示される。
相当する区道台帳図のメッシュ番号をキーボード１５か
ら入力して、地図データベース３内の地図データをディ
スプレイ９上に表示させる（ステップ３０３）。図７
は、このとき指定される地図データベース３内の地図デ
ータ４１を示す。そして図６に示されるように、この地
図データ４１がディスプレイ９上に表示される。
【００１６】オペレータはディスプレイ９上で点Ａに相
当する点Ａ´を入力する（ステップ３０４）。すなわち
指示具１３を動かしてディスプレイ９上のカーソルを動
かし、カーソルが点Ａ´上に位置した状態で指示具１３
をクリックして点Ａ´の座標（ＸA ´，ＹA ´）を入力
する。
当する点Ａ´を入力する（ステップ３０４）。すなわち
指示具１３を動かしてディスプレイ９上のカーソルを動
かし、カーソルが点Ａ´上に位置した状態で指示具１３
をクリックして点Ａ´の座標（ＸA ´，ＹA ´）を入力
する。
【００１７】ホストコンピュータ１はワークファイルＷ
２に格納されたデータおよび点Ａ´のデータを用いて、
ワークファイルＷ１のデータの変換係数を次式に従って
求める（ステップ３０５）。
２に格納されたデータおよび点Ａ´のデータを用いて、
ワークファイルＷ１のデータの変換係数を次式に従って
求める（ステップ３０５）。
【００１８】θ＝π／２−tan -1［（ｙB −ｙC ）／
（ｘB −ｘC ）］ ここで、θは直線Ｂ、Ｃの傾きを表わしている。
（ｘB −ｘC ）］ ここで、θは直線Ｂ、Ｃの傾きを表わしている。
【００１９】つぎにホストコンピュータ１はワークファ
イルＷ１に格納された座標値（ｘ，ｙ）を次式に従って
変換し、ワークファイルＷ３に保管する（ステップ３０
６） 。 Ｘ´＝cos θ・（ｘ−ｘA ）＋sin θ・（ｙ−ｙ
A ）＋ＸA ´ Ｙ´＝−sin θ・（ｘ−ｘA ）＋cos θ・（ｙ−ｙA ）
＋ＹA ´ ここでＸ´、Ｙ´は変換された座標である。
イルＷ１に格納された座標値（ｘ，ｙ）を次式に従って
変換し、ワークファイルＷ３に保管する（ステップ３０
６） 。 Ｘ´＝cos θ・（ｘ−ｘA ）＋sin θ・（ｙ−ｙ
A ）＋ＸA ´ Ｙ´＝−sin θ・（ｘ−ｘA ）＋cos θ・（ｙ−ｙ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 ・Ｙ3 とすると Ｌ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 ´）を入
力する。
当する点Ａ´を入力する（ステップ１７０４）。すなわ
ち指示具１３を動かしてディスプレイ９上のカーソルを
動かし、カーソルが点Ａ´上に位置した状態で指示具１
３をクリックして点Ａ´の座標（ＸA ´，ＹA ´）を入
力する。
【００６４】ホストコンピュータ１はワークファイルＷ
１２に格納されたデータおよび点Ａ´のデータを用い
て、ワークファイルＷ１１、Ｗ１１´のデータの変換係
数を次式に従って求める（ステップ１７０５）。
１２に格納されたデータおよび点Ａ´のデータを用い
て、ワークファイルＷ１１、Ｗ１１´のデータの変換係
数を次式に従って求める（ステップ１７０５）。
【００６５】θ＝π／２−tan -1［（ｙB −ｙC ）／
（ｘB −ｘC ）］ すなわちθは直線Ｂ、Ｃの傾きを表わしている。
（ｘB −ｘC ）］ すなわちθは直線Ｂ、Ｃの傾きを表わしている。
【００６６】つぎにホストコンピュータ１はワークファ
イルＷ１１、Ｗ１１´に格納された座標値（ｘ，ｙ）を
次式に従って変換し、ワークファイルＷ１３、Ｗ１３´
に保管する（ステップ１７０６）。
イルＷ１１、Ｗ１１´に格納された座標値（ｘ，ｙ）を
次式に従って変換し、ワークファイルＷ１３、Ｗ１３´
に保管する（ステップ１７０６）。
【００６７】Ｘ´＝cos θ・（ｘ−ｘA ）＋sin θ・
（ｙ−ｙA ）＋ＸA ´ Ｙ´＝−sin θ・（ｘ−ｘA ）＋cos θ・（ｙ−ｙA ）
＋ＹA ´ ここでＸ´、Ｙ´は変換された座標である。
（ｙ−ｙA ）＋ＸA ´ Ｙ´＝−sin θ・（ｘ−ｘA ）＋cos θ・（ｙ−ｙA ）
＋ＹA ´ ここでＸ´、Ｙ´は変換された座標である。
【００６８】つぎにワークファイルＷ１３、Ｗ１３´の
データをステップ１７０３の表示内容と異なる色でディ
スプレイ９上に表示させる（ステップ１７０７）。
データをステップ１７０３の表示内容と異なる色でディ
スプレイ９上に表示させる（ステップ１７０７）。
【００６９】ホストコンピュータは組数を示すパラメー
タｉを初期化し（ステップ１７０８）、デジタイザ１１
の指示具１３を用いてワークファイルＷ１３´の線分ｇ
１と、これに対応する線分ｇ２を指定する（ステップ１
７０９）。このように対応する線分を指定する際、ワー
クファイルＷ１３´のデータを用いることにより、無修
正部分のデータを用いて、以下に述べるアフィン変換用
パラメータを求める。つぎにホストコンピュータ１は、
線分ｇ１と線分ｇ２との距離Ｌi と角度θiを計算し、
線分ｇ２上の点Ｈの座標（Ｘi ，Ｙi ）と共に格納する
（ステップ１７１０）。
タｉを初期化し（ステップ１７０８）、デジタイザ１１
の指示具１３を用いてワークファイルＷ１３´の線分ｇ
１と、これに対応する線分ｇ２を指定する（ステップ１
７０９）。このように対応する線分を指定する際、ワー
クファイルＷ１３´のデータを用いることにより、無修
正部分のデータを用いて、以下に述べるアフィン変換用
パラメータを求める。つぎにホストコンピュータ１は、
線分ｇ１と線分ｇ２との距離Ｌi と角度θiを計算し、
線分ｇ２上の点Ｈの座標（Ｘi ，Ｙi ）と共に格納する
（ステップ１７１０）。
【００７０】つぎにパラメータｉを「１」増加させ（ス
テップ１７１１）、オペレータが終了の指示を行うまで
ステップ１７０９からステップ１７１１までの処理が繰
り返される（ステップ１７１２）。すなわち異なる線分
ｇ１、線分ｇ２の組が複数組指定され、各組ごとに、そ
の距離Ｌi と角度θi が計算される。
テップ１７１１）、オペレータが終了の指示を行うまで
ステップ１７０９からステップ１７１１までの処理が繰
り返される（ステップ１７１２）。すなわち異なる線分
ｇ１、線分ｇ２の組が複数組指定され、各組ごとに、そ
の距離Ｌi と角度θi が計算される。
【００７１】このときの組の数をＮに代入し（ステップ
１７１３）、アフィン変換用パラメータを計算する（ス
テップ１７１４）。
１７１３）、アフィン変換用パラメータを計算する（ス
テップ１７１４）。
【００７２】つぎにオペレータはディスプレイ９上から
ワークファイルＷ１３、Ｗ１３´のデータを消去し（ス
テップ１７１５）、ワークファイルＷ１３、Ｗ１３´内
のデータをステップ３１４で求められたアフィン変換用
パラメータを用いて変換し、ワークファイルＷ１４、Ｗ
１４´に格納すると共にディスプレイ９上に表示する
（ステップ１７１６）。
ワークファイルＷ１３、Ｗ１３´のデータを消去し（ス
テップ１７１５）、ワークファイルＷ１３、Ｗ１３´内
のデータをステップ３１４で求められたアフィン変換用
パラメータを用いて変換し、ワークファイルＷ１４、Ｗ
１４´に格納すると共にディスプレイ９上に表示する
（ステップ１７１６）。
【００７３】すなわちワークファイルＷ１３、Ｗ１３´
内のデータの座標をＸ、Ｙとし、アフィン変換用パラメ
ータを用いて変換した座標をＸ＊、Ｙ＊とすると、 Ｘ＊＝Ｘ＋（ＡＸ＋ＢＹ＋Ｃ） Ｙ＊＝Ｙ＋（ＤＸ＋ＥＹ＋Ｆ）となる。この座標変換を
ワークファイルＷ１３、Ｗ１３´内にある地図データに
ついて行う。
内のデータの座標をＸ、Ｙとし、アフィン変換用パラメ
ータを用いて変換した座標をＸ＊、Ｙ＊とすると、 Ｘ＊＝Ｘ＋（ＡＸ＋ＢＹ＋Ｃ） Ｙ＊＝Ｙ＋（ＤＸ＋ＥＹ＋Ｆ）となる。この座標変換を
ワークファイルＷ１３、Ｗ１３´内にある地図データに
ついて行う。
【００７４】そしてワークファイルＷ１４と地図データ
ベース３のデータについて不要な線等を消去し、線デー
タ等の追加、補正等の編集を行い（ステップ１７１
７）、編集されたデータを地図データベース３に格納す
る（ステップ１７１８）。
ベース３のデータについて不要な線等を消去し、線デー
タ等の追加、補正等の編集を行い（ステップ１７１
７）、編集されたデータを地図データベース３に格納す
る（ステップ１７１８）。
【００７５】なお第７実施例では道路の出来型図を対象
にしたが、ガス管、水道管等の設備データの出来型図を
第７実施例と同様に処理することができる。また、この
第７実施例についても、第２実施例から第６実施例を応
用することができる。
にしたが、ガス管、水道管等の設備データの出来型図を
第７実施例と同様に処理することができる。また、この
第７実施例についても、第２実施例から第６実施例を応
用することができる。
【００７６】
【発明の効果】以上、詳細に説明したように本発明によ
れば、一方の地図を高精度でかつ効率よく他の地図に整
合させることができる。
れば、一方の地図を高精度でかつ効率よく他の地図に整
合させることができる。
【図１】 本発明の原理説明図
【図２】 本発明の１実施例に係る図形処理方法を実行
するための装置構成図
するための装置構成図
【図３】 実施例の動作を示すフローチャート
【図４】 実施例の動作を示すフローチャート
【図５】 実施例の動作を示すフローチャート
【図６】 ディスプレイ９の表示内容とデジタイザ１１
上に貼られた国道台帳図を示す図
上に貼られた国道台帳図を示す図
【図７】 区道台帳図４１を示す図
【図８】 区道台帳図４１の領域４３を示す図
【図９】 線分ｇと点Ｈとの位置関係を示す図
【図１０】 アフィン変換を説明する図
【図１１】 目視によった場合と、本実施例による方法
を用いて国道台帳図を区道台帳図に整合させた場合の誤
差の比較を示す図
を用いて国道台帳図を区道台帳図に整合させた場合の誤
差の比較を示す図
【図１２】 アフィン変換用パラメータの決定要素を自
動的に除外する場合の動作を示すフローチャート
動的に除外する場合の動作を示すフローチャート
【図１３】 アフィン変換用パラメータの決定要素を半
自動的に除外する場合の動作を示すフローチャート
自動的に除外する場合の動作を示すフローチャート
【図１４】 区域別にアフィン変換用パラメータを決定
する場合の動作を示すフローチャート
する場合の動作を示すフローチャート
【図１５】 各地図上の基準点が、それぞれ指定された
場合の処理を説明する図
場合の処理を説明する図
【図１６】 第７実施例の動作を示すフローチャート
【図１７】 第７実施例の動作を示すフローチャート
【図１８】 第７実施例の動作を示すフローチャート
１………ホストコンピュータ ３………地図データベース ５………アプリケーションプログラム ７………ワーク領域 ９………ディスプレイ １１………デジタイザ １５………キーボード １７………プリンタ
Claims (10)
- 【請求項１】 デジタイザと、ディスプレイと、記憶装
置とを備えたシステムにおいて、 （ａ）第１の地図を前記デジタイザから入力して、前記
記憶装置に記憶させる工程と、 （ｂ）記憶された前記第１の地図を方位変換し、前記記
憶装置に記憶させる工程と、 （ｃ）前記記憶装置に記憶された方位変換された前記第
１の地図を前記ディスプレイ上に表示させる工程と、 （ｄ）前記ディスプレイ上の第１の線分を指定する工程
と、 （ｅ）前記記憶装置に記憶され、前記第１の地図に対応
する部分を有する第２の地図データを前記ディスプレイ
に表示させる工程と、 （ｆ）前記ディスプレイに表示された前記第１の地図上
の第１の線分と、この第１の線分に対応する前記第２の
地図上の第２の線分とを指定する工程と、 （ｇ）前記第１の線分と、前記第２の線分との距離と角
度を算出する工程と、 （ｈ）工程（ｄ）〜（ｇ）を繰り返して、複数組の第１
の線分と第２の線分との距離と角度を用いて、アフィン
変換用パラメータを算出する工程と、 （ｉ）前記記憶装置に記憶された前記第１の地図のデー
タを前記アフィン変換用パラメータを用いてアフィン変
換する工程と、 を具備する図形処理方法。 - 【請求項２】 デジタイザと、ディスプレイと、記憶装
置とを備えたシステムにおいて、 （ａ）第１の地図を前記デジタイザから入力して、前記
記憶装置に記憶させる工程と、 （ｂ）記憶された前記第１の地図を方位変換し、前記記
憶装置に記憶させる工程と、 （ｃ）前記記憶装置に記憶された方位変換された前記第
１の地図を前記ディスプレイ上に表示させる工程と、 （ｄ）前記ディスプレイ上の第１の地図上の第１の点を
指定する工程と、 （ｅ）前記記憶装置に記憶され、前記第１の地図に対応
する部分を有する第２の地図データを前記ディスプレイ
に表示させる工程と、 （ｆ）前記ディスプレイ上で、前記第１の点に対応する
前記第２の地図上の第２の点を指定する工程と、 （ｇ）前記第１の点と前記第２の点を対角線とする矩形
の長辺および短辺のそれぞれの長さおよび傾きを距離と
角度として算出する工程と、 （ｈ）工程（ｄ）〜（ｇ）を繰り返して、複数組の第１
の点と第２の点との距離と角度を用いて、アフィン変換
用パラメータを算出する工程と、 （ｉ）前記記憶装置に記憶された前記第１の地図のデー
タを前記アフィン変換用パラメータを用いてアフィン変
換する工程と、 を具備する図形処理方法。 - 【請求項３】 デジタイザと、ディスプレイと、記憶装
置とを備えたシステムにおいて、 （ａ）第１の地図を前記デジタイザから入力して、前記
記憶装置に記憶させる工程と、 （ｂ）記憶された前記第１の地図を方位変換し、前記記
憶装置に記憶させる工程と、 （ｃ）前記記憶装置に記憶された方位変換された前記第
１の地図を前記ディスプレイ上に表示させる工程と、 （ｄ）前記ディスプレイ上の第１の線分または、第３の
点を指定する工程と、 （ｅ）前記記憶装置に記憶され、前記第１の地図に対応
する部分を有する第２の地図データを前記ディスプレイ
に表示させる工程と、 （ｆ）前記ディスプレイ上で第２の地図上の前記第１の
線分に対応する第２の線分または、前記第３の点に対応
する第４の点を指定する工程と、 （ｇ）前記第１の線分と前記第２の線分との距離と角度
を算出し、または前記第３の点と前記第４の点を対角線
とする矩形の長辺および短辺のそれぞれの長さおよび傾
きを距離と角度として算出する工程と、 （ｈ）工程（ｄ）〜（ｇ）を繰り返して、複数組の第１
の線分と参照点との距離と角度または、第３の点と第４
の点との距離と角度を用いて、アフィン変換用パラメー
タを算出する工程と、 （ｉ）前記記憶装置に記憶された前記第１の地図のデー
タを前記アフィン変換用パラメータを用いてアフィン変
換する工程と、 を具備する図形処理方法。 - 【請求項４】 前記工程（ｇ）は、前記第１の線分から
前記第２の線分におろされた垂線の長さを距離とし、前
記垂線の水平軸に対する傾きを角度とすることを特徴と
する請求項１および請求項３記載の図形処理方法。 - 【請求項５】 前記工程（ｈ）は、距離が標準偏差から
大きくずれた組の距離と角度を除外して、アフィン変換
用パラメータを算出することを特徴とする請求項１から
３のいずれかに記載された図形処理方法。 - 【請求項６】 前記工程（ｈ）は、指定された区域別に
アフィンは変換用パラメータを算出することを特徴とす
る請求項１から３のいずれかに記載された図形処理方
法。 - 【請求項７】 前記工程（ｈ）は、指定された線の種類
に対して重みづけを行ってアフィン変換用パラメータを
算出することを特徴とする請求項１から３のいずれかに
記載された図形処理方法。 - 【請求項８】 前記第１の地図は任意図郭の道路台帳図
であり、前記第２の地図は国家基準点を基準として作成
された道路台帳図である請求項１から７のいずれかに記
載された図形処理方法。 - 【請求項９】 前記第１の地図は道路の出来型図であ
り、修正部分と無修正部分とを有しており、無修正部分
を基にして、アフィン変換用パラメータを求めることを
特徴とする請求項１から７のいずれかに記載された図形
処理方法。 - 【請求項１０】 前記第１の地図は設備の出来型図であ
り、道路データを基にして、アフィン変換用パラメータ
を求めることを特徴とする請求項１から７のいずれかに
記載された図形処理方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP06292793A JP3298973B2 (ja) | 1993-02-28 | 1993-02-28 | 図形処理方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP06292793A JP3298973B2 (ja) | 1993-02-28 | 1993-02-28 | 図形処理方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH06251121A JPH06251121A (ja) | 1994-09-09 |
| JP3298973B2 true JP3298973B2 (ja) | 2002-07-08 |
Family
ID=13214408
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP06292793A Expired - Lifetime JP3298973B2 (ja) | 1993-02-28 | 1993-02-28 | 図形処理方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3298973B2 (ja) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3367548B2 (ja) * | 1995-03-06 | 2003-01-14 | 日本電信電話株式会社 | 道路自動認識方法及びその装置並びに地図情報処理方法及びその装置 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2859004B2 (ja) * | 1991-08-30 | 1999-02-17 | 東京瓦斯株式会社 | 図形処理装置 |
-
1993
- 1993-02-28 JP JP06292793A patent/JP3298973B2/ja not_active Expired - Lifetime
Also Published As
| Publication number | Publication date |
|---|---|
| JPH06251121A (ja) | 1994-09-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4773973B2 (ja) | 単一像から、３ｄ物体の測定、モデリングおよびマッピングのためのシステム、コンピュータプログラムおよび方法 | |
| JP2642070B2 (ja) | 四角形メッシュの生成方法及びシステム | |
| EP0619554A2 (en) | A method for geo-registration of imported bit-mapped spatial data | |
| EP2251833A2 (en) | Pitch determination systems and methods for aerial roof estimation | |
| JP6591131B1 (ja) | 構造物計測装置および構造物計測方法 | |
| JP3803113B2 (ja) | データのモザイクを構築するための装置および方法 | |
| JP2859004B2 (ja) | 図形処理装置 | |
| JP3298972B2 (ja) | 図形処理方法 | |
| JP3298973B2 (ja) | 図形処理方法 | |
| JP3535595B2 (ja) | 図形処理方法 | |
| JP3402651B2 (ja) | 地図データベース検定方法 | |
| JP3367548B2 (ja) | 道路自動認識方法及びその装置並びに地図情報処理方法及びその装置 | |
| JP2006293809A (ja) | 図形データ処理装置 | |
| JP5199860B2 (ja) | 地図データの編集方法 | |
| JP3535594B2 (ja) | 図形処理方法 | |
| JP4643507B2 (ja) | 地図位置合せ装置および方法 | |
| JP4164745B2 (ja) | ３次元測量支援方法及びその装置 | |
| JP3535597B2 (ja) | 図形処理装置 | |
| JP3617939B2 (ja) | デジタル地図における登録情報の更新方法 | |
| JP3259981B2 (ja) | 図形データの一部修正システムおよび図形データの一部修正方法 | |
| JPH11283004A (ja) | 地籍図データ管理方法 | |
| JP2002024336A (ja) | 送電線図面管理方法及びシステム | |
| CN108614552B (zh) | 一种基于离散位姿的路径规划方法 | |
| JP3087992B2 (ja) | 建物図作成方法及び装置 | |
| JP2001282801A (ja) | 設備データ変換システム |
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 |