(57)【要約】【課題】多くの種類の家屋異動を精度よく推定できる家屋異動推定装置及びプログラムを提供する。【解決手段】家屋オブジェクト抽出部５４は、地物を含んだ地上表層の標高のデータに基づいて、家屋異動推定領域内に家屋が表された家屋領域を抽出する。特定部５６は、第１の時期の標高のデータに基づいて抽出される家屋領域について、第２の時期における標高のデータが示す標高、第１の時期と第２の時期との間における標高のデータが示す標高の比較結果、又は、第１の時期と第２の時期に撮影された画像に表された色情報の比較結果の少なくとも１つを特定する。家屋異動推定部５８は、特定される少なくとも１つに基づいて、第１の時期と第２の時期との間における家屋についての家屋異動を推定する。【選択図】図１１
2020016923.tif 000002
【特許請求の範囲】
【請求項１】
地物を含んだ地上表層の標高のデータに基づいて、家屋異動推定領域内に家屋が表された家屋領域を抽出する抽出手段と、
第１の時期の前記標高のデータに基づいて抽出される前記家屋領域について、第２の時期における前記標高のデータが示す標高、前記第１の時期と前記第２の時期との間における前記標高のデータが示す標高の比較結果、又は、前記第１の時期と前記第２の時期に撮影された画像に表された色情報の比較結果の少なくとも１つを特定する特定手段と、
特定される前記少なくとも１つに基づいて、前記第１の時期と前記第２の時期との間における前記家屋についての家屋異動を推定する家屋異動推定手段と、
を含むことを特徴とする家屋異動推定装置。
【請求項２】
前記家屋異動推定手段は、前記第１の時期と前記第２の時期との間における前記標高のデータが示す標高の比較結果に基づいて、前記家屋についての建替の発生の有無を推定する、
ことを特徴とする請求項１に記載の家屋異動推定装置。
【請求項３】
前記家屋異動推定手段は、前記第１の時期と前記第２の時期との間における前記標高のデータが示す標高の比較結果、及び、前記色情報の比較結果に基づいて、前記家屋についての建替の発生の有無を推定する、
ことを特徴とする請求項１に記載の家屋異動推定装置。
【請求項４】
前記家屋異動推定手段は、ＨＳＩ表色系により表現された前記色情報における色相成分の比較結果に基づいて、前記家屋についての家屋異動を推定する、
ことを特徴とする請求項１から３のいずれか一項に記載の家屋異動推定装置。
【請求項５】
前記第１の時期に撮影された画像又は前記第２の時期に撮影された画像の少なくとも一方に対して色調整を実行する色調整手段、をさらに含み、
前記家屋異動推定手段は、前記色調整が実行された画像に表された色情報の比較結果に基づいて、前記家屋についての家屋異動を推定する、
ことを特徴とする請求項１から４のいずれか一項に記載の家屋異動推定装置。
【請求項６】
地物を含んだ地上表層の標高のデータに基づいて、家屋異動推定領域内に家屋が表された家屋領域を抽出する手順、
第１の時期の前記標高のデータに基づいて抽出される前記家屋領域について、第２の時期における前記標高のデータが示す標高、前記第１の時期と前記第２の時期との間における前記標高のデータが示す標高の比較結果、又は、前記第１の時期と前記第２の時期に撮影された画像に表された色情報の比較結果の少なくとも１つを特定する手順、
特定される前記少なくとも１つに基づいて、前記第１の時期と前記第２の時期との間における前記家屋についての家屋異動を推定する手順、
をコンピュータに実行させることを特徴とするプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、家屋異動推定装置及びプログラムに関する。
【背景技術】
【０００２】
２時期における家屋等の固定資産の異動を推定する技術が存在する。
【０００３】
このような技術の一例として、特許文献１には、ＤＳＭ（Digital Surface Model）画像データやオルソフォト画像データを用いてメッシュ単位での比較結果に基づいて異動を判別する技術が記載されている。
【０００４】
また特許文献２には、異なる２時期についてのＤＳＭデータに基づいて特定される家屋オブジェクトを対比して家屋異動を判定する技術が記載されている。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２０１６−９９３１６号公報
【特許文献２】国際公開第２００９／５１２５８号
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかし特許文献１に記載の技術では、画像全体に対してメッシュ単位でＤＳＭ画像データやオルソフォト画像データが比較される。そのため、ＤＳＭデータのノイズ、位置合わせ（レジストレーション）におけるずれ、地形、樹木などの地物や移動体の影響によって、家屋の異動を精度よく推定できないことがあった。
【０００７】
また特許文献２に記載の技術では、家屋オブジェクトを対比して家屋異動を判定する。そのため例えば家屋の高さ変化が小さな建替や増改築などといった、家屋異動推定がうまくできない家屋異動の種類があった。
【０００８】
本発明は上記課題に鑑みてなされたものであって、その目的の１つは、多くの種類の家屋異動を精度よく推定できる家屋異動推定装置及びプログラムを提供することにある。
【課題を解決するための手段】
【０００９】
（１）本発明に係る家屋異動推定装置は、地物を含んだ地上表層の標高のデータに基づいて、家屋異動推定領域内に家屋が表された家屋領域を抽出する抽出手段と、第１の時期の前記標高のデータに基づいて抽出される前記家屋領域について、第２の時期における前記標高のデータが示す標高、前記第１の時期と前記第２の時期との間における前記標高のデータが示す標高の比較結果、又は、前記第１の時期と前記第２の時期に撮影された画像に表された色情報の比較結果の少なくとも１つを特定する特定手段と、特定される前記少なくとも１つに基づいて、前記第１の時期と前記第２の時期との間における前記家屋についての家屋異動を推定する家屋異動推定手段と、を含む。
【００１０】
（２）上記（１）に記載の家屋異動推定装置において、前記家屋異動推定手段は、前記第１の時期と前記第２の時期との間における前記標高のデータが示す標高の比較結果に基づいて、前記家屋についての建替の発生の有無を推定する構成とすることができる。
【００１１】
（３）上記（１）に記載の家屋異動推定装置において、前記家屋異動推定手段は、前記第１の時期と前記第２の時期との間における前記標高のデータが示す標高の比較結果、及び、前記色情報の比較結果に基づいて、前記家屋についての建替の発生の有無を推定する構成とすることができる。
【００１２】
（４）上記（１）〜（３）に記載の家屋異動推定装置において、前記家屋異動推定手段は、ＨＳＩ表色系により表現された前記色情報における色相成分の比較結果に基づいて、前記家屋についての家屋異動を推定する構成とすることができる。
【００１３】
（５）上記（１）〜（４）に記載の家屋異動推定装置において、前記第１の時期に撮影された画像又は前記第２の時期に撮影された画像の少なくとも一方に対して色調整を実行する色調整手段、をさらに含み、前記家屋異動推定手段は、前記色調整が実行された画像に表された色情報の比較結果に基づいて、前記家屋についての家屋異動を推定する構成とすることができる。
【００１４】
（６）本発明に係るプログラムは、地物を含んだ地上表層の標高のデータに基づいて、家屋異動推定領域内に家屋が表された家屋領域を抽出する手順、第１の時期の前記標高のデータに基づいて抽出される前記家屋領域について、第２の時期における前記標高のデータが示す標高、前記第１の時期と前記第２の時期との間における前記標高のデータが示す標高の比較結果、又は、前記第１の時期と前記第２の時期に撮影された画像に表された色情報の比較結果の少なくとも１つを特定する手順、特定される前記少なくとも１つに基づいて、前記第１の時期と前記第２の時期との間における前記家屋についての家屋異動を推定する手順、をコンピュータに実行させる。
【発明の効果】
【００１５】
本発明によれば、多くの種類の家屋異動を精度よく推定できる。
【図面の簡単な説明】
【００１６】
【図１】本発明の一実施形態に係る家屋異動推定装置の構成図である。
【図２Ａ】地上の模式的な垂直断面図の一例を示す図である。
【図２Ｂ】ＮＤＳＭの一例を示す図である。
【図３】時期ｔ１のＤＳＭデータに基づいて生成される植生除去処理後の対象領域におけるＮＤＳＭの凹凸の一例を、便宜的に二値化して表した図である。
【図４】時期ｔ２のＤＳＭデータに基づいて生成される植生除去処理後の対象領域におけるＮＤＳＭの凹凸の一例を、便宜的に二値化して表した図である。
【図５】図３に示すＮＤＳＭに対する領域分割処理の結果である複数の凸部領域の一例を示す図である。
【図６】図４に示すＮＤＳＭに対する領域分割処理の結果である複数の凸部領域の一例を示す図である。
【図７】時期ｔ１における対象領域内の家屋オブジェクトの一例を示す図である。
【図８】時期ｔ２における対象領域内の家屋オブジェクトの一例を示す図である。
【図９】時期ｔ１に対象領域を撮影したオルソ画像の一例を模式的に示す図である。
【図１０】時期ｔ２に対象領域を撮影したオルソ画像の一例を模式的に示す図である。
【図１１】本発明の一実施形態に係る家屋異動推定装置の機能の一例を示す機能ブロック図である。
【図１２】家屋異動の推定結果を表す画面の一例を示す図である。
【図１３Ａ】本発明の一実施形態に係る家屋異動推定装置で行われる処理の流れの一例を示すフロー図である。
【図１３Ｂ】本発明の一実施形態に係る家屋異動推定装置で行われる処理の流れの一例を示すフロー図である。
【図１３Ｃ】本発明の一実施形態に係る家屋異動推定装置で行われる処理の流れの一例を示すフロー図である。
【発明を実施するための形態】
【００１７】
以下、本発明の一実施形態について図面に基づき詳細に説明する。
【００１８】
図１は、本発明の一実施形態に係る家屋異動推定装置１０の構成図である。本実施形態に係る家屋異動推定装置１０は、パーソナルコンピュータなどのコンピュータである。図１に示すように家屋異動推定装置１０は、例えば、プロセッサ１２、記憶部１４、表示部１６、操作部１８を含んでいる。
【００１９】
プロセッサ１２は、例えば家屋異動推定装置１０にインストールされるプログラムに従って動作するＣＰＵ等のプログラム制御デバイスである。
【００２０】
記憶部１４は、ＲＯＭやＲＡＭ等の記憶素子やハードディスクドライブなどである。記憶部１４には、プロセッサ１２によって実行されるプログラムなどが記憶される。
【００２１】
表示部１６は、液晶ディスプレイ等の表示デバイスであって、プロセッサ１２の指示に従って各種の画像を表示する。
【００２２】
操作部１８は、キーボードやマウスなどといったユーザインタフェースであって、ユーザの操作入力を受け付けて、その内容を示す信号をプロセッサ１２に出力する。
【００２３】
なお、家屋異動推定装置１０は、ネットワークボードなどの通信インタフェース、ＤＶＤ−ＲＯＭやＢｌｕ−ｒａｙ（登録商標）ディスクなどの光ディスクを読み取る光ディスクドライブ、ＵＳＢ（Universal Serial Bus）ポートなどを含んでいてもよい。
【００２４】
本実施形態に係る家屋異動推定装置１０は、以下で説明するようにして、所定の対象領域についての時期ｔ１と時期ｔ２との間における家屋異動を推定する。以下の説明では、時期ｔ２は時期ｔ１よりも後の時期であることとする。
【００２５】
本実施形態では例えば、まず、図２Ａ及び図２Ｂに示すようにして、時期ｔ１と時期ｔ２とのＤＳＭ（Digital Surface Model）データが表すＤＳＭ２０に含まれるＤＴＭ（Digital Terrain Model：数値地形モデル）２２の影響を除去する正規化処理が実行される。
【００２６】
図２Ａには、地上の模式的な垂直断面図の一例が示されている。ＤＳＭ２０は、家屋２４、樹木２６等の地物表面の水平座標・標高と、地物の間に露出する地表面２８の水平座標・標高とを表す。ここで、ＤＳＭ２０における家屋の標高は、その下の地表面の標高と家屋の高さとの合計となる。
【００２７】
正規化処理では、ＤＳＭ２０の値からＤＴＭ２２の値を減算することで、図２Ｂに示すＮＤＳＭ（Normalized Digital Surface Model）３０が生成される。図２Ｂには、図２Ａの地上に対応して得られるＮＤＳＭ３０が模式的に示されている。ＮＤＳＭ３０では、地物（家屋３２、樹木３４）は地表からの高さ情報を有する一方、地表３６は基本的に、高さ０に応じた高さとなる。例えばＤＴＭ２２として、ＤＳＭ２０からモフォロジカルフィルタリング処理等の種々のフィルタ処理により推定されたものを用いることができる。あるいは対象領域に関する既存のＤＴＭをＤＴＭ２２として利用してもよい。
【００２８】
ＮＤＳＭ３０では、家屋３２以外の地物も有意な高さを有し、特に、樹木３４等の植生は家屋３２に相当する高さとなり得るため、家屋３２の判定に際してノイズ要因となる。そこで本実施形態では例えば、ＮＤＳＭ３０から植生に関わる部分を除去する植生除去処理が実行される。植生除去処理では、植生の領域が判別され、当該領域のＮＤＳＭデータに対し０への書き換えなどとったマスキングが実行される。
【００２９】
家屋の屋根や屋上は比較的大きな水平距離にわたって平滑な面が広がるのに対し、植生領域では、比較的小さい空間周期で高さの変動が生じる。例えば、この違いを利用して植生領域を判別し、ＮＤＳＭ３０から当該植生領域を除去することができる。具体的には、ＮＤＳＭデータの隣接点間での勾配の分散が所定の閾値を超える領域を植生領域と判断することができる。
【００３０】
また、高さの変動により航空機等から撮影した画像に生じるテクスチャに基づいて植生領域を判別することも可能であり、ＡＩ技術の画像解析から植生領域を判別することも可能である。
【００３１】
さらに、マルチスペクトル画像を使用し、樹木の葉が近赤外光を強く反射する性質から植生領域を検知することもできる。例えば、赤（Ｒ）成分の画素値Ｄ（Ｒ）及び近赤外（ＩＲ）成分の画素値Ｄ（ＩＲ）を用いて算出されるＮＤＶＩ（Normalized Difference Vegetation Index：正規化植生指標）が所定の閾値以上の領域を植生領域と判定することができる。ＮＤＶＩは例えば、ＮＤＶＩ＝（Ｄ（ＩＲ）−Ｄ（Ｒ））/（Ｄ（ＩＲ）＋Ｄ（Ｒ））との数式で算出される。
【００３２】
そして図３〜図６に示すように、植生除去処理が実行されたＮＤＳＭで表される対象領域を、当該ＮＤＳＭに現れる凸部４０をそれぞれ１つずつ含む凸部領域４２に分割する領域分割処理が実行される。上述のＤＳＭからＮＤＳＭへの変換処理によって、基本的に家屋の広がりより大きなスケールで現れる地表の凹凸は除去され、家屋に応じた大きさの凸部が残る。よって、家屋毎に凸部領域４２が生成され得る。
【００３３】
図３は、時期ｔ１のＤＳＭデータに基づいて生成される植生除去処理後の対象領域におけるＮＤＳＭの凹凸の一例を、便宜的に二値化して表した図である。図４は、時期ｔ２のＤＳＭデータに基づいて生成される植生除去処理後の対象領域におけるＮＤＳＭの凹凸を、便宜的に二値化して表した図である。図３及び図４では、凸部４０の位置が斜線によって示されている。
【００３４】
そして領域分割処理では、図５及び図６に示すように、複数の凸部４０がそれぞれ互いに分離されるように対象領域が複数の凸部領域４２に分割される。図５には、図３に示すＮＤＳＭに対する領域分割処理の結果である複数の凸部領域４２の一例が示されている。図６には、図４に示すＮＤＳＭに対する領域分割処理の結果である複数の凸部領域４２が示されている。この領域分割の手法として、例えば、ウォーターシェッド（Watershed）法を用いることができる。
【００３５】
そして以上のようにして、正規化処理、植生除去処理、及び、領域分割処理により抽出された凸部領域４２に対して、家屋の輪郭を抽出する輪郭抽出処理が実行される。輪郭抽出処理では、例えば、航空機等から取得された地上表層の標高の変化（例えば、勾配データ）に基づいて、エッジを構成する点群が抽出される。
【００３６】
また、マルチスペクトル画像を併用して、航空機等から撮影された地上表層の画像を用いた輪郭抽出を行うこともできる。
【００３７】
具体的には例えば、輝度や色で表示した画像の変化に基づいて、エッジを構成する点群が抽出される。エッジの抽出は各種エッジフィルタを用いて行うことができる。
【００３８】
そしてエッジを構成する点群に対して、線分を抽出し線分画像を生成する処理が行われる。例えば、トラッキング及び分割処理を行うフィルタリングを実行することで、当該点群をそれぞれひとまとまりとして捉えることができる単位の点の集合に分割し、また十分な大きさのまとまりを構成せずノイズと把握し得る点を除去する。そして例えばハフ変換により、線分が抽出される。そして線分画像から家屋の輪郭を抽出する処理が実行される。例えば、線分画像にPerceptual Grouping技術を適用することにより家屋の輪郭が抽出される。
【００３９】
そして本実施形態では例えば、抽出された輪郭により囲まれるポリゴンが家屋オブジェクト４４として特定される。本実施形態では、対象領域において家屋オブジェクト４４が占める領域が、当該家屋オブジェクト４４に対応付けられる家屋が表された領域であることとなる。以下、家屋オブジェクト４４が占める領域を家屋領域と呼ぶこととする。
【００４０】
図７には、図５に示す凸部領域４２に基づいて特定される、時期ｔ１における対象領域内の家屋オブジェクト４４の一例が示されている。図８には、図６に示す凸部領域４２に基づいて特定される、時期ｔ２における対象領域内の家屋オブジェクト４４の一例が示されている。
【００４１】
そして本実施形態では、時期ｔ１における対象領域内の家屋オブジェクト４４のそれぞれについて、時期ｔ１と時期ｔ２との間における、当該家屋オブジェクト４４に対応付けられる家屋の家屋異動が推定される。また本実施形態では、時期ｔ２における対象領域内の家屋オブジェクト４４のそれぞれについて、時期ｔ１と時期ｔ２との間における、当該家屋オブジェクト４４に対応付けられる家屋の家屋異動が推定される。
【００４２】
以下、時期ｔ１における対象領域内の１つの家屋オブジェクト４４に着目して、当該家屋オブジェクト４４に対応付けられる家屋の家屋異動の推定の一例について説明する。
【００４３】
まず、当該１つの家屋オブジェクト４４に対応付けられる家屋領域の時期ｔ１における標高が特定される。ここで例えば、時期ｔ１におけるＤＳＭが示す当該家屋領域内における複数の位置のそれぞれにおける標高の平均値が、時期ｔ１における当該家屋領域の標高として特定されてもよい。
【００４４】
そして、時期ｔ２における当該家屋領域の標高が特定される。ここで例えば、時期ｔ２におけるＤＳＭが示す当該家屋領域内における複数の位置のそれぞれにおける標高の平均値が、時期ｔ２における当該家屋領域の標高として特定されてもよい。
【００４５】
そして例えば時期ｔ１における当該家屋領域の標高と時期ｔ２における当該家屋領域の標高とが比較される。
【００４６】
ここで例えば時期ｔ１における当該家屋領域内の各画素の標高情報と時期ｔ２における当該家屋領域内の各画素の標高情報とが比較されてもよい。標高情報の比較結果の一例としては、例えば時期ｔ１における当該家屋領域内の各画素の標高と時期ｔ２における当該家屋領域内の各画素の標高との差を示す標高差値などが挙げられる。
【００４７】
また本実施形態では例えば、当該家屋領域について、時期ｔ１に撮影された画像に表された色情報が特定される。図９は、時期ｔ１に対象領域を撮影したオルソ画像の一例を模式的に示す図である。例えば図９に例示するオルソ画像において、当該家屋領域内の各画素の色情報が特定される。図９に示すオルソ画像には、複数の家屋４６が表されている。
【００４８】
また本実施形態では例えば、当該家屋領域について、時期ｔ２に撮影された画像に表された色情報が特定される。図１０は、時期ｔ２に対象領域を撮影したオルソ画像の一例を模式的に示す図である。例えば図１０に例示するオルソ画像において、当該家屋領域内の各画素の色情報が特定される。図１０に示すオルソ画像にも、複数の家屋４６が表されている。
【００４９】
そして本実施形態では例えば、時期ｔ１における当該家屋領域内の各画素の色情報と時期ｔ２における当該家屋領域内の各画素の色情報とが比較される。色情報の比較結果の一例としては、例えば時期ｔ１における当該家屋領域内の各画素の色と時期ｔ２における当該家屋領域内の各画素の色との差を示す色差値などが挙げられる。
【００５０】
ここで例えば、当該家屋領域内のすべての画素についての、時期ｔ１における当該画素と時期ｔ２における当該画素との、Ｒ値の差の二乗とＧ値の差の二乗とＢ値の差の二乗との和の平方根の平均値である平均変化ベクトルが色差値として特定されてもよい。ここで、Ｒ値、Ｇ値、Ｂ値は、それぞれ、赤成分の画素値、緑成分の画素値、青成分の画素値を指すこととする。
【００５１】
また例えば、時期ｔ１における当該家屋領域内の画素のＲ値、Ｇ値、Ｂ値のそれぞれの平均値と、時期ｔ２における当該家屋領域内の画素のＲ値、Ｇ値、Ｂ値のそれぞれの平均値と、が特定されてもよい。そしてＲ値の平均値の差の二乗とＧ値の平均値の差の二乗とＢ値の平均値の差の二乗との和の平方根が色差値として特定されてもよい。
【００５２】
また例えばＲＧＢ値による画素値の表現がＨＳＩ値による画素値の表現（ＨＳＩ表色系による画素値の表現）に変換されてもよい。そして当該家屋領域内のすべての画素について算出される色相成分（Ｈ成分）の差の平均値や色相成分（Ｈ成分）の差の絶対値の平均値が色差値として特定されてもよい。
【００５３】
また例えば家屋領域に対応するオルソ画像内の領域についてエッジ抽出が実行されてもよい。そして家屋領域に含まれる全画素数に対するエッジとして特定される画素数の割合の差が色差値として特定されてもよい。
【００５４】
そして例えば、当該家屋オブジェクト４４に対応する家屋領域についての、時期ｔ２におけるＮＤＳＭ、標高の比較結果、及び、色情報の比較結果のうちの少なくとも１つに基づいて、当該家屋オブジェクト４４に対応付けられる家屋の家屋異動が推定される。
【００５５】
ここで例えば時期ｔ２における標高を示す値が、所定の閾値よりも小さい場合（ほぼ０である場合）は、当該家屋について滅失が発生したと推定されてもよい。
【００５６】
また例えば時期ｔ１における標高を示す値と時期ｔ２における標高を示す値との差が、所定の閾値以上（例えば２メートル以上）である場合は、当該家屋について建替が発生したと推定されてもよい。
【００５７】
また例えば、時期ｔ１における標高を示す値と時期ｔ２における標高を示す値との差が所定の閾値未満（例えば２メートル未満）であり、色情報の差を示す値が所定値以上である場合に、当該家屋について建替が発生したと推定されてもよい。
【００５８】
また本実施形態において、当該家屋オブジェクト４４に対応する家屋領域についての時期ｔ２におけるＮＤＳＭ、標高の比較結果、及び、色情報の比較結果に基づいて、当該家屋オブジェクト４４に対応付けられる家屋の面積の変化が推定されてもよい。
【００５９】
そして時期ｔ１における面積を示す値よりも時期ｔ２における面積を示す値が、所定の閾値以上大きい（例えば１０平方メートル以上大きい）場合は、当該家屋について大きな増改築が発生したと推定されてもよい。
【００６０】
そして時期ｔ１における面積を示す値よりも時期ｔ２における面積を示す値の方が大きいが、その差が所定の閾値未満（例えば１０平方メートル未満）である場合は、当該家屋について小さな増改築が発生したと推定されてもよい。
【００６１】
また時期ｔ１における面積を示す値よりも時期ｔ２における面積を示す値が小さい場合は、当該家屋について減築が発生したと推定されてもよい。
【００６２】
また本実施形態において、当該家屋オブジェクト４４に対応する家屋領域についての、時期ｔ２におけるＮＤＳＭ、標高の比較結果、及び、色情報の比較結果に基づいて、当該家屋領域に配置されている家屋の数の変化が推定されてもよい。そして例えば時期ｔ１の家屋オブジェクト４４に対応する家屋領域について、時期ｔ２において複数の家屋オブジェクト４４が配置されていることが推定された場合は、当該家屋オブジェクト４４に対応付けられる家屋について建替が発生したと推定されてもよい。
【００６３】
本実施形態では、時期ｔ１における対象領域内のすべての家屋オブジェクト４４について、以上のようにして、時期ｔ１と時期ｔ２との間における、当該家屋オブジェクト４４に対応付けられる家屋の家屋異動が推定される。
【００６４】
また、時期ｔ２における対象領域内のすべての家屋オブジェクト４４についても同様に、時期ｔ１と時期ｔ２との間における、当該家屋オブジェクト４４に対応付けられる家屋の家屋異動が推定される。
【００６５】
この場合、例えば時期ｔ１における標高を示す値が、所定の閾値よりも小さい場合（ほぼ０である場合）は、当該家屋について新築が発生したと推定されてもよい。
【００６６】
また例えば、時期ｔ２における家屋オブジェクト４４に対応する家屋領域に、時期ｔ１において複数の家屋オブジェクト４４が配置されていることが推定された場合は、当該家屋について建替が発生したと推定されてもよい。
【００６７】
また本実施形態において、オルソ画像の色調整が実行されてもよい。例えばオルソ画像のＲ値、Ｇ値、Ｂ値のそれぞれについてのヒストグラムが生成されてもよい。そして例えば、時期ｔ１におけるオルソ画像のヒストグラムが表す平均値や標準偏差を、時期ｔ２におけるオルソ画像のヒストグラムが表す平均値や標準偏差に一致させるあるいは近づける色調整が時期ｔ１におけるオルソ画像に対して実行されてもよい。また例えば、時期ｔ２におけるオルソ画像のヒストグラムが表す平均値や標準偏差を、時期ｔ１におけるオルソ画像のヒストグラムが表す平均値や標準偏差に一致させるあるいは近づける色調整が時期ｔ２におけるオルソ画像に対して実行されてもよい。また時期ｔ１におけるオルソ画像と時期ｔ２におけるオルソ画像の両方に対して色調整が実行されてもよい。
【００６８】
そして色調整が実行されたオルソ画像に表された色情報の比較結果に基づいて、家屋オブジェクト４４に対応付けられる家屋の家屋異動が推定されてもよい。
【００６９】
なお以上で説明した閾値は、例えば、既存の家屋異動判読結果が大量にあれば大津の自動判別手法を利用して決定されるようにしてもよい。あるいは以上で説明した閾値は、サンプルを用いた目視判読により人手で設定されるようにしてもよい。
【００７０】
以上のようにして本実施形態によれば、多くの種類の家屋異動を精度よく推定できることとなる。
【００７１】
以下、本実施形態に係る家屋異動推定装置１０の機能並びに本実施形態に係る家屋異動推定装置１０で実行される処理についてさらに説明する。
【００７２】
図１１は、本実施形態に係る家屋異動推定装置１０で実装される機能の一例を示す機能ブロック図である。なお、本実施形態に係る家屋異動推定装置１０で、図１１に示す機能のすべてが実装される必要はなく、また、図１１に示す機能以外の機能が実装されていても構わない。
【００７３】
図１１に示すように、本実施形態に係る家屋異動推定装置１０には、機能的には例えば、ＤＳＭデータ記憶部５０、画像記憶部５２、家屋オブジェクト抽出部５４、特定部５６、家屋異動推定部５８、が含まれる。そして特定部５６には、標高特定部５６ａ、色情報特定部５６ｂ、が含まれる。ＤＳＭデータ記憶部５０、画像記憶部５２は、記憶部１４を主として実装される。家屋オブジェクト抽出部５４、特定部５６、家屋異動推定部５８は、プロセッサ１２を主として実装される。
【００７４】
以上の機能は、コンピュータである家屋異動推定装置１０にインストールされた、以上の機能に対応する指令を含むプログラムをプロセッサ１２で実行することにより実装されてもよい。このプログラムは、例えば、光ディスク、磁気ディスク、磁気テープ、光磁気ディスク、フラッシュメモリ等のコンピュータ読み取り可能な情報記憶媒体を介して、あるいは、インターネットなどを介して家屋異動推定装置１０に供給されてもよい。
【００７５】
ＤＳＭデータ記憶部５０は、本実施形態では例えば、時期ｔ１及び時期ｔ２を含む複数の時期における対象領域のＤＳＭデータを記憶する。
【００７６】
画像記憶部５２は、本実施形態では例えば、時期ｔ１及び時期ｔ２を含む複数の時期における対象領域の画像を記憶する。画像記憶部５２が記憶する画像は、例えば、航空写真画像、高分解能衛星画像などといった、対象領域の撮影画像であってもよい。また画像記憶部５２が記憶する画像は、例えば、オルソ画像であってもよい。
【００７７】
家屋オブジェクト抽出部５４は、本実施形態では例えば、家屋異動の推定の対象となる家屋異動推定領域についての地物を含んだ地上表層の標高のデータに基づいて、当該家屋異動推定領域内に家屋が表された家屋領域を抽出する。上述の例では、家屋オブジェクト抽出部５４は例えば時期ｔ１と時期ｔ２における対象領域内の家屋オブジェクト４４に対応付けられる家屋領域を抽出する。
【００７８】
特定部５６は、本実施形態では例えば、第１の時期と第２の時期との間における、家屋領域に表された家屋についての家屋異動の推定に用いられる、標高、標高の比較結果、又は、色情報の比較結果の少なくとも１つを特定する。
【００７９】
標高特定部５６ａは、本実施形態では例えば、第１の時期の標高のデータに基づいて抽出される家屋領域について、第１の時期における標高のデータが示す標高と第２の時期における標高のデータが示す標高とを特定する。標高特定部５６ａは、第１の時期における標高のデータが示す標高と第２の時期における標高のデータが示す標高との比較結果を特定してもよい。上述の例において、時期ｔ１が第１の時期に対応する場合は、時期ｔ２が第２の時期に対応する。また、時期ｔ２が第１の時期に対応する場合は、時期ｔ１が第２の時期に対応する。
【００８０】
色情報特定部５６ｂは、本実施形態では例えば、第１の時期の標高のデータに基づいて抽出される家屋領域について、第１の時期に撮影された画像に表された色情報と第２の時期に撮影された画像に表された色情報とを特定する。ここでも上述のように、上述の例において、時期ｔ１が第１の時期に対応する場合は、時期ｔ２が第２の時期に対応する。また、時期ｔ２が第１の時期に対応する場合は、時期ｔ１が第２の時期に対応する。
【００８１】
色情報特定部５６ｂは、第１の時期に撮影された画像又は第２の時期に撮影された画像の少なくとも一方に対して色調整を実行してもよい。この場合、色情報特定部５６ｂは、色調整が実行された画像に表された色情報を特定してもよい。
【００８２】
家屋異動推定部５８は、本実施形態では例えば、第１の時期と第２の時期との間における、家屋領域に表された家屋についての家屋異動を推定する。
【００８３】
ここで家屋異動推定部５８は、上述のように、第１の時期の標高のデータに基づいて抽出される家屋領域について、第２の時期における標高のデータが示す標高に基づいて当該家屋領域に対応付けられる家屋の家屋異動を推定してもよい。また家屋異動推定部５８は、上述のように、第１の時期の標高のデータに基づいて抽出される家屋領域に対応付けられる家屋の家屋異動を第１の時期と第２の時期との間における標高のデータが示す標高の比較結果に基づいて推定してもよい。
【００８４】
また家屋異動推定部５８は、上述のように、第１の時期の標高のデータに基づいて抽出される家屋領域に対応付けられる家屋の家屋異動を第１の時期と第２の時期に撮影された画像に表された色情報の比較結果に基づいて推定してもよい。ここで上述のように上述の色差値が色情報の比較結果として用いられてもよい。
【００８５】
また標高のデータが示す標高、標高のデータが示す標高の比較結果、及び、撮影された画像に表された色情報の比較結果のうちの複数の組み合わせに基づいて、家屋異動が推定されてもよい。
【００８６】
また家屋異動推定部５８は、第１の時期と第２の時期との間における標高のデータが示す標高の比較結果に基づいて、家屋領域に表された家屋についての建替の発生の有無を推定してもよい。
【００８７】
また家屋異動推定部５８は、第１の時期と第２の時期との間における標高のデータが示す標高の比較結果、及び、色情報の比較結果に基づいて、家屋領域に表された家屋についての建替の発生の有無を推定してもよい。
【００８８】
また家屋異動推定部５８は、ＨＳＩ表色系により表現された色情報における色相成分の比較結果に基づいて、家屋領域に表された家屋についての家屋異動を推定してもよい。
【００８９】
また家屋異動推定部５８は、色調整が実行された画像に表された色情報の比較結果に基づいて、家屋領域に表された家屋についての家屋異動を推定してもよい。
【００９０】
また家屋異動推定部５８は、家屋領域に表された家屋についての家屋異動の推定結果を出力してもよい。図１２に、家屋異動の推定結果を表す画面の一例を示す。図１２の例では、図１０に例示するオルソ画像に時期ｔ１と時期ｔ２との間における対象領域内の家屋についての家屋異動の推定結果を表す画像が重畳配置されている。
【００９１】
以下、本実施形態に係る家屋異動推定装置１０において行われる処理の流れの一例を、図１３Ａ〜図１３Ｃに例示するフロー図を参照しながら説明する。
【００９２】
まず、家屋オブジェクト抽出部５４が、ＤＳＭデータ記憶部５０から、時期ｔ１における対象領域のＤＳＭ及び時期ｔ２における対象領域のＤＳＭを示すＤＳＭデータを取得する（Ｓ１０１）。
【００９３】
そして、色情報特定部５６ｂが、画像記憶部５２から、時期ｔ１における対象領域の画像及び時期ｔ２における対象領域の画像を取得する（Ｓ１０２）。ここで色情報特定部５６ｂは、上述のように、時期ｔ１における対象領域の画像及び時期ｔ２における対象領域の画像の少なくとも一方に対して色調整を実行してもよい。
【００９４】
そして、家屋オブジェクト抽出部５４が、時期ｔ１における対象領域のＤＳＭと時期ｔ２における対象領域のＤＳＭのそれぞれに対して上述の正規化処理を実行する（Ｓ１０３）。Ｓ１０３に示す処理によって、時期ｔ１における対象領域のＮＤＳＭと時期ｔ２における対象領域のＮＤＳＭとが生成される。
【００９５】
そして、家屋オブジェクト抽出部５４が、Ｓ１０３に示す処理で生成されたＮＤＳＭに対して植生除去処理を実行する（Ｓ１０４）。
【００９６】
そして、家屋オブジェクト抽出部５４が、Ｓ１０４に示す処理が実行されたＮＤＳＭに対して領域分割処理を実行する（Ｓ１０５）。Ｓ１０５に示す処理によって、時期ｔ１における対象領域、及び、時期ｔ２における対象領域は複数の凸部領域４２に分割される。
【００９７】
そして、家屋オブジェクト抽出部５４が、時期ｔ１における複数の凸部領域４２、及び、時期ｔ２における複数の凸部領域４２のそれぞれに対して輪郭抽出処理を実行する（Ｓ１０６）。Ｓ１０６に示す処理によって、時期ｔ１における対象領域内の家屋オブジェクト４４、及び、時期ｔ２における対象領域内の家屋オブジェクト４４が特定される。以下、時期ｔ１における対象領域について特定された家屋オブジェクト４４を第１家屋オブジェクトと呼び、時期ｔ２における対象領域について特定された家屋オブジェクト４４を第２家屋オブジェクトと呼ぶこととする。
【００９８】
そして、家屋異動推定部５８が、Ｓ１０６に示す処理で特定された第１家屋オブジェクトのそれぞれに１から始まる順序番号を設定し、Ｓ１０６に示す処理で特定された第２家屋オブジェクトのそれぞれに１から始まる順序番号を設定する（Ｓ１０７）。以下、特定された第１家屋オブジェクトの数がＭであり、特定された第２家屋オブジェクトの数がＮであることとする。そして特定された第１家屋オブジェクトのそれぞれには１以上Ｍ以下の順序番号が設定され、特定された第２家屋オブジェクトのそれぞれには１以上Ｎ以下の順序番号が設定されることとする。
【００９９】
そして、家屋異動推定部５８が、変数ｉの値として１を設定する（Ｓ１０８）。
【０１００】
そして、標高特定部５６ａが、時期ｔ１のＤＳＭが示す、順序番号がｉである第１家屋オブジェクトに対応する家屋領域に含まれる複数の位置についての時期ｔ１の標高に基づいて、当該家屋領域の時期ｔ１における標高を特定する（Ｓ１０９）。ここで例えば、当該家屋領域に含まれる複数の位置についての時期ｔ１の標高の平均値が当該家屋領域の時期ｔ１における標高として特定されてもよい。
【０１０１】
そして、標高特定部５６ａが、時期ｔ２のＤＳＭが示す、順序番号がｉである第１家屋オブジェクトに対応する家屋領域に含まれる複数の位置についての時期ｔ２の標高に基づいて、当該家屋領域の時期ｔ２における標高を特定する（Ｓ１１０）。ここで例えば、当該家屋領域に含まれる複数の位置についての時期ｔ２の標高の平均値が当該家屋領域の時期ｔ２における標高として特定されてもよい。
【０１０２】
そして、色情報特定部５６ｂが、Ｓ１０２に示す処理で取得された時期ｔ１における対象領域の画像における、順序番号がｉである第１家屋オブジェクトに対応する家屋領域に含まれる各画素の、時期ｔ１における色情報を特定する（Ｓ１１１）。ここで色調整が実行された画像における色情報が特定されてもよい。
【０１０３】
そして、色情報特定部５６ｂが、Ｓ１０２に示す処理で取得された時期ｔ２における対象領域の画像における、順序番号がｉである第１家屋オブジェクトに対応する家屋領域に含まれる各画素の、時期ｔ２における色情報を特定する（Ｓ１１２）。ここで色調整が実行された画像における色情報が特定されてもよい。
【０１０４】
そして、家屋異動推定部５８が、Ｓ１０９〜Ｓ１１２に示す処理で特定された標高及び色情報に基づいて、順序番号がｉである第１家屋オブジェクトに対応する家屋についての時期ｔ１と時期ｔ２との間における家屋異動を推定する（Ｓ１１３）。ここでＳ１０９〜Ｓ１１２に示す処理で特定された標高及び色情報の一部に基づいて家屋異動が推定されてもよい。
【０１０５】
そして、家屋異動推定部５８が、変数ｉの値がＭであるか否かを確認する（Ｓ１１４）。
【０１０６】
変数ｉの値がＭでないことが確認された場合は（Ｓ１１４：Ｎ）、家屋異動推定部５８が、変数ｉの値を１だけ増加させて（Ｓ１１５）、Ｓ１０９に示す処理に戻る。
【０１０７】
変数ｉの値がＭであることが確認された場合は（Ｓ１１４：Ｙ）、家屋異動推定部５８が、変数ｉの値を１に変更する（Ｓ１１６）。
【０１０８】
そして、標高特定部５６ａが、時期ｔ２のＤＳＭが示す、順序番号がｉである第２家屋オブジェクトに対応する家屋領域に含まれる複数の位置についての時期ｔ２の標高に基づいて、当該家屋領域の時期ｔ２における標高を特定する（Ｓ１１７）。ここで例えば、当該家屋領域に含まれる複数の位置についての時期ｔ２の標高の平均値が当該家屋領域の時期ｔ２における標高として特定されてもよい。
【０１０９】
そして、標高特定部５６ａが、時期ｔ１のＤＳＭが示す、順序番号がｉである第２家屋オブジェクトに対応する家屋領域に含まれる複数の位置についての時期ｔ１の標高に基づいて、当該家屋領域の時期ｔ１における標高を特定する（Ｓ１１８）。ここで例えば、当該家屋領域に含まれる複数の位置についての時期ｔ１の標高の平均値が当該家屋領域の時期ｔ１における標高として特定されてもよい。
【０１１０】
そして、色情報特定部５６ｂが、Ｓ１０２に示す処理で取得された時期ｔ２における対象領域の画像における、順序番号がｉである第２家屋オブジェクトに対応する家屋領域に含まれる各画素の、時期ｔ２における色情報を特定する（Ｓ１１９）。ここで色調整が実行された画像における色情報が特定されてもよい。
【０１１１】
そして、色情報特定部５６ｂが、Ｓ１０２に示す処理で取得された時期ｔ１における対象領域の画像における、順序番号がｉである第２家屋オブジェクトに対応する家屋領域に含まれる各画素の、時期ｔ１における色情報を特定する（Ｓ１２０）。ここで色調整が実行された画像における色情報が特定されてもよい。
【０１１２】
そして、家屋異動推定部５８が、Ｓ１１７〜Ｓ１２０に示す処理で特定された標高及び色情報に基づいて、順序番号がｉである第２家屋オブジェクトに対応する家屋についての時期ｔ１と時期ｔ２との間における家屋異動を推定する（Ｓ１２１）。ここでＳ１１７〜Ｓ１２０に示す処理で特定された標高及び色情報の一部に基づいて家屋異動が推定されてもよい。
【０１１３】
そして、家屋異動推定部５８が、変数ｉの値がＮであるか否かを確認する（Ｓ１２２）。
【０１１４】
変数ｉの値がＮでないことが確認された場合は（Ｓ１２２：Ｎ）、家屋異動推定部５８が、変数ｉの値を１だけ増加させて（Ｓ１２３）、Ｓ１１７に示す処理に戻る。
【０１１５】
変数ｉの値がＮであることが確認された場合は（Ｓ１２２：Ｙ）、本処理例に示す処理は終了される。
【０１１６】
なお、本発明は上述の実施形態に限定されるものではない。
【０１１７】
また、上述の具体的な文字列や数値、並びに、図面中の具体的な文字列は例示であり、これらの文字列や数値には限定されない。
【符号の説明】
【０１１８】
１０  家屋異動推定装置、１２  プロセッサ、１４  記憶部、１６  表示部、１８  操作部、２０  ＤＳＭ、２２  ＤＴＭ、２４  家屋、２６  樹木、２８  地表面、３０  ＮＤＳＭ、３２  家屋、３４  樹木、３６  地表、４０  凸部、４２  凸部領域、４４  家屋オブジェクト、４６  家屋、５０  ＤＳＭデータ記憶部、５２  画像記憶部、５４  家屋オブジェクト抽出部、５６  特定部、５６ａ  標高特定部、５６ｂ  色情報特定部、５８  家屋異動推定部。
【図１】
2020016923.tif 000003
【図２Ａ】
2020016923.tif 000004
【図２Ｂ】
2020016923.tif 000005
【図３】
2020016923.tif 000006
【図４】
2020016923.tif 000007
【図５】
2020016923.tif 000008
【図６】
2020016923.tif 000009
【図７】
2020016923.tif 000010
【図８】
2020016923.tif 000011
【図９】
2020016923.tif 000012
【図１０】
2020016923.tif 000013
【図１１】
2020016923.tif 000014
【図１２】
2020016923.tif 000015
【図１３Ａ】
2020016923.tif 000016
【図１３Ｂ】
2020016923.tif 000017
【図１３Ｃ】
2020016923.tif 000018