(57)【要約】【課題】  人工衛星等から撮像した画像データを用いて、自然災害や人的災害による被害域を自動抽出する被害域自動抽出システムを提供することを目的とする。【解決手段】  被害前画像データと被害後画像データを読み込む画像データ読込手段と、被害後画像データを構成する画素の輝度により被害域候補画素を判定する被害域候補判定手段と、被害後画像データに於いて、画素の輝度が不連続な部分をエッジとし、エッジの連続性と方向性を判定することにより非被害域候補画素を判定する非被害域候補判定手段と、被害前画像データを構成する画素の輝度により第１誤認候補画素を判定する第１誤認候補判定手段と、被害後画像データと被害域候補画素と非被害域候補画素と第１誤認候補画素とにより、被害域を判定する被害域判定手段と、を有する被害域自動抽出システム１である。【選択図】  図１
2005284539.tif 000002
【特許請求の範囲】
【請求項１】
撮像された画像データから自然災害や人的災害の被害域を自動抽出する被害域自動抽出システムであって、
被害前画像データと被害後画像データを読み込む画像データ読込手段と、
前記被害後画像データの画素の輝度により被害域候補画素を判定する被害域候補判定手段と、
前記被害後画像データに於いて、前記画素の輝度が不連続な部分をエッジとし、前記エッジの連続性及び／又は方向性を判定することにより非被害域候補画素を判定する非被害域候補判定手段と、
前記被害前画像データの画素の輝度により第１誤認候補画素を判定する第１誤認候補判定手段と、
前記被害後画像データと前記被害域候補画素と前記非被害域候補画素と第１誤認候補画素とにより、被害域を判定する被害域判定手段と、
を有することを特徴とする被害域自動抽出システム。
【請求項２】
撮像された画像データから自然災害や人的災害の被害域を自動抽出する被害域自動抽出システムであって、
被害前画像データと被害後画像データを読み込む画像データ読込手段と、
前記被害後画像データの画素の輝度により被害域候補画素を判定する被害域候補判定手段と、
前記被害後画像データに於いて、前記画素の輝度が不連続な部分をエッジとし、前記エッジの連続性及び／又は方向性を判定することにより非被害域候補画素を判定する非被害域候補判定手段と、
前記被害前画像データに於いて、前記画素の輝度が不連続な部分をエッジとし、前記エッジの連続性及び／又は方向性を判定することにより第２誤認候補画素を判定する第２誤認候補判定手段と、
前記画像データと前記被害域候補画素と前記非被害域候補画素と第２誤認候補画素とにより、被害域を判定する被害域判定手段と、
を有することを特徴とする被害域自動抽出システム。
【請求項３】
撮像された画像データから自然災害や人的災害の被害域を自動抽出する被害域自動抽出システムであって、
被害前画像データと被害後画像データを読み込む画像データ読込手段と、
前記被害後画像データの画素の輝度により被害域候補画素を判定する被害域候補判定手段と、
前記被害後画像データに於いて、前記画素の輝度が不連続な部分をエッジとし、前記エッジの連続性及び／又は方向性を判定することにより非被害域候補画素を判定する非被害域候補判定手段と、
前記被害前画像データの画素の輝度により第１誤認候補画素を判定する第１誤認候補判定手段と、
前記被害前画像データに於いて、前記画素の輝度が不連続な部分をエッジとし、前記エッジの連続性及び／又は方向性を判定することにより第２誤認候補画素を判定する第２誤認候補判定手段と、
前記被害後画像データと前記被害域候補画素と前記非被害域候補画素と第１誤認候補画素と第２誤認候補画素とにより、被害域を判定する被害域判定手段と、
を有することを特徴とする被害域自動抽出システム。
【請求項４】
前記被害域候補判定手段は、
前記被害後画像データに於いて、前記画素を中心とする所定範囲の部分画像を設定し、前記各部分画像に於いて、前記部分画像内で、各輝度毎にその輝度を有する画素数の対応関係を作成し、前記対応関係から最頻度の画素と、前記輝度の最小値と最大値の幅を輝度幅として算出し、前記最頻度の画素数が所定の閾値以下であり、且つ前記輝度幅が所定の範囲内にある前記部分画像の一又は複数の画素を被害域候補画素として判定する、
ことを特徴とする請求項１から請求項３のいずれかに記載の被害域自動抽出システム。
【請求項５】
前記非被害域候補判定手段は、
前記被害後画像データに於いて、前記輝度が不連続な画素をエッジとして抽出し、前記エッジの強度が所定の閾値以上のエッジに対して細線化し、前記被害後画像データに於ける画素を中心とする所定範囲の部分画像を設定し、前記各部分画像に於いて、前記部分画像内で、前記細線化したエッジが所定の長さ以上連続しており、且つ一定の方向性を有している画素数が所定の閾値以上の、前記部分画像の一又は複数の画素を非被害域候補画素として判定する、
ことを特徴とする請求項１から請求項３のいずれかに記載の被害域自動抽出システム。
【請求項６】
前記第１誤認候補判定手段は、
前記被害前画像データの画素を中心とする所定範囲の部分画像を設定し、前記各部分画像に於いて、前記部分画像内で、各輝度毎にその輝度を有する画素数の対応関係を作成し、前記対応関係から最頻度の画素と、前記輝度の最小値と最大値の幅を輝度幅として算出し、前記最頻度の画素数が所定の閾値以下であり、且つ前記輝度幅が所定の範囲内にある前記部分画像の一又は複数の画素を第１誤認候補画素として判定する、
ことを特徴とする請求項１又は請求項３に記載の被害域自動抽出システム。
【請求項７】
前記第２誤認候補判定手段は、
前記被害前画像データに於いて、前記輝度が不連続な画素をエッジとして抽出し、前記エッジの強度が所定の閾値以上のエッジに対して細線化し、前記被害前画像データに於ける画素を中心とする所定範囲の部分画像を設定し、前記各部分画像に於いて、前記部分画像内で、前記細線化したエッジが所定の長さ以上連続しており、且つ一定の方向性を有している画素数が所定の閾値以上の、前記部分画像の一又は複数の画素を第２誤認候補画素として判定する、
ことを特徴とする請求項２又は請求項３に記載の被害域自動抽出システム。
【請求項８】
前記被害域判定手段は、
前記被害後画像データの画素を中心とする所定範囲の部分画像を設定し、前記各部分画像に於いて、前記部分画像内の被害域候補画素数が所定の閾値以上で、且つ非被害域候補画素数が所定の閾値以下であって、更に前記第１誤認候補画素が所定の閾値以下の、前記部分画像の一又は複数の画素の一定範囲を被害域として判定する、
ことを特徴とする請求項１又は請求項３に記載の被害域自動抽出システム。
【請求項９】
前記被害域判定手段は、
前記被害後画像データの画素を中心とする所定範囲の部分画像を設定し、前記各部分画像に於いて、前記部分画像内の被害域候補画素数が所定の閾値以上で、且つ非被害域候補画素数が所定の閾値以下であって、更に前記第２誤認候補画素が所定の閾値以上の、前記部分画像の一又は複数の画素の一定範囲を被害域として判定する、
ことを特徴とする請求項２又は請求項３に記載の被害域自動抽出システム。
【請求項１０】
前記被害域自動抽出システムは、
前記被害域判定手段で判定した被害域を前記読み込んだ画像データに重ね合わせ、表示装置に出力する、
ことを特徴とする請求項１から請求項９のいずれかに記載の被害域自動抽出システム。
【請求項１１】
データを記憶する記憶部とプログラムの処理演算を実行する処理演算部とを有しているコンピュータ端末で実行する被害域自動抽出プログラムであって、前記被害域自動抽出プログラムは、前記処理演算部に読み込まれる画像データ読込機能と被害域候補判定機能と非被害域候補判定機能と第１誤認候補判定機能と被害域判定機能とを有しており、
前記画像データ読込機能は、
被害前画像データと被害後画像データとを読み込み前記記憶部に記憶し、
前記被害域候補判定機能は、
前記記憶部に記憶した被害後画像データの画素の輝度により被害域候補画素を判定し、それを前記記憶部に記憶し、
前記非被害域候補判定機能は、
前記被害後画像データに於いて、前記画素の輝度が不連続な部分をエッジとし、前記エッジの連続性及び／又は方向性を判定することにより非被害域候補画素を判定し、それを前記記憶部に記憶し、
前記第１誤認候補機能は、
前記被害前画像データの画素の輝度により第１誤認候補画素を判定し、それを前記記憶部に記憶し、
前記被害域判定機能は、
前記記憶部に記憶した被害後画像データと前記被害域候補画素と前記非被害域候補画素と前記第１誤認候補画素とにより、被害域を判定する、
ことを特徴とする被害域自動抽出プログラム。
【請求項１２】
データを記憶する記憶部とプログラムの処理演算を実行する処理演算部とを有しているコンピュータ端末で実行する被害域自動抽出プログラムであって、前記被害域自動抽出プログラムは、前記処理演算部に読み込まれる画像データ読込機能と被害域候補判定機能と非被害域候補判定機能と第２誤認候補判定機能と被害域判定機能とを有しており、
前記画像データ読込機能は、
被害前画像データと被害後画像データとを読み込み前記記憶部に記憶し、
前記被害域候補判定機能は、
前記記憶部に記憶した被害後画像データの画素の輝度により被害域候補画素を判定し、それを前記記憶部に記憶し、
前記非被害域候補判定機能は、
前記被害後画像データに於いて、前記画素の輝度が不連続な部分をエッジとし、前記エッジの連続性及び／又は方向性を判定することにより非被害域候補画素を判定し、それを前記記憶部に記憶し、
前記第２誤認候補判定機能は、
前記被害前画像データに於いて、前記画素の輝度が不連続な部分をエッジとし、前記エッジの連続性及び／又は方向性を判定することにより第２誤認候補画素を判定し、それを前記記憶部に記憶し、
前記記憶部に記憶した被害後画像データと前記被害域候補画素と前記非被害域候補画素と前記第２誤認候補画素とにより、被害域を判定する、
ことを特徴とする被害域自動抽出プログラム。
【請求項１３】
データを記憶する記憶部とプログラムの処理演算を実行する処理演算部とを有しているコンピュータ端末で実行する被害域自動抽出プログラムであって、前記被害域自動抽出プログラムは、前記処理演算部に読み込まれる画像データ読込機能と被害域候補判定機能と非被害域候補判定機能と第１誤認候補判定機能と第２誤認候補判定機能と被害域判定機能とを有しており、
前記画像データ読込機能は、
被害前画像データと被害後画像データとを読み込み前記記憶部に記憶し、
前記被害域候補判定機能は、
前記記憶部に記憶した被害後画像データの画素の輝度により被害域候補画素を判定し、それを前記記憶部に記憶し、
前記非被害域候補判定機能は、
前記被害後画像データに於いて、前記画素の輝度が不連続な部分をエッジとし、前記エッジの連続性及び／又は方向性を判定することにより非被害域候補画素を判定し、それを前記記憶部に記憶し、
前記第１誤認候補機能は、
前記被害前画像データの画素の輝度により第１誤認候補画素を判定し、それを前記記憶部に記憶し、
前記第２誤認候補機能は、
前記被害前画像データに於いて、前記画素の輝度が不連続な部分をエッジとし、前記エッジの連続性及び／又は方向性を判定することにより第２誤認候補画素を判定し、それを前記記憶部に記憶し、
前記被害域判定機能は、
前記記憶部に記憶した被害後画像データと前記被害域候補画素と前記非被害域候補画素と前記第１誤認候補画素と前記第２誤認候補画素とにより、被害域を判定する、
ことを特徴とする被害域自動抽出プログラム。
【請求項１４】
データを記憶する記憶部とプログラムの処理演算を実行する処理演算部とを有しているコンピュータ端末で実行する被害域自動抽出プログラムを記録した記録媒体であって、前記被害域自動抽出プログラムは、前記処理演算部に読み込まれる画像データ読込機能と被害域候補判定機能と非被害域候補判定機能と第１誤認候補判定機能と被害域判定機能とを有しており、
前記画像データ読込機能は、
前記画像データを読み込み前記記憶部に記憶し、
前記被害域候補判定機能は、
前記記憶部に記憶した画像データの画素の輝度により被害域候補画素を判定し、それを前記記憶部に記憶し、
前記非被害域候補判定機能は、
前記画像データに於いて、前記画素の輝度が不連続な部分をエッジとし、前記エッジの連続性及び／又は方向性を判定することにより非被害域候補画素を判定し、それを前記記憶部に記憶し、
前記第１誤認候補機能は、
前記被害前画像データの画素の輝度により第１誤認候補画素を判定し、それを前記記憶部に記憶し、
前記被害域判定機能は、
前記記憶部に記憶した被害後画像データと前記被害域候補画素と前記非被害域候補画素と前記第１誤認候補画素とにより、被害域を判定する、
ことを特徴とする被害域自動抽出プログラムを記録した記録媒体。
【請求項１５】
データを記憶する記憶部とプログラムの処理演算を実行する処理演算部とを有しているコンピュータ端末で実行する被害域自動抽出プログラムを記録した記録媒体であって、前記被害域自動抽出プログラムは、前記処理演算部に読み込まれる画像データ読込機能と被害域候補判定機能と非被害域候補判定機能と第２誤認候補判定機能と被害域判定機能とを有しており、
前記画像データ読込機能は、
前記画像データを読み込み前記記憶部に記憶し、
前記被害域候補判定機能は、
前記記憶部に記憶した画像データの画素の輝度により被害域候補画素を判定し、それを前記記憶部に記憶し、
前記非被害域候補判定機能は、
前記画像データに於いて、前記画素の輝度が不連続な部分をエッジとし、前記エッジの連続性及び／又は方向性を判定することにより非被害域候補画素を判定し、それを前記記憶部に記憶し、
前記第２誤認候補判定機能は、
前記被害前画像データに於いて、前記画素の輝度が不連続な部分をエッジとし、前記エッジの連続性及び／又は方向性を判定することにより第２誤認候補画素を判定し、それを前記記憶部に記憶し、
前記被害域判定機能は、
前記記憶部に記憶した被害後画像データと前記被害域候補画素と前記非被害域候補画素と前記第２誤認候補画素とにより、被害域を判定する、
ことを特徴とする被害域自動抽出プログラムを記録した記録媒体。
【請求項１６】
データを記憶する記憶部とプログラムの処理演算を実行する処理演算部とを有しているコンピュータ端末で実行する被害域自動抽出プログラムを記録した記録媒体であって、前記被害域自動抽出プログラムは、前記処理演算部に読み込まれる画像データ読込機能と被害域候補判定機能と非被害域候補判定機能と第１誤認候補判定機能と第２誤認候補判定機能と被害域判定機能とを有しており、
前記画像データ読込機能は、
前記画像データを読み込み前記記憶部に記憶し、
前記被害域候補判定機能は、
前記記憶部に記憶した画像データの画素の輝度により被害域候補画素を判定し、それを前記記憶部に記憶し、
前記非被害域候補判定機能は、
前記画像データに於いて、前記画素の輝度が不連続な部分をエッジとし、前記エッジの連続性及び／又は方向性を判定することにより非被害域候補画素を判定し、それを前記記憶部に記憶し、
前記第１誤認候補機能は、
前記被害前画像データの画素の輝度により第１誤認候補画素を判定し、それを前記記憶部に記憶し、
前記第２誤認候補判定機能は、
前記被害前画像データに於いて、前記画素の輝度が不連続な部分をエッジとし、前記エッジの連続性及び／又は方向性を判定することにより第２誤認候補画素を判定し、それを前記記憶部に記憶し、
前記被害域判定機能は、
前記記憶部に記憶した被害後画像データと前記被害域候補画素と前記非被害域候補画素と前記第１誤認候補画素と前記第２誤認候補判定機能とにより、被害域を判定する、
ことを特徴とする被害域自動抽出プログラムを記録した記録媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、人工衛星等から撮像した画像データを用いて、自然災害や人的災害による被害領域（被害域）を自動抽出する被害域自動抽出システムに関する。
【背景技術】
【０００２】
人工衛星や航空機によるリモートセンシングは、広範囲での地上の状況把握が可能であり、火山災害、土砂災害、洪水、林野火災等の自然火災や重油流出等の人的災害の状況把握に利用されている。
【０００３】
このような災害に関しては被害状況をいち早く把握することが非常に重要であり、被災地の上空から人工衛星や航空機を用いて被害状況を撮像している。しかしこのような撮像が行われても、従来は、下記非特許文献１に示すように、目視に依存している面が大きく、被害状況の把握に多大な日時を要しており、迅速な救援活動に対する障害となっている。
【０００４】
その為、目視による被害状況の把握ではなく、撮像した画像データから被害域を自動的に抽出するシステムが下記特許文献１及び特許文献２に開示されている。
【０００５】
【特許文献１】特開２００３−１２９４８７号公報
【特許文献２】特開２０００−２２２５６５号公報
【非特許文献１】小川直樹、山崎文雄著  、”航空写真を用いた兵庫県南部地震における建物被害の目視判読”、地域安全学会論文集、Ｖｏｌ．２、ｐ．１１９−１２８、２０００年
【発明の開示】
【発明が解決しようとする課題】
【０００６】
上述した特許文献１では画像データから災害危険箇所を自動的に抽出する方法が開示されているが、この方法を用いる場合には、ディジタルマッピングデータや数値標高データ等を利用する為、これらのデータを予め用意しておく必要がある。そうすると、いつどこで発生するかの予測が困難な自然災害や人的災害の場合、災害の発生地域に応じたこれらのデータを全ての地域について予め用意しておく必要があるが、現実的ではなく、事実上それは困難である。従って、広範囲で且つ任意の場所での被害域を自動抽出することが難しい。
【０００７】
特許文献２には、２つの画像の変化を抽出することによって被害域を自動的に判定する方法及び装置が開示されている。しかし２つの画像の変化を被害域として判定しているので、被害とは関係なく変化している箇所も被害域として抽出する可能性があり、誤認率が高い。
【課題を解決するための手段】
【０００８】
そこで本発明者は、人工衛星や航空機等から撮像した被害後の画像データに、更に被害前の画像データから誤認する可能性のある画素を抽出することによって、自動的に被害域を抽出し、且つその精度も高い被害域自動抽出システムを発明した。
【０００９】
請求項１の発明は、撮像された画像データから自然災害や人的災害の被害域を自動抽出する被害域自動抽出システムであって、被害前画像データと被害後画像データを読み込む画像データ読込手段と、前記被害後画像データの画素の輝度により被害域候補画素を判定する被害域候補判定手段と、前記被害後画像データに於いて、前記画素の輝度が不連続な部分をエッジとし、前記エッジの連続性及び／又は方向性を判定することにより非被害域候補画素を判定する非被害域候補判定手段と、前記被害前画像データの画素の輝度により第１誤認候補画素を判定する第１誤認候補判定手段と、前記被害後画像データと前記被害域候補画素と前記非被害域候補画素と第１誤認候補画素とにより、被害域を判定する被害域判定手段と、を有する被害域自動抽出システムである。
【００１０】
請求項２の発明は、撮像された画像データから自然災害や人的災害の被害域を自動抽出する被害域自動抽出システムであって、被害前画像データと被害後画像データを読み込む画像データ読込手段と、前記被害後画像データの画素の輝度により被害域候補画素を判定する被害域候補判定手段と、前記被害後画像データに於いて、前記画素の輝度が不連続な部分をエッジとし、前記エッジの連続性及び／又は方向性を判定することにより非被害域候補画素を判定する非被害域候補判定手段と、前記被害前画像データに於いて、前記画素の輝度が不連続な部分をエッジとし、前記エッジの連続性及び／又は方向性を判定することにより第２誤認候補画素を判定する第２誤認候補判定手段と、前記画像データと前記被害域候補画素と前記非被害域候補画素と第２誤認候補画素とにより、被害域を判定する被害域判定手段と、を有する被害域自動抽出システムである。
【００１１】
請求項３の発明は、撮像された画像データから自然災害や人的災害の被害域を自動抽出する被害域自動抽出システムであって、被害前画像データと被害後画像データを読み込む画像データ読込手段と、前記被害後画像データの画素の輝度により被害域候補画素を判定する被害域候補判定手段と、前記被害後画像データに於いて、前記画素の輝度が不連続な部分をエッジとし、前記エッジの連続性及び／又は方向性を判定することにより非被害域候補画素を判定する非被害域候補判定手段と、前記被害前画像データの画素の輝度により第１誤認候補画素を判定する第１誤認候補判定手段と、前記被害前画像データに於いて、前記画素の輝度が不連続な部分をエッジとし、前記エッジの連続性及び／又は方向性を判定することにより第２誤認候補画素を判定する第２誤認候補判定手段と、前記被害後画像データと前記被害域候補画素と前記非被害域候補画素と第１誤認候補画素と第２誤認候補画素とにより、被害域を判定する被害域判定手段と、を有する被害域自動抽出システムである。
【００１２】
画像データに於ける輝度は、明るさ又は濃度（以下、これらを総称して明るさとする）に相当するものである。明るさが変化することはそこに被害が発生している可能性があることを示唆するものである。一方、画像データから抽出できたエッジは、建造物等の輪郭を識別するものである。そうするとそのエッジが連続していたり、一定方向に方向性を有していればその建造物等の輪郭ははっきりしている、即ち建造物等に被害は発生していないことを示唆するものである。このような観点から、輝度とエッジの２つの要素により被害域を判定することで、従来よりも精度を向上させた被害域の自動抽出システムが可能となる。
【００１３】
又、本発明では更に被害前の画像データに対しても同様の分析を行うことで、その誤認率を低減させている。即ち被害前の画像データに対して輝度による分析を行い、被害域候補画素（第１誤認候補画素）を判定している。しかし被害前の画像データでは、実際には被害域があるはずがないが、被害域として誤認する可能性がある。しかしその誤認した第１誤認候補画素の数は、被害後の画像データの被害域候補画素の数よりは少ないと考えられる。そうすると、逆に言えばこの関係が成立してなければそれは誤認として排除することが出来るので、これにより誤認率を低減させている。又、被害前の画像データに対してエッジによる分析を行い、非被害域候補画素（第２誤認候補画素）を判定している。被害前の画像データであっても非被害域として抽出できる部分があるが、被害前の画像データの非被害域候補画素（第２誤認候補画素）の方が、被害後画像データの被害域候補よりも多いはずである。そうすると、逆に言えばこの関係が成立してなければそれは誤認として排除することが出来るので、これにより誤認率を低減させている。
【００１４】
請求項４の発明は、前記被害域候補判定手段は、前記被害後画像データに於いて、前記画素を中心とする所定範囲の部分画像を設定し、前記各部分画像に於いて、前記部分画像内で、各輝度毎にその輝度を有する画素数の対応関係を作成し、前記対応関係から最頻度の画素と、前記輝度の最小値と最大値の幅を輝度幅として算出し、前記最頻度の画素数が所定の閾値以下であり、且つ前記輝度幅が所定の範囲内にある前記部分画像の一又は複数の画素を被害域候補画素として判定する、被害域自動抽出システムである。
【００１５】
輝度による被害域候補画素は、本発明のように判定すると良い。輝度の最頻度の画素数が所定の閾値以下であることは、部分画像に於いて特定の輝度が突出していないことを示しており、上述のように、これが被害域の候補であることの一つの条件となる。輝度幅が所定の範囲内にあることは、極端に輝度の変化がないことを示している。つまり、被害が発生していない建造物等であれば、部分画像内に於いて輝度幅に変化は少ない（同一の建造物等であるので、その屋根や壁等は同一の明るさ、それに類似する明るさであることが想定される）。逆に言えば、部分画像内に於いて輝度幅に変化が大きい場合、それは屋根や壁等が被害を受けたことによって、部分画像内に様々な明るさが含まれていると考えられる。そうすると、その部分画像の画素は、被害域の候補の可能性があることを示唆することとなる。そしてこの２つの条件を満たした場合に、その部分画像に於ける画素を被害域の候補画素として判定することが出来る。
【００１６】
請求項５の発明は、前記非被害域候補判定手段は、前記被害後画像データに於いて、前記輝度が不連続な画素をエッジとして抽出し、前記エッジの強度が所定の閾値以上のエッジに対して細線化し、前記被害後画像データに於ける画素を中心とする所定範囲の部分画像を設定し、前記各部分画像に於いて、前記部分画像内で、前記細線化したエッジが所定の長さ以上連続しており、且つ一定の方向性を有している画素数が所定の閾値以上の、前記部分画像の一又は複数の画素を非被害域候補画素として判定する、被害域自動抽出システムである。
【００１７】
エッジによる非被害域候補画素は、本発明のように判定すると良い。建造物等の輪郭は、その建造物等の輝度と、建造物等外の輝度では異なる。例えば道路に隣接している建物を上空から撮像した場合、建物の屋根の明るさと道路の明るさは異なる。従って建造物等の輝度が不連続（輝度の境目）をエッジとして抽出することによって、その輪郭を把握することが出来る。上述のようにエッジを抽出する際には、ノイズが含まれている。従ってその不連続性が一定の強度で示されていないと、ノイズの可能性がある。逆に言えば一定の強度で示されていればそれはノイズではないと考えられる。従ってそのエッジが所定の長さ以上であって、一定の方向性を有していれば、建物のような建造物等の輪郭を示している可能性があり、そしてそれが所定以上の長さと一定の方向性を有していれば、建造物等が被害を受けていないことが想定される。被害を受けていれば、輪郭が短くなる、或いは一定方向に方向性を有さない可能性があるからである。
【００１８】
請求項６の発明は、前記第１誤認候補判定手段は、前記被害前画像データの画素を中心とする所定範囲の部分画像を設定し、前記各部分画像に於いて、前記部分画像内で、各輝度毎にその輝度を有する画素数の対応関係を作成し、前記対応関係から最頻度の画素と、前記輝度の最小値と最大値の幅を輝度幅として算出し、前記最頻度の画素数が所定の閾値以下であり、且つ前記輝度幅が所定の範囲内にある前記部分画像の一又は複数の画素を第１誤認候補画素として判定する、被害域自動抽出システムである。
【００１９】
請求項７の発明は、前記第２誤認候補判定手段は、前記被害前画像データに於いて、前記輝度が不連続な画素をエッジとして抽出し、前記エッジの強度が所定の閾値以上のエッジに対して細線化し、前記被害前画像データに於ける画素を中心とする所定範囲の部分画像を設定し、前記各部分画像に於いて、前記部分画像内で、前記細線化したエッジが所定の長さ以上連続しており、且つ一定の方向性を有している画素数が所定の閾値以上の、前記部分画像の一又は複数の画素を第２誤認候補画素として判定する、被害域自動抽出システムである。
【００２０】
被害前画像データから誤認候補とする画素は、請求項６及び請求項７の発明のように判定すると良い。第１誤認候補判定手段で、被害域候補判定手段と同様の処理をすることで、誤認候補を判定できるのは、以下の理由による。つまり画像データの撮影対象物の特性（画像データの解像度に対する対象物の大きさや密集の仕方等）によって、被害後画像データから被害域として誤認される場合、同様の処理を被害前画像データを用いても被害域として判定される部分がある。しかし被害前の画像データでは、実際には被害域があるはずがないので、それは誤認と考えられる。そうすると、誤認した第１誤認候補画素の数は、被害後の画像データの被害域候補画素の数よりは少ないと考えられる。逆に言えばこの関係が成立してなければそれは誤認として排除することが出来るので、これにより誤認率を低減させている。
【００２１】
又、第２誤認候補判定手段で、非被害域候補判定手段と同様の処理をすることで、誤認候補を判定できるのは、以下の理由による。つまり画像データの撮影対象物の特性（画像データの解像度に対する対象物の大きさや密集の仕方等）によって、被害後画像データから非被害域として誤認される場合、同様の処理を被害前画像データを用いると非被害域として判定される部分がある。そして、被害前画像データの非被害域（第２誤認候補）の方が、被害後画像データの非被害域候補よりも多いはずである。そうすると、逆に言えばこの関係が成立していなければそれは誤認であると考えられる。
【００２２】
請求項８の発明は、前記被害域判定手段は、前記被害後画像データの画素を中心とする所定範囲の部分画像を設定し、前記各部分画像に於いて、前記部分画像内の被害域候補画素数が所定の閾値以上で、且つ非被害域候補画素数が所定の閾値以下であって、更に前記第１誤認候補画素が所定の閾値以下の、前記部分画像の一又は複数の画素の一定範囲を被害域として判定する、被害域自動抽出システムである。
【００２３】
請求項９の発明は、前記被害域判定手段は、前記被害後画像データの画素を中心とする所定範囲の部分画像を設定し、前記各部分画像に於いて、前記部分画像内の被害域候補画素数が所定の閾値以上で、且つ非被害域候補画素数が所定の閾値以下であって、更に前記第２誤認候補画素が所定の閾値以下の、前記部分画像の一又は複数の画素の一定範囲を被害域として判定する、被害域自動抽出システムである。
【００２４】
被害域候補画素と非被害域候補画素と第１誤認候補又は第２誤認候補画素とに基づいて、本発明のように被害域判定を行うことで、被害後画像データのみを用いる場合よりも精度を向上させることが出来る。
【００２５】
請求項１０の発明は、前記被害域自動抽出システムは、前記被害域判定手段で判定した被害域を前記読み込んだ画像データに重ね合わせ、表示装置に出力する、被害域自動抽出システムである。
【００２６】
判定した被害域は、読み込んだ被害後画像データに重ね合わせて表示することによって、どこの部分（地域）で被害が発生しているかを認識することが出来る。
【００２７】
請求項１１の発明は、データを記憶する記憶部とプログラムの処理演算を実行する処理演算部とを有しているコンピュータ端末で実行する被害域自動抽出プログラムであって、前記被害域自動抽出プログラムは、前記処理演算部に読み込まれる画像データ読込機能と被害域候補判定機能と非被害域候補判定機能と第１誤認候補判定機能と被害域判定機能とを有しており、前記画像データ読込機能は、被害前画像データと被害後画像データとを読み込み前記記憶部に記憶し、前記被害域候補判定機能は、前記記憶部に記憶した被害後画像データの画素の輝度により被害域候補画素を判定し、それを前記記憶部に記憶し、前記非被害域候補判定機能は、前記被害後画像データに於いて、前記画素の輝度が不連続な部分をエッジとし、前記エッジの連続性及び／又は方向性を判定することにより非被害域候補画素を判定し、それを前記記憶部に記憶し、前記第１誤認候補機能は、前記被害前画像データの画素の輝度により第１誤認候補画素を判定し、それを前記記憶部に記憶し、前記被害域判定機能は、前記記憶部に記憶した被害後画像データと前記被害域候補画素と前記非被害域候補画素と前記第１誤認候補画素とにより、被害域を判定する、被害域自動抽出プログラムである。
【００２８】
請求項１２の発明は、データを記憶する記憶部とプログラムの処理演算を実行する処理演算部とを有しているコンピュータ端末で実行する被害域自動抽出プログラムであって、前記被害域自動抽出プログラムは、前記処理演算部に読み込まれる画像データ読込機能と被害域候補判定機能と非被害域候補判定機能と第２誤認候補判定機能と被害域判定機能とを有しており、前記画像データ読込機能は、被害前画像データと被害後画像データとを読み込み前記記憶部に記憶し、前記被害域候補判定機能は、前記記憶部に記憶した被害後画像データの画素の輝度により被害域候補画素を判定し、それを前記記憶部に記憶し、前記非被害域候補判定機能は、前記被害後画像データに於いて、前記画素の輝度が不連続な部分をエッジとし、前記エッジの連続性及び／又は方向性を判定することにより非被害域候補画素を判定し、それを前記記憶部に記憶し、前記第２誤認候補判定機能は、前記被害前画像データに於いて、前記画素の輝度が不連続な部分をエッジとし、前記エッジの連続性及び／又は方向性を判定することにより第２誤認候補画素を判定し、それを前記記憶部に記憶し、前記記憶部に記憶した被害後画像データと前記被害域候補画素と前記非被害域候補画素と前記第２誤認候補画素とにより、被害域を判定する、被害域自動抽出プログラムである。
【００２９】
請求項１３の発明は、データを記憶する記憶部とプログラムの処理演算を実行する処理演算部とを有しているコンピュータ端末で実行する被害域自動抽出プログラムであって、前記被害域自動抽出プログラムは、前記処理演算部に読み込まれる画像データ読込機能と被害域候補判定機能と非被害域候補判定機能と第１誤認候補判定機能と第２誤認候補判定機能と被害域判定機能とを有しており、前記画像データ読込機能は、被害前画像データと被害後画像データとを読み込み前記記憶部に記憶し、前記被害域候補判定機能は、前記記憶部に記憶した被害後画像データの画素の輝度により被害域候補画素を判定し、それを前記記憶部に記憶し、前記非被害域候補判定機能は、前記被害後画像データに於いて、前記画素の輝度が不連続な部分をエッジとし、前記エッジの連続性及び／又は方向性を判定することにより非被害域候補画素を判定し、それを前記記憶部に記憶し、前記第１誤認候補機能は、前記被害前画像データの画素の輝度により第１誤認候補画素を判定し、それを前記記憶部に記憶し、前記第２誤認候補機能は、前記被害前画像データに於いて、前記画素の輝度が不連続な部分をエッジとし、前記エッジの連続性及び／又は方向性を判定することにより第２誤認候補画素を判定し、それを前記記憶部に記憶し、前記被害域判定機能は、前記記憶部に記憶した被害後画像データと前記被害域候補画素と前記非被害域候補画素と前記第１誤認候補画素と前記第２誤認候補画素とにより、被害域を判定する、被害域自動抽出プログラムである。
【００３０】
上述の被害域自動抽出システムは、プログラムとしてネットワークからダウンロードされ、それがコンピュータ端末に読み込まれることによって実現されても良い。
【００３１】
請求項１４の発明は、データを記憶する記憶部とプログラムの処理演算を実行する処理演算部とを有しているコンピュータ端末で実行する被害域自動抽出プログラムを記録した記録媒体であって、前記被害域自動抽出プログラムは、前記処理演算部に読み込まれる画像データ読込機能と被害域候補判定機能と非被害域候補判定機能と第１誤認候補判定機能と被害域判定機能とを有しており、前記画像データ読込機能は、前記画像データを読み込み前記記憶部に記憶し、前記被害域候補判定機能は、前記記憶部に記憶した画像データの画素の輝度により被害域候補画素を判定し、それを前記記憶部に記憶し、前記非被害域候補判定機能は、前記画像データに於いて、前記画素の輝度が不連続な部分をエッジとし、前記エッジの連続性及び／又は方向性を判定することにより非被害域候補画素を判定し、それを前記記憶部に記憶し、前記第１誤認候補機能は、前記被害前画像データの画素の輝度により第１誤認候補画素を判定し、それを前記記憶部に記憶し、前記被害域判定機能は、前記記憶部に記憶した被害後画像データと前記被害域候補画素と前記非被害域候補画素と前記第１誤認候補画素とにより、被害域を判定する、被害域自動抽出プログラムを記録した記録媒体である。
【００３２】
請求項１５の発明は、データを記憶する記憶部とプログラムの処理演算を実行する処理演算部とを有しているコンピュータ端末で実行する被害域自動抽出プログラムを記録した記録媒体であって、前記被害域自動抽出プログラムは、前記処理演算部に読み込まれる画像データ読込機能と被害域候補判定機能と非被害域候補判定機能と第２誤認候補判定機能と被害域判定機能とを有しており、前記画像データ読込機能は、前記画像データを読み込み前記記憶部に記憶し、前記被害域候補判定機能は、前記記憶部に記憶した画像データの画素の輝度により被害域候補画素を判定し、それを前記記憶部に記憶し、前記非被害域候補判定機能は、前記画像データに於いて、前記画素の輝度が不連続な部分をエッジとし、前記エッジの連続性及び／又は方向性を判定することにより非被害域候補画素を判定し、それを前記記憶部に記憶し、前記第２誤認候補判定機能は、前記被害前画像データに於いて、前記画素の輝度が不連続な部分をエッジとし、前記エッジの連続性及び／又は方向性を判定することにより第２誤認候補画素を判定し、それを前記記憶部に記憶し、前記被害域判定機能は、前記記憶部に記憶した被害後画像データと前記被害域候補画素と前記非被害域候補画素と前記第２誤認候補画素とにより、被害域を判定する、被害域自動抽出プログラムを記録した記録媒体である。
【００３３】
請求項１６の発明は、データを記憶する記憶部とプログラムの処理演算を実行する処理演算部とを有しているコンピュータ端末で実行する被害域自動抽出プログラムを記録した記録媒体であって、前記被害域自動抽出プログラムは、前記処理演算部に読み込まれる画像データ読込機能と被害域候補判定機能と非被害域候補判定機能と第１誤認候補判定機能と第２誤認候補判定機能と被害域判定機能とを有しており、前記画像データ読込機能は、前記画像データを読み込み前記記憶部に記憶し、前記被害域候補判定機能は、前記記憶部に記憶した画像データの画素の輝度により被害域候補画素を判定し、それを前記記憶部に記憶し、前記非被害域候補判定機能は、前記画像データに於いて、前記画素の輝度が不連続な部分をエッジとし、前記エッジの連続性及び／又は方向性を判定することにより非被害域候補画素を判定し、それを前記記憶部に記憶し、前記第１誤認候補機能は、前記被害前画像データの画素の輝度により第１誤認候補画素を判定し、それを前記記憶部に記憶し、前記第２誤認候補判定機能は、前記被害前画像データに於いて、前記画素の輝度が不連続な部分をエッジとし、前記エッジの連続性及び／又は方向性を判定することにより第２誤認候補画素を判定し、それを前記記憶部に記憶し、前記被害域判定機能は、前記記憶部に記憶した被害後画像データと前記被害域候補画素と前記非被害域候補画素と前記第１誤認候補画素と前記第２誤認候補判定機能とにより、被害域を判定する、被害域自動抽出プログラムを記録した記録媒体である。
【００３４】
上述の被害域自動抽出システムは、そのプログラムを記録した記録媒体がコンピュータ端末に読み込まれることによって実現されても良い。
【発明の効果】
【００３５】
本発明により、被害後の画像データから自動的に被害域を抽出し、又、被害前の画像データから誤認発生箇所を抽出することで、精度が高い被害域自動抽出システムが可能となる。
【発明を実施するための最良の形態】
【００３６】
本発明の被害域自動抽出システム１のシステム構成の一例を示すシステム構成図を図１に示す。被害域自動抽出システム１は、画像データ読込手段２と被害域候補判定手段３と非被害域候補判定手段４と第１誤認候補判定手段５と第２誤認候補判定手段６と被害域判定手段７と画像データ出力手段８とを有している。
【００３７】
画像データ読込手段２は、人工衛星や航空機等により撮像した被害後の画像データを被害後画像データとして、被害前の画像データを被害前画像データとして、被害域自動抽出システム１に読み込む手段である。尚、画像データ読込手段２で読み込む画像データは、人工衛星や航空機等により撮像した画像データをそのまま読み込んでも良いし、システム補正、オルソ補正、精密幾何補正を施した後の画像データであっても良い。人工衛星や航空機等により撮像した場合、地球の自転、地球の湾曲、衛星の姿勢変化等によって画像データに歪みが発生しているので、それを補正した方が画像データとして扱いやすいからである。又、撮像した対象物の高さについても、歪みが発生しているので、それを正射投影に対する投射位置誤差を補正する補正処理（オルソ補正）を行うことによって、画像データとしてより扱いやすいものとなる。
【００３８】
被害域候補判定手段３は、画像データ読込手段２で読み込んだ被害後画像データから、被害域（被害が発生している領域）の画素候補を判定する手段である。
【００３９】
具体的には、被害後画像データ内の画素を中心として、所定範囲（例えば１５ｍ四方相当部分の画像範囲）を部分画像として選択し、部分画像内の輝度（単位面積あたりの光束の量。一定の広さ部分の明るさのこと。アスキーデジタル用語辞典、http://yougo.ascii24.com）をレベルスライス（各輝度のレベルに応じて振り分ける）し、そのヒストグラムを作成する。図７にヒストグラムの一例を示す。このヒストグラムから最頻度値、輝度の最大値と最小値の差（輝度幅）を算出する。そして最頻度値が所定の閾値（例えば部分画像内の画素数の１２％）より小さく、且つ輝度幅が所定の２つの閾値の範囲内（例えば５００以上、２０４８以下）にあれば、その中心画素を被害域候補画素として判定する。尚、以下の実施例では部分画像の中心画素を被害域候補画素とする場合を説明するが、部分画像内の一又は複数の任意の画素であっても良い。
【００４０】
輝度の最頻度の画素数が所定の閾値以下であることは、部分画像に於いて特定の輝度が突出していないことを示しており、上述のように、これが被害域の候補であることの一つの条件となる。輝度幅が所定の範囲内にあることは、極端に輝度の変化がないことを示している。つまり、被害が発生していない建造物等であれば、部分画像内に於いて輝度幅に変化は少ない（同一の建造物等であるので、その屋根や壁等は同一の明るさ、それに類似する明るさであることが想定される）。逆に言えば、部分画像内に於いて輝度幅に変化が大きい場合、それは屋根や壁等が被害を受けたことによって、部分画像内に様々な明るさが含まれていると考えられる。そうすると、その部分画像の画素は、被害域の候補の可能性があることを示唆することとなる。そしてこの２つの条件を満たした場合に、その部分画像に於ける画素を被害域の候補画素として判定することが出来る。
【００４１】
被害域候補の概念を図１６に示す。図１６（ａ）は被害を受けていない建造物であり、図１６（ｂ）は被害を受けている建造物である。そして図１６（ｃ）は図１６（ａ）の被害を受けていない建造物の各輝度毎のヒストグラムであり、図１６（ｄ）は図１６（ｂ）の被害を受けた建造物の各輝度毎のヒストグラムである。被害を受けていない建造物の場合、上述のように、所定の輝度の頻度が高くなり、且つ輝度の分布のばらつきが少ない。一方、被害を受けた建造物の場合、上述のように、突出した輝度はなくなり、各輝度にばらつきが見られる。その為、このように輝度に基づいて被害域候補画素の判定が可能となる。
【００４２】
非被害域候補判定手段４は、画像データ読込手段２で読み込んだ被害後画像データから、非被害域（被害が発生していない領域）の画素候補を判定する手段である。
【００４３】
具体的には、被害後画像データから輝度の不連続な部分をエッジとして抽出し、エッジ強度（エッジ強度とは、隣接する画素の輝度差を言う）が所定の閾値以上（例えば０．０８以上）のエッジだけを残し、それに対して細線化処理を行う。尚、細線化処理とは、２値画像に含まれる各種の連結図形に対し、連結性を失うことなく線幅が１の線図形とする処理である。細線化処理を行った画像（細線化処理画像）に対して、細線化されたエッジが連続している画素数を調べ、所定の閾値以上（例えば１０画素以上）に連続しているエッジのみを抽出し、それ以外はエッジから除去する。そして被害後画像データ内の画素を中心として、所定範囲（例えば１５ｍ四方相当部分の画像範囲）を部分画像として選択し、その部分画像内にエッジとして存在する画素の数を調べる。この際に、エッジの各点に於ける方向が、ある一定の方向に向いている画素の数のみを数え、各点の方向が分散しているエッジは除外する。このようにして数えたエッジの数が所定の閾値以上（部分画像内の画素数の１％以上）であれば、部分画像の中心画素を非被害域候補画素として判定する。尚、以下の実施例では部分画像の中心画素を非被害域候補画素とする場合を説明するが、部分画像内の一又は複数の任意の画素であっても良い。
【００４４】
尚、上述の方向性の判定は、以下のように行えば良い。まず、部分画像内の所定の長さ以上を有するエッジを構成する全部又はその連続する一部の画素を抽出する。そしてある画素とその隣り合う画素の位置関係、即ち、上方向、右上方向、右方向、右下方向、下方向、左下方向、左方向、左上方向のいずれか、を判定する。その位置関係の判定を、抽出した隣り合う各画素同士について行う。そして隣り合う画素同士の位置関係のうち、どの位置関係が最も多いか（最頻度方向）を算出し、それが所定の閾値以上であれば、一定の方向性を有すると判定する。つまり、隣り合う各画素同士の位置関係を調べれば、隣り合う画素間の方向が分かる。そしてエッジは隣り合う各画素同士の連続性で構成されていることから、エッジに於ける最頻度方向を算出し、且つ最頻度方向が所定の閾値以上であれば、エッジ全体がどの方向性を有しているかが分かる。
【００４５】
非被害域候補の概念を図１７に示す。図１７（ａ）は被害を受けていない建造物であり、図１７（ｂ）は被害を受けている建造物である。図１７（ａ）の被害を受けていない建造物のエッジは、その建造物の輪郭がはっきりしているので、所定の長さ、方向性を有しているが、図１７（ｂ）の被害を受けた建造物のエッジは、その建造物が倒壊したことにより、輪郭がはっきりせず、所定の長さ、方向性を有していない。その為、このようにエッジに基づいて非被害域候補画素の判定が可能となる。
【００４６】
第１誤認候補判定手段５は、画像データ読込手段２で読み込んだ被害前画像データから、画素の輝度に基づいて、第１誤認候補画素を判定する手段である。
【００４７】
具体的には、被害域候補判定手段３と同様な処理を行う。即ち、被害前画像データ内の画素を中心として、所定範囲（例えば１５ｍ四方相当部分の画像範囲）を部分画像として選択し、部分画像内の輝度をレベルスライス（各輝度のレベルに応じて振り分ける）し、そのヒストグラム（図６）を作成する。このヒストグラムから最頻度値、輝度の最大値と最小値の差（輝度幅）を算出する。そして最頻度値が所定の閾値（例えば部分画像内の画素数の１２％）より小さく、且つ輝度幅が所定の２つの閾値の範囲内（例えば５００以上、２０４８以下）にあれば、その画素を第１誤認候補画素として判定する。尚、以下の実施例では部分画像の中心画素を第１誤認候補画素とする場合を説明するが、部分画像内の一又は複数の任意の画素であっても良い。
【００４８】
このように第１誤認候補判定手段５で、被害域候補判定手段３と同様の処理をすることで、誤認候補を判定できるのは、以下の理由による。つまり画像データの撮影対象物の特性（画像データの解像度に対する対象物の大きさや密集の仕方等）によって、被害後画像データから被害域として誤認される場合、同様の処理を被害前画像データを用いても被害域として判定される部分がある。しかし被害前の画像データでは、実際には被害域があるはずがないので、それは誤認と考えられる。そうすると、誤認した第１誤認候補画素の数は、被害後の画像データの被害域候補画素の数よりは少ないと考えられる。逆に言えばこの関係が成立してなければそれは誤認として排除することが出来るので、これにより誤認率を低減させている。
【００４９】
第２誤認候補判定手段６は、画像データ読込手段２で読み込んだ被害前画像データから、画素のエッジに基づいて、第２誤認候補画素を判定する手段である。
【００５０】
具体的には、非被害域候補判定手段４と同様な処理を行う。即ち、被害前画像データから輝度の不連続な部分をエッジとして抽出し、エッジ強度が所定の閾値以上（例えば０．０８以上）のエッジだけを残し、それに対して細線化処理を行う。細線化処理を行った画像（細線化処理画像）に対して、細線化されたエッジが連続している画素数を調べ、所定の閾値以上（例えば１０画素以上）に連続しているエッジのみを抽出し、それ以外はエッジから除去する。そして被害前画像データ内の画素を中心として、所定範囲（例えば１５ｍ四方相当部分の画像範囲）を部分画像として選択し、その部分画像内にエッジとして存在する画素の数を調べる。この際に、エッジの各点に於ける方向が、ある一定の方向に向いている画素の数のみを数え、各点の方向が分散しているエッジは除外する。このようにして数えたエッジの数が所定の閾値以上（例えば部分画像内の画素数の１０％以上）であれば、部分画像の中心画素を第２誤認候補画素として判定する。尚、以下の実施例では部分画像の中心画素を第２誤認候補画素とする場合を説明するが、部分画像内の一又は複数の任意の画素であっても良い。
【００５１】
このように第２誤認候補判定手段６で、非被害域候補判定手段４と同様の処理をすることで、誤認候補を判定できるのは、以下の理由による。つまり画像データの撮影対象物の特性（画像データの解像度に対する対象物の大きさや密集の仕方等）によって、被害後画像データから非被害域として誤認される場合、同様の処理を被害前画像データを用いると非被害域として判定される部分がある。そして、被害前画像データの非被害域（第２誤認候補）の方が、被害後画像データの非被害域候補よりも多いはずである。そうすると、逆に言えばこの関係が成立していなければそれは誤認であると考えられる。
【００５２】
尚、上述の第１誤認候補判定手段５及び第２誤認候補判定手段６は、被害域自動抽出システム１にいずれかが備えられていればよいが、双方あることによって、精度を向上させることが出来る。
【００５３】
被害域判定手段７は、被害域候補判定手段３で判定した被害域候補画素と、非被害域候補判定手段４で判定した非被害域候補画素と、第１誤認候補判定手段５で判定した第１誤認候補と、第２誤認候補判定手段６で判定した第２誤認候補とに基づいて、被害域を判定する手段である。
【００５４】
具体的には、被害後画像データ内の画素を中心として、所定範囲（例えば１５ｍ四方相当部分の画像範囲）を部分画像として選択し、被害域候補画素の数が所定の閾値以上（例えば部分画像内の画素数の２０％以上）で、且つ非被害域候補画素の数が所定の閾値以下（例えば部分画像内の画素数の２５％以下）であって、更に第１誤認候補画素の数が所定の閾値以下（例えば部分画像内の画素数の１０％以下）であり、第２誤認候補がその数が所定の閾値以上（例えば部分画像内の画素数の２５％以上）である場合、その画素を中心とした所定の範囲を、被害域として判定する。尚、以下の実施例では部分画像の中心画素を中心とした所定範囲を被害域とする場合を説明するが、部分画像内の一又は複数の任意の画素を中心とした所定範囲を被害域としても良い。又、被害域の判定には、第１誤認候補画素、第２誤認候補画素の双方を利用しているが、いずれか一方のみであっても良い。
【００５５】
又、被害域候補画素の閾値を閾値Ａ、非被害域候補画素の閾値を閾値Ｂ、第１誤認候補画素の閾値を閾値Ｃ、第２誤認候補画素の閾値を閾値Ｄとする。被害前画像データでは、実際は被害が発生しているはずがないので、本来ならばそれに対する輝度の分析の結果抽出される第１誤認候補画素は判定されない。しかし実際はノイズや撮像された対象物の影響等の様々な理由から第１誤認候補画素が複数判定される。しかし、その数は被害域候補画素よりは少なくなると考えられる。即ち被害域として判定される部分画像では、「被害域候補画素の数＞第１誤認候補画素の数」の関係が成立すると考えられる。
【００５６】
一方、同様に被害前画像データに対するエッジの分析の結果抽出される第２誤認候補画素は、建造物等が倒壊する（被害を受ける）とその数は減るはずである。従って、被害前画像データから判定した第２誤認候補画素の数より、被害後画像データから判定した非被害域候補画素の数の方が少なくなると考えられる。即ち、被害域として判定される部分画像では、「非被害域候補画素の数＜第２誤認候補画素の数」の関係が成立すると考えられる。
【００５７】
従って少なくとも、閾値Ａは当該部分画像に於ける第１誤認候補画素の数、閾値Ｂは当該部分画像に於ける第２誤認候補画素の数、閾値Ｃは当該部分画像に於ける被害域候補画素の数、閾値Ｄは当該部分画像に於ける非被害域候補画素の数を満たしていればよい。そうすると、被害域判定手段７での被害域の判定は、「被害域候補画素の数＞第１誤認候補画素の数」と「非被害域候補画素の数＜第２誤認候補画素の数」の２つを満たす場合に、部分画像内の一又は複数の任意の画素を中心とした所定範囲を被害域とすることが出来る。
【００５８】
より精度を上げる為には、上述のように、「被害域候補画素の数≧閾値Ａ、非被害域候補画素の数≦閾値Ｂ、第１誤認候補画素の数≦閾値Ｃ、第２誤認候補の数≧閾値Ｄ（但し、閾値Ａ≦閾値Ｃ、閾値Ｄ≦閾値Ｂ）」とすると良い。
【００５９】
画像データ出力手段８は、画像データ読込手段２で読み込んだ被害後画像データに、被害域判定手段７で判定した被害域を重ね合わせ、表示装置等に出力する手段である。
【００６０】
次に本発明の被害域自動抽出システム１の処理の流れの一例を図２から図６のフローチャート図と図１のシステム構成図とを用いて説明する。
【００６１】
人工衛星や航空機で撮像した被害前の画像データ（被害前画像データ）と被害後の画像データ（被害後画像データ）とを画像データ読込手段２で被害域自動抽出システム１に取り込む。（Ｓ１００）。この取り込んだ被害後画像データは、コンピュータ端末のメモリ、キャッシュ、ハードディスク等の記憶部に記憶される。この被害後画像データの一例を図８に示す。
【００６２】
次に記憶部に記憶した被害後画像データから、輝度情報に基づいて被害域候補画素の判定処理を被害域候補判定手段３が実行する（Ｓ２００）。
【００６３】
被害域候補画素の判定処理（Ｓ２００）は、まず、被害後画像データ内の画素を中心として、所定範囲（例えば１５ｍ四方相当部分の画像範囲）を部分画像として選択し（Ｓ２１０）、部分画像内の輝度をレベルスライス（各輝度のレベルに応じて振り分ける）する（Ｓ２３０）。部分画像内の各画素の輝度をレベルスライスした後、そのヒストグラムを図７に示すように作成する（Ｓ２３０）。尚、実際にヒストグラムとしてグラフ化されていなくても、各輝度毎にその輝度を有する画素がいくつあるかの、対応関係を示すものをヒストグラムとする。
【００６４】
そして輝度とその輝度を有する画素の数の対応関係から、画素数が最も多い輝度を最頻度値とし、又、輝度の最大値と最小値の差を輝度幅として算出する（Ｓ２４０）。
【００６５】
Ｓ２４０で算出した最頻度値と輝度幅とに基づいて、最頻度値が所定の閾値以下（例えば部分画像内の画素数の１２％以下）であり、且つ輝度幅が所定の範囲内（例えば５００以上、２０４８以下）にあれば（Ｓ２５０）、その部分画像の中心画素を被害域候補画素として判定する（Ｓ２６０）。判定した被害域候補画素は、上述の記憶部に記憶する。
【００６６】
被害域候補画素の判定を行った場合、又はＳ２５０の条件を満たさなければ、次の部分画像に対して、Ｓ２３０からＳ２６０の処理を実行する（Ｓ２７０）。これを、中心とする画素が一定間隔毎（例えば４画素毎）の全ての部分画像に対して実行するまで反復する（Ｓ２２０）。
【００６７】
被害域候補画素を判定後、非被害域候補判定手段４は、記憶部に記憶した被害後画像データから、エッジ情報に基づいて、非被害域候補画素の判定処理を実行する（Ｓ３００）。
【００６８】
非被害域候補画素の判定処理（Ｓ３００）は、まず、被害後画像データの画素に対して、その輝度が不連続な部分をエッジとして抽出する（Ｓ３１０）。図８の被害後画像データからエッジを抽出した場合の画像データを図９に示す。そしてＳ３１０に於いてエッジを抽出した画像データから、エッジ強度が所定の閾値以上（例えば０．０８以上）のエッジだけを抽出する（Ｓ３２０）。この場合の画像データを図１０に示す。
【００６９】
Ｓ３２０で抽出した強度の大きなエッジを抽出した画像データに対して、細線化処理を行う（Ｓ３３０）。図１０の画像データに対して細線化処理を施した細線化処理画像データを図１１に示す。
【００７０】
細線化処理画像データに対して、細線化されたエッジが連続している画素数を調べ、その数が所定の閾値以上（例えば１０画素以上）に連続しているエッジのみを抽出し、それ以外のエッジは除去する（Ｓ３４０）。つまり細線化処理画像データから、所定の閾値以上（例えば１０画素以上）に長いエッジのみを抽出することとなる。この場合の画像データを図１２に示す。
【００７１】
次に、被害後画像データ内の画素を中心として、所定範囲（例えば１５ｍ四方相当部分の画像範囲）を部分画像として選択し、Ｓ３４０で作成した画像データと対比させ、その部分画像内にエッジとして存在する画素の数をカウントする。この際に、エッジの各点に於ける方向が、ある一定の方向に向いている画素の数のみを数え、各点の方向が分散しているエッジは除去する（Ｓ３６０）。
【００７２】
このようにして数えたエッジの数が所定の閾値以上（例えば部分画像内の画素数の１％以上）であれば（Ｓ３７０）、部分画像の中心画素を非被害域候補画素として判定する（Ｓ３８０）。判定した非被害域候補画素は、上述の記憶部に記憶する。
【００７３】
非被害域候補画素の判定を行った場合、又はＳ３７０の条件を満たさなければ、次の部分画像に対して、Ｓ３６０からＳ３８０の処理を実行する（Ｓ３９０）。これを、中心とする画素が一定間隔毎（例えば１画素毎）の全ての部分画像に対して実行するまで反復する（Ｓ３５０）。図１１の画像データから抽出した、一定方向に向いたエッジを抽出した画像データを図１２に示す。
【００７４】
このように、単にエッジを抽出しただけの画像データでは、様々なノイズや建造物倒壊による瓦礫によるエッジも抽出されてしまうが、エッジ強度、エッジを構成するがその連続する長さ、エッジの方向によって選別することによって、上述のノイズや瓦礫等によるエッジは除去され、建造物としての構造がはっきり残っているものだけを非被害候補として抽出することが出来る。
【００７５】
非被害域候補判定手段４によるＳ３００の判定処理後、第１誤認候補判定手段５は、被害前画像データに基づいて、第１誤認候補画素を判定する（Ｓ４００）。
【００７６】
具体的には、被害域候補判定手段３と同様な処理を行う。即ち、被害前画像データ内の画素を中心として、所定範囲（例えば１５ｍ四方相当部分の画像範囲）を部分画像として選択し（Ｓ４１０）、部分画像内の輝度をレベルスライス（各輝度のレベルに応じて振り分ける）し、そのヒストグラム（図７）を作成する（Ｓ４３０）。このヒストグラムから最頻度値、輝度の最大値と最小値の差（輝度幅）を算出する（Ｓ４４０）。そして最頻度値が所定の閾値（例えば部分画像内の画素数の１２％）より小さく、且つ輝度幅が所定の２つの閾値の範囲内（例えば５００以上、２０４８以下）にあれば（Ｓ４５０）、その画素を第１誤認候補画素として判定する（Ｓ４６０）。判定した第１誤認候補画素は、上述の記憶部に記憶する。
【００７７】
第１誤認候補画素の判定を行った場合、又はＳ４５０の条件を満たさなければ、次の部分画像に対して、Ｓ４３０からＳ４６０の処理を実行する（Ｓ４７０）。これを、中心とする画素が一定間隔毎（例えば４画素毎）の全ての部分画像に対して実行するまで反復する（Ｓ４２０）。
【００７８】
第１誤認候補判定手段５によるＳ４００の判定処理後、第２誤認候補判定手段６は、被害前画像データに基づいて、第２誤認候補画素を判定する（Ｓ５００）。
【００７９】
具体的には、非被害域候補判定手段４と同様な処理を行う。即ち、被害前画像データから輝度の不連続な部分をエッジとして抽出し（Ｓ５１０）、エッジ強度が所定の閾値以上（例えば０．０８以上）のエッジだけを残し（Ｓ５２０）、それに対して細線化処理を行う（Ｓ５３０）。細線化処理を行った画像（細線化処理画像）に対して、細線化されたエッジが連続している画素数を調べ、所定の閾値以上（例えば１０画素以上）に連続しているエッジのみを抽出し、それ以外はエッジから除去する（Ｓ５４０）。そして被害前画像データ内の画素を中心として、所定範囲（例えば１５ｍ四方相当部分の画像範囲）を部分画像として選択し、その部分画像内にエッジとして存在する画素の数を調べる。この際に、エッジの各点に於ける方向が、ある一定の方向に向いている画素の数のみを数え、各点の方向が分散しているエッジは除外する（Ｓ５６０）。このようにして数えたエッジの数が所定の閾値以上（例えば部分画像内の画素数の１％以上）であれば（Ｓ５７０）、部分画像の中心画素を第２誤認候補画素として判定する（Ｓ５８０）。判定した第２誤認候補画素は、上述の記憶部に記憶する。
【００８０】
第２誤認候補画素の判定を行った場合、又はＳ５７０の条件を満たさなければ、次の部分画像に対して、Ｓ５６０からＳ５８０の処理を実行する（Ｓ５９０）。これを、中心とする画素が一定間隔毎（例えば１画素毎）の全ての部分画像に対して実行するまで反復する（Ｓ５５０）。
【００８１】
以上のような処理プロセスによって、被害域候補画素、非被害域候補画素、第１誤認候補画素、第２誤認候補画素を判定後、被害域判定手段７は、記憶部に記憶した被害後画像データと被害域候補画素と非被害域候補画素と第１誤認候補画素と第２誤認候補画素とに基づいて、被害域判定処理を実行する（Ｓ６００）。
【００８２】
具体的には、被害後画像データ内の画素を中心として、所定範囲（例えば１５ｍ四方相当部分の画像範囲）を部分画像として選択し、被害域候補画素の数が所定の閾値以上（例えば部分画像内の画素数の２０％以上）で、且つ非被害域候補画素の数が所定の閾値以下（例えば部分画像内の画素数の２５％以下）であって、更に第１誤認候補画素の数が所定の閾値以下（例えば部分画像内の画素数の１０％以下）であり、第２誤認候補がその数が所定の閾値以上（例えば部分画像内の画素数の２５％以上）である場合、その画素を中心とした所定の範囲を、被害域として判定する。
【００８３】
このような処理プロセスを経ることによって、被害域自動抽出システム１は、人工衛星や航空機により撮像した画像データ（被害後画像データ）から、自動的に被害域を判定することが可能となり、又、更に被害前画像データから誤認可能性のある画素を判定していることによって、誤認率も下げることが可能となる。
【００８４】
そして画像データ出力手段８は、被害域判定手段７で判定した被害域を被害後画像データに重ね合わせることによって、ディスプレイ等の表示装置に出力する（Ｓ６１０）。
【００８５】
尚、被害域自動抽出システム１に於いて、予め被害前画像データ、被害後画像データに対して、植生域、水域、雲域、影域と見なすことが出来る画素については、事前にマスク処理を行い、被害域判定を行わないこととしても良い。
【００８６】
輝度のみによる被害域の画素候補の判定の場合では、輝度による特徴量が同じ値であっても、実際に被害域である場合と、例えば輝度の異なる複数の建造物が密集しているような場合があり、このような場合には被害域の誤認が生ずる可能性がある。又、エッジによる特徴量が同じ値であっても、例えば被害域に影が出来ていて、その影のエッジが特徴量に反映される場合と、非被害域で建造物によるエッジが反映される場合とがあり、このような場合にも非被害域の誤認が生ずる可能性がある。
【００８７】
一方、上述したように本発明では、輝度による被害域の画素候補の判定とエッジによる非被害域の画素候補の判定とを併用している。これによって、例えば異なる複数の建造物が密集している場合では、建造物のエッジによってエッジの特徴量が非被害を示したり、又、被害域に於いて影がある場合では輝度の特徴量が被害であることを示したりするので、輝度のみ、エッジのみによる判定よりも、非被害域の判定精度の向上を図ることが出来る。この概念図を図１４に示す。
【００８８】
本発明の被害域自動抽出システム１を用いた実例を図１５に示す。図１５（ａ）は被害前画像データであり、図１５（ｂ）は被害前画像データを用いずに被害後画像データのみで被害域を自動抽出した結果であり、図１５（ｃ）は被害後画像データと被害前画像データを用いて被害域を自動抽出した結果（即ち本発明を用いた場合）である。
【００８９】
このように被害前の画像データを用いない場合では誤判定されていた領域を、被害前の画像データを用いることによって除外することが出来、精度の向上を図ることが出来る。
【００９０】
尚、上述の実施例では、被害前画像データから誤認候補画素を判定する際に、画素の輝度に基づく第１誤認候補と、エッジに基づく第２誤認候補の２つを用いた場合を説明したが、そのうちのいずれかのみを用いても良い。
【００９１】
又、上述した被害域自動抽出システム１の処理の内容を記録した被害域自動抽出プログラムがＣＤ−ＲＯＭやＤＶＤ等の記録媒体に記録され、それがコンピュータ端末に読み込まれることによって、被害域自動抽出システム１を実現しても良い。又、当該被害域自動抽出プログラムが所定のサーバからネットワークを介してコンピュータ端末にダウンロードされることで、コンピュータ端末に読み込まれ、被害域自動抽出システム１を実現しても良い。
【００９２】
この場合被害域自動抽出プログラムは、画像データ読込手段２をコンピュータ端末で実行する機能（モジュールやプログラム。以下同様）と、被害域候補判定手段３をコンピュータ端末で実行する機能と、非被害域候補判定手段４をコンピュータ端末で実行する機能と、被害域判定手段７をコンピュータ端末で実行する機能と、第１誤認候補判定手段５をコンピュータ端末で実行する機能と、第２誤認候補判定手段６をコンピュータ端末で実行する機能と、画像データ出力手段８をコンピュータ端末で実行する機能とを有している。
【００９３】
この被害域自動抽出プログラムを記録した記録媒体がコンピュータ端末にセットされ、それがコンピュータ端末に読み込まれることによって、上述の各機能がコンピュータ端末の処理演算部（ＣＰＵやレジスタ等の演算装置）に読み込まれる。そして、処理演算部に読み込まれた画像データ読込手段２は、被害域自動抽出システム１で読み込んだ被害前画像データと被害後画像データを、コンピュータ端末の記憶部（メモリやキャッシュ、ハードディスク等）に記憶する。そして処理演算部に読み込まれた被害域候補判定手段３と非被害域候補判定手段４は、記憶部に記憶した被害後画像データに基づいて、被害域候補画素と非被害域候補画素とを判定し、それらを記憶部に記憶する。又、処理演算部に読み込まれた第１誤認候補判定手段５と第２誤認候補判定手段６は、記憶部に記憶した被害前画像データに基づいて、第１誤認候補画素と第２誤認候補画素とを判定し、それらを記憶部に記憶する。そして、処理演算部に読み込まれた被害域判定手段７が、記憶部に記憶した被害後画像データと被害域候補画素と非被害域候補画素と第１誤認候補画素と第２誤認候補画素とに基づいて被害域の判定を行い、それを画像データ出力手段８が、被害後画像データに重ね合わせて出力することによって、コンピュータ端末のディスプレイ等の表示装置に於いて、被害域の自動判定の結果を確認することが出来る。
【００９４】
当該被害域自動抽出プログラムがネットワークからコンピュータ端末にダウンロードされた場合も、被害域自動抽出プログラムが記録された記録媒体をコンピュータ端末で実行する場合と同様である。
【００９５】
本発明に於ける各手段は、その機能が論理的に区別されているのみであって、物理上あるいは事実上は同一の領域を為していても良い。
【００９６】
尚、本発明を実施するにあたり本実施態様の機能を実現するソフトウェアのプログラムを記録した記憶媒体をシステムに供給し、そのシステムのコンピュータが記憶媒体に格納されたプログラムを読み出し実行することによって実現されることは当然である。
【００９７】
この場合、記憶媒体から読み出されたプログラム自体が前記した実施態様の機能を実現することとなり、そのプログラムを記憶した記憶媒体は本発明を当然のことながら構成することになる。
【００９８】
プログラムを供給する為の記憶媒体としては、例えば磁気ディスク、ハードディスク、光ディスク、光磁気ディスク、磁気テープ、不揮発性のメモリカード等を使用することができる。
【００９９】
又、コンピュータが読み出したプログラムを実行することにより、上述した実施態様の機能が実現されるだけではなく、そのプログラムの指示に基づき、コンピュータ上で稼働しているオペレーティングシステムなどが実際の処理の一部又は全部を行い、その処理によって前記した実施態様の機能が実現される場合も含まれることは言うまでもない。
【０１００】
更に、記憶媒体から読み出されたプログラムが、コンピュータに挿入された機能拡張ボードやコンピュータに接続された機能拡張ユニットに備わる不揮発性あるいは揮発性の記憶手段に書き込まれた後、そのプログラムの指示に基づき、機能拡張ボードあるいは機能拡張ユニットに備わる演算処理装置などが実際の処理の一部あるいは全部を行い、その処理により前記した実施態様の機能が実現される場合も含まれることは当然である。 
【産業上の利用可能性】
【０１０１】
上述したように、被害後画像データだけでは除去できない誤認とされる領域を、被害前画像データを併用することによって除去でき、被害後画像データのみを用いた場合よりもより精度の良い抽出結果を得ることが出来る。
【図面の簡単な説明】
【０１０２】
【図１】本発明のシステム構成の一例を示すシステム構成図である。
【図２】本発明の全体処理のプロセスの一例を示すフローチャート図である。
【図３】本発明の被害域候補画素の判定のプロセスの一例を示すフローチャート図である。
【図４】本発明の非被害域候補画素の判定のプロセスの一例を示すフローチャート図である。
【図５】本発明の第１誤認候補画素の判定のプロセスの一例を示すフローチャート図である。
【図６】本発明の第２誤認候補画素の判定のプロセスの一例を示すフローチャート図である。
【図７】ヒストグラムの一例である。
【図８】被害後画像データの一例である。
【図９】被害後画像データからエッジを抽出した場合の画像データ一例である。
【図１０】エッジ強度が所定の閾値以上のエッジだけを抽出した場合の画像データの一例である。
【図１１】細線化処理を施した場合の画像データの一例である。
【図１２】所定の閾値以上に長いエッジのみを抽出するした場合の画像データの一例である。
【図１３】一定方向に向いたエッジを抽出した場合の画像データの一例である。
【図１４】本発明の概念図である。
【図１５】本発明の被害域自動抽出システムを用いた実例である。
【図１６】被害域候補の概念を示す概念図である。
【図１７】非被害域候補の概念を示す概念図である。
【符号の説明】
【０１０３】
１：被害域自動抽出システム
２：画像データ読込手段
３：被害域候補判定手段
４：非被害域候補判定手段
５：第１誤認候補判定手段
６：第２誤認候補判定手段
７：被害域判定手段
８：画像データ出力手段
【図１】
2005284539.tif 000003
【図２】
2005284539.tif 000004
【図３】
2005284539.tif 000005
【図４】
2005284539.tif 000006
【図５】
2005284539.tif 000007
【図６】
2005284539.tif 000008
【図７】
2005284539.tif 000009
【図８】
2005284539.tif 000010
【図９】
2005284539.tif 000011
【図１０】
2005284539.tif 000012
【図１１】
2005284539.tif 000013
【図１２】
2005284539.tif 000014
【図１３】
2005284539.tif 000015
【図１４】
2005284539.tif 000016
【図１５】
2005284539.tif 000017
【図１６】
2005284539.tif 000018
【図１７】
2005284539.tif 000019