(57)【要約】      （修正有）【課題】農地区画データの作成処理を簡単にするとともに精度を良くする。【解決手段】区画データに対応する位置の詳細エッジを選んだ特徴ベクトルの垂直及び平行な方向に投影し、その方向で見たエッジの量を調べ、集計値によりピークを見つける。ピークが見つかった場合は、ピークの位置には畦道と推定されるエッジの集合があると考えられる。ピークの位置Ｐ１、Ｐ２で区画データを分割し、それぞれを一つの区画として詳細区画データに登録する。複数の投影方向でピークが見つかった場合には、得られたピークの強さを比較し、最大の一方向のみで区画を分割する。１つの投影方向で複数のピークが見つかった場合は、ピークで分割する。【選択図】図８
2008234163.tif 000002
【特許請求の範囲】
【請求項１】
地理画像データを解析し、土地区画データを作成する土地区画データ作成システムであって、
地理画像データからエッジデータを抽出するエッジ抽出処理部と、
該エッジ抽出処理部において第１のエッジ強調度でエッジ抽出した概要エッジデータに基づいて多角形を形成するエッジを区画として抽出する概要区画抽出処理部と、
該概要区画抽出処理部で得られた各区画から概要区画特徴データを抽出する概要区画解析処理部と、
前記エッジ抽出処理部において前記第１のエッジ強調度よりも高い第２のエッジ強調度でエッジ抽出した詳細エッジデータと、該詳細エッジデータと同じ区画の概要区画データと、該概要区画データに基づいて前記概要区画解析処理部により得られた概要区画特徴データと、に基づいて、詳細区画データを作成する詳細区画抽出処理部と
を有することを特徴とする土地区画データ作成システム。
【請求項２】
前記詳細区画抽出処理部は、
前記詳細エッジデータと対応する概要区画データにおいて概要区画特徴データとして特徴ベクトルが存在する場合に、該特徴ベクトルに水平な方向及び垂直な方向において区画内の前記詳細エッジデータを集計してエッジデータの集計値のピークを検出し、ピークが検出された場合に、該ピークの位置において前記概要区画データを分割して前記詳細区画データとして登録することを特徴とする請求項１に記載の土地データ作成システム。
【請求項３】
前記概要区画解析処理部において抽出する特徴データとして、区画を構成する主たるベクトルの数と方向とを用いることを特徴とする請求項１又は２に記載の区画データ作成システム。
【請求項４】
前記概要区画解析処理部は、
概要エッジデータに基づいて多角形を形成するエッジを区画として抽出する処理であって、隣接する頂点を結ぶベクトルのうち方向が近いベクトルを統合し、ある値以上の長さを有するベクトルを集計しヒストグラムを作成し、該ヒストクラムのピークが２又は３である場合に、該ピークを前記概要区画特徴データとすることを特徴とする請求項１から３までのいずれか１項に記載の区画データ作成システム。
【請求項５】
前記エッジ抽出処理部は、入力によりエッジ強調度を調整するエッジ強調度調整部を有し、該エッジ抽出処理部におけるエッジ抽出結果を表示部に表示可能に構成されていることを特徴とする請求項１から４までのいずれか１項に記載の区画データ作成システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、衛星や航空機等を用いて高度上空から地表を撮影した観測画像を解析する画像認識システムに関し、特に、画像から道路、河川、畦道などの土地の境界線を抽出し、農地区画等を解析してその地区の区画データを作成するシステムに関する。
【背景技術】
【０００２】
近年、衛星画像や航空写真の利用が広まってきており、特に地図の作成や更新における広大な領域の探索作業を効率的に行うために利用されている。
【０００３】
衛星画像や航空写真の利点は、従来のような現地で測量する際のコストを省けること、広い範囲のデータを一度に得ることができることである。さらに、撮影頻度を上げることにより、常に新しいデータを用意することも可能となる。
【０００４】
衛星画像を用いて地図の作成や更新をする作業は、従来は、画像を目視して地物や境界線を確認し、人手によりそれらを入力するという手法を用いていた。目視による地図作成作業は、作業者への負担が大きい上、作業者の熟練度により作業の進み具合や結果として得られるデータの精度が大きく異なるという問題点があった。現在では、安定した結果やコストの削減のため、コンピュータを利用して画像を解析する様々な手法が研究され、提案されている。
【０００５】
画像から道路や河川などの境界線を検出し、その境界線のつながりを解析して土地区画を自動的に抽出する手法の一例として、下記特許文献１では、まず、地表を撮影した画像からエッジ(輪郭線)を抽出し、そのエッジ画像上の任意の１点から放射状に伸ばした直線とエッジとの交点を求め、複数の交点を連結させてできる多角形の区画領域を抽出する処理技術が提案されている。
【０００６】
【特許文献１】特開２００３−２５６８０７号公報
【非特許文献１】J.Canny,「A Computational Approach to Edge Detection」IEEE Transactions on Pattern Analysis and Machine Intelligence,Vol.PAMI-8,No6,1986年11月
【非特許文献２】尾上守,「画像処理ハンドブック」,昭晃堂出版,1987年
【発明の開示】
【発明が解決しようとする課題】
【０００７】
一般的に、画像の解析を行なう際は、画像から形や模様、境界線等の特徴を抽出するフェイズと、そこで取れた特徴を解析して必要な情報に加工するフェイズと、に分かれる。上記特許文献１に記載の手法も、画像からエッジを抽出するフェイズと、そのつながりを補完したり、閉曲線を検出し、区画の情報に加工するフェイズとから成っている。
【０００８】
このような手法による画像の解析結果は、前半のフェイズで抽出された特徴の種類と精度とに大きく影響される。従って、画像から特徴を適切に抽出することが重要である。
【０００９】
衛星画像や航空写真から地図を作る際は、地表に見られる各種の境界線が特に重要な特徴である。これらの境界線は多くの場合、建物や道路、水域等のエッジとして検出することができる。画像からエッジを抽出する方法は、これまでにも多くの研究がなされており、実用的な種々の提案もなされている(上記非特許文献１，２)。
【００１０】
画像からエッジを抽出する方法を用いた場合の理想的な結果は、後半の解析過程において必要な情報を全て備え、その他の不必要な情報を含まないものである。しかしながら、実際に利用される方法では理想的な結果が得られることは少ない。より多くのエッジを抽出した場合には、不必要な情報（以下「ノイズ」と称する。）の量も多くなる。一方、ノイズの量を減らそうとすると、実際には必要なエッジに関する情報量も減ってしまうため、エッジ抽出の手法やその中で用いるパラメータを調整する必要があるという問題がある。
【００１１】
抽出したいエッジの特徴が例えばエッジの周方向に向けてほぼ均等である場合には、従来の手法でもパラメータの調整を行なうことにより、良い結果を得ることができるが、抽出したいエッジの特徴が周方向で異なる場合には調整が困難になる。
【００１２】
例えば、農地区画データを作成する場合は、農地の境界となるエッジ、つまり、例えば、道路や河川、畦道などを抽出する必要がある。しかしながら、それらのエッジを抽出する場合に、道路や河川のエッジの特徴と畦道のエッジの特徴とには差があるため、道路のエッジだけを取るとエッジの量が不十分となり、反面、畦道のエッジまで取るとノイズの量が多くなりすぎるという問題があった。
【００１３】
このような場合に、従来の方法では自動的な解析が難しく、特徴を抽出した結果を改めて目視で確認し、修正する必要があった。従って、人手で入力する手法と比べると安定した結果が得られず、また、安定した結果を得ようとして人手をかけるとコスト高になるという問題があった。
【００１４】
本発明は、農地区画データの作成処理を簡単にするとともに精度を良くすることを目的とする。
【課題を解決するための手段】
【００１５】
本発明は、前記の課題に対応するため、衛星画像や航空写真画像などの地理画像から特徴の異なるエッジを複数の段階に分けて抽出することを特徴とする。
【００１６】
画像から特徴を抽出するフェイズを複数に分割し、前段階の特徴抽出で得られたデータを分析した情報を後段階の特徴抽出で利用する段階型の処理を行なうことにより特徴に差があるエッジを高い精度で抽出することを可能にした。農地区画データを作る際には、まず特徴がはっきりしていてノイズを含まずに取ることができる道路や河川のエッジを抽出する。その後に道路や河川のエッジから大まかな区画を抽出し、畦道は道路や河川に平行または垂直なものが多いという特徴を利用し、畦道のエッジとノイズを含むエッジ画像からノイズを分離して畦道のエッジのみを抽出する。そして両方の手順で得られたエッジを総合して農地区画を抽出する。
【００１７】
具体的には、衛星や航空機から農地地区の地表を撮影した画像を入力とし、道路や河川、畦道のエッジを、複数の段階に分けて抽出するエッジ抽出処理と、エッジ抽出処理で得られた複数のエッジ画像のうちエッジの量が少ないものから多角形を形成するエッジを区画として抽出する概要区画抽出処理と、概要区画抽出処理で得られた各区画から特徴を抽出する概要区画解析処理と、概要区画解析処理で得られた特徴を利用し、エッジ抽出処理で得られたエッジが多いエッジ画像からより詳細な区画を抽出する詳細区画抽出処理を経て、撮影された地区の農地の境界を解析し、正確な農地区画データを作成する。  
【００１８】
さらなる本発明の特徴は、以下の発明を実施する形態及び添付図面によって詳細を示す。
【発明の効果】
【００１９】
以上、説明したように、本発明の土地区画データ作成システムによれば、農地区画等の土地区画の抽出及び区画データの作成を効率化的に行うことができる。作業工程の自動化が可能となるため、作業コストの削減に加え、確認、修正作業者の熟練度に左右されない均質な区画データを得ることができるという利点がある。
【発明を実施するための最良の形態】
【００２０】
以下、図面を参照しながら、本発明の一実施の形態による土地区画データ作成システムについて詳細に説明する。図１から図９までは、本実施の形態による土地区画データ作成システムを例示する図である。これらの図において、同一の符号を付した部分は同一物を表し、基本的な構成及び動作は同様であるものとする。
【００２１】
本実施の形態による土地区画データ作成システムは、高空から農地地区の地表を撮影した画像を入力とし、各画像解析処理過程を経て、撮影された地区の農地の管理に必要な詳細区画データを作成し、出力するものである。
【００２２】
＜システム構成＞
図１は、本実施の形態による土地区画データ作成システムの概略構成例を示す機能ブロック図である。本実施の形態による土地区画データ作成システムは、入出力装置（１００）、パソコン、ワークステーション等からなる処理装置（２００）、及びメインメモリとして使用されるＲＡＭ(ランダムアクセスメモリ)及び磁気ディスク記憶装置等の補助記憶装置とを含む記憶装置（３００）を備える。
【００２３】
入出力装置（１００）は、キーボード及びマウス等のポインティングデバイスを含む入力装置（１０１）と、ＣＲＴディスプレイ装置等の表示装置（１０２）と、プリンタ等の出力装置（１０３）と、を備える。入力装置（１０１）は、ユーザによるコマンドの起動やパラメータの入力等に用いられる。表示装置（１０２）及びプリンタ等の出力装置（１０３）は、本実施の形態による土地区画データ作成システムで処理される画像や区画データをユーザに提示するために用いられる。
【００２４】
処理装置（２００）は、区画抽出プログラム（２０１）を含む。区画抽出プログラム（２０１）は、プログラムのモジュールとして、エッジ抽出処理部（２０２）、概要区画抽出処理部（２０３）、概要区画解析処理部（２０４）、及び、詳細区画抽出処理部（２０５）を含む。区画抽出プログラムの詳細な処理の内容については図３以降を参照して説明する。エッジ抽出処理部（２０２）には、さらに、エッジの強調度を調整するためのエッジ強調度調整部２０２ａが設けられている。
【００２５】
記憶装置（３００）は、画像データ（３０１）、概要エッジデータ（３０２）、概要区画データ（３０３）、概要区画特徴データ（３０４）、詳細エッジデータ（３０５）、詳細区画データ（３０６）、を格納する。これらのデータのうち、画像データ（３０１）は、衛星等で撮影された画像データであり、システム外部から入力される。それ以外の各種データ（３０２）〜（３０６）は、画像解析処理によって作成されるデータであり、詳細区画データ（３０６）を作成することがシステムの最終目的である。それぞれのデータの詳細例については図２に示される。
【００２６】
＜取り扱うデータの種別と内容＞
図２Ａは、本実施の形態による土地区画データ作成システムによる処理中に記憶装置（３００）に格納される各データ例を示す図である。図２Ａ(a)は、衛星から農地を撮影した画像の一例を示す図であり、画像データ（３０１）に基づく画像に相当する。このデータは、画像解析処理の実行前に予め記憶装置３００内に格納されている必要がある。
【００２７】
図２Ａ(b)は、画像データ（３０１）に基づいて作成された第１のエッジ画像(画像の輪郭を抽出した画像)の一例であり、概要エッジデータ（３０２）に相当する。この画像には、道路や河川等のはっきりしたエッジのみが含まれており、はっきりしないエッジやノイズ等をほとんど含まない。
【００２８】
図２Ａ(c)は、画像データ（３０１）を元に作られた第２のエッジ画像の一例であり、詳細エッジデータ（３０５）に相当する画像である。この画像は、上記のようなはっきりしたエッジに加えて、畦道のようなあいまいなエッジも抽出されているが、ノイズ成分も相対的に多く含まれる。
【００２９】
図２Ｂ(d)は、土地の区画を示すデータの一例であり、概要区画データ（３０３）、詳細区画データ（３０６）に相当するデータの例を示す図である。区画ＩＤによってユニークに特定される１区画が１エントリとして扱われ、その内容は区画ＩＤと区画に含まれる頂点の数と、頂点１から頂点Ｎまでの各頂点の座標（ｘ，ｙ）の配列と、により表現されている。
【００３０】
図２Ｂ(e)は、図２Ｂ(d)の区画データを画像として表現した一例を示す図である。その中の一部を示す右側の拡大部分において、区画は多数の頂点を結んだ多角形で表されている。この区画に割り振られたＩＤ,頂点数と座標の配列とが図２Ｂ(d)のテーブルにおいて示されるデータの１エントリに対応する。すなわち、図２Ｂ（ｄ）は、画像中に含まれる区画毎にその位置・形状等を特定することができるデータとなっている。
【００３１】
図２Ｂ(f)は、区画データの特徴を示すデータ構成を示す図であり、概要区画特徴データ(３０４)に相当するデータの例を示す図である。このシステムでは、特徴として区画を構成するベクトルを集計した場合にピークとなったベクトルの数とその方向(単位ベクトル)とを使用している。対応する区画ＩＤは区画データと対応付けられている。例えば区画ＩＤ＝１の区画は、特徴ベクトル数が２であり、特徴ベクトル１（１，０）と特徴ベクトル２（０，１）とによって構成される区画であることを示している。一方、区画ＩＤ＝２の区画は、特徴ベクトル数が０であり、集計した場合にピークが無いため、図形を構成するベクトルが本システムで扱う特徴に一致しないことを示す。これらの区画は(後述する)詳細区画抽出処理の処理対象とならない。これらのデータと区画抽出プログラムとの関係の詳細を図３に示す。
【００３２】
＜地理画像処理の概要＞
図３は、本実施の形態による土地区画データ作成システムにより行なわれる画像解析の概要を示す図である。本実施の形態による土地区画データ作成システムは、まず、予め用意された画像データ（３０１）を入力とし、それに対してエッジ抽出処理（２０２）を行う。エッジ抽出処理の出力として、エッジ強調度調整部２０２ａによりエッジ強調度を相対的に低くした概要エッジデータ（３０２）と、エッジ強調度を相対的に高くした詳細エッジデータ（３０５）とを得る。
【００３３】
概要エッジデータと詳細エッジデータとの違いは、図２Ａ(b)と図２Ａ(c)とで対比させて示したように、エッジ強調度を低くした概要エッジデータは、道路や河川等のはっきりしたエッジのみが含まれ、ノイズをほとんど含まない。詳細エッジデータはエッジ強調度を高めて、はっきりしたエッジに加えて、畦道のようなあいまいなエッジも抽出されるようにしているが、反面、ノイズも多く含んでいる。尚、入力装置（１０１）から、任意のエッジ強調度を入力してエッジ抽出処理を行わせることができ、そのエッジ抽出処理結果は、表示装置（１０２）において適宜参照することができる。
【００３４】
次に、概要エッジデータ（３０２）を入力として概要区画抽出処理部による処理（２０３）を行い、概要区画データ（３０３）を得る。概要区画データ（３０３）は、はっきりとしたエッジのみで構成される区画データであり、実際の農地区画としては正確ではなく、本来畦道などで区切られるべき複数の区画が一つの区画としてまとめられたデータとなっているケースが多い。区画抽出処理の詳細な内容は図４に示す（後述する）。
【００３５】
次に、概要区画データ（３０３）を入力として概要区画解析処理部２０４による処理を行い、概要区画特徴データ（３０４）を得ることができる。概要区画特徴データは、概要区画データに含まれる各区画から、特徴として区画を構成するベクトルの主な方向とその数を抽出したデータである。
【００３６】
最後に、概要区画データ（３０３）と概要区画特徴データ（３０４）、詳細エッジデータ（３０５）、を入力として詳細区画抽出処理部（２０５）による処理を行い、詳細区画データ（３０６）を得る。詳細区画データ（３０６）は、本実施の形態によるシステムの目的とするデータであり、畦道等で正確に区切られた農地区画のデータが含まれ、農地の管理や解析に利用することができるデータである。
【００３７】
＜地理画像処理の詳細＞
以下、本実施の形態による地理解析システムの処理装置２００におけるエッジ抽出処理、概要区画抽出処理、概要区画解析処理、及び詳細区画抽出処理のそれぞれについて詳細を説明する。
【００３８】
(1)エッジ抽出処理（２０２）
エッジを抽出する処理としては種々の周知技術が適用可能である。例えば、Canny法（上記非特許文献１参照）や、Sobelフィルタ（上記非特許文献２参照）などの各種エッジフィルタを地理画像データに適用することで、畦道などのエッジを抽出することができる。これらの処理は、画像中の濃度の変化が大きい部分をエッジとして抽出する。どこまでをエッジとみなすかは処理のパラメータに依存するため、そのパラメータを調整することにより抽出するエッジを限定することが出来る。ここでは、エッジの強調度を調整するエッジ強調度調整部により異なる２種類のエッジ強調度を用い、これら２種類のエッジ強調度に基づくエッジ抽出処理の結果、すなわちエッジ強調度を相対的に低くし道路や河川等のはっきりしたエッジのみが含まれノイズをほとんど含まない概要エッジデータと、エッジ強調度を高くすることによってはっきりしたエッジに加えて畦道のようなあいまいなエッジも抽出されているがノイズも多く含んでいる詳細エッジデータと、を出力し、その後の処理に利用する。
【００３９】
(2)概要区画抽出処理（２０３）
概要エッジ画像に対して、大まかな区画を抽出する概要区画抽出処理には、例えば直線交点法（特許文献１参照）を用いる。
【００４０】
図４は、この概要区画抽出処理の流れを画像で示した図である。この処理は、表示４０１において、エッジ抽出処理（２０２）を経て与えられたエッジ画像と、シード点と呼ばれる画像内の一点の座標から、シード点を含む区画を抽出する。シード点とは、区画を抽出する際の出発点となる座標のことを指す。シード点は、格子状に点在させたり、ランダムに配置したり、ユーザが予め入力しておくなどの方法で決定することができる。
【００４１】
すなわち、図４(a)に示すように、エッジ画像とシード点４０３とが与えられている場合に、まず、シード点４０３から例えば等角度間隔に全方位に直線４０５を伸ばす（図４(b)）。そして、それぞれの直線４０５とエッジ画像におけるエッジ画素との交点４０７を求め、得られた交点４０７を時計回り(あるいは反時計回り)に連結する(図４(c)：４１１)。この処理によって抽出された多角形をひとつの区画とみなし、この区画に対してＩＤを与えて区画データに登録する(図４(ｄ))。区画の頂点の数は、シード点から伸ばす直線数に依存する。例えば、図４では直線の数を８本としている。この直線数を増やすことにより、より精密な区画領域の抽出を行うことができる。図４(ｄ)に示す処理においては、台形に近い多角形の形状が農地の区画として抽出されている。以上の処理を、全てのシード点に対して行い、画像内の全ての区画を抽出する。
【００４２】
(3)概要区画解析処理（２０４）
概要区画解析処理（２０４）では、概要区画抽出処理（２０３）により得られた各区画データから詳細区画抽出処理で利用する特徴量を抽出し、概要区画特徴データを作成する。本実施の形態による土地区画データ作成システムでは、農地の区画が畦道を挟んでいくつか並んで正方形や長方形、または台形をしていることが多いことに注目し、畦道は農地を区画する辺に平行または垂直であることを利用するために、区画の各辺の方向とその分布を調べ、区画を構成する特徴ベクトルを抽出する。
【００４３】
図５は、概要区画解析処理（２０４）の流れを示すフローチャート図であり、図６は、その内容を画像で示した図である。
【００４４】
処理を開始すると、まず入力として概要区画データを読み込み(ステップ５０１)、格納されている区画データのうちでステップ５０３以下の処理を行なっていないエントリを１つ選択する(ステップ５０２)。そのエントリに含まれる頂点(座標)を配列の順に繋げたそれぞれの線をベクトルとみなして以下の処理を行う。
【００４５】
図６(a)に示す表示６０１は区画を含むエッジデータ６０７において、頂点６０５を配列順に線６０３で繋げた概要区画データ３０３（図１）を示した図である。ベクトルは、図６（ａ）中で区画を構成する各辺を意味する。まず、隣り合ったベクトルであって方向が近いベクトルを１つのベクトルに統合する(ステップ５０３)。すなわち、隣り合ったベクトルの成す角度θを求め、その角度θが一定値以下である場合は、ベクトルを繋ぐ点（頂点）を取り除き、２つのベクトルを１つのベクトルに統合する。図６（ｂ）は図６（ａ）の区画データに対してベクトルの統合処理を行い、頂点６１１のみを残し、いくつかのベクトルを統合した結果を示す図である。頂点が取り除かれて頂点６１１の数が減少していることがわかる。
【００４６】
次に、統合後のベクトルであって一定以上の長さのベクトルを集計し、ヒストグラムを作成する(ステップ５０４)。すなわち、ベクトルの角度で階級を分類し、ベクトルの長さを度数とみなして集計する。図６（ｃ）中の細い矢印６２１・６３１は、図６（ｂ）の区画データのベクトルを集計したものを示している。実際の処理では、テーブル形式で処理されるが、方向の偏りを分かりやすくするため、図６（ｃ）において円形に表現している。
【００４７】
次に、集計されたヒストグラムの分布を調べ、ピークの数を調べる(ステップ５０５)。すなわち、ある範囲の角度のベクトルが他の角度のものに比べて特に突出して集計されているかどうか、また突出している階級がいくつあるかを調べる。これにより、調べている区画が大まかにどのような形状をしているか判定することが出来る。ピークが２つの場合は区画が正方形、長方形又は平行四辺形をしていることを示し、三つの場合は三角形か台形をしていることを示す。それ以上の場合は、五角形以上の多角形、一つ以下の場合は円形に近い形状と推定できる。
【００４８】
図６（ｃ）中の太い矢印６２３・６３３は、集計されたベクトルのピークを示しており、この例では集計結果のピークは二つでありほぼ直交する２つの特徴ベクトルが得られることから、区画は長方形をしていると推定できる。
【００４９】
ステップ５０５における処理の結果として、ピークの数が２つ又は３つの場合、その区画に対して特徴データを登録する(ステップ５０６)。すなわち、調べた区画ＩＤの特徴データとして、ピークの数を特徴ベクトル数として、ピークの方向を特徴ベクトルとして、それぞれ登録する。ピークが、上記に該当しない他の数の場合は特徴データを入力せず、次のエントリのチェックに移る。そして、ステップ５０３〜５０６の処理を、概要区画データに含まれる全てのエントリに対してチェックが終了するまで繰り返し、この処理を終了する(ステップ５０７)。
【００５０】
(4)詳細区画抽出処理（２０５）
詳細区画抽出処理（２０５）では、エッジ抽出処理で得られた詳細エッジデータと、概要区画抽出処理で得られた概要区画データと、概要区画解析処理（２０４）で得られた概要区画特徴データと、から、より正確な区画の境界を抽出し、詳細区画データを作成する。
【００５１】
従来の手法では、詳細エッジデータのようなノイズを多く含むエッジデータは自動的な解析への利用には適していなかったが、本実施の形態によるシステムでは、農地区画の境界の多くを占める畦道が、道路や河川に平行または垂直なものが多いという特徴と概要区画特徴データを利用して、エッジ画像からノイズを分離して畦道のエッジのみを抽出する。
【００５２】
図７は、詳細区画抽出処理の流れを示すフローチャート図であり、図８は、その内容を画像で示した図である。この処理では、まず入力として概要区画データ（３０３）、概要区画特徴データ（３０４）、及び詳細エッジデータ（３０５）を読み込み(ステップ７０１)、概要区画データの中で７０３〜７０９までの処理を行っていないエントリを１つ選択する(ステップ７０２)。次いで、概要区画特徴データ（３０４）の中にそのエントリの区画ＩＤに対応するデータがあるか否かを調べ、対応するデータがない場合には、詳細区画データ（３０６）に概要区画データ（３０３）をそのまま登録し(ステップ７０９)、他のエントリのチェックに移る。
【００５３】
対応する概要区画特徴データ（３０４）が存在する場合には、その区画データに対してステップ７０４〜７０９までの処理を行う。
【００５４】
図８（ａ）は、概要区画データの一つのエントリ８０１と、対応する特徴データ８０３・８０５を示す図である。この図では、区画データに対応する２本の特徴ベクトル８１１・８１３があり、太い矢印で示されている方向（短辺と長辺に沿った２つの白抜き矢印）がその方向である。次に、特徴データに登録されている特徴ベクトルのうちステップ７０５の処理を行っていないベクトルを１つ選び(ステップ７０４)、そのベクトルに垂直及び平行な方向で区画内の詳細エッジデータを集計する(ステップ７０５)。すなわち、区画データに対応する位置の詳細エッジを選んだ特徴ベクトルの垂直及び平行な方向に投影し、その方向で見たエッジの量を調べる。この集計値をグラフ化し、ピークを見つけることによって、ノイズが多い場合でも、特徴ベクトルに垂直または平行な方向にエッジが集中している場合にそれを見つけ出すことが出来る。図８（ｂ）に示す表示８２１は、図８（ａ）で示した区画８２３を、１つの特徴ベクトルに垂直及び平行な方向でエッジを集計している状況を図示したものである。図示した例では、縦方向に集計したエッジでピークＰ２・Ｐ３が見られ、畦道と思われるエッジの集合を発見している。そのたのピークＰ１・Ｐ４・Ｐ６・Ｐ５は、区画に対応すると考えられる。
【００５５】
ステップ７０５の処理を全ての特徴ベクトルに対して繰り返し(ステップ７０６)、それぞれの集計したデータでピークが存在するか否かを調べる(ステップ７０７)。ある位置においてピークが存在するか否かの判断は、集計した軸に対するエッジの量の平均値との差や偏差値をピークの強さとみなし、それが適切に定めた閾値との比較によって真にピークと言えるかどうかで判定することが出来る。
【００５６】
この処理においてピークが見つかった場合は、ピークの位置には畦道と推定されるエッジの集合があると考えられる。そこで、ピークの位置Ｐ２、Ｐ３で区画データを分割し、それぞれを一つの区画として詳細区画データに登録する(ステップ７０８)。もし、複数の投影方向でピークが見つかった場合には、得られたピークの強さを比較し、最大の一方向のみで区画を分割する。１つの投影方向で複数のピークが見つかった場合は、それぞれのピークで分割する。図８（ｃ）は、図８（ａ）で示した概要区画データが分割されて詳細区画データに登録される様子を示した図である（８４１）。
【００５７】
太い矢印で示した特徴ベクトル８４３に平行な方向に二本のピークＰ２・Ｐ３が見つかったため、それぞれで１つの区画８２３を確定する区画線８５３を線８４５と線８４７とにより３分割して３つの区画が得られる。ピークが見つからない場合には、詳細区画データに概要区画データをそのまま登録する(ステップ７０９)。そして、ステップ７０３〜７０９の処理を概要区画データに含まれる全てのエントリに対してチェックするまで繰り返し、処理を終了する(ステップ７１０)。
【００５８】
＜まとめ＞
以上のように、本実施の形態によれば、地表を撮影して得られた地理画像データを解析し、農地区画データを作成する農地区画データ作成システムにおいて、地理画像データからエッジを抽出するエッジ抽出処理部と、エッジ抽出結果から大まかな区画を抽出する概要区画抽出処理部と、得られた大まかな区画にみられる特徴ベクトルを抽出する概要区画解析部と、それまでの処理で得られたデータを元に詳細な区画を抽出する詳細区画抽出部と、を備える。
【００５９】
これにより、農地判定の処理を自動化できると共に、主観的な判定及びそれによる判定のばらつきを防止でき、客観的な農地判定を実現することができるようになる。
【００６０】
但し、本実施の形態に限定されるものではない。当業者であれば、本発明の要旨を逸脱しない範囲内において、実施形態にかかる発明の構成及び機能に様々な変更、改良を加えることが可能である。
【００６１】
尚、本実施の形態において説明した上記の機能は、ソフトウェアプログラムによっても実現できる。この場合、プログラムコードを記録した記憶媒体をシステム或は装置に提供し、そのシステム或は装置のコンピュータ（又はＣＰＵやＭＰＵ）が記憶媒体に格納されたプログラムコードを読み出すことができるようにすれば良い。この場合、記憶媒体から読み出されたプログラムコード自体が前述した実施の形態による機能を実現することになり、そのプログラムコード自体及びそれを記憶した記憶媒体は本発明によるシステムを構成することになる。このようなプログラムコードを供給するための記憶媒体としては、例えば、フロッピィ（登録商標）ディスク、ＣＤ−ＲＯＭ、ＤＶＤ−ＲＯＭ、ハードディスク、光ディスク、光磁気ディスク、ＣＤ−Ｒ、磁気テープ、不揮発性のメモリカード、ＲＯＭなどが用いられる。
【００６２】
また、プログラムコードの指示に基づき、コンピュータ上で稼動しているＯＳ（オペレーティングシステム）などが実際の処理の一部又は全部を行い、その処理によって前述した実施の形態の機能が実現されるようにしてもよい。さらに、記憶媒体から読み出されたプログラムコードが、コンピュータ上のメモリに書きこまれた後、そのプログラムコードの指示に基づき、コンピュータのＣＰＵなどが実際の処理の一部又は全部を行い、その処理によって前述した実施の形態の機能が実現されるようにしてもよい。
【００６３】
また、実施の形態の機能を実現するソフトウェアのプログラムコードがネットワークを介して配信されることにより、システム又は装置のハードディスクやメモリ等の記憶手段又はＣＤ-ＲＷ、ＣＤ-Ｒ等の記憶媒体に格納され、そのシステム又は装置のコンピュータ(又はＣＰＵやＭＰＵ)が当該記憶手段や当該記憶媒体に格納されたプログラムコードを読み出して実行することによっても、達成されるようにしてもよい。
【産業上の利用可能性】
【００６４】
本発明は、農地区画データ作成装置として利用可能である。
【図面の簡単な説明】
【００６５】
【図１】本発明の一実施の形態による土地区画データの画像解析システムの一構成例を示す機能ブロック図である。
【図２Ａ】図１に示す画像解析システムにより利用される各種データの一例を示す図であり、図２Ａ（ａ）は画像データの一例、図２Ａ（ｂ）は、概要区画データに相当するエッジ画像の一例を、図２Ａ（ｃ）は詳細区画データに相当するエッジ画像の一例を示す図である。
【図２Ｂ】図１に示す画像解析システムにより利用される各種データの一例を示す図であり、図２Ｂ（ｄ）は、概要区画データの一例を示す図であり、図２Ｂ（ｅ）は図２Ｂ(ｄ)の区画データを画像として表現した一例を示す図であり、図２Ｂ（ｆ）は、詳細区画データの一例を示す図である。
【図３】図１に示す画像解析システムの区画抽出プログラムの各処理部の流れと、その入出力データを概略的に示す図である。
【図４】図４（ａ）〜（ｄ）までは、図１に示す地理解析システムの区画抽出処理部によるデータの処理例を示す図である。
【図５】図１に示す地理解析システムの概要区画解析処理部による処理の流れを示すフローチャート図である。
【図６】図１に示す概要区画解析処理部によるデータの処理の概要を示す図である。
【図７】図１に示す画像解析システムの詳細区画抽出処理部による処理の流れを示すフローチャート図である。
【図８】図１に示す詳細区画抽出処理部によるデータの処理の概要を示す図である。
【符号の説明】
【００６６】
１００…入出力部、１０１…入力装置、１０２…表示装置、１０３…出力装置、２００…処理装置、２０１…区画抽出プログラム、２０２…エッジ抽出処理部、２０３…概要区画抽出処理部、２０４…概要区画解析処理部、２０５…詳細区画抽出処理部、３００…記憶装置、３０１…画像データ、３０２…概要エッジデータ、３０３…概要区画データ、３０４…概要区画特徴データ、３０５…詳細エッジデータ、３０６…詳細区画データ。
【図１】
2008234163.tif 000003
【図２Ａ】
2008234163.tif 000004
【図２Ｂ】
2008234163.tif 000005
【図３】
2008234163.tif 000006
【図４】
2008234163.tif 000007
【図５】
2008234163.tif 000008
【図６】
2008234163.tif 000009
【図７】
2008234163.tif 000010
【図８】
2008234163.tif 000011