(57)【要約】【課題】  市街地の形態を定量的に解析し、分類する。【解決手段】  所定の形態の市街地が撮影されている航空写真に基づいて生成された画像データが格納されている比較画像データ格納部１０と、任意の市街地が撮影されている航空写真に基づいて生成された入力画像データの全体のフラクタルシーケンスを算出し、比較画像データ格納部１０に格納されている格納画像データの全体のフラクタルシーケンスを算出するフラクタルシーケンス演算部１１と、フラクタルシーケンス演算部１１により算出された、入力画像データの全体のフラクタルシーケンスと、格納画像データの全体のフラクタルシーケンスとに基づいて、画像特徴距離ｄを算出する画像特徴距離演算部１２と、画像特徴距離演算部１２により算出された画像特徴距離ｄに基づいて、任意の市街地の形態を分類する形態分類部１３とを備える。【選択図】  図２
2007139421.tif 000002
【特許請求の範囲】
【請求項１】
所定の形態の市街地が撮影されている航空写真に基づいて生成された画像データが格納されている画像データ格納手段と、
任意の市街地が撮影されている航空写真に基づいて生成された入力画像データのフラクタルシーケンスを算出し、上記画像データ格納手段に格納されている上記画像データのフラクタルシーケンスを算出するフラクタルシーケンス演算手段と、
上記フラクタルシーケンス演算手段により算出された、上記入力画像データのフラクタルシーケンスと、上記画像データのフラクタルシーケンスとに基づいて、画像特徴距離を算出する画像特徴距離演算手段と、
上記画像特徴距離演算手段により算出された上記画像特徴距離に基づいて、上記任意の市街地の形態を分類する形態分類手段とを備えることを特徴とする形態分類装置。
【請求項２】
任意の市街地を撮影した航空写真に基づいて生成された入力画像データのフラクタルシーケンスを算出し、
画像データ格納部に格納されており、所定の形態の市街地を撮影した航空写真に基づいて生成された格納画像データのフラクタルシーケンスを算出し、
上記入力画像データのフラクタルシーケンスと、上記格納画像データのフラクタルシーケンスとに基づいて画像特徴距離を算出し、
上記画像特徴距離に基づいて、上記任意の市街地の形態を分類することを特徴とする形態分類方法。
【請求項３】
所定の形態の市街地を撮影した航空写真に基づいて生成された画像データが格納されている画像データ格納手段と、
上記航空写真が撮影された時期とは異なる時期に撮影された上記市街地の航空写真に基づいて生成された入力画像データを所定のサイズの小領域に分割し、上記画像データ格納手段に格納されている上記画像データを上記所定のサイズと同等サイズの小領域に分割する小領域分割手段と、
上記小領域分割手段により小領域に分割された上記入力画像データを、小領域ごとにフラクタルシーケンスを算出し、上記小領域分割手段により小領域に分割された上記画像データを、小領域ごとにフラクタルシーケンスを算出するフラクタルシーケンス演算手段と、
上記フラクタルシーケンス演算手段により算出された、上記入力画像データの小領域ごとのフラクタルシーケンスと、上記画像データの小領域ごとのフラクタルシーケンスとに基づいて、対応する小領域ごとに画像特徴距離を算出する画像特徴距離演算手段と、
上記画像特徴距離演算手段により算出された上記画像特徴距離に基づいて、上記市街地の変化領域を抽出する変化領域抽出手段とを備えることを特徴とする変化領域抽出装置。
【請求項４】
上記変化抽出手段は、変化領域の抽出を画像データの光強度分布を走査することにより行うことを特徴とする請求項３記載の変化領域抽出装置。
【請求項５】
画像データ格納部に格納されており、所定の形態の市街地が撮影されている航空写真に基づいて生成された格納画像データを所定のサイズの小領域に分割し、
上記航空写真が撮影された時期とは異なる時期に撮影された上記市街地の航空写真に基づいて生成された入力画像データを所定のサイズの小領域に分割し、
小領域に分割された上記入力画像データを、小領域ごとにフラクタルシーケンスを算出し、
小領域に分割された上記格納画像データを、小領域ごとにフラクタルシーケンスを算出し、
上記入力画像データの小領域ごとのフラクタルシーケンスと、上記格納画像データの小領域ごとのフラクタルシーケンスとに基づいて、対応する小領域ごとに画像特徴距離を算出し、
算出された上記画像特徴距離に基づいて、上記市街地の変化領域を抽出することを特徴とする変化領域抽出方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、任意の市街地を撮影した航空写真に基づき、当該市街地の形態を分類する形態分類装置及び方法、並びに当該市街地の変化領域の抽出を行う変化領域抽出装置及び方法に関する。
【背景技術】
【０００２】
都市計画を実施する際は、周辺の環境への影響を考慮した事前の評価が必要であり、そのためには複雑な都市の構造を定量的に解析することが必要である。また、該当地域における建物や道路の工事、災害による建物の倒壊などによる環境の変化も迅速に把握する必要がある。そのような場合、従来から広く用いられているのは、市街地を撮影した航空写真を人間の目視により判別する人的手法である。この場合、莫大な時間と労力が必要となるため効率が悪く、多大なコストが必要となる。そこで最近では、パターン認識を応用した手法が提案されており、計算機による自動化の試みがなされている（例えば、特許文献１を参照）。
【０００３】
【特許文献１】特開平１１−８３４７８号公報
【発明の開示】
【発明が解決しようとする課題】
【０００４】
しかしながら、パターン認識を利用した手法の場合、一般にアルゴリズムが複雑であり、装置に実装した場合にはコストがかかり、莫大な処理時間を要する。さらに、判定結果は必ずしも人間の感覚とは一致しないという問題点がある。
【０００５】
そこで、本発明は、上述した点に鑑みてなされたものであり、人的手法により作業を行った場合と同等の信頼性の高い解析を行うことができ、かつ迅速・容易に市街地の形態を分類する形態分類装置及び方法、並びに変化領域抽出装置及び方法を提供することを目的とする。
【課題を解決するための手段】
【０００６】
本願発明に係る形態分類装置は、上述したような課題を解決するために、所定の形態の市街地が撮影されている航空写真に基づいて生成された画像データが格納されている画像データ格納手段と、任意の市街地が撮影されている航空写真に基づいて生成された入力画像データのフラクタルシーケンスを算出し、上記画像データ格納手段に格納されている上記画像データのフラクタルシーケンスを算出するフラクタルシーケンス演算手段と、上記フラクタルシーケンス演算手段により算出された、上記入力画像データのフラクタルシーケンスと、上記画像データのフラクタルシーケンスとに基づいて、画像特徴距離を算出する画像特徴距離演算手段と、上記画像特徴距離演算手段により算出された上記画像特徴距離に基づいて、上記任意の市街地の形態を分類する形態分類手段とを備える。
【０００７】
また、本発明に係る形態分類方法は、上述したような課題を解決するために、任意の市街地を撮影した航空写真に基づいて生成された入力画像データのフラクタルシーケンスを算出し、画像データ格納部に格納されており、所定の形態の市街地を撮影した航空写真に基づいて生成された格納画像データのフラクタルシーケンスを算出し、上記入力画像データのフラクタルシーケンスと、上記格納画像データのフラクタルシーケンスとに基づいて画像特徴距離を算出し、上記画像特徴距離に基づいて、上記任意の市街地の形態を分類する。
【０００８】
また、本発明に係る変化領域抽出装置は、上述したような課題を解決するために、所定の形態の市街地を撮影した航空写真に基づいて生成された画像データが格納されている画像データ格納手段と、上記航空写真が撮影された時期とは異なる時期に撮影された上記市街地の航空写真に基づいて生成された入力画像データを所定のサイズの小領域に分割し、上記画像データ格納手段に格納されている上記画像データを上記所定のサイズと同等サイズの小領域に分割する小領域分割手段と、上記小領域分割手段により小領域に分割された上記入力画像データを、小領域ごとにフラクタルシーケンスを算出し、上記小領域分割手段により小領域に分割された上記画像データを、小領域ごとにフラクタルシーケンスを算出するフラクタルシーケンス演算手段と、上記フラクタルシーケンス演算手段により算出された、上記入力画像データの小領域ごとのフラクタルシーケンスと、上記画像データの小領域ごとのフラクタルシーケンスとに基づいて、対応する小領域ごとに画像特徴距離を算出する画像特徴距離演算手段と、上記画像特徴距離演算手段により算出された上記画像特徴距離に基づいて、上記市街地の変化領域を抽出する変化領域抽出手段とを備える。
【０００９】
また、本発明に係る変化領域抽出方法は、上述したような課題を解決するために、画像データ格納部に格納されており、所定の形態の市街地が撮影されている航空写真に基づいて生成された格納画像データを所定のサイズの小領域に分割し、上記航空写真が撮影された時期とは異なる時期に撮影された上記市街地の航空写真に基づいて生成された入力画像データを所定のサイズの小領域に分割し、小領域に分割された上記入力画像データを、小領域ごとにフラクタルシーケンスを算出し、小領域に分割された上記格納画像データを、小領域ごとにフラクタルシーケンスを算出し、上記入力画像データの小領域ごとのフラクタルシーケンスと、上記格納画像データの小領域ごとのフラクタルシーケンスとに基づいて、対応する小領域ごとに画像特徴距離を算出し、算出された上記画像特徴距離に基づいて、上記市街地の変化領域を抽出する。
【発明の効果】
【００１０】
本発明によれば、人手により作業を行った場合に匹敵する、信頼性の高い市街地の形態に関する解析を迅速かつ容易に行うことができ、都市計画における事前調査や、地図情報の作成、及び災害時の被災地域の把握等に役立てることができる。
【発明を実施するための最良の形態】
【００１１】
本発明は、市街地の形態を分類する形態分類装置及び方法、並びに市街地の変化領域の抽出を行う変化領域抽出装置及び方法に適用される。
【００１２】
本発明の発明者等は、フラクタル理論に基づいて、任意の市街地の航空写真のフラクタルシーケンスと、サンプルとなる市街地の航空写真のフラクタルシーケンスとに基づいて、画像特徴距離ｄを求め、当該画像特徴距離ｄにより任意の市街地の形態を定量化できることを解明し、本発明を創作した。
【００１３】
以下に本発明の実施の形態について図面を参照しながら詳細に説明する。
【００１４】
まず、フラクタルシーケンスの算出方法と、画像特徴距離ｄの算出方法について説明する。なお、以下では、画像データの光強度分布（航空写真の画素値）をf(x,y)で表すこととする。
【００１５】
まず、ｒ×ｒ[Pixel]の任意の小領域を考え、その四隅における光強度｛f(xi,yj),f(xi+r,yj),f(xi,yj+r),f(xi+r,yj+r)｝（但し、i,j=0,1,2,・・・,N-1）の最小値をf0とし、その四隅の光強度からf0を差し引き、それぞれを頂点とする小立方体の仮想的な体積をVr(xi,yj)とし、このときの体積の平均V(r)を（１）式に基づいて求める（図１）。
【００１６】
【数１】
2007139421.tif 000003
【００１７】
ここで、小領域｛r0,r1,・・・,rM-1｝に対する体積の平均｛V(r0),V(r1),・・・,V(rM-1)｝を求めると、フラクタルシーケンスは、｛logV(r0),logV(r1),・・・,logV(rM-1)｝dで求められる。
【００１８】
つぎに、解像度の等しい二つの画像データのフラクタルシーケンスをそれぞれS=｛s0,s1,・・・,sM-1｝,T=｛t0,t1,・・・,tM-1｝とすると、（２）式に基づいて画像特徴距離ｄが算出される。
【００１９】
【数２】
2007139421.tif 000004
【００２０】
画像特徴距離ｄは、二つの画像の類似度を表し、例えば、二つの画像が同一であればｄ＝０となる。一方、二つの画像が異なるほど画像特徴距離ｄの値は大きくなり、その結果は、人間の視覚的感覚に近いものとなる。
【００２１】
＜形態分類装置に関して＞
形態分類装置１は、図２に示すように、所定の形態の市街地が撮影されている航空写真に基づいて生成された画像データが格納されている比較画像データ格納部１０と、任意の市街地が撮影されている航空写真に基づいて生成された入力画像データの全体のフラクタルシーケンスを算出し、比較画像データ格納部１０に格納されている格納画像データの全体のフラクタルシーケンスを算出するフラクタルシーケンス演算部１１と、フラクタルシーケンス演算部１１により算出された、入力画像データの全体のフラクタルシーケンスと、格納画像データの全体のフラクタルシーケンスとに基づいて、画像特徴距離ｄを算出する画像特徴距離演算部１２と、画像特徴距離演算部１２により算出された画像特徴距離ｄに基づいて、任意の市街地の形態を分類する形態分類部１３とを備える。
【００２２】
比較画像データ格納部１０は、例えば、形態が格子状である市街地の航空写真に基づいて生成されたグレースケールの格子状画像データが格納されている第１の形状の比較画像データ格納部２０と、形態が同心円状である市街地の航空写真に基づいて生成されたグレースケールの同心円状画像データが格納されている第２の形状の比較画像データ格納部２１と、形態が放射状である市街地の航空写真に基づいて生成されたグレースケールの放射状画像データが格納されている第３の形状の比較画像データ格納部２２とを備える。なお、比較画像データ格納部１０は、格子状、同心円状及び放射状以外の形態の市街地の航空写真に基づいて生成されたグレースケールの比較画像データが格納されている他の格納部を有していても良い。なお、フラクタルシーケンス演算部１１に入力される入力画像データは、フラクタル性のある画像データである。
【００２３】
フラクタルシーケンス演算部１１は、フラクタル性を有する入力画像データの全体のフラクタルシーケンスを算出する。また、フラクタルシーケンス演算部１１は、第１の形状の画像データ格納部２０に格納されている格子状画像データを読み出し、当該格子状画像データの全体のフラクタルシーケンスを算出する。また、フラクタルシーケンス演算部１１は、第２の形状の画像データ格納部２１に格納されている同心円状画像データを読み出し、当該同心円状画像データの全体のフラクタルシーケンスを算出する。また、フラクタルシーケンス演算部１１は、第３の形状の画像データ格納部２２に格納されている放射状画像データを読み出し、当該放射状画像データの全体のフラクタルシーケンスを算出する。
【００２４】
つぎに、画像特徴距離演算部１２の構成について説明する。画像特徴距離演算部１２は、フラクタルシーケンス演算部１１により算出された、入力画像データの全体のフラクタルシーケンスと、格子状画像データの全体のフラクタルシーケンスとに基づいて、画像特徴距離ｄ１を算出する。また、画像特徴距離演算部１２は、フラクタルシーケンス演算部１１により算出された、入力画像データの全体のフラクタルシーケンスと、同心円状画像データの全体のフラクタルシーケンスとに基づいて、画像特徴距離ｄ２を算出する。また、画像特徴距離演算部１２は、フラクタルシーケンス演算部１１により算出された、入力画像データの全体のフラクタルシーケンスと、放射状画像データの全体のフラクタルシーケンスとに基づいて、画像特徴距離ｄ３を算出する。
【００２５】
また、画像特徴距離ｄは、二つの画像の類似度を表し、例えば、二つの画像が同一であればｄ＝０となる。一方、二つの画像が異なるほど画像特徴距離ｄの値は大きくなり、その結果は、人間の視覚的感覚に近いものとなる。
【００２６】
また、画像特徴距離ｄは、上述したように二つの画像間の変化の度合いを表しており、主に、建物や道路等の３次元的な配置によって変化するが、車や人等のように極端に小さい物体は誤差（雑音分）となり、変化に寄与しない特徴を有している。
【００２７】
形態分類部１３は、画像特徴距離演算部１２により得られた画像特徴距離ｄ１と、画像特徴距離ｄ２と、画像特徴距離ｄ３とを比較し、当該比較結果に基づいて、入力画像データの元になった任意の市街地の形態を、例えば、格子状、同心円状又は放射状と分類する。
【００２８】
また、形態分類部１３により分類された形態に基づいて、入力画像データを入力画像データ格納部３０に登録する登録部３１を備える構成であっても良い。また、入力画像データ格納部３０は、形態分類部１３に分類された結果に応じて、例えば、格子状の形態の入力画像データを格納する第１の形状の入力画像データ格納部３０Ａと、同心円状の形態の入力画像データを格納する第２の形状の入力画像データ格納部３０Ｂと、放射状の形態の入力画像データを格納する第３の形状の入力画像データ格納部３０Ｃ等を備える。
【００２９】
＜具体例＞
ここで、図３に実際の航空写真に基づく画像データから画像特徴距離ｄを算出した結果を示す。また、図３では、航空写真の一辺の大きさは、約４００［ｍ］の正方形であり、当該航空写真から生成されたグレースケールの画像データの解像度は、５１２×５１２［pixel］である。また、撮影条件の違いによる影響を軽減するため、いずれの航空写真も１０月中旬の正午頃に撮影された航空写真を用いた。
【００３０】
図３は、「京都府中京区二条城町」を基準として、複数の市街地の画像特徴距離ｄを求めた結果である。「京都府中京区二条城町」は、「東京都北区小山町」（格子状の市街地）との画像特徴距離ｄが約０．０１となり、「東京都大田区田園調布」及び「神奈川県港北区日吉」（同心円状の市街地）との画像特徴距離ｄが約０．０３となり、「東京都豊島区南大塚」及び「東京都豊島区北大塚」（放射状の市街地）との画像特徴距離ｄが約０．１１となっている。したがって、「京都府中京区二条城町」の市街地の形態は、格子状の市街地である「東京都北区小山町」に近似していることが分かり、また、同心円状の市街地に比べ、放射状の市街地の方が画像特徴距離ｄが大きいことが分かる。このことは、人間の視覚的感覚とも良く一致している。
【００３１】
また、本発明の発明者等は、上述と同様に、他の形態の市街地を基準として画像特徴距離ｄを求めた結果、市街地の形態が相互に類似している場合には、画像特徴距離ｄが約０．０１となり、また、格子状の市街地と同心円状の市街地の場合には、画像特徴距離ｄが約０．０３となり、同心円状の市街地と放射状の市街地の場合には、画像特徴距離ｄが約０．０９となり、格子状の市街地と放射状の市街地の場合には、画像特徴距離ｄが約０．１１となることを解明した。
【００３２】
したがって、市街地の形態の違いが画像特徴距離ｄに十分に反映されているので、画像特徴距離ｄを用いることにより、市街地の形態分類が可能となるのである。
【００３３】
このようにして、本発明に係る形態分類装置１は、基準とする市街地と異なる複数の市街地との間の画像特徴距離ｄを求め、求めた画像特徴距離ｄの大小関係を比較し、その結果、画像特徴距離ｄがほぼ等しい場合には、同じ形態の市街地として分類することができ、画像特徴距離ｄが異なる場合には、別の形態の市街地として分類することができる。
【００３４】
＜変化領域抽出装置に関して＞
また、画像特徴距離ｄは、上述したように二つの画像間の変化の度合いを表している。したがって、本願発明によれば、航空写真に基づく二つの画像データをそれぞれ同等サイズの小領域に分割し、それぞれ対応する小領域ごとの画像特徴距離ｄを算出することにより、精密に変化領域の抽出を行うことができる。以下に、変化領域抽出装置の構成について説明する。
【００３５】
変化領域抽出装置２は、図４に示すように、ある年代において、任意の市街地を撮影した航空写真に基づいて生成された画像データ（以下、格納画像データという。）が格納されている画像データ格納部４０と、当該航空写真が撮影された時期とは異なる時期に撮影された同一の市街地の航空写真に基づいて生成された入力画像データを所定のサイズの小領域に分割し、また、画像データ格納部４０に格納されている格納画像データを所定のサイズと同等サイズの小領域に分割する小領域分割部４１と、小領域分割部４１により小領域に分割された入力画像データを、小領域ごとにフラクタルシーケンスを算出し、また、小領域分割部４１により小領域に分割された格納画像データを、小領域ごとにフラクタルシーケンスを算出するフラクタルシーケンス演算部４２と、フラクタルシーケンス演算部４２により算出された、入力画像データの小領域ごとのフラクタルシーケンスと、格納画像データの小領域ごとのフラクタルシーケンスとに基づいて、対応する小領域ごとに画像特徴距離を算出する画像特徴距離演算部４３と、画像特徴距離演算部４３により算出された画像特徴距離に基づいて、市街地の変化領域を抽出する変化領域抽出部４４とを備える。なお、小領域分割部４１に入力される入力画像データは、フラクタル性のある画像データである。
【００３６】
画像データ格納部４０には、撮影年代の異なる市街地の航空写真に基づいて生成された画像データが格納されている。
【００３７】
小領域分割部４１は、フラクタル性を有する入力画像データを所定のサイズの小領域に分割する。また、小領域分割部４１は、画像データ格納部４０に格納されている格納画像データを読み出し、当該格納画像データを所定のサイズの小領域に分割する。
【００３８】
フラクタルシーケンス演算部４２は、小領域に分割された入力画像データを、小領域ごとにフラクタルシーケンスを算出する。また、フラクタルシーケンス演算部４２は、小領域に分割された格納画像データを、小領域ごとにフラクタルシーケンスを算出する。
【００３９】
画像特徴距離演算部４３は、フラクタルシーケンス演算部４２により算出された、入力画像データの小領域ごとのフラクタルシーケンスと、格納画像データの小領域ごとのフラクタルシーケンスとに基づいて、それぞれ対応する小領域ごとの画像特徴距離を算出する。
【００４０】
ここで、フラクタルシーケンス検算部４２と、画像特徴距離演算部４３による一連の動作について説明する。
【００４１】
フラクタルシーケンス演算部４２は、一の画像データを図５（Ａ）に示すように、大きさｕ×ｕ［pixel］の小領域{pn}（但し、n=1,2,・・・L2）に分割し、それぞれのフラクタルシーケンスを求める。また、フラクタルシーケンス演算部４２は、他の画像データを図５（Ｂ）に示すように、大きさｕ×ｕ［pixel］の小領域{qn}（但し、n=1,2,・・・L2）に分割し、それぞれのフラクタルシーケンスを求める。
【００４２】
そして、画像特徴距離演算部４３は、一の画像データの小領域{p1}のフラクタルシーケンスと、小領域{p1}に対応する他の画像データの小領域{q1}のフラクタルシーケンスとに基づいて画像特徴距離ｄ１を算出し、順次、対応する位置にある小領域のフラクタルシーケンスに基づいて画像特徴距離ｄｎ（但し、n=1,2,・・・L2）を算出する。
【００４３】
変化領域抽出部４４は、画像特徴距離演算部４３により算出された小領域ごとの画像特徴距離ｄに基づいて、一定の値以上の画像特徴距離ｄになっている小領域のみを抽出し、抽出した小領域の光強度分布を算出する。変化領域抽出部４４は、算出した値が一定値以上の小領域を変化領域として抽出する。
【００４４】
また、小領域分割部４１により、画像データを小領域に分割する際のサイズを設定するサイズ設定部４５を備える構成であっても良い。なお、サイズ設定部４５で設定されるサイズが小さい場合には、より精密にフラクタルシーケンスを求められるメリットがあるが、演算に時間が掛かってしまうデメリットがあり、一方、サイズを大きくすれば、演算時間は少なくて済むが、フラクタルシーケンスの値が大まかになってしまう。したがって、設定サイズは、どの程度の精密さが要求されるのかによって適宜決定される。
【００４５】
＜変化領域の抽出＞
ここで、図６に示すように、過去の画像データの小領域{pＬ}の内部にある大きさ［wp×hp］の建物が、現在の画像データの小領域{qＬ}においては、大きさ［wｑ×hｑ］の建物に変化している場合を考える。
【００４６】
一般に、地面と建物上部の光強度分布は、かなり異なるため、例えば、小領域の中央部から所定量だけｘ方向、ｙ方向に沿って光強度分布を求めると、図６（Ａ）及び図６（Ｂ）のように、地面と建物上部の境界で光強度分布は変化する。
【００４７】
また、図６（Ｃ）には、過去の画像データの小領域{pＬ}及び現在の画像データの小領域{qＬ}について、ｘ方向に沿って走査して求めた光強度分布の比較例をそれぞれ示す。
【００４８】
したがって、光強度分布が異なる範囲（ｗｐ≦ｘ≦ｗｑ）を調べることにより、変化している小領域を抽出することができる。
【００４９】
このようにして、本発明に係る変化領域抽出装置２は、画像データにおいて、目視では判別しにくい小領域の変化を抽出することができるので、例えば、ある市街地において何らかの災害が発生した場合に、災害前の航空写真による画像データと、災害後の航空写真による画像データとに基づいて、小領域ごとの画像特徴距離ｄ及び光強度分布を算出することにより、市街地の変化を迅速に検出することができ、災害場所の把握に利用することができる。
【００５０】
また、本発明を応用すれば、例えば、ある場所の過去の航空写真による画像データと、同一場所の現在の航空写真による画像データとに基づいて、小領域ごとの画像特徴距離ｄ及び光強度分布を算出し、変化があった小領域を抽出し、抽出された小領域に対応する画像データのみを更新することにより、地図の作成に役立てることができる。
【００５１】
また、本発明は、図面を参照して説明した上述の実施例に限定されるものではなく、添付の請求の範囲及びその主旨を逸脱することなく、様々な変更、置換又はその同等のものを行うことができることは勿論である。
【図面の簡単な説明】
【００５２】
【図１】小立方体の仮想体積を示す図である。
【図２】本発明を適用した形態分類装置の構造を示すブロック図である。
【図３】実際の航空写真に基づく画像データから画像特徴距離ｄを算出した結果を示す図である。
【図４】本発明を適用した変化領域抽出装置の構造を示すブロック図である。
【図５】フラクタルシーケンス演算部により行われる画像の分割の説明に供する図である。
【図６】変化領域の抽出を行う説明に供する図である。
【符号の説明】
【００５３】
１  形態分類装置、２  変化領域抽出装置、１０  比較画像データ格納部、１１  フラクタルシーケンス演算部、１２  画像特徴距離演算部、１３  形態分類部、２０  第１の形状の比較画像データ格納部、２１  第２の形状の比較画像データ格納部、２２  第３の形状の比較画像データ格納部、３０  入力画像データ格納部、３１  登録部、４０  画像データ格納部、４１  小領域分割部、４２  フラクタルシーケンス演算部、４３  画像特徴距離演算部、４４  変化領域抽出部、４５  サイズ設定部
【図１】
2007139421.tif 000005
【図２】
2007139421.tif 000006
【図３】
2007139421.tif 000007
【図４】
2007139421.tif 000008
【図５】
2007139421.tif 000009
【図６】
2007139421.tif 000010