JP5384429B2 - コンクリート構造物画像のひび割れ検知装置、ひび割れ検知方法及びそのプログラム - Google Patents

Description

本発明は、画像を用いた構造物診断法に関するものであり、特に、建物や電柱等のコンクリート構造物を撮像した画像から、コンクリート構造物の表面に生じたひび割れを画像処理によって自動検知するひび割れ検知装置、ひび割れ検知方法及びそのプログラムに関するものである。

戦後の高度成長期に建設された高速道路などの老朽化したコンクリート構造物に対する不安が社会問題化している。コンクリート構造物の耐久性を検査する際の最重要項目のひとつとして、コンクリート表面のひび割れ検査が挙げられる。ひび割れは、地震等によって構造物に応力が加わるとコンクリートが破断して生じる。コンクリートの表面でひび割れが発生すると、そこから雨水等が入り込み、内部の鉄筋に到達すると腐食が始まり、鉄筋の強度が落ちて構造物の耐久性低下を引き起こす。このような構造物の耐久性低下を未然に防ぐには、表面ひび割れの早期発見と補修の実施が重要となる。

補修が必要となるひび割れの程度を示すには、コンクリート表面に現れたひび割れの「幅」を指標として用いる場合が多い。鉄筋への到達と直接的に関わる「深さ」も重要な指標であるが、ひび割れ幅から深さをある程度予測できるため、表面のひび割れ幅を用いて補修実施の要否を判断するのが一般的となっている。さて、構造物の耐久性に関する許容ひび割れ幅に関しては、各国の提案基準によって若干のばらつきがある。おおよそ、乾燥空気中であれば幅０．３〜０．４ｍｍ、湿空中であればそれ以下の値を許容ひび割れ幅と規定しているものが多く、幅が０．２ｍｍ以上に達したひび割れについては、充填材や被覆による補修、あるいは鋼製アンカー等による補強を実施する対象としている（例えば、非特許文献１参照）。

このような理由によって現在行われているひび割れの診断方法は目視が主流である。しかしながら、熟練した作業員であっても長時間の作業においてサブミリメートル幅のひび割れを見落とさない集中力を維持することは困難である。ゆえに、機械によるひび割れ検査の自動化が望まれているが、機材の入手し易さや現場作業の簡易さという観点から、デジタルカメラによって構造物の壁面を遠隔撮影し、その画像から画像処理ソフトウエアによってひび割れの自動検知を試みるというアプローチが有力視されている。

さて、コンクリート壁面の画像からひび割れの自動検知を行うアルゴリズムは既にいくつか提案されているが、これら既存のアルゴリズムは次の理由から実用化課題を抱えている。

コンクリート構造物の表面には、ひび割れ以外にも汚れや傷等が多数ある。ゆえに、画像処理によってひび割れ検知を試みる場合、ノイズを含んだ入力画像からひび割れパターンを認識する処理が必要となる。これらのパターン認識を行うには、一般的にニューラルネットワークが最も有力なアルゴリズムのひとつであると知られている。しかしながら、従来のニューラルネットワークによるひび割れ検知アルゴリズムには、以下に挙げる二つの問題点が指摘されている（例えば、非特許文献２参照）。

一つの問題点は、ひび割れ抽出精度の低さによる劣化見落としの危険性である。非特許文献２の方法によれば、一次元のエッジフィルタによって抽出された特徴量からひび割れの探索を行っているが、この方法では、汚れ、傷、或いはコンクリート表面に浮き出ている骨材等をひび割れとして誤検知する場合がしばしば生じる。なぜならば、一次元のエッジ情報には二次元の空間情報が含まれておらず、局在点として現れる汚れ、傷、或いは骨材と、線分として現れるひび割れの間にて特徴量の差が少ないためである。それゆえ、一次元のエッジフィルタによる特徴量を用いてニューラルネットワークでひび割れ収束を行った場合、ひび割れ以外に収束することがある。

もう一つの問題点は、ひび割れの自動判定部が無いことによる作業時間の長さである。ひび割れ自動判定部が無い現状のアルゴリズムにおいては、ひび割れ収束結果を人間が眼で見て判断しないとならないため、大量にひび割れ画像を処理する場合において膨大な作業時間が必要となる。

以上のように、ひび割れの抽出精度が低く、かつ自動判定部が無い従来アルゴリズムでは、現場運用に供することができない。

この問題点を克服するには、ひび割れの抽出精度の向上と、ひび割れ自動判定機能の追加が必要となるが、ひび割れの形状、ひび割れ以外の汚れ、傷、骨材等の形状が千差万別であるゆえ、例えばテンプレートマッチングのような既存のパターン認識アルゴリズムを用いて簡易にひび割れ抽出精度を向上させることは容易ではない。

また、実際のコンクリート構造物に生じるひび割れは直線とは限らず湾曲している場合があるため、例えばハフ変換のような既存の直線検出アルゴリズムを用いて簡易に自動判定することは容易ではない（例えば、非特許文献３参照）。

以上のような理由から、コンクリート構造物の画像からひび割れを高精度で抽出し自動判定するシステムが望まれている反面、これまでの技術では実現困難である。

「コンクリートのひび割れ調査、補修・補強指針」、社団法人日本コンクリート工学協会、２００３年６月２０日発行、ｐ８５−９４ S. Oka, S. Mochizuki, H. Togo, and N. Kukutsu, "A Neural Network Algorithm for Detecting Invisible Concrete Surface Cracks in Near-field Millimeter-wave Images," Proc. of IEEE International Conference on Systems, Man, and Cybernetics, Neural Networks and Their Applications II, W2, Oct. 11-14, 2009, in Texas, p3803-3804 酒井幸市、「Ｖｉｓｕａｌ Ｂａｓｉｃ＆Ｖｉｓｕａｌ Ｃ＋＋によるデジタル画像処理入門」、ＣＱ出版株式会社、２００８年２月１日第８版発行、ｐ７２−７３

上述したように、コンクリート構造物の画像からひび割れを高精度で抽出し自動判定するシステムが望まれている反面、これまでの技術では実現困難であった。

本発明の目的は、上記の課題に鑑みてなされたものであり、建物や電柱等のコンクリート構造物を撮像した画像から、コンクリート構造物の表面に生じたひび割れを画像処理によって、ひび割れの抽出精度を向上させて自動検知するコンクリート構造物画像のひび割れ検知装置、ひび割れ検知方法及びそのプログラムを提供することにある。

本発明のひび割れ検知装置は、コンクリート構造物を撮像した画像から、コンクリート構造物の表面に生じたひび割れを画像処理によって自動検知するひび割れ検知装置であって、コンクリート構造物を撮像した画像を入力してコンクリート構造物の表面に生じたひび割れの候補画像を生成するひび割れ探査部と、該ひび割れの候補画像からコンクリート構造物の表面に生じたひび割れを抽出するひび割れ判定部とを備え、前記ひび割れ探査部は、前記コンクリート構造物を撮像した画像に対して方向及び幅の異なる複数種類の２次元エッジフィルタを個別に適用し、画素ごとに、複数種類のエッジフィルタ適用後の画素値をそれぞれ比較し、最大の画素値をひび割れ特徴量として選定するひび割れ特徴量抽出部と、前記ひび割れ特徴量からなる画素値を有する処理対象画像に対してニューラルネットの動作式を用いてニューラル収束演算を実行し、ひび割れ候補画像を生成するひび割れ収束部とを備え、前記ひび割れ判定部は、前記ひび割れ候補画像に対して、前記ひび割れ収束部で発火したニューロンの対応画素をひび割れ候補画素とし、前記ひび割れ候補画素中の判定対象画素に対して所定の近傍範囲内にある画素についてひび割れと判定する近傍連結判定部と、前記ひび割れ候補画素の特徴量を合計した値に基づいてひび割れを判定する特徴量判定部と、前記ひび割れ候補画素中の前記判定対象画素に対して隣接する画素のユークリッド距離に基づいてひび割れと判定する飛び地判定部と、前記近傍連結判定部、前記特徴量判定部及び前記飛び地判定部にて判定された結果の論理積を行う論理積演算部とを備えることを特徴とする。

また、本発明のひび割れ検知装置において、前記ひび割れ収束部は、２次元のニューロン配列を有するＨｏｐｆｉｅｌｄ型ニューラルネットワークとして構成されていることを特徴とする。

また、本発明のひび割れ検知装置は、コンクリート構造物を撮像した画像から、コンクリート構造物の表面に生じたひび割れを画像処理によって自動検知するひび割れ検知装置であって、コンクリート構造物を撮像した画像を入力してコンクリート構造物の表面に生じたひび割れの候補画像を生成するひび割れ探査部と、該ひび割れの候補画像からコンクリート構造物の表面に生じたひび割れを抽出するひび割れ判定部とを備え、前記ひび割れ探査部は、前記コンクリート構造物を撮像した画像に対して方向及び幅の異なる複数種類の２次元エッジフィルタを個別に適用し、画素ごとに、複数種類のエッジフィルタ適用後の画素値をそれぞれ比較し、最大の画素値をひび割れ特徴量として選定するひび割れ特徴量抽出部と、前記ひび割れ特徴量からなる画素値を有する処理対象画像に対してニューラルネットの動作式を用いてニューラル収束演算を実行し、ひび割れ候補画像を生成するひび割れ収束部とを備え、前記ひび割れ特徴量抽出部は、前記２次元エッジフィルタの各対応画素についてｍ種類のエッジ方向とｎ種類のエッジ幅の組み合わせを持つｍ×ｎ個の２次元エッジフィルタ処理を行い、そのｍ×ｎ個の出力値の中から最大値となるものをひび割れ特徴量ξとして選択し、前記対応画素から一定範囲に存在し、識別情報がひび割れを撮像した領域に該当することを示した周辺ニューロンに対応する周辺画素の位置と、前記対応画素の位置との間のユークリッド距離に基づいて、前記周辺画素に対する前記対応画素の連結度ηを算出し、前記ひび割れ特徴量ξと前記連結度ηの各々に動的係数を乗算したものを線形結合して前記ひび割れ収束部（Ｈｏｐｆｉｅｌｄ型ニューラルネットワーク）の入力値とすることによって、前記ひび割れ収束部における前記２次元ニューロン配列の発火パターンをひび割れパターンへと収束させる手段を有することを特徴とする。

また、本発明のひび割れ検知装置において、前記近傍連結判定部は、前記ひび割れ候補画像を短冊状に分割したセグメント内で、前記ひび割れ候補画素中の判定対象画素に対して所定の近傍範囲内にある前記ひび割れ候補画素以外の各画素とのユークリッド距離を合計した値が一定値以下の場合に近傍連結カウント値ｃｎｔとしてカウントし、前記近傍連結カウント値ｃｎｔが一定値以上の場合にひび割れセグメントと判定し、さらに前記ひび割れセグメントの数が一定数以上の場合にひび割れと判定する手段を有することを特徴とする。

また、本発明のひび割れ検知装置において、前記特徴量判定部は、前記ひび割れ候補画素の特徴量ξを合計した値が当該ひび割れ候補画像全体の平均輝度から算出した閾値θを超えた場合にひび割れと判定する手段を有することを特徴とする。

また、本発明のひび割れ検知装置において、前記飛び地判定部は、前記ひび割れ候補画素中の前記判定対象画素と前記判定対象画素に隣接した画素とのユークリッド距離が一定値以上であり、かつ前記判定対象画素の前後の所定数画素において隣接する画素間のユークリッド距離が全て一定値内である場合、飛び地でありひび割れではないと判定する手段を有することを特徴とする。

さらに、本発明のひび割れ検知方法は、ひび割れ検知装置によってコンクリート構造物を撮像した画像から、コンクリート構造物の表面に生じたひび割れを画像処理によって自動検知するひび割れ検知方法であって、前記ひび割れ検知装置の処理手順は、（Ａ）前記コンクリート構造物を撮像した画像に対して方向及び幅の異なる複数種類の２次元エッジフィルタを個別に適用し、画素ごとに、複数種類のエッジフィルタ適用後の画素値をそれぞれ比較し、最大の画素値をひび割れ特徴量として選定するステップと、（Ｂ）前記ひび割れ特徴量からなる画素値を有する処理対象画像に対してニューラルネットの動作式を用いてニューラル収束演算を実行し、ひび割れ候補画像を生成するステップと、（Ｃ）前記ひび割れ候補画像に対して、前記ステップにて発火したニューロンの対応画素をひび割れ候補画素とし、前記ひび割れ候補画素中の判定対象画素に対して所定の近傍範囲内にある画素についてひび割れと判定するステップと、（Ｄ）前記ひび割れ候補画素の特徴量を合計した値に基づいてひび割れを判定するステップと、（Ｅ）前記ひび割れ候補画素中の前記判定対象画素に対して隣接する画素のユークリッド距離に基づいてひび割れと判定するステップと、前記ステップ（Ｃ）、前記ステップ（Ｄ）及び前記ステップ（Ｅ）にて判定された結果の論理積を行うステップと、を含むことを特徴とする。

さらに、本発明のプログラムは、コンクリート構造物を撮像した画像から、コンクリート構造物の表面に生じたひび割れを画像処理によって自動検知するひび割れ検知装置として構成するコンピュータに、（Ａ）前記コンクリート構造物を撮像した画像に対して方向及び幅の異なる複数種類の２次元エッジフィルタを個別に適用し、画素ごとに、複数種類のエッジフィルタ適用後の画素値をそれぞれ比較し、最大の画素値をひび割れ特徴量として選定するステップと、（Ｂ）前記ひび割れ特徴量からなる画素値を有する処理対象画像に対してニューラルネットの動作式を用いてニューラル収束演算を実行し、ひび割れ候補画像を生成するステップと、（Ｃ）前記ひび割れ候補画像に対して、前記ステップにて発火したニューロンの対応画素をひび割れ候補画素とし、前記ひび割れ候補画素中の判定対象画素に対して所定の近傍範囲内にある画素についてひび割れと判定するステップと、（Ｄ）前記ひび割れ候補画素の特徴量を合計した値に基づいてひび割れを判定するステップと、（Ｅ）前記ひび割れ候補画素中の前記判定対象画素に対して隣接する画素のユークリッド距離に基づいてひび割れと判定するステップと、前記ステップ（Ｃ）、前記ステップ（Ｄ）及び前記ステップ（Ｅ）にて判定された結果の論理積を行うステップと、を実行させるためのプログラムとして構成される。

本発明によれば、コンクリート構造物の画像から精度良くひび割れを抽出し、かつ、自動判定することが可能となる。

本発明による一実施例のひび割れ検知装置のブロック図である。 本発明による一実施例のひび割れ検知装置の動作説明図である。 従来のニューラルネットワークを用いたひび割れ検知ソフトによる実行結果を示す図である。 本発明による一実施例のひび割れ検知装置におけるひび割れ特徴量抽出部の処理の詳細を説明する図である。 本発明による一実施例のひび割れ検知装置におけるひび割れ特徴量抽出部の処理の詳細をさらに説明する図である。 本発明による一実施例のひび割れ検知装置におけるひび割れ特徴量抽出部の特徴量抽出例を示す図である。 本発明による一実施例のひび割れ検知装置におけるひび割れ収束部の処理の詳細を説明する図である。 本発明による一実施例のひび割れ検知装置におけるひび割れ収束部の収束速度を高めるための処理を説明する図である。 本発明による一実施例のひび割れ検知装置におけるひび割れ判定部のデジタル画像上に現れうるひび割れ形状の定義を模式的に表した図である。 本発明による一実施例のひび割れ検知装置における近傍連結判定部の動作を示すフローチャートである。 本発明による一実施例のひび割れ検知装置における近傍連結判定部の計算例を説明する図である。 本発明による一実施例のひび割れ検知装置における特徴量判定部の動作を示すフローチャートである。 本発明による一実施例のひび割れ検知装置における飛び地判定部における処理の詳細を説明する図である。 本発明による一実施例のひび割れ検知装置における飛び地判定部の動作を示すフローチャートである。 本発明による一実施例のひび割れ検知装置における飛び地判定部の処理例を示す図である。 本発明による一実施例のひび割れ検知装置におけるひび割れ検知結果例を示す図である。

以下、本発明による一実施例のひび割れ検知装置について説明する。本発明に係るひび割れ検知方法及びプログラムは、本実施例のひび割れ検知装置の説明から明らかになる。

〔装置構成〕
図１は、本発明による一実施例のひび割れ検知装置のブロック図である。本実施例のひび割れ検知装置１００は、制御部１０及び記憶部４０を備えており、コンピュータとして構成できるものである。制御部１０の各機能を実現する処理内容を記述したプログラムを、当該コンピュータの記憶部４０に格納しておき、当該コンピュータの中央演算処理装置（ＣＰＵ）によってこのプログラムを読み出して実行させることで制御部１０の各機能を実現することができる。

制御部１０は、コンクリート構造物を撮像した画像を入力してコンクリート構造物の表面に生じたひび割れの候補画像を生成するひび割れ探査部２０と、ひび割れの候補画像からコンクリート構造物の表面に生じたひび割れを抽出するひび割れ判定部３０とを備え、ひび割れ判定部３０は、ひび割れ抽出画像を生成して表示装置（図示せず）に表示させることができる。尚、表示装置は、例えばパーソナルコンピュータのモニタを用いて実現できる。

ひび割れ探査部２０は、コンクリート構造物が映った画像からひび割れと思しきパターンの検知を行うことができ、より具体的には、ひび割れ特徴量抽出部２０１と、ひび割れ収束部２０２とを備える。ひび割れ判定部３０は、ひび割れ探査部２０で検知されたひび割れ候補のパターンを自動判定することができ、近傍連結判定部３０１と、特徴量判定部３０２と、飛び地判定部３０３と、論理積演算部３０４とを備える。

ひび割れ特徴量抽出部２０１は、コンクリート構造物を撮像した画像に対して方向及び幅の異なる複数種類の２次元エッジフィルタを個別に適用し、画素ごとに、複数種類のエッジフィルタ適用後の画素値をそれぞれ比較し、最大の画素値をひび割れ特徴量として選定する機能を有する。より具体的には、ひび割れ特徴量抽出部２０１は、２次元エッジフィルタの各対応画素についてｍ種類のエッジ方向とｎ種類のエッジ幅の組み合わせを持つｍ×ｎ個の２次元エッジフィルタ処理を行い、そのｍ×ｎ個の出力値の中から最大値となるものをひび割れ特徴量ξとして選択し、２次元エッジフィルタの対応画素から一定範囲に存在し、識別情報がひび割れを撮像した領域に該当することを示した周辺ニューロンに対応する周辺画素の位置と、２次元エッジフィルタの対応画素の位置との間のユークリッド距離に基づいて、当該周辺画素に対する前記対応画素の連結度ηを算出し、ひび割れ特徴量ξと連結度ηの各々に動的係数を乗算したものを線形結合してひび割れ収束部２０２（Ｈｏｐｆｉｅｌｄ型ニューラルネットワーク）の入力値とすることによって、ひび割れ収束部２０２における２次元ニューロン配列の発火パターンをひび割れパターンへと収束させる機能を有する。

ひび割れ収束部２０２は、ひび割れ特徴量からなる画素値を有する処理対象画像に対してニューラルネットの動作式を用いてニューラル収束演算を実行し、ひび割れ候補画像を生成する機能を有する。より具体的には、ひび割れ収束部２０２は、ひび割れ特徴量抽出部２０１からの、ひび割れ特徴量ξと連結度ηの各々に動的係数を乗算したものを線形結合して得られる２次元ニューロン配列の発火パターンのニューロンを入力値とする２次元のニューロン配列を有するＨｏｐｆｉｅｌｄ型ニューラルネットワークとして構成することができる。

近傍連結判定部３０１は、ひび割れ候補画像に対して、ひび割れ収束部２０２で発火したニューロンの対応画素をひび割れ候補画素とし、ひび割れ候補画素中の判定対象画素に対して所定の近傍範囲内にある画素についてひび割れと判定する機能を有する。より具体的には、近傍連結判定部３０１は、ひび割れ候補画像を短冊状に分割したセグメント内で、ひび割れ候補画素中の判定対象画素に対して所定の近傍範囲内にあるひび割れ候補画素以外の各画素とのユークリッド距離を合計した値が一定値以下の場合に近傍連結カウント値ｃｎｔとしてカウントし、近傍連結カウント値ｃｎｔが一定値以上の場合にひび割れセグメントと判定し、さらにひび割れセグメントの数が一定数以上の場合にひび割れと判定する機能を有する。

特徴量判定部３０２は、ひび割れ候補画素の特徴量を合計した値に基づいてひび割れを判定する機能を有する。より具体的には、特徴量判定部３０２は、ひび割れ候補画素の特徴量ξを合計した値が当該入力画像全体（ひび割れ候補画像全体）の平均輝度から算出した閾値θを超えた場合にひび割れと判定する機能を有する。

飛び地判定部３０３は、ひび割れ候補画素中の判定対象画素に対して隣接する画素のユークリッド距離に基づいてひび割れと判定する機能を有する。より具体的には、飛び地判定部３０３は、ひび割れ候補画素中の判定対象画素と判定対象画素に隣接した画素とのユークリッド距離が一定値以上であり、かつ判定対象画素の前後の所定数画素において隣接する画素間のユークリッド距離が全て一定値内である場合、飛び地でありひび割れではないと判定する機能を有する。

論理積演算部３０４は、近傍連結判定部３０１、特徴量判定部３０２及び飛び地判定部３０３にて判定された結果の論理積を行う機能を有する。

また、本実施例のひび割れ検知装置１００の動作について図２を参照して説明する。

〔装置動作〕
ひび割れ検知装置１００は、まず、コンクリート構造物を撮像した画像を入力して（ステップＳ１００）、ひび割れ特徴量抽出部２０１によって、コンクリート構造物を撮像した画像に対して方向及び幅の異なる複数種類の２次元エッジフィルタを個別に適用し、画素ごとに、複数種類のエッジフィルタ適用後の画素値をそれぞれ比較し、最大の画素値をひび割れ特徴量として選択する（ステップＳ２００）。

次に、ひび割れ検知装置１００は、ひび割れ収束部２０２によって、ひび割れ特徴量からなる画素値を有する処理対象画像に対してニューラルネットの動作式を用いてニューラル収束演算を実行し、ひび割れ候補画像を生成する（ステップＳ３００）。

次に、ひび割れ検知装置１００は、近傍連結判定部３０１によって、ひび割れ候補画像に対して、ひび割れ収束部２０２で発火したニューロンの対応画素をひび割れ候補画素とし、ひび割れ候補画素中の判定対象画素に対して所定の近傍範囲内にある画素についてひび割れと判定するとともに、特徴量判定部３０２によって、ひび割れ候補画素の特徴量を合計した値に基づいてひび割れを判定し、さらに、飛び地判定部３０３によって、ひび割れ候補画素中の前記判定対象画素に対して隣接する画素のユークリッド距離に基づいてひび割れと判定し、論理積演算部３０４によって、近傍連結判定部３０１、特徴量判定部３０２及び飛び地判定部３０３にて判定された結果の論理積演算を実行し、最終的なひび割れ画像を抽出して所定の表示装置（図示せず）に表示する（ステップＳ４００）。

これにより、ひび割れ検知装置１００は、建物や電柱等のコンクリート構造物を撮像した画像から、コンクリート構造物の表面に生じたひび割れを画像処理によって、ひび割れの抽出精度を向上させて自動検知することが可能となる。

以下より具体的に、ひび割れ検知装置１００の処理について、図３（Ａ）に示すひび割れ画像を例に挙げて各機能の詳細を説明する。図３（Ｂ）は、図３（Ａ）に示す処理対象の入力画像を従来のニューラルネットワークによるひび割れ検知アルゴリズム（前述した非特許文献２参照）によって処理した結果であり、実際のひび割れ部分へと収束せず、正確なひび割れ形状の抽出に成功していない。なぜならば、この非特許文献２のアルゴリズムでは、一次元のエッジフィルタによって抽出された特徴量をニューラルネットワークの入力情報としてひび割れ探索を行うことになるが、図３（Ｃ）に示すように、この一次元のエッジフィルタではひび割れパターンとひび割れ以外のノイズパターンの間に特徴量の差が出難く、ノイズパターンへ収束してしまうケースがよくあるからである。ゆえに、より高い確率でひび割れに収束させるには、ひび割れの特徴のみを優位に捉えてニューラルネットワークの入力情報として設定する必要がある。尚、ニューラルネットワークの詳細は、以下のＵＲＬを参考にされたい。

＜参考ＵＲＬ＞
http://msdn.microsoft.com/ja-jp/academic/cc998609.aspx
http://ja.wikipedia.org/wiki/%E3%83%8B%E3%83%A5%E3%83%BC%E3%83%A9%E3%83%AB%E3%83%8D%E3%83%83%E3%83%88%E3%83%AF%E3%83%BC%E3%82%AF

図４及び図５は、ひび割れ検知装置１００におけるひび割れ特徴量抽出部２０１の処理の詳細を示している。まず、図４及び図５に示すような方向および幅の異なる複数の２次元エッジフィルタを処理対象の入力画像全体に適用し、式（１）に示すように各画素毎にその中から最大値となるものをひび割れ特徴量ξとして選択する。

これらのエッジフィルタのサイズは、入力画像の解像度と探索するひび割れの平均的な幅の両方を加味して設計しなければならない。例えば、図３（Ａ）の画像は、約１４００万画素のデジタルカメラを用いて、このデジタルカメラの８メートル先にある幅０．１ｍｍのひび割れを有するコンクリート構造物の表面を撮像したものであるが、これを１８ｘ２６、３ｘ５、９９ｘ１０１の三種類のエッジフィルタで処理した結果を、それぞれ図６（Ａ）、図６（Ｂ）、図６（Ｃ）に示す。この中では、図６（Ａ）においてひび割れの特徴量が最も優位に抽出されていることが分かるが、図６（Ｂ）ではひび割れと汚れ・傷との間に特徴量の差が少なく、前者のエッジフィルタの方が抽出効果として高くなることが分かる。しかしながら、図６（Ｃ）に示すようにエッジフィルタのサイズを大きくし過ぎても、かえってひび割れの特徴量が埋没してしまい逆効果であることが分かる。従って、ターゲットまでの距離、探索するひび割れ幅、撮像に用いるデジタルカメラの画素数等のパラメータから、予め最適フィルタサイズを調整する必要がある。

図６（Ａ）における２次元エッジフィルタのサイズ（１８×２６）は、前述した「カメラの８メートル先にある幅０．１ｍｍのひび割れを約１４００万画素のデジタルカメラを用いて」撮影した数１０枚のひび割れ画像をひび割れ抽出・判定し、成功が最も多かったサイズである。したがって、ターゲットまでの距離、探索するひび割れ幅、カメラ画素数等のパラメータから、予め最適フィルタサイズを調べておくのが好適である。

図７は、本発明に係るひび割れ収束部２０２における処理の詳細として、前述した非特許文献２に記載されている処理を引用している。ただし、本発明に係るひび割れ検知装置１００では、（前述した非特許文献２に記載される）式（２）におけるｇ（ｉ，ｊ）を、ひび割れ特徴量抽出部２０１で用いる２次元エッジフィルタの出力ξで置き換えることに留意する。

ここで、実用上の課題としてこのニューラルネットワークの処理速度について言及しておく。一般的にニューラルネットは繰り返し演算を用いており、入力画像が大きくなるにつれて計算コストが大きくなるという問題は本発明に係るひび割れ収束部２０２のアルゴリズムにおいても例外ではない。この問題を回避するための第１の工夫を図８（Ａ）に示す。

従来のアルゴリズムにおいては式（２）の係数Ａ（ｔ）と係数Ｂ（ｔ）は定数であるが、本発明に係るひび割れ収束部２０２では係数Ａ（ｔ）と係数Ｂ（ｔ）とを動的に変化させることで高速化を行う。ニューラルネットにおける高速化の手段としては繰り返し回数を減らすことが有効手段であり、繰り返し回数を減らすためには式（２）におけるひび割れパターンへ収束する力を向上させる必要がある。式（２）において、ひび割れ特有の線形状を持つパターンへ収束する力は係数Ａ（ｔ）の大きさに反比例し、係数Ｂ（ｔ）の大きさに比例するという性質がある。しかしながら、ひび割れパターンへ収束する力を大きくするために単純に係数Ａ（ｔ）を小さな、また係数Ｂ（ｔ）を大きな定数として固定してしまうと、ひび割れ以外の画像において線形状に収束してしまうエラーが起きる。この問題を回避するため、係数Ａ（ｔ），Ｂ（ｔ）を繰り返し回数ｔに対して動的に変化させる。

例えば、図８（Ａ）に示すように始めの数回τは、係数Ａ（ｔ）を大きく、また係数Ｂ（ｔ）を小さく設定してエッジ特徴量ξが大きい画素に収束し易くする。そして、繰り返し回数がτ＋１以降は、徐々にＡ（ｔ）を小さくＢ（ｔ）を大きくして画素連結の影響を強くしていくことで、少ない繰り返し回数でひび割れに収束させることができ、ひび割れを抽出する力とひび割れ以外を排他する力を共存させることが可能となる。

次に、収束速度を高めるための第２の工夫を図８（Ｂ）に示す。従来のアルゴリズムでは、全ての繰り返しステップにおいて全てのニューロンについて計算を行っていたが、本発明に係るひび割れ収束部２０２では繰り返しステップが進むにつれてニューロンの計算量を減らすことでさらなる高速化を行う。

図８（Ｂ）に示すように、発火ニューロンまでの距離合計が一番小さくなるような重心列を求め、そこから一定距離Ｒ以内で発火しているニューロンのみ計算を行うようにする。その範囲外で発火しているニューロンについては消火し、当該範囲内のどこかでランダム発火させる。但し、一定距離Ｒを小さくすればするほど計算量が減って高速化される反面、ひび割れが斜め方向に生じている場合には発見し辛くなるというジレンマがあるので、時間と精度のどちらを優先すべきかを用途に応じて考慮し、一定距離Ｒを慎重に決定しなければならない。本実施例では一定距離Ｒをウィンドウサイズの６分の１とした。ここで、「発火」とは、式（２）におけるＶ_ｉ，ｊが１の値を持つ状態を云う。一方、「消火」とは、式（２）におけるＶ_ｉ，ｊの値を０にすることを云う。同時にＵ_ｉ，ｊも二度と発火しないように、小さな値（例：‐１０００００００）を代入する。

図９は、ひび割れ判定部におけるデジタル画像上に現れたひび割れ形状の定義を模式的に表した図である。一般的に、ひび割れの形状を局所的な画素の連なりとして見た場合、図９（Ａ）に示すように連続かつ適度に蛇行した分布となる。なぜならば、コンクリート構造物のひび割れというものは、内部に混在した大きさの異なる骨材に沿って連続的に破断しているため、図９（Ｂ）に示すような断絶分布となるものや、或いは、図９（Ｃ）に示すように規則的な直線分布となるものは、ひび割れではないからである。例えば、ニューラルネットが表面に浮き出た小さな骨材に収束した場合などは図９（Ｂ）のような分布になり易く、或いは、人工的な型枠跡に収束した場合などは図９（Ｃ）のような分布になり易く、この定義を用いればいずれもひび割れではないと判定することができる。

図１０は、ひび割れ検知装置１００における近傍連結判定部３０１におけるアルゴリズムのフローチャートを示している。前述したように、本発明に係る近傍連結判定部３０１では、連続かつ適度に蛇行した分布であるひび割れは局所的に連結しているという特徴を利用してひび割れを判定する。

まず、近傍連結判定部３０１は、ひび割れカウント変数ｃｒａｃｋ＿ｃｎｔを０、近傍連結探索範囲ｋａｚｕを任意数、式（３）を用いてひび割れ候補画素の各ｉ行のｙ座標を座標変数ｐｓ（ｉ）で初期化する。入力画像を複数のセグメントに分割し、一番上のセグメントから計算を行う（ステップＳ１）。

次に、近傍連結判定部３０１は、セグメントごとの初期化として、近傍連結カウント変数ｃｎｔを０で初期化し、各セグメントの最上部から近傍連結の計算を行う（ステップＳ２）。

続いて、近傍連結判定部３０１は、式（４）を用いて近傍連結値ｋｙｏｒｉ（ｉ，ｐｓ（ｉ））を計算する（ステップＳ３）。

近傍連結判定部３０１は、発火ニューロン位置のばらつきを考慮し設定した閾値よりｋｙｏｒｉ（ｉ，ｐｓ（ｉ））が小さければ（ステップＳ４）、近傍で連結していると判定し近傍連結カウント変数ｃｎｔに１を加算する（ステップＳ５）。

近傍連結判定部３０１は、セグメント内の最後の画素であれば次の処理に進み、最後の画素でなければ次の行に移動し（ステップＳ６）、ｋｙｏｒｉ（ｉ，ｐｓ（ｉ））の計算を繰り返す（ステップＳ７）。

近傍連結判定部３０１は、近傍連結カウント変数ｃｎｔをセグメントサイズで除算した値が、近傍連結カウント閾値ｃｎｔ＿ｔｈｒｅｓｈより大きければひび割れセグメントと判定し（ステップＳ８）、ひび割れカウント変数ｃｒａｃｋ＿ｃｎｔに１を加算する（ステップＳ９）。

近傍連結判定部３０１は、全セグメントが計算を終了していれば次の処理に進み（ステップＳ１０）、終了していなければ次のセグメントに移動し、セグメントごとの初期化から同様の計算を行う（ステップＳ１１）。

最後に、近傍連結判定部３０１は、ひび割れカウント変数ｃｒａｃｋ＿ｃｎｔがひび割れカウント閾値ｃｒａｃｋ＿ｃｎｔ＿ｔｈｒｅｓｈより大きければ（ステップＳ１２）、ひび割れと判定し（ステップＳ１３）、小さければひび割れではないと判定する（ステップＳ１４）。

図１１は、近傍連結判定部３０１における計算例を示している。図１１（図示左）は、ひび割れに収束した画像を複数のセグメントに分割したものである。各セグメントにおいて、図１１（図示右）のようにひび割れ候補画素に対して近傍連結の計算を行う。太線の四角が近傍連結の計算対象画素であり、近傍連結探索範囲内のひび割れ候補画素とのユークリッド距離を計算する。仮に近傍連結閾値を３０と設定し、計算対象画素の上下３画素を近傍連結探索範囲とする。

ひび割れの例では、計算対象画素から１行上の画素との距離は１であり、２行上は３であり、３行上は５であり、１行下は１であり、２行下は３であり、３行下は６であることから、近傍連結値ｋｙｏｒｉ（ｉ，ｐｓ（ｉ））は１９となる。この場合、近傍連結判定部３０１は、近傍連結値が近傍連結閾値より小さいことから、連結していると判定し近傍連結カウント変数に１を加算する。

同様にひび割れではない例では、計算対象画素から１行上の画素との距離は４０であり、２行上は２２であり、３行上は５であり、１行下は１５であり、２行下は２０であり、３行下は３であることから、近傍連結値ｋｙｏｒｉ（ｉ，ｐｓ（ｉ））は１０５となる。この場合、近傍連結判定部３０１は、近傍連結値が近傍連結閾値より大きいことから、連結していないと判定し近傍連結カウント変数にはカウントしない。

図１２は、ひび割れ検知装置１００における特徴量判定部３０２におけるアルゴリズムのフローチャートを示している。近傍連結判定部３０１のみでは、ひび割れではない入力画像においてひび割れの形状にひび割れ候補画素が収束した場合、ひび割れと誤って判定してしまう可能性がある。そこで、特徴量判定部３０２は、ひび割れ候補画素の２次元エッジ特徴量の大きさからひび割れを判定する。ひび割れがあると２次元エッジ特徴量は大きくなり、ひび割れが無いと２次元エッジ特徴量は小さくなる。また、通常ひび割れは黒いため、影などで暗い画像はひび割れがあったとしても２次元エッジ特徴量が小さくなってしまう。それゆえ、暗いほどひび割れを判定する閾値が小さくなるように入力画像全体の輝度に応じて閾値を動的に設定し、ひび割れ抽出力を高める。

まず、特徴量判定部３０２は、処理対象の入力画像全体の初期化として、特徴量カウント変数ｆｅａｔｕｒｅ＿ｃｎｔに０を代入する。入力画像を複数のセグメントに分割し、一番上のセグメントから計算を行う（ステップＳ２１）。

次に、特徴量判定部３０２は、式（５）に示すように当該入力画像全体の輝度から特徴量閾値ｆｅａｔｕｒｅ＿ｓｕｍ＿ｔｈｒｅｓｈを計算する（ステップＳ２２）。

続いて、特徴量判定部３０２は、セグメントごとの初期化として、２次元エッジ合計特徴量ｆｅａｔｕｒｅ＿ｓｕｍに０を代入する（ステップＳ２３）。

特徴量判定部３０２は、式（６）に示すように、セグメント内の全ひび割れ候補画素の２次元エッジ合計特徴量ｆｅａｔｕｒｅ＿ｓｕｍを計算する（ステップＳ２４）。

特徴量判定部３０２は、もし、２次元エッジ合計特徴量ｆｅａｔｕｒｅ＿ｓｕｍが特徴量閾値ｆｅａｔｕｒｅ＿ｓｕｍ＿ｔｈｒｅｓｈより大きければ（ステップＳ２５）、ひび割れセグメントと判定し特徴量カウント変数ｆｅａｔｕｒｅ＿ｃｎｔに１を加算する（ステップＳ２６）。

特徴量判定部３０２は、全セグメントが計算を終了していれば次の処理に進み（ステップＳ２７）、終了していなければ次のセグメントに移動し、セグメントごとの初期化から同様の計算を行う（ステップＳ２８）。

最後に、特徴量判定部３０２は、特徴量カウント変数ｆｅａｔｕｒｅ＿ｃｎｔが特徴量カウント閾値ｆｅａｔｕｒｅ＿ｃｎｔ＿ｔｈｒｅｓｈより大きければ（ステップＳ２９）、ひび割れと判定し、小さければひび割れではないと判定する。

図１３は、ひび割れ検知装置１００における飛び地判定部３０３の処理の詳細を示しており、図１４は、ひび割れ検知装置１００における飛び地判定部３０３のアルゴリズムのフローチャートを示している。図１３に示すように、近傍連結判定部３０２と特徴量判定部３０３のみではセグメントごとのひび割れ判定であるため、汚れをひび割れと判定とすることがある。それゆえ、図１３に示すように、ひび割れは入力画像全体で連続しているとの観点から、飛び地判定部３０３は、飛び地がある場合はひび割れではないとする判定を行う。

まず、飛び地判定部３０３は、飛び地カウント変数ｔｏｂｉ＿ｃｎｔを０、島サイズｓｈｉｍａ＿ｓｉｚｅを任意の数、式（３）を用いてひび割れ候補画素の各ｉ行のｙ座標を座標変数ｐｓ（ｉ）で初期化する。入力画像の最上部の行から計算を開始する（ステップＳ３１）。

次に、飛び地判定部３０３は、隣接距離カウント変数ｋｙｏｒｉ＿ｃｎｔと隣接距離計算用変数ｉｉを０で初期化する（ステップＳ３２）。

続いて、飛び地判定部３０３は、ｉ行のひび割れ候補画素に対して式（７）を用いて上部隣接距離ｔｏｎａｒｉ＿ｋｙｏｒｉ＿ｕを計算し、式（８）を用いて下部隣接距離ｔｏｎａｒｉ＿ｋｙｏｒｉ＿ｓを計算する（ステップＳ３３）。

飛び地判定部３０３は、上部隣接距離ｔｏｎａｒｉ＿ｋｙｏｒｉ＿ｕが隣接距離閾値ｔｏｎａｒｉ＿ｋｙｏｒｉ＿ｔｈｒｅｓｈより小さければ（ステップＳ３４）隣接距離カウント変数ｋｙｏｒｉ＿ｃｎｔに１を加算する（ステップＳ３５）。同様に、飛び地判定部３０３は、下部隣接距離ｔｏｎａｒｉ＿ｋｙｏｒｉ＿sが隣接距離閾値ｔｏｎａｒｉ＿ｋｙｏｒｉ＿ｔｈｒｅｓｈより小さければ（ステップＳ３６）、隣接距離カウント変数ｋｙｏｒｉ＿ｃｎｔに１を加算する（ステップＳ３７）。

飛び地判定部３０３は、隣接距離計算用変数ｉｉが島サイズｓｈｉｍａ＿ｓｉｚｅであれば次の処理に進み（ステップＳ３８）、島サイズ未満であれば隣接距離計算用変数ｉｉに１を加算する（ステップＳ３９）。

飛び地判定部３０３は、式（９）を用いて飛び距離ｔｏｂｉ＿ｋｙｏｒｉを計算する（ステップＳ４０）。

ここで、飛び地判定部３０３は、隣接距離カウント変数ｋｙｏｒｉ＿ｃｎｔが０であり（ステップＳ４１）かつ飛び距離ｔｏｂｉ＿ｋｙｏｒｉが飛び距離閾値ｔｏｂｉ＿ｋｙｏｒｉ＿ｔｈｒｅｓｈより大きければ（ステップＳ４２）飛び地カウント変数ｔｏｂｉ＿ｃｎｔに１を加算する（ステップＳ４３）。

飛び地判定部３０３は、最終行の計算が終了していれば次の処理に進み（Ｓ１４)，終了していなければ次の行に移動する（ステップＳ４５）。

最後に、飛び地判定部３０３は、飛び地カウント変数ｔｏｂｉ＿ｃｎｔが飛び地カウント閾値ｔｏｂｉ＿ｃｎｔ＿ｔｈｒｅｓｈより大きければ（ステップＳ４６）、ひび割れではないと判定し、それ以外はひび割れと判定する。

図１５は、飛び地判定部３０３の処理例を示している。飛び地判定部３０３は、判定対象画素と判定対象画素に隣接したひび割れ候補画素間のユークリッド距離が一定値以上であり、かつ判定対象画素の前後の所定数画素において隣接するひび割れ候補画素間のユークリッド距離が全て一定値内である場合、飛び地でありひび割れではないと判定する。

図１６に、本発明による一実施例のひび割れ検知装置１００を用いてひび割れを抽出した結果を示す。従来手法では抽出困難であったひび割れ形状が正確に抽出され、かつ、ひび割れと自動判定されている。

以上に説明した通り、本発明によれば、コンクリート構造物の画像から、ひび割れを高精度に抽出し、かつ、自動判定することが可能となる。

本発明に係るひび割れ検知装置をコンピュータで構成した場合、各機能を実現する処理内容を記述したプログラムを、当該コンピュータの内部又は外部の記憶部に格納しておき、当該コンピュータの中央演算処理装置（ＣＰＵ）によってこのプログラムを読み出して実行させることで実現することができる。また、このようなプログラムは、例えばＤＶＤ又はＣＤ−ＲＯＭなどの可搬型記録媒体の販売、譲渡、貸与等により流通させることができるほか、そのようなプログラムを、例えばネットワーク上にあるサーバの記憶部に記憶しておき、ネットワークを介してサーバから他のコンピュータにそのプログラムを転送することにより、流通させることができる。また、そのようなプログラムを実行するコンピュータは、例えば、可搬型記録媒体に記録されたプログラム又はサーバから転送されたプログラムを、一旦、自己の記憶部に記憶することができる。また、このプログラムの別の実施態様として、コンピュータが可搬型記録媒体から直接プログラムを読み取り、そのプログラムに従った処理を実行することとしてもよく、更に、このコンピュータにサーバからプログラムが転送される度に、逐次、受け取ったプログラムに従った処理を実行することとしてもよい。従って、本発明は、前述した実施例に限定されるものではなく、その主旨を逸脱しない範囲において種々変更可能である。

本発明によれば、コンクリート構造物の画像から、ひび割れを高精度に抽出し、かつ、自動判定することが可能となるので、ひび割れの検査を要する用途に有用である。

１０ 制御部
２０ ひび割れ探査部
３０ ひび割れ判定部
４０ 記憶部
１００ ひび割れ検知装置
２０１ ひび割れ特徴量抽出部
２０２ ひび割れ収束部
３０１ 近傍連結判定部
３０２ 特徴量判定部
３０３ 飛び地判定部
３０４ 論理積演算部

Claims (8)

1. コンクリート構造物を撮像した画像から、コンクリート構造物の表面に生じたひび割れを画像処理によって自動検知するひび割れ検知装置であって、
   コンクリート構造物を撮像した画像を入力してコンクリート構造物の表面に生じたひび割れの候補画像を生成するひび割れ探査部と、該ひび割れの候補画像からコンクリート構造物の表面に生じたひび割れを抽出するひび割れ判定部とを備え、
   前記ひび割れ探査部は、
   前記コンクリート構造物を撮像した画像に対して方向及び幅の異なる複数種類の２次元エッジフィルタを個別に適用し、画素ごとに、複数種類のエッジフィルタ適用後の画素値をそれぞれ比較し、最大の画素値をひび割れ特徴量として選定するひび割れ特徴量抽出部と、
   前記ひび割れ特徴量からなる画素値を有する処理対象画像に対してニューラルネットの動作式を用いてニューラル収束演算を実行し、ひび割れ候補画像を生成するひび割れ収束部とを備え、
   前記ひび割れ判定部は、
   前記ひび割れ候補画像に対して、前記ひび割れ収束部で発火したニューロンの対応画素をひび割れ候補画素とし、前記ひび割れ候補画素中の判定対象画素に対して所定の近傍範囲内にある画素についてひび割れと判定する近傍連結判定部と、
   前記ひび割れ候補画素の特徴量を合計した値に基づいてひび割れを判定する特徴量判定部と、
   前記ひび割れ候補画素中の前記判定対象画素に対して隣接する画素のユークリッド距離に基づいてひび割れと判定する飛び地判定部と、
   前記近傍連結判定部、前記特徴量判定部及び前記飛び地判定部にて判定された結果の論理積を行う論理積演算部とを備えることを特徴とする、ひび割れ検知装置。
2. 前記ひび割れ収束部は、２次元のニューロン配列を有するＨｏｐｆｉｅｌｄ型ニューラルネットワークとして構成されていることを特徴とする、請求項１に記載のひび割れ検知装置。
3. 前記近傍連結判定部は、前記ひび割れ候補画像を短冊状に分割したセグメント内で、前記ひび割れ候補画素中の判定対象画素に対して所定の近傍範囲内にある前記ひび割れ候補画素以外の各画素とのユークリッド距離を合計した値が一定値以下の場合に近傍連結カウント値としてカウントし、前記近傍連結カウント値が一定値以上の場合にひび割れセグメントと判定し、さらに前記ひび割れセグメントの数が一定数以上の場合にひび割れと判定する手段を有することを特徴とする、請求項１又は２に記載のひび割れ検知装置。
4. 前記特徴量判定部は、前記ひび割れ候補画素の特徴量を合計した値が当該ひび割れ候補画像全体の平均輝度から算出した閾値を超えた場合にひび割れと判定する手段を有することを特徴とする、請求項１〜３のいずれか一項に記載のひび割れ検知装置。
5. 前記飛び地判定部は、前記ひび割れ候補画素中の前記判定対象画素と前記判定対象画素に隣接した画素とのユークリッド距離が一定値以上であり、かつ前記判定対象画素の前後の所定数画素において隣接する画素間のユークリッド距離が全て一定値内である場合、飛び地でありひび割れではないと判定する手段を有することを特徴とする、請求項１〜４のいずれか一項に記載のひび割れ検知装置。
6. コンクリート構造物を撮像した画像から、コンクリート構造物の表面に生じたひび割れを画像処理によって自動検知するひび割れ検知装置であって、
   コンクリート構造物を撮像した画像を入力してコンクリート構造物の表面に生じたひび割れの候補画像を生成するひび割れ探査部と、該ひび割れの候補画像からコンクリート構造物の表面に生じたひび割れを抽出するひび割れ判定部とを備え、
   前記ひび割れ探査部は、
   前記コンクリート構造物を撮像した画像に対して方向及び幅の異なる複数種類の２次元エッジフィルタを個別に適用し、画素ごとに、複数種類のエッジフィルタ適用後の画素値をそれぞれ比較し、最大の画素値をひび割れ特徴量として選定するひび割れ特徴量抽出部と、
   前記ひび割れ特徴量からなる画素値を有する処理対象画像に対してニューラルネットの動作式を用いてニューラル収束演算を実行し、ひび割れ候補画像を生成するひび割れ収束部とを備え、
   前記ひび割れ特徴量抽出部は、前記２次元エッジフィルタの各対応画素についてｍ種類のエッジ方向とｎ種類のエッジ幅の組み合わせを持つｍ×ｎ個の２次元エッジフィルタ処理を行い、そのｍ×ｎ個の出力値の中から最大値となるものをひび割れ特徴量として選択し、前記対応画素から一定範囲に存在し、識別情報がひび割れを撮像した領域に該当することを示した周辺ニューロンに対応する周辺画素の位置と、前記対応画素の位置との間のユークリッド距離に基づいて、前記周辺画素に対する前記対応画素の連結度を算出し、前記ひび割れ特徴量と前記連結度の各々に動的係数を乗算したものを線形結合して前記ひび割れ収束部の入力値とすることによって、前記ひび割れ収束部における前記２次元ニューロン配列の発火パターンをひび割れパターンへと収束させる手段を有することを特徴とする、ひび割れ検知装置。
7. ひび割れ検知装置によってコンクリート構造物を撮像した画像から、コンクリート構造物の表面に生じたひび割れを画像処理によって自動検知するひび割れ検知方法であって、
   前記ひび割れ検知装置の処理手順は、
   （Ａ）前記コンクリート構造物を撮像した画像に対して方向及び幅の異なる複数種類の２次元エッジフィルタを個別に適用し、画素ごとに、複数種類のエッジフィルタ適用後の画素値をそれぞれ比較し、最大の画素値をひび割れ特徴量として選定するステップと、
   （Ｂ）前記ひび割れ特徴量からなる画素値を有する処理対象画像に対してニューラルネットの動作式を用いてニューラル収束演算を実行し、ひび割れ候補画像を生成するステップと、
   （Ｃ）前記ひび割れ候補画像に対して、前記ステップにて発火したニューロンの対応画素をひび割れ候補画素とし、前記ひび割れ候補画素中の判定対象画素に対して所定の近傍範囲内にある画素についてひび割れと判定するステップと、
   （Ｄ）前記ひび割れ候補画素の特徴量を合計した値に基づいてひび割れを判定するステップと、
   （Ｅ）前記ひび割れ候補画素中の前記判定対象画素に対して隣接する画素のユークリッド距離に基づいてひび割れと判定するステップと、
   前記ステップ（Ｃ）、前記ステップ（Ｄ）及び前記ステップ（Ｅ）にて判定された結果の論理積を行うステップと、
   を含むことを特徴とする、ひび割れ検知方法。
8. コンクリート構造物を撮像した画像から、コンクリート構造物の表面に生じたひび割れを画像処理によって自動検知するひび割れ検知装置として構成するコンピュータに、
   （Ａ）前記コンクリート構造物を撮像した画像に対して方向及び幅の異なる複数種類の２次元エッジフィルタを個別に適用し、画素ごとに、複数種類のエッジフィルタ適用後の画素値をそれぞれ比較し、最大の画素値をひび割れ特徴量として選定するステップと、
   （Ｂ）前記ひび割れ特徴量からなる画素値を有する処理対象画像に対してニューラルネットの動作式を用いてニューラル収束演算を実行し、ひび割れ候補画像を生成するステップと、
   （Ｃ）前記ひび割れ候補画像に対して、前記ステップにて発火したニューロンの対応画素をひび割れ候補画素とし、前記ひび割れ候補画素中の判定対象画素に対して所定の近傍範囲内にある画素についてひび割れと判定するステップと、
   （Ｄ）前記ひび割れ候補画素の特徴量を合計した値に基づいてひび割れを判定するステップと、
   （Ｅ）前記ひび割れ候補画素中の前記判定対象画素に対して隣接する画素のユークリッド距離に基づいてひび割れと判定するステップと、
   前記ステップ（Ｃ）、前記ステップ（Ｄ）及び前記ステップ（Ｅ）にて判定された結果の論理積を行うステップと、
   を実行させるためのプログラム。