(57)【特許請求の範囲】
【請求項１】
地理画像を解析し、道路地図データを作成する道路地図データ作成装置であって、
前記地理画像データから抽出された前記地理画像の領域において複数の開始点をまんべんなく与え、与えられた該複数の開始点のそれぞれから道路を追跡していくことにより得られる追跡線に沿って、各開始点における座標及びそれと結合しているノードの座標とにより定義される道路のベクトルデータを抽出する道路追跡処理部と、
前記地理画像内の各道路領域において、前記道路のベクトルデータのうち重複データと推定されるベクトルデータを統合する道路統合処理部と、
前記地理画像データ内の各道路領域の幅員を統合された前記ベクトルデータの数に基づいて推定する幅員推定処理部と、
を備えることを特徴とする道路地図データ作成装置。
【請求項２】
前記幅員推定処理部は、統合された前記ベクトルデータの数が多いほど幅員を大きくし、統合した道路の前記ベクトルデータの長さが短いほど幅員を大きくする方向に推定することを特徴とする請求項１に記載の道路地図データ作成装置。
【請求項３】
前記道路統合処理部は、抽出された第１の道路ベクトルと第２の道路ベクトルとのノードの中から最も近接しているノードを探索してノードペアとし、該ノードペアにおいて前記追跡線方向に沿った前後のノードについてノード間の距離が一定範囲内であるノードペアを探索し、該探索されたノードペアを統合対象のノードとし、該統合対象のノードペアの間に新しいノードを生成することを特徴とする請求項１又は２に記載の道路地図データ作成装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、道路地図データ作成技術に関し、特に、衛星や航空機等を用いて高度・上空から地表を撮影して得られた観測画像を解析し、特に、地理画像に含まれる道路の抽出を容易にするための道路地図データ作成装置に関する。
【背景技術】
【０００２】
近年、衛星画像や航空写真画像（以下、「地理画像」と称する。）の利用が広まってきており、それらの画像を使った地図の作成や地表の分析等に関する技術が開発されている。その中で、例えばカーナビゲーション（以下、「カーナビ」と略す。）用の道路地図を作成するにためには、画像内の道路領域を抽出することが必要である。
【０００３】
従来、地理画像からの道路領域の抽出処理は、目視作業で行なわれていた。そのため、多大な人的コストがかかるという問題があった。そこで、画像内の道路領域を抽出する処理をコンピュータにより自動化する方法が提案されている。その一例として、下記非特許文献１では、まず、地理画像に対して被覆分類処理を適用し、道路らしい色調特徴を持つ領域を抽出する。また、地理画像に対してエッジ抽出処理及びノイズ除去処理を適用し、道路らしいエッジ特徴を持つ領域を抽出する。そして、それら２つの領域の情報をもとに、道路中心線の追跡を行なう。具体的には、上記２つの領域を重ね合わせた画像に対して、ユーザが、道路１本につき１点の開始点を指定する。すると、プログラムが開始点から自動的に道路中心線を追跡していき道路１本の中心線を抽出することができる。その後、ユーザが画像内のすべての道路に対して開始点を指定していくことにより、画像内の全ての道路中心線を抽出することができる。このような処理により、道路中心線の抽出を半自動的に行なうことができる。その後、目視作業により、道路の幅員情報を付加し、最終的な道路領域データが作成される。
【０００４】
【非特許文献１】熊谷潤,中川雅史,趙卉菁,柴崎亮介,「高解像度衛星画像からの道路抽出」 日本写真測量学会平成１４年度年次学術講演会発表論文集,155-158,2002年。
【非特許文献２】J.Canny,「A Computational Approach to Edge Detection」IEEE Transactions on Pattern Analysis and Machine Intelligence,Vol.PAMI-8,No6,1986年11月．
【発明の開示】
【発明が解決しようとする課題】
【０００５】
しかしながら、非特許文献１の手法においては、道路１本毎に開始点を指定する必要があり、この作業が煩雑であり、コストがかかるという問題がある。また、この手法は、道路中心線を追跡する際に、同一道路の中心線を複数回抽出する処理であり、抽出後のデータが重複してしまうことがある。さらに、この手法で半自動的に抽出された道路情報は、道路の中心線のベクトルデータのみであり、道路の幅員情報を持たない。地図情報において、道路の幅員情報は重要度が高い情報であるため、この手法で抽出した道路情報をもとに、後に幅員情報を付加する作業を行なう必要がある。また、広範囲の画像に対して限られた時間で作業を行なうためには、一般に複数の作業者が必要となるが、この場合、各作業者の主観的判断基準の差によって、幅員の精度にばらつきが生じてしまう可能性が高いという問題もある。
【０００６】
本発明はこのような状況に鑑みてなされたものであり、衛星画像や航空写真画像などの地理画像から、幅員情報を含む道路領域データを容易かつ均質に作成することができる道路地図データ作成技術を提供することを目的とする。
【課題を解決するための手段】
【０００７】
上記解決課題に鑑みて鋭意研究の結果、本発明者は、非特許文献１の手法によって抽出された道路の候補領域と道路の候補エッジを重ね合わせた画像データに対して、開始点を自動的に与える処理と、抽出した道路データの重複部分を統合する処理と、道路の幅員を推定する処理を加えることによって、上記問題が解決されることに想到した。
【０００８】
本発明の一観点によれば、地理画像を解析し、道路地図データを作成する道路地図データ作成装置であって、前記地理画像データから抽出された前記地理画像の領域において複数の開始点をまんべんなく与え、与えられた該複数の開始点のそれぞれから道路を追跡していくことにより得られる追跡線に沿って、各開始点における座標及びそれと結合しているノードの座標とにより定義される道路のベクトルデータを抽出する道路追跡処理部と、前記地理画像内の各道路領域において、前記道路のベクトルデータのうち重複データと推定されるベクトルデータを統合する道路統合処理部と、前記地理画像データ内の各道路領域の幅員を統合された前記ベクトルデータの数に基づいて推定する幅員推定処理部と、を備えることを特徴とする道路地図データ作成装置が提供される。
【０００９】
前記幅員推定処理部は、統合された前記ベクトルデータの数が多いほど幅員を大きくし、統合した道路の前記ベクトルデータの長さが短いほど幅員を大きくする方向に推定することが好ましい。
【００１０】
前記道路統合処理部は、抽出された第１の道路ベクトルと第２の道路ベクトルとのノードの中から最も近接しているノードを探索してノードペアとし、該ノードペアにおいて前記追跡線方向に沿った前後のノードについてノード間の距離が一定範囲内であるノードペアを探索し、該探索されたノードペアを統合対象のノードとし、該統合対象のノードペアの間に新しいノードを生成することが好ましい。
【発明の効果】
【００１１】
以上に説明したように、本発明による地理画像処理装置では、衛星画像や航空写真画像などの地理画像から道路の中心線および幅員を抽出する作業を効率化することができる。また、作業者の技量に頼っていた作業工程を自動化することができるため、作業者の熟練度による差の影響を低減できる。
【発明を実施するための最良の形態】
【００１２】
以下、本発明の一実施の形態による地理画像処理技術について、図面を参照しながら詳細に説明する。図１から図１３までにおいて、同一の符号を付した部分は同一の構成を表し、基本的な構成及び動作は同様であるものとする。尚、本発明は、本実施の形態に記載した範囲に限定されるものではない。当業者であれば、本発明の要旨を逸脱しない範囲内において、実施形態にかかる発明の構成及び機能に様々な変更・改良を加えることが可能である。
【００１３】
＜システム構成＞
図１は、本実施の形態による地理画像処理装置の一構成例を示す機能ブロック図である。図１に示すように、地理画像処理装置は、パーソナルコンピュータ、ワークステーション等から構成される処理装置１０と、メインメモリとして使用されるＲＡＭ（ランダムアクセスメモリ）及び磁気ディスク記憶装置等の補助記憶装置とを含む記憶装置２０と、入出力装置３０と、を備えている。
【００１４】
入出力装置３０は、例えば、キーボード及びマウス等のポインティングデバイスを含む入力装置３１と、液晶ディスプレイ装置等の表示装置３２と、プリンタ３３と、を備えている。入力装置３１は、ユーザによるパラメータの入力やコマンドの起動等、並びに、地理画像処理装置により処理を施した地理画像を用いた土地区画データ生成処理に用いられる。表示装置３２及びプリンタ３３は、地理画像や土地区画データをユーザに提示するために用いられる。尚、表示装置としては、表示装置３２及びプリンタ３３のいずれか一方のみを備える構成としてもよい。
【００１５】
処理装置１０は、地理画像処理プログラム４０を含んでいる。地理画像処理プログラム４０は、プログラムモジュールとして、地理画像に含まれる色調情報をもとにして道路らしい領域（道路候補）を抽出する道路候補領域抽出部１０１と、地理画像のエッジ情報をもとにして道路らしいエッジを抽出する道路候補エッジ抽出部１０２と、道路らしい領域および道路らしいエッジをもとに道路の中心線を抽出する道路追跡処理部１０３と、抽出した道路中心線の重複部分を統合する道路統合処理部１０４と、抽出した道路の幅員情報を推定する幅員推定処理部１０５と、を含む。
【００１６】
記憶装置２０は、地理画像データ記憶領域２１と、道路候補領域データ記憶領域２２と、道路候補エッジデータ記憶領域２３と、道路領域データ記憶領域２４と、を有する。図では、領域という用語を除いて示している。これらの記憶領域のうち、地理画像データ記憶領域２１は、人工衛星又は航空写真画像から得られるデータを記憶する領域であり、地理画像処理装置による処理の実行前に予め記憶されているデータである。一方、道路候補領域データ記憶領域２２、道路候補エッジデータ記憶領域２３、道路領域データ記憶領域２４は、地理画像処理装置において地理画像データ記憶領域２１データに基づいて生成される。
【００１７】
＜地理画像処理の概要＞
図２は、画像処理の概要を示す図である。図２（ａ）は、予め記憶装置２０に記憶されている地理画像データに基づく画像を示している。この地理画像データに対して、被覆分類処理を行ない、さらに、その結果をもとに、道路らしい領域を抽出することにより、図２（ｂ）に示す道路候補領域画像が得られる。この道路候補領域画像は、道路候補領域データ記憶領域２２に格納される。また、地理画像データに対してエッジ抽出処理を行ない、さらに抽出したエッジから道路らしくないノイズを除去することにより、図２（ｃ）に示す道路候補エッジ画像が得られる。
【００１８】
この道路候補エッジ画像は、各画素でエッジの方向を表すデータを保持している。このデータは道路候補エッジデータとして道路候補エッジデータ記憶領域２３に格納される。そして、道路候補領域データと道路候補エッジデータとをもとに、道路の中心線を抽出するため道路追跡処理を行なう。道路追跡処理では、自動的に与えられた開始点から道路を追跡していく。各開始点から道路を追跡すると、図２（ｄ）に示す道路の追跡線を表示した画像が得られる。この道路追跡線（図中の白く細い線）は、各ノードにおける座標および結合しているノード情報を持つベクトルデータである。このベクトルデータは、道路領域データとして道路領域データ記憶領域２４に格納される。但し、右上の領域には、ノイズのような表示が存在する。道路領域データにおける道路抽出線は、同一道路上に複数存在する可能性があり、その場合、推定されるデータが冗長となる。そのため重複したデータを統合する処理（道路統合処理）を行なうことでノイズを除去する。
【００１９】
図２（ｅ）は、図２（ｄ）における道路データを統合した結果を示す図である。統合後のベクトルデータは、道路領域データ記憶領域２４に上書きされて格納される。統合データでは、ノイズはほぼ除去されてきれいな白い線として表されていることがわかる。次いで、幅員推定処理を行なう。道路領域データにおける各道路の道路ベクトル長と、統合したベクトルの本数から、道路の幅員を推定する。図２（ｆ）は、道路データにおける各道路のベクトルデータの幅員を推定した結果を示す図である。明るい領域が道路として抽出された領域であり、幅員が示されている。幅員データは道路領域データ記憶部２４に追加される。
【００２０】
＜地理画像処理の詳細＞
以下に、上記の地理画像処理における各処理の詳細について説明を行う。本実施の形態による地理画像処理装置において、地理画像処理プログラム４０が起動されると、道路候補領域抽出部１０１、道路候補エッジ抽出部１０２、道路追跡部１０３、道路統合処理部１０４、幅員推定処理部１０５、が順に起動される。以下、それぞれの処理部における各処理について詳細に説明する。
【００２１】
（１）道路候補領域抽出処理
図３（ａ）は、地理画像処理プログラム４０の道路候補領域抽出部１０１による処理の流れを示すフローチャート図である。図３（ａ）において、図１に示す〈以下同様〉道路候補領域抽出部１０１は、記憶装置２０の地理画像データ記憶領域２１から地理画像データを読み込み（Ｓ３０１）、被覆分類（上記非特許文献１参照。）を行なう（Ｓ３０２）。被覆分類は、地理画像内に含まれる各画素を、予め決められた様々な物体のクラス（道路、各色の建物、樹木、水域、土壌など）に分類する手法である。近年では、リモートセンシング技術の発達により、１ｍ以下の高分解能衛星画像データが利用可能になった。従来の低分解能の衛星画像と異なり、人工構造物、樹木、道路などの地表物について、コンピュータによる被覆分類でもある程度の認識が可能である。また、航空写真画像も同様の特徴を持っている。
【００２２】
図３（ｂ）は、図２（ａ）の画像に対する被覆分類処理の結果を示す図である。この画像においては、道路、各色の建物、樹木など数クラスの地表物に分類されている。しかしながら、現実的には「道路」と分類された領域が正確に実際の道路を表しているとは限らず、通常は、抽出漏れや、誤抽出を起こしている場合もある。そこで、この被覆分類結果をもとに、抽出漏れした領域を補い、また誤抽出の影響を緩和するため、道路候補領域抽出処理を行なう（Ｓ３０３）。道路候補領域を抽出する際には、被覆分類結果の中から、道路クラス以外に道路と誤認しやすいクラス（例えば、白色系や黒色系の色調特徴を持つ建造物のクラス）を、誤認しやすさの重みを考慮して抽出する。また、道路と誤認する可能性が低いクラス（例えば、樹木や土壌などのクラス）はこの処理の対象から除外する。そのような処理を加えることにより、被覆分類結果画像から、道路らしい領域（道路候補領域）を抽出することができる。図３（ｃ）は、道路候補領域を示す図であり、明るく示されている領域が道路候補領域である。図３（ｃ）に示す道路候補領域を特定するデータは、道路候補領域データとして道路候補領域データ記憶領域２２に出力され格納される（Ｓ３０４）。
【００２３】
（２）道路候補エッジ抽出処理
図４（ａ）は、図１に示す（以下同様）地理画像処理プログラム４０の道路候補エッジ抽出部１０２による処理の流れを示すフローチャート図である。図４（ａ）において道路候補エッジ抽出部１０２は、記憶装置２０の地理画像データ記憶領域２１から地図画像データを読み込み（Ｓ４０１）、エッジ抽出を行なう（Ｓ４０２）。尚、画像におけるエッジ抽出処理には様々な手法があるが、特にＣａｎｎｙ法（上記非特許文献２参照。）が適している。図４（ｂ）は、図２（ａ）の地理画像に対してＣａｎｎｙ法によりエッジ抽出を行なった結果として得られた画像である。この画像には、道路のエッジが含まれているが、同時にノイズも多数含まれている可能性がある。そこで、次に、この画像から道路以外のノイズを除去する（Ｓ４０３）。ノイズを除去する際には、短いエッジや、道路らしくない急激に方向が変化しているエッジを除去する（上記非特許文献１参照。）。図４（ｃ）は、図４（ｂ）の画像に対してノイズ除去処理を行なった結果として得られた画像である。ノイズ除去後の各画素上のエッジデータは道路候補エッジデータとして道路候補エッジデータ記憶領域２３に格納される（Ｓ４０４）。
【００２４】
（３）道路追跡処理
図５は、地理画像処理プログラム４０の道路追跡処理部１０３における処理の流れを示すフローチャート図である。図５において、図１に示す（以下同様）道路追跡処理部１０３は、記憶装置２０の道路候補領域データ記憶領域２１と道路候補エッジデータ記憶領域２２からそれぞれのデータ（図３（ｃ）の道路候補領域データと図４（ｃ）の道路候補エッジデータ）を読み込み、その２つのデータを画像上で重ねあわせる処理を行う（Ｓ５０１）。図６（ａ）は、これら２つのデータを重ね合わせた様子を示す図である。道路追跡処理部１０３（図１）では、この図において道路追跡を行い道路の中心線のデータを抽出する。道路追跡を行なう際には、まず追跡時の開始点を設定する（Ｓ５０２）。開始点は、画像全体に対してまんべんなく多数点を設定するのが好ましい。例えば、図６（ｂ）に示すように、画像における互いに直交する座標軸であるｘ軸とｙ軸とに対して平行で、かつ、等間隔に開始点を多数設定する方法が好ましい。このような方法によれば、ユーザが開始点を設定する必要がないという利点がある。
【００２５】
次に、各開始点から道路の追跡を行なう。道路の追跡処理としては、ラインテンプレートマッチング法により行なうのが好ましい（上記非特許文献１参照。）。まず、未処理の開始点を１つ選択し（Ｓ５０３）、スタート方向を決定する（Ｓ５０４）。スタート方向を決定する際には、図６（ｃ）に示すように、長方形のテンプレートＴ１を画像に重ね、道路候補領域データと道路候補エッジデータとを参照する。そしてテンプレートＴ１の角度を変えて各角度での、両データから算出される道路らしさを測定し、より道路らしい方向をスタート方向とする。図６（ｃ）の場合で説明すると、左上の方向がより道路らしい（道路候補領域データおよび道路候補エッジデータの明るい画素が多いことにより判定することができる）ため、左上の方向をスタート方向として決定する（スタート方向は、図中の矢印の方向である）。スタート方向を決定した後に、延長方向を探索する（Ｓ５０５）。延長方向を決定する際にも、ラインテンプレートの角度を変えながら、より道路らしい方向（を延長方向とする処理を継続する（図６（ｄ）の明るい矩形と破線で示される角度の異なる矩形）。
【００２６】
但し、スタート方向の決定時と異なり、延長方向を決定する際には、変化させる角度範囲を狭くすることが好ましい。それは、一般に道路は交差点や山間部などの特殊な場合を除けば、それほど急激に曲がることは少なく、比較的緩やかな曲線を描くことが多いためである。そして、延長方向を決定すると（道路候補エッジデータの明るい画素が多い方向）、ラインテンプレートの長さの分だけ道路候補エッジデータを延長し、延長先の点を、道路ベクトルを構成するノード（起点）とする。
【００２７】
ここで、ラインテンプレートの角度を変えて延長方向を決定する際に、どの方向でも道路らしさが低いと判定される場合には、延長処理を中止し、道路の追跡を終了しフローチャートのＳ５０７の反対側の処理を行なう（Ｓ５０６）。延長方向が決定した場合は、延長先で再び延長方向を決定する（Ｓ５０５）。この処理を、延長方向が決定できなくなるまで、すなわち道路が途切れるまで行なう。このような道路の追跡を、Ｓ５０４で決定したスタート方向の反対側でも行なう。Ｓ５０７は反対側を追跡したかどうかの判定処理である。Ｓ５０３で設定した開始点に対して両側で道路の追跡を完了した後に、追跡時に設定されたノードを結合させることで生成された線分の連結データを１本の道路領域データとする（図６（ｅ）のＬ１）。以上のＳ５０３からＳ５０７までの処理を、Ｓ５０２で設定した各開始点に関して行なう。Ｓ５０８はその判定処理である。図６（ｆ）は道路追跡処理後の結果を示した図である。図６（ｆ）に示すように、多数の開始点に基づいて、道路の追跡を行うため、多数の細い線が抽出されていることがわかる。各道路ベクトルのデータは、図１０に示すような形式で記憶装置２０〈図１〉に格納される。但し、“統合数”および“幅員”の項目は、この段階では求められていないため、これらの項目のデータは格納されていない。“道路ベクトル長”の項目は、後述する幅員推定処理で算出するため、この段階ではデータが格納されていないが、この段階でも算出可能であるためデータを格納しておいてもよい。
【００２８】
（４）道路統合処理
図７は、地理画像処理プログラム４０の図１に示す（以下同様）道路追跡処理部１０３における処理の流れを示すフローチャート図である。図７において、道路統合処理部１０４は、記憶装置２０の道路領域データ記録領域２４から道路領域データを読み込み（Ｓ７０１）、道路領域データの中から道路ベクトルを１本だけ読み込む（Ｓ７０２）。各道路ベクトルに対して近傍の道路ベクトルと統合処理を行なう。まず、注目している道路ベクトル（Ａ）の近傍に、距離と方向とが近い道路ベクトル（Ｂ）があるか否かをチェックする（Ｓ７０３）。距離と方向が近い道路ベクトルが近くに見つからなかった場合はＳ７０２に戻り、別の道路ベクトルを読み込む。距離と方向が近い道路ベクトルが近くに見つかった場合はＳ７０５に進み、２つの道路ベクトル（Ａ）と道路ベクトル（Ｂ）とが統合すべきデータと推定されるためこれらを統合する。図８（ｂ）および（ｃ）は、図８（ａ）に示すような道路に対して抽出された２つの道路ベクトル（Ａ）と（Ｂ）とが統合される様子を示すステップ図である。道路ベクトル（Ａ）と（Ｂ）とを統合する際には、まず２つの道路ベクトルのノードの中から最も近接しているノードを探索し、このノードペアＰ１、Ｐ２、…を基準とする。この基準のノードペアＰ１、Ｐ２、…において前後のノードを調べ、それぞれのノードペアＰ１、Ｐ２、…におけるノード間の距離が一定範囲内であるノードペアを探索する。
【００２９】
ここで探索されたノードペアを統合対象のノードとする。この際、ノード間の距離が一定以上となり近接するノードがなくなったと判定されると探索処理を終了する。次に、統合対象のノードペアの中間位置に統合後のノードを設定する。図８（ｂ）では、統合対象のノードペアが６ペアあるため、各ノードペアの中間に６個の新しいノードＰＡ１、…が設定される。そしてこれらの新しいノードＰＡ１、…を近い順に結合し、結合させて生成される連結した線分を新しい道路ベクトルとする。この際、統合比で、新たなノード座標を求めても良い、図８（ｃ）では、統合比が２：５であるため、下に引き出した図に示すように、ノード座標が統合数の大きい方に近い位置になる。
【００３０】
また、前の処理で近接したノードがなく切断された部分では、切断された部分のノードと新しい道路ベクトルのノードとを結合させる。この際、枝分かれした部分のノードは３つの道路ベクトルで共有される。また、統合の際には、これまでに統合した道路ベクトルの数（統合数）を記憶する。この統合数の情報は、後述する幅員推定処理で幅員を推定する際に使用することができる。道路ベクトル（Ａ）と（Ｂ）との統合後に、道路領域データを更新する（Ｓ７０６）。この際、統合によって道路ＩＤに増減が生じることがあるため、道路ＩＤを付与し直す。Ｓ７０２からＳ７０６までの処理をすべての道路ベクトルに対して行なう（Ｓ７０７）。最後に、Ｓ７０２からＳ７０７までの処理の結果、道路領域データが更新されたか否かをチェックする（Ｓ７０８）。Ｓ７０２からＳ７０７までの処理の中で少なくとも１つの道路ベクトルが更新されていれば、Ｓ７０２に戻り再帰的に処理を繰り返す。更新がなければ処理を終了する。
【００３１】
図９は道路統合処理を示す概略図である。二又に分かれた道路上に４本の道路ベクトル（ＩＤ：０００１〜０００４）が存在しており、道路統合処理を適用した結果、統合数４の道路ベクトル（ＩＤ：０００１）１本と、統合数２の道路ベクトル（ＩＤ：０００２、０００３）が２本生成されている様子が示されている。
【００３２】
（５）幅員推定処理
図１１は、地理画像処理プログラム４０の幅員推定処理部１０５（図１参照）における処理の流れを示すフローチャート図である。図１１において、幅員推定処理部１０５は、記憶装置２０の道路領域データ記憶領域２４から道路領域データを読み込む（Ｓ１１０１）。次いで、各道路ベクトルに対してＳ１１０２以降の処理を行なう。まず、各道路ベクトルに対し、道路ベクトル長（各ノード間のショートベクトルの長さの合計）を計算する（Ｓ１１０２）。次に、各道路ベクトルの道路ベクトル長と統合数とから幅員を推定する。本実施の形態では、道路追跡処理部１０３において道路ベクトルを抽出する際に、画像全体に対して均等に開始点を設定しているため、道路の面積が大きいほど多くの道路ベクトルが抽出される。そのため、道路統合処理部１０４において重複した道路ベクトルを統合する際には、道路の面積が大きいほど統合される道路ベクトルの本数も多くなる。面積は道路の長さと幅員とから概算することができるため、道路ベクトル長と道路ベクトルの統合数とから幅員を推定することができる。幅員の概算の例は、道路ベクトルの統合数／道路ベクトル（比）の長さである。幅員の推定は、基準となる幅員に上記の比を乗算することにより行う。より具体的には、下記の式（１）により求めることができる。
【００３３】
（道路ベクトルの統合数／道路ベクトルの長さ）×ａ（比例定数）    （１）
比例定数は、経験的に求めても良いし、同じ国、地域であれば、同じと推定して予め決めておいても良い。後に、実測値と比較して修正しても良い。
【００３４】
図１２は、本実施の形態によって道路のベクトル及び幅員を抽出される様子を示す図である。図１２（ａ）は、Ｙ字型の道路の画像と、その上に設定された開始点を表している。図１２（ｂ）は、図１２（ａ）の画像に対して道路追跡処理を適用した結果の画像である。道路上の複数の開始点から道路を追跡したため同一道路上に複数の道路ベクトルが存在する。図１２（ｃ）は、図１２（ｂ）の結果に対して道路統合処理を適用した結果画像である。道路統合処理により、同一道路内に複数存在した道路ベクトルが統合されている様子が示されている。図１２（ｄ）は、図１２（ｃ）の結果に対して幅員推定処理を適用した結果画像である。道路ベクトル長データと統合数データとをもとに、幅員を推定した結果を表している。
【００３５】
図１３は、図１０に対応する図であり、道路領域データ２４の最終的な処理結果を例示した図である。各道路ベクトルのＩＤごとに各データが格納される。道路ＩＤ００００１では、統合数が６０、道路ベクトル長が１５０であり、道路ＩＤ００００２では、統合数が４０、道路ベクトル長が１００である。Ｓ１１０３では、例えば、幅員は、統合数／ベクトル長で概算され、より具体的には上記の式（１）により求めると、以下のようになる。
【００３６】
すなわち、上記の式（１）において、ａ＝１２５（この値は任意）とすると、ＩＤ００００１の道路の幅員は（６０／１５０）×１２５＝５０となり、ＩＤ００００２の道路の幅員は（４０／１００）×１２５＝５０となり、共に５０となる。
【００３７】
＜まとめ＞
以上のように本発明の実施形態は、地表を撮影して得られた地理画像データを解析し、地表上の道路領域データを容易に作成する道路地図データ作成システムに関するものである。そしてそのシステムは、地理画像データから色調情報をもとにして道路の候補となる領域を抽出する道路候補領域抽出部（モジュール）と、地理画像データからエッジ情報をもとにして道路の候補となる領域を抽出する道路候補エッジ抽出部（モジュール）と、道路領域を追跡して道路のベクトル情報を抽出する道路追跡処理部（モジュール）と、重複した道路のベクトルを統合する道路統合処理部（モジュール）と、道路の幅員情報を推定する幅員推定処理部（モジュール）と、を備える。これにより幅員情報を含む重複の無い道路領域の抽出を自動化できる。
【００３８】
尚、本実施の形態の機能を実現するソフトウェアのプログラムコードによっても本発明は実現できる。この場合、プログラムコードを記録した記憶媒体をシステム或は装置に提供し、そのシステム或は装置のコンピュータ（又はＣＰＵやＭＰＵ）が記憶媒体に格納されたプログラムコードを読み出す。この場合、記憶媒体から読み出されたプログラムコード自体が前述した実施形態の機能を実現することになり、そのプログラムコード自体、及びそれを記憶した記憶媒体は本発明を構成することになる。このようなプログラムコードを供給するための記憶媒体としては、例えば、プロッピィ（登録商標）ディスク、ＣＤ−ＲＯＭ、ＤＶＤ−ＲＯＭ、ハードディスク、光ディスク、光磁気ディスク、ＣＤ−Ｒ、磁気テープ、不揮発性のメモリカード、ＲＯＭなどが用いられる。
【００３９】
また、プログラムコードの指示に基づき、コンピュータ上で稼動しているＯＳ（オペレーティングシステム）などが実際の処理の一部又は全部を行い、その処理によって前述した実施の形態の機能が実現されるようにしてもよい。さらに、記憶媒体から読み出されたプログラムコードが、コンピュータ上のメモリに書き込まれた後、そのプログラムコードの指示に基づき、コンピュータのＣＰＵなどが実際の処理の一部又は全部を行い、その処理によって前述した実施の形態の機能が実現されるようにしてもよい。
【００４０】
また、実施の形態の機能を実現するソフトウェアのプログラムコードがディスクやメモリ等の記憶手段又はＣＤ−ＲＷ、ＣＤ−Ｒ等の記憶媒体に格納され、そのシステム又は装置のコンピュータ（又はＣＰＵやＭＰＵ）が当該記憶手段や当該記憶媒体に格納されたプログラムコードを読み出して実行することによっても、達成されるようにしてもよい。
【産業上の利用可能性】
【００４１】
本発明は、道路の地図データ作成装置に利用可能である。
【図面の簡単な説明】
【００４２】
【図１】本発明の一実施の形態による地理画像処理装置の一構成例を示す機能ブロック図である。
【図２】本実施の形態による地理画像処理装置ムにより行われる画像処理の概要を例示する図である。
【図３】道路候補領域抽出部の処理の流れを示すフローチャート図である。
【図４】道路候補エッジ抽出部の処理の流れを示すフローチャート図である。
【図５】道路追跡処理部の処理の流れを示すフローチャート図である。
【図６】道路追跡処理部の処理の概要を例示する図である。
【図７】道路統合処理部の処理の流れを示すフローチャート図である。
【図８】道路統合処理部の処理の概要を例示する図である。
【図９】道路統合処理部の処理の結果を例示する図である。
【図１０】道路領域データの内容を例示する図である。
【図１１】幅員推定処理部の処理の流れを示すフローチャート図である。
【図１２】幅員推定処理部の処理の概要を例示する図である。
【図１３】最終的な処理結果を例示する図である。
【符号の説明】
【００４３】
１０  処理装置
２０  記憶装置
２１  地理画像データ記憶領域
２２  道路候補領域データ記憶領域
２３  道路候補エッジデータ記憶領域
２４  道路領域データ記憶領域
３０  入出力装置
３１  入力装置
３２  表示装置
３３  プリンタ
４０  地理画像処理プログラム
１０１  道路候補領域抽出部
１０２  道路候補エッジ抽出部
１０３  道路追跡処理部
１０４  道路統合処理部
１０５  幅員推定処理部
【図１】
0004850161.tif 000002
【図２】
0004850161.tif 000003
【図３】
0004850161.tif 000004
【図４】
0004850161.tif 000005
【図５】
0004850161.tif 000006
【図６】
0004850161.tif 000007
【図７】
0004850161.tif 000008
【図８】
0004850161.tif 000009
【図９】
0004850161.tif 000010
【図１０】
0004850161.tif 000011
【図１１】
0004850161.tif 000012
【図１２】
0004850161.tif 000013
【図１３】
0004850161.tif 000014