建物屋根の劣化判定方法
【課題】建物屋根の錆の進行による劣化を精度良く把握し、定量的、効率的に評価する。
【解決手段】撮影時期の異なる２つの画像ＰＡ，ＰＢの差分画像ＰＧを求める。ＰＧの画素値ベクトルＧの大きさに相当する変化強度|Ｇ|に基づいて錆が発生した領域の画素か否かを判定する（Ｓ２５）。|Ｇ|が閾値αＲを超えた画素は錆ピクセルとして抽出する（Ｓ４０）。また、αＳ＜|Ｇ|≦αＲの画素については、さらに、３時期以上の画像を用いて、|Ｇ|の変化傾向や、周辺領域の|Ｇ|＞αＳとなる画素の数の変化傾向を調べ（Ｓ４５，Ｓ５０）、それらが単調増加である場合には当該画素を錆ピクセルとして抽出する（Ｓ４０）。錆ピクセルは、その|Ｇ|の大きさに応じてランク分けして、建物の平面図上にてランク毎に劣化部分が図示される。
【解決手段】撮影時期の異なる２つの画像ＰＡ，ＰＢの差分画像ＰＧを求める。ＰＧの画素値ベクトルＧの大きさに相当する変化強度|Ｇ|に基づいて錆が発生した領域の画素か否かを判定する（Ｓ２５）。|Ｇ|が閾値αＲを超えた画素は錆ピクセルとして抽出する（Ｓ４０）。また、αＳ＜|Ｇ|≦αＲの画素については、さらに、３時期以上の画像を用いて、|Ｇ|の変化傾向や、周辺領域の|Ｇ|＞αＳとなる画素の数の変化傾向を調べ（Ｓ４５，Ｓ５０）、それらが単調増加である場合には当該画素を錆ピクセルとして抽出する（Ｓ４０）。錆ピクセルは、その|Ｇ|の大きさに応じてランク分けして、建物の平面図上にてランク毎に劣化部分が図示される。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、工場、倉庫等の大規模の建物の屋根の劣化を判定する方法に関する。
【背景技術】
【０００２】
工場、倉庫等の大面積の屋根を有する建物では、鋼板やめっき鋼板を屋根材として用いることが多い。これら屋根材は塗装され、これにより錆の発生を抑制可能である。しかし、太陽光線や雨・風などに晒され、塗装も劣化し、それと共に急速に錆が発生し得る。これを放置すると腐食が進行し、その後、大規模な修繕工事が必要となり、費用は塗装作業の何倍もかかってしまう。このようなことが起こる前に、屋根の劣化部分を発見し、塗替え等の適切な処置を行なうことが重要である。
【０００３】
従来、劣化状況の把握は、建物屋根に作業者が上がって目視観察することにより行われている。しかし、当該作業には、高所作業に伴う安全性の確保が必要であり、また専門の技術を要すると共に大面積の検査には人手や時間がかかるという問題点が存在する。そのため、下記特許文献１に示されるように、屋根を撮影した画像を処理して、劣化度を評価することが考えられている。
【特許文献１】特開平７−３１８５１０号公報
【発明の開示】
【発明が解決しようとする課題】
【０００４】
錆が発生している部分等の劣化部分をより精度良く把握し、定量的、効率的に評価することが望まれている。本発明は、屋根を上方より撮影した画像を用いて、上述の要求に応える、建物屋根の劣化判定方法を提供することを目的とする。
【課題を解決するための手段】
【０００５】
本発明に係る建物屋根の劣化判定方法は、判定対象の建物屋根を撮影した異なる時期の複数の屋根画像間での画素値の変化に基づいて、当該建物屋根の劣化部分を判定する方法である。
【０００６】
本発明の好適な態様は、前記画素値の変化量が所定の劣化判定閾値を超えたことに基づいて前記劣化部分を判定する方法である。
【０００７】
本発明の他の好適な態様は、前記画素値の変化量に応じた複数の劣化等級を設定し、前記劣化部分を前記劣化等級に応じて区分する方法である。
【０００８】
また、本発明に係る建物屋根の劣化判定方法は、前記画素値の変化量が前記劣化判定閾値以下であるが所定の予備閾値を超えている劣化部分候補画素について、撮影時期が異なる３つ以上の前記屋根画像間にて、前記画素値が前記建物屋根の劣化進行に応じた変化傾向を有することに基づき、前記劣化部分を構成する劣化部分画素と判定する方法とすることもできる。
【０００９】
また、本発明に係る建物屋根の劣化判定方法は、前記画素値の変化量が前記劣化判定閾値以下であるが所定の予備閾値を超えている劣化部分候補画素について、撮影時期が異なる３つ以上の前記屋根画像間にて、当該劣化部分候補画素の所定周辺領域内に存在し前記劣化部分を構成する劣化部分画素の数が前記建物屋根の劣化進行に応じた変化傾向を有することに基づき、前記劣化部分画素と判定する方法とすることもできる。
【００１０】
本発明に係る建物屋根の劣化判定方法は、例えば、錆の発生部分を前記劣化部分として判定する。
【発明の効果】
【００１１】
塗装のはがれや錆の発生等による劣化は時間と共に進行し、また基本的に不可逆であるという特性を有する。本発明はこの性質を利用することで、劣化部分を自動的に精度よく判定し、定量的な評価を効率的に行うことができる。
【発明を実施するための最良の形態】
【００１２】
以下、本発明の実施の形態（以下実施形態という）について、図面に基づいて説明する。
【００１３】
図１は、本発明の実施形態である建物屋根の劣化判定方法を説明する概略の処理フロー図である。
【００１４】
判定対象の建物の屋根の画像を、例えば、数ヶ月、数年といった時間間隔をおいた複数の時期にて撮影する（Ｓ１０）。撮影は、航空機等の飛行体から行うことができる。航空機からマルチスペクトルで撮影された地上の高分解能画像は、屋根の形状や色等の情報を与え、屋根の錆びた部分の抽出に適している。ただし、撮影された航空画像は、航空機に搭載したカメラのレンズ中心から平面に投影した画像であり、建物等が倒れるといった歪みが生じる。そこで、各時期の航空画像の当該歪みを補正して、地表面を真上から見た正射影画像であるトゥルーオルソ画像に変換する（Ｓ１５）。
【００１５】
次に、２つの時期のトゥルーオルソ画像の差分画像を生成し、変化強度を求める（Ｓ２０）。例えば、直近に撮影された航空画像に基づく画像ＰＡと、画像ＰＡより期間Ｔ１だけ前に撮影された航空画像に基づく画像ＰＢとから、差分画像ＰＧを生成する。画像ＰＡ，ＰＢ，ＰＧそれぞれの互いに対応する任意の画素をＡ，Ｂ，Ｇと表す。画像ＰＡ，ＰＢはマルチスペクトル画像であり、Ａ，Ｂの画素値は、各スペクトル毎の成分からなるベクトルで表すことができる。これを次のように表す。ここで、ｎはマルチスペクトル画像のバンド数であり、Ｔは転置行列を表す。
Ａ≡（ａ１，ａ２，…，ａｎ）Ｔ ………（１）
Ｂ≡（ｂ１，ｂ２，…，ｂｎ）Ｔ ………（２）
【００１６】
例えば、Ａ，Ｂとして、Ｒ（赤色光），Ｇ（緑色光），Ｂ（青色光）の各成分で構成されるマルチスペクトル画像や、さらにＩＲ（赤外光）成分を含んだマルチスペクトル画像を用いることができる。
【００１７】
差分画像の各画素Ｇは、次式で与えられる。
Ｇ≡Ａ−Ｂ＝（ａ１−ｂ１，ａ２−ｂ２，…，ａｎ−ｂｎ）Ｔ ………（３）
【００１８】
このようにして得られた差分画像ＰＧを用いて、錆が生じている部分に対応する画素（錆ピクセル）を抽出する。錆ピクセルは、画像ＰＡ，ＰＢが得られた２時期の間での画素値Ａ，Ｂの変化量に基づいて判定され、具体的には、差分画像ＰＧの画素値ベクトルＧから次式、
|Ｇ|≡｛(ａ１−ｂ１)２＋(ａ２−ｂ２)２＋…＋(ａｎ−ｂｎ)２｝１／２ ………（４）
で算出される変化強度|Ｇ|に基づいて判定される（Ｓ２５）。このとき、予め得られている家屋ポリゴン情報３０を参照して、建物に対応する画素だけを抽出し、判定処理Ｓ２５を行う。これにより、建物以外の地物に対応する画素が錆ピクセルとして抽出されることを防止できる。例えば、２時期の間にて成長した樹木等は大きな変化強度|Ｇ|を有し得るが、家屋ポリゴン情報を参照することで、そのようなものを錆ピクセルとして誤検出することが防止される。
【００１９】
なお、|Ｇ|は、各成分に重み付けをして算出してもよい。例えば、劣化が進むと当初の塗装の色が失われ、錆の色が強まると考えられることから、塗装や錆の色に対応する成分の重みを大きくし、それ以外の成分の重みを小さくすることにより、劣化による変化とそうでない変化との弁別度を向上できる。
【００２０】
錆ピクセルの判定は、各画素毎に|Ｇ|と劣化判定閾値αＲ及び予備閾値αＳ（αＳ＜αＲ）を比較することにより行われる。例えば、劣化判定閾値及び予備閾値は以下のように設定することができる。劣化の有無を判定する閾値αは、例えば錆画像の分析結果より経験的に求めるが、一般に変化強度から錆であるか否かの判定を高精度に行うことは困難である。そこで、閾値αに±Δαの幅の許容範囲を設定する。α＋Δαを顕著な劣化と判断できる下限値とし、この値を超える|Ｇ|を生じる画素を錆ピクセルと断定する。またα−Δαを劣化の進行が無視できる上限値とし、この値以下である|Ｇ|を生じる画素を錆ピクセルでないと断定する。これら下限値・上限値をそれぞれ劣化判定閾値αＲ及び予備閾値αＳとして設定する。なお、この比較に際して、|Ｇ|を画像の撮影時期の間隔Ｔ１で規格化したり、Ｔ１に応じてαＲ，αＳを変化させてもよい。
【００２１】
建物を構成する各画素は、２つの閾値αＲ，αＳとの比較結果に応じてインデックスＦを付与され、当該インデックスに応じて後続の処理が選択される（Ｓ３５）。具体的には、|Ｇ|がαＲより大きい場合には、Ｆ＝１を付与する。一方、|Ｇ|がαＳ以下の場合には、Ｆ＝２を付与する。それらの中間の状態、すなわち、αＳ＜|Ｇ|≦αＲの場合は、Ｆ＝０を付与する。
【００２２】
Ｆ＝１を付与された画素は上述のように錆ピクセルと断定され、錆領域（劣化部分）を構成する画素として抽出される（Ｓ４０）。
【００２３】
Ｆ＝０を付与された画素は、|Ｇ|に基づいては錆ピクセルか否かを判定せずに錆ピクセルの候補である錆候補ピクセルとし、さらに３つ以上の時期に亘る画像が利用できる場合には、これらの画素値を利用して、その変化傾向についての分析を行う。例えば、画像ＰＡ，ＰＢを含め、ｎ時期（ｎ≧３）に撮影された画像ＰＣｉ（ｉは０〜ｎ−１なる自然数）の互いに対応する画素Ｃｉの画素値の間の変化傾向を調べる（Ｓ４５）。例えば、各ＰＣｉの撮影時期ｔｉをｔ０＜ｔ１＜…＜ｔｎ−１とし、最初の画像ＰＣ０を基準として、後続の画像ＰＣ１〜ＰＣ（ｎ−１）それぞれと画像ＰＣ０との差分画像ＰＧ１〜ＰＧ（ｎ−１）を求める。そして各ＰＧｊ（ｊは１〜ｎ−１なる自然数）における互いに対応する画素Ｇｊの変化強度|Ｇｊ|の時期を追った変化を分析する。錆による劣化は不可逆に進行するので、基本的に、|Ｇ１|＜|Ｇ２|＜…＜|Ｇｎ−１|なる変化傾向が存在する場合には、錆による変化であると判断して、当該画素のインデックスＦをＦ＝１に変更する。一方、|Ｇ１|＜|Ｇ２|＜…＜|Ｇｎ−１|なる|Ｇｊ|の単調増加傾向が認められない場合には、Ｆ＝０を維持する。
【００２４】
ちなみに、|Ｇｊ|の単調増加傾向が存在することは、１≦ｋ≦ｎ−２である任意の自然数ｋについて|Ｇｋ|＜|Ｇｋ＋１|が成立することに基づいて判定することができる。ここで、塗装が維持され劣化が抑制されている場合には、|Ｇｋ|は基本的には単調には増加せず、撮影条件等に応じて増減して|Ｇｋ|≧|Ｇｋ＋１|となる場合が起こる。よって、この場合はＦ＝０を維持することとする。
【００２５】
また、錆が進行していても、撮影条件の相違等の影響が勝って、|Ｇｋ|≧|Ｇｋ＋１|となることも考えられる。そのような影響を除去するために、例えば、回帰分析等を行い、その結果、変化強度の単調増加が認められる場合にはＦ＝１とし、有意な増加傾向が認められない場合にはＦ＝０を維持する処理とすることもできる。
【００２６】
以上、変化強度|Ｇｊ|の変化傾向の分析について説明したが、さらに他の変化傾向分析を併用することができる。例えば、注目する画素の周辺領域の変化傾向を考慮することができる（Ｓ５０）。例えば、各画像ＰＣｉにおいて注目画素Ｄの所定周辺領域内に存在する|Ｇ|＞αＳとなる画素、すなわち錆ピクセルまたは錆候補ピクセルに相当する画素の数Ｍｉを求める。そしてＭｉの変化を分析する。処理Ｓ２５での|Ｇ|に基づく判定では錆ピクセルとは断定されなかった画素に関して、Ｍ１＜Ｍ２＜…＜Ｍｎなる変化傾向が存在する場合には、当該画素の周辺領域での劣化の進行により、当該画素も錆により|Ｇ|が変化していると判断して、当該画素のインデックスＦをＦ＝１に変更する。一方、Ｍ１＜Ｍ２＜…＜Ｍｎなる変化傾向が認められない場合には、Ｆ＝０を維持する。
【００２７】
例えば、処理Ｓ２５にてＦ＝０とされた各画素について、上述の処理Ｓ４５の判定及び処理Ｓ５０の判定を行い、両方判定にて錆の進行に対応した変化傾向が認められた場合に、当該画素を錆ピクセルと判断することができる。また、いずれか一方の判定で錆の進行に対応した変化傾向が認められた場合に、当該画素を錆ピクセルと判断してもよい。なお、処理Ｓ４５と処理Ｓ５０は順序を入れ替えても構わない。
【００２８】
時系列の画像に基づく分析処理Ｓ４５，Ｓ５０にてＦ＝１を付与された画素は（Ｓ５５）、錆領域（劣化部分）を構成する画素として抽出される（Ｓ４０）。
【００２９】
錆ピクセルは、例えば|Ｇ|に応じてランク分けされる。例えば、閾値β１，β２，…，βｍ−１（β１＞β２＞…＞βｍ−１）を用いて、|Ｇ|≧β１，β１＞|Ｇ|≧β２，…βｍ−１＞|Ｇ|のｍランクを設定する。そして、例えば、|Ｇ|が大きなランクから順に、ランク値Ｒ＝１〜ｍを定義し、各ランクに属する錆ピクセルに当該ランクのランク値Ｒを付与する（Ｓ６０）。なお、２時期の画像しか利用できない場合には、Ｓ４５、Ｓ５０の各処理は行わず、Ｆ＝０とする。あるいは、Ｓ３５にてαＲまたはαＳの一方を用いてＦ＝１またはＦ＝２と判断してもよい。
【００３０】
以上の処理により画素に対応して求められたＦ値やＲ値はファイルに出力され（Ｓ６５）、当該画素についての処理を完了する。一方、Ｆ＝２を付与された画素や、処理Ｓ４５，Ｓ５０にてＦ＝０を維持された画素は、錆ピクセルとは判定されずに処理を終了する。
【００３１】
各画素について上述の処理が繰り返され（Ｓ７０）、画像を構成する全ての画素について当該処理が完了すると、建物屋根の劣化判定処理を終了する。
【００３２】
例えば、上述の劣化判定処理でファイルに出力された結果を用いて、建物屋根に生じた錆（又は劣化部分）の面積や程度（ランク）を定量的に把握でき、効率的に屋根の点検や補修の計画や作業実行のための資料を作成することができる。本方法によれば、劣化部分が自動的に判定され、熟練者の経験を要せずに客観的な判定が可能となる。例えば、建物屋根の劣化判定結果は、アセットマネジメントに利用するといった用途も考えられる。
【００３３】
上述の実施形態の劣化判定では航空機画像を用いる例を示したが、これに限られるものではなく、高分解能衛星画像を用いて同様の処理を行うこともできる。
【図面の簡単な説明】
【００３４】
【図１】本発明の実施形態である建物屋根の劣化判定方法を説明する概略の処理フロー図である。
【技術分野】
【０００１】
本発明は、工場、倉庫等の大規模の建物の屋根の劣化を判定する方法に関する。
【背景技術】
【０００２】
工場、倉庫等の大面積の屋根を有する建物では、鋼板やめっき鋼板を屋根材として用いることが多い。これら屋根材は塗装され、これにより錆の発生を抑制可能である。しかし、太陽光線や雨・風などに晒され、塗装も劣化し、それと共に急速に錆が発生し得る。これを放置すると腐食が進行し、その後、大規模な修繕工事が必要となり、費用は塗装作業の何倍もかかってしまう。このようなことが起こる前に、屋根の劣化部分を発見し、塗替え等の適切な処置を行なうことが重要である。
【０００３】
従来、劣化状況の把握は、建物屋根に作業者が上がって目視観察することにより行われている。しかし、当該作業には、高所作業に伴う安全性の確保が必要であり、また専門の技術を要すると共に大面積の検査には人手や時間がかかるという問題点が存在する。そのため、下記特許文献１に示されるように、屋根を撮影した画像を処理して、劣化度を評価することが考えられている。
【特許文献１】特開平７−３１８５１０号公報
【発明の開示】
【発明が解決しようとする課題】
【０００４】
錆が発生している部分等の劣化部分をより精度良く把握し、定量的、効率的に評価することが望まれている。本発明は、屋根を上方より撮影した画像を用いて、上述の要求に応える、建物屋根の劣化判定方法を提供することを目的とする。
【課題を解決するための手段】
【０００５】
本発明に係る建物屋根の劣化判定方法は、判定対象の建物屋根を撮影した異なる時期の複数の屋根画像間での画素値の変化に基づいて、当該建物屋根の劣化部分を判定する方法である。
【０００６】
本発明の好適な態様は、前記画素値の変化量が所定の劣化判定閾値を超えたことに基づいて前記劣化部分を判定する方法である。
【０００７】
本発明の他の好適な態様は、前記画素値の変化量に応じた複数の劣化等級を設定し、前記劣化部分を前記劣化等級に応じて区分する方法である。
【０００８】
また、本発明に係る建物屋根の劣化判定方法は、前記画素値の変化量が前記劣化判定閾値以下であるが所定の予備閾値を超えている劣化部分候補画素について、撮影時期が異なる３つ以上の前記屋根画像間にて、前記画素値が前記建物屋根の劣化進行に応じた変化傾向を有することに基づき、前記劣化部分を構成する劣化部分画素と判定する方法とすることもできる。
【０００９】
また、本発明に係る建物屋根の劣化判定方法は、前記画素値の変化量が前記劣化判定閾値以下であるが所定の予備閾値を超えている劣化部分候補画素について、撮影時期が異なる３つ以上の前記屋根画像間にて、当該劣化部分候補画素の所定周辺領域内に存在し前記劣化部分を構成する劣化部分画素の数が前記建物屋根の劣化進行に応じた変化傾向を有することに基づき、前記劣化部分画素と判定する方法とすることもできる。
【００１０】
本発明に係る建物屋根の劣化判定方法は、例えば、錆の発生部分を前記劣化部分として判定する。
【発明の効果】
【００１１】
塗装のはがれや錆の発生等による劣化は時間と共に進行し、また基本的に不可逆であるという特性を有する。本発明はこの性質を利用することで、劣化部分を自動的に精度よく判定し、定量的な評価を効率的に行うことができる。
【発明を実施するための最良の形態】
【００１２】
以下、本発明の実施の形態（以下実施形態という）について、図面に基づいて説明する。
【００１３】
図１は、本発明の実施形態である建物屋根の劣化判定方法を説明する概略の処理フロー図である。
【００１４】
判定対象の建物の屋根の画像を、例えば、数ヶ月、数年といった時間間隔をおいた複数の時期にて撮影する（Ｓ１０）。撮影は、航空機等の飛行体から行うことができる。航空機からマルチスペクトルで撮影された地上の高分解能画像は、屋根の形状や色等の情報を与え、屋根の錆びた部分の抽出に適している。ただし、撮影された航空画像は、航空機に搭載したカメラのレンズ中心から平面に投影した画像であり、建物等が倒れるといった歪みが生じる。そこで、各時期の航空画像の当該歪みを補正して、地表面を真上から見た正射影画像であるトゥルーオルソ画像に変換する（Ｓ１５）。
【００１５】
次に、２つの時期のトゥルーオルソ画像の差分画像を生成し、変化強度を求める（Ｓ２０）。例えば、直近に撮影された航空画像に基づく画像ＰＡと、画像ＰＡより期間Ｔ１だけ前に撮影された航空画像に基づく画像ＰＢとから、差分画像ＰＧを生成する。画像ＰＡ，ＰＢ，ＰＧそれぞれの互いに対応する任意の画素をＡ，Ｂ，Ｇと表す。画像ＰＡ，ＰＢはマルチスペクトル画像であり、Ａ，Ｂの画素値は、各スペクトル毎の成分からなるベクトルで表すことができる。これを次のように表す。ここで、ｎはマルチスペクトル画像のバンド数であり、Ｔは転置行列を表す。
Ａ≡（ａ１，ａ２，…，ａｎ）Ｔ ………（１）
Ｂ≡（ｂ１，ｂ２，…，ｂｎ）Ｔ ………（２）
【００１６】
例えば、Ａ，Ｂとして、Ｒ（赤色光），Ｇ（緑色光），Ｂ（青色光）の各成分で構成されるマルチスペクトル画像や、さらにＩＲ（赤外光）成分を含んだマルチスペクトル画像を用いることができる。
【００１７】
差分画像の各画素Ｇは、次式で与えられる。
Ｇ≡Ａ−Ｂ＝（ａ１−ｂ１，ａ２−ｂ２，…，ａｎ−ｂｎ）Ｔ ………（３）
【００１８】
このようにして得られた差分画像ＰＧを用いて、錆が生じている部分に対応する画素（錆ピクセル）を抽出する。錆ピクセルは、画像ＰＡ，ＰＢが得られた２時期の間での画素値Ａ，Ｂの変化量に基づいて判定され、具体的には、差分画像ＰＧの画素値ベクトルＧから次式、
|Ｇ|≡｛(ａ１−ｂ１)２＋(ａ２−ｂ２)２＋…＋(ａｎ−ｂｎ)２｝１／２ ………（４）
で算出される変化強度|Ｇ|に基づいて判定される（Ｓ２５）。このとき、予め得られている家屋ポリゴン情報３０を参照して、建物に対応する画素だけを抽出し、判定処理Ｓ２５を行う。これにより、建物以外の地物に対応する画素が錆ピクセルとして抽出されることを防止できる。例えば、２時期の間にて成長した樹木等は大きな変化強度|Ｇ|を有し得るが、家屋ポリゴン情報を参照することで、そのようなものを錆ピクセルとして誤検出することが防止される。
【００１９】
なお、|Ｇ|は、各成分に重み付けをして算出してもよい。例えば、劣化が進むと当初の塗装の色が失われ、錆の色が強まると考えられることから、塗装や錆の色に対応する成分の重みを大きくし、それ以外の成分の重みを小さくすることにより、劣化による変化とそうでない変化との弁別度を向上できる。
【００２０】
錆ピクセルの判定は、各画素毎に|Ｇ|と劣化判定閾値αＲ及び予備閾値αＳ（αＳ＜αＲ）を比較することにより行われる。例えば、劣化判定閾値及び予備閾値は以下のように設定することができる。劣化の有無を判定する閾値αは、例えば錆画像の分析結果より経験的に求めるが、一般に変化強度から錆であるか否かの判定を高精度に行うことは困難である。そこで、閾値αに±Δαの幅の許容範囲を設定する。α＋Δαを顕著な劣化と判断できる下限値とし、この値を超える|Ｇ|を生じる画素を錆ピクセルと断定する。またα−Δαを劣化の進行が無視できる上限値とし、この値以下である|Ｇ|を生じる画素を錆ピクセルでないと断定する。これら下限値・上限値をそれぞれ劣化判定閾値αＲ及び予備閾値αＳとして設定する。なお、この比較に際して、|Ｇ|を画像の撮影時期の間隔Ｔ１で規格化したり、Ｔ１に応じてαＲ，αＳを変化させてもよい。
【００２１】
建物を構成する各画素は、２つの閾値αＲ，αＳとの比較結果に応じてインデックスＦを付与され、当該インデックスに応じて後続の処理が選択される（Ｓ３５）。具体的には、|Ｇ|がαＲより大きい場合には、Ｆ＝１を付与する。一方、|Ｇ|がαＳ以下の場合には、Ｆ＝２を付与する。それらの中間の状態、すなわち、αＳ＜|Ｇ|≦αＲの場合は、Ｆ＝０を付与する。
【００２２】
Ｆ＝１を付与された画素は上述のように錆ピクセルと断定され、錆領域（劣化部分）を構成する画素として抽出される（Ｓ４０）。
【００２３】
Ｆ＝０を付与された画素は、|Ｇ|に基づいては錆ピクセルか否かを判定せずに錆ピクセルの候補である錆候補ピクセルとし、さらに３つ以上の時期に亘る画像が利用できる場合には、これらの画素値を利用して、その変化傾向についての分析を行う。例えば、画像ＰＡ，ＰＢを含め、ｎ時期（ｎ≧３）に撮影された画像ＰＣｉ（ｉは０〜ｎ−１なる自然数）の互いに対応する画素Ｃｉの画素値の間の変化傾向を調べる（Ｓ４５）。例えば、各ＰＣｉの撮影時期ｔｉをｔ０＜ｔ１＜…＜ｔｎ−１とし、最初の画像ＰＣ０を基準として、後続の画像ＰＣ１〜ＰＣ（ｎ−１）それぞれと画像ＰＣ０との差分画像ＰＧ１〜ＰＧ（ｎ−１）を求める。そして各ＰＧｊ（ｊは１〜ｎ−１なる自然数）における互いに対応する画素Ｇｊの変化強度|Ｇｊ|の時期を追った変化を分析する。錆による劣化は不可逆に進行するので、基本的に、|Ｇ１|＜|Ｇ２|＜…＜|Ｇｎ−１|なる変化傾向が存在する場合には、錆による変化であると判断して、当該画素のインデックスＦをＦ＝１に変更する。一方、|Ｇ１|＜|Ｇ２|＜…＜|Ｇｎ−１|なる|Ｇｊ|の単調増加傾向が認められない場合には、Ｆ＝０を維持する。
【００２４】
ちなみに、|Ｇｊ|の単調増加傾向が存在することは、１≦ｋ≦ｎ−２である任意の自然数ｋについて|Ｇｋ|＜|Ｇｋ＋１|が成立することに基づいて判定することができる。ここで、塗装が維持され劣化が抑制されている場合には、|Ｇｋ|は基本的には単調には増加せず、撮影条件等に応じて増減して|Ｇｋ|≧|Ｇｋ＋１|となる場合が起こる。よって、この場合はＦ＝０を維持することとする。
【００２５】
また、錆が進行していても、撮影条件の相違等の影響が勝って、|Ｇｋ|≧|Ｇｋ＋１|となることも考えられる。そのような影響を除去するために、例えば、回帰分析等を行い、その結果、変化強度の単調増加が認められる場合にはＦ＝１とし、有意な増加傾向が認められない場合にはＦ＝０を維持する処理とすることもできる。
【００２６】
以上、変化強度|Ｇｊ|の変化傾向の分析について説明したが、さらに他の変化傾向分析を併用することができる。例えば、注目する画素の周辺領域の変化傾向を考慮することができる（Ｓ５０）。例えば、各画像ＰＣｉにおいて注目画素Ｄの所定周辺領域内に存在する|Ｇ|＞αＳとなる画素、すなわち錆ピクセルまたは錆候補ピクセルに相当する画素の数Ｍｉを求める。そしてＭｉの変化を分析する。処理Ｓ２５での|Ｇ|に基づく判定では錆ピクセルとは断定されなかった画素に関して、Ｍ１＜Ｍ２＜…＜Ｍｎなる変化傾向が存在する場合には、当該画素の周辺領域での劣化の進行により、当該画素も錆により|Ｇ|が変化していると判断して、当該画素のインデックスＦをＦ＝１に変更する。一方、Ｍ１＜Ｍ２＜…＜Ｍｎなる変化傾向が認められない場合には、Ｆ＝０を維持する。
【００２７】
例えば、処理Ｓ２５にてＦ＝０とされた各画素について、上述の処理Ｓ４５の判定及び処理Ｓ５０の判定を行い、両方判定にて錆の進行に対応した変化傾向が認められた場合に、当該画素を錆ピクセルと判断することができる。また、いずれか一方の判定で錆の進行に対応した変化傾向が認められた場合に、当該画素を錆ピクセルと判断してもよい。なお、処理Ｓ４５と処理Ｓ５０は順序を入れ替えても構わない。
【００２８】
時系列の画像に基づく分析処理Ｓ４５，Ｓ５０にてＦ＝１を付与された画素は（Ｓ５５）、錆領域（劣化部分）を構成する画素として抽出される（Ｓ４０）。
【００２９】
錆ピクセルは、例えば|Ｇ|に応じてランク分けされる。例えば、閾値β１，β２，…，βｍ−１（β１＞β２＞…＞βｍ−１）を用いて、|Ｇ|≧β１，β１＞|Ｇ|≧β２，…βｍ−１＞|Ｇ|のｍランクを設定する。そして、例えば、|Ｇ|が大きなランクから順に、ランク値Ｒ＝１〜ｍを定義し、各ランクに属する錆ピクセルに当該ランクのランク値Ｒを付与する（Ｓ６０）。なお、２時期の画像しか利用できない場合には、Ｓ４５、Ｓ５０の各処理は行わず、Ｆ＝０とする。あるいは、Ｓ３５にてαＲまたはαＳの一方を用いてＦ＝１またはＦ＝２と判断してもよい。
【００３０】
以上の処理により画素に対応して求められたＦ値やＲ値はファイルに出力され（Ｓ６５）、当該画素についての処理を完了する。一方、Ｆ＝２を付与された画素や、処理Ｓ４５，Ｓ５０にてＦ＝０を維持された画素は、錆ピクセルとは判定されずに処理を終了する。
【００３１】
各画素について上述の処理が繰り返され（Ｓ７０）、画像を構成する全ての画素について当該処理が完了すると、建物屋根の劣化判定処理を終了する。
【００３２】
例えば、上述の劣化判定処理でファイルに出力された結果を用いて、建物屋根に生じた錆（又は劣化部分）の面積や程度（ランク）を定量的に把握でき、効率的に屋根の点検や補修の計画や作業実行のための資料を作成することができる。本方法によれば、劣化部分が自動的に判定され、熟練者の経験を要せずに客観的な判定が可能となる。例えば、建物屋根の劣化判定結果は、アセットマネジメントに利用するといった用途も考えられる。
【００３３】
上述の実施形態の劣化判定では航空機画像を用いる例を示したが、これに限られるものではなく、高分解能衛星画像を用いて同様の処理を行うこともできる。
【図面の簡単な説明】
【００３４】
【図１】本発明の実施形態である建物屋根の劣化判定方法を説明する概略の処理フロー図である。
【特許請求の範囲】
【請求項１】
判定対象の建物屋根を撮影した異なる時期の複数の屋根画像間での画素値の変化に基づいて、当該建物屋根の劣化部分を判定することを特徴とする建物屋根の劣化判定方法。
【請求項２】
前記画素値の変化量が所定の劣化判定閾値を超えたことに基づいて前記劣化部分を判定すること、を特徴とする請求項１に記載の建物屋根の劣化判定方法。
【請求項３】
前記画素値の変化量に応じた複数の劣化等級を設定し、
前記劣化部分を前記劣化等級に応じて区分すること、
を特徴とする請求項２に記載の建物屋根の劣化判定方法。
【請求項４】
前記画素値の変化量が前記劣化判定閾値以下であるが所定の予備閾値を超えている劣化部分候補画素について、撮影時期が異なる３つ以上の前記屋根画像間にて、前記画素値が前記建物屋根の劣化進行に応じた変化傾向を有することに基づき、前記劣化部分を構成する劣化部分画素と判定すること、を特徴とする請求項２又は請求項３に記載の建物屋根の劣化判定方法。
【請求項５】
前記画素値の変化量が前記劣化判定閾値以下であるが所定の予備閾値を超えている劣化部分候補画素について、撮影時期が異なる３つ以上の前記屋根画像間にて、当該劣化部分候補画素の所定周辺領域内に存在し前記劣化部分を構成する劣化部分画素の数が前記建物屋根の劣化進行に応じた変化傾向を有することに基づき、前記劣化部分画素と判定すること、を特徴とする請求項２から請求項４のいずれか１つに記載の建物屋根の劣化判定方法。
【請求項６】
錆の発生部分を前記劣化部分として判定することを特徴とする請求項１から請求項５のいずれか１つに記載の建物屋根の劣化判定方法。
【請求項１】
判定対象の建物屋根を撮影した異なる時期の複数の屋根画像間での画素値の変化に基づいて、当該建物屋根の劣化部分を判定することを特徴とする建物屋根の劣化判定方法。
【請求項２】
前記画素値の変化量が所定の劣化判定閾値を超えたことに基づいて前記劣化部分を判定すること、を特徴とする請求項１に記載の建物屋根の劣化判定方法。
【請求項３】
前記画素値の変化量に応じた複数の劣化等級を設定し、
前記劣化部分を前記劣化等級に応じて区分すること、
を特徴とする請求項２に記載の建物屋根の劣化判定方法。
【請求項４】
前記画素値の変化量が前記劣化判定閾値以下であるが所定の予備閾値を超えている劣化部分候補画素について、撮影時期が異なる３つ以上の前記屋根画像間にて、前記画素値が前記建物屋根の劣化進行に応じた変化傾向を有することに基づき、前記劣化部分を構成する劣化部分画素と判定すること、を特徴とする請求項２又は請求項３に記載の建物屋根の劣化判定方法。
【請求項５】
前記画素値の変化量が前記劣化判定閾値以下であるが所定の予備閾値を超えている劣化部分候補画素について、撮影時期が異なる３つ以上の前記屋根画像間にて、当該劣化部分候補画素の所定周辺領域内に存在し前記劣化部分を構成する劣化部分画素の数が前記建物屋根の劣化進行に応じた変化傾向を有することに基づき、前記劣化部分画素と判定すること、を特徴とする請求項２から請求項４のいずれか１つに記載の建物屋根の劣化判定方法。
【請求項６】
錆の発生部分を前記劣化部分として判定することを特徴とする請求項１から請求項５のいずれか１つに記載の建物屋根の劣化判定方法。
【図１】
【公開番号】特開２００８−２９２２６３（Ｐ２００８−２９２２６３Ａ）
【公開日】平成２０年１２月４日（２００８．１２．４）
【国際特許分類】
【出願番号】特願２００７−１３７３６６（Ｐ２００７−１３７３６６）
【出願日】平成１９年５月２４日（２００７．５．２４）
【出願人】（０００１３５７７１）株式会社パスコ (102)
【Ｆターム（参考）】
【公開日】平成２０年１２月４日（２００８．１２．４）
【国際特許分類】
【出願日】平成１９年５月２４日（２００７．５．２４）
【出願人】（０００１３５７７１）株式会社パスコ (102)
【Ｆターム（参考）】
[ Back to top ]