JP7008408B2 - ひび割れ検出方法 - Google Patents

ひび割れ検出方法 Download PDF

Info

Publication number
JP7008408B2
JP7008408B2 JP2017020499A JP2017020499A JP7008408B2 JP 7008408 B2 JP7008408 B2 JP 7008408B2 JP 2017020499 A JP2017020499 A JP 2017020499A JP 2017020499 A JP2017020499 A JP 2017020499A JP 7008408 B2 JP7008408 B2 JP 7008408B2
Authority
JP
Japan
Prior art keywords
crack
image
pixel
wavelet coefficient
path
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2017020499A
Other languages
English (en)
Other versions
JP2018128315A (ja
Inventor
昌美 本澤
三馨 鈴木
賢一 堀口
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Taisei Corp
Original Assignee
Taisei Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Taisei Corp filed Critical Taisei Corp
Priority to JP2017020499A priority Critical patent/JP7008408B2/ja
Publication of JP2018128315A publication Critical patent/JP2018128315A/ja
Application granted granted Critical
Publication of JP7008408B2 publication Critical patent/JP7008408B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Investigating Materials By The Use Of Optical Means Adapted For Particular Applications (AREA)
  • Image Processing (AREA)

Description

本発明は、コンクリート表面に生じているひび割れの検出をおこなうひび割れ検出方法に係り、特に、ひび割れ検出の解析時間の短縮を図りながら、高精度のひび割れ検出をおこなうことのできるひび割れ検出方法に関するものである。
コンクリート表面上のひび割れを検出する方法としては、従来、調査員が目視観察をおこない、クラックスケールやメジャー等を使用しながらひび割れの幅や長さを測定する方法が一般的であった。しかし、この目視観察による方法は調査員の測定技量などによって精度のばらつきが大きくなることや、ひび割れが大量に存在する場合においては大量の情報を正確に処理するために莫大な労力および時間を要するといった問題があった。
上記の問題に対して、コンクリート表面の撮影画像をコンピュータに取り込み、画像をひび割れ領域とそれ以外の領域とに2値化処理する画像処理手法が適用されている。画像の2値化処理とは、ある濃度値に対して画像の濃度を0または1に表現することであり、例えば、入力画像f(i,j)に対して2値化処理で得られる2値化画像b(i,j)はb(i,j)=1(f(i,j)>k)、0(f(i,j)≦k)となる。ここで、kは2値化する際の閾値であり、したがって2値化画像の良し悪しは閾値kの選定によって決まるといってよい。
従来の閾値を求める手法としては、固定閾値または可変閾値による処理方法がある。固定閾値による処理方法には、Pタイル法やモード法、相関比を用いた方法などが挙げられる。固定閾値による処理方法は、対象画像の濃度ヒストグラムを作成し、画像の背景(コンクリート表面)の濃度値とひび割れの濃度値との間に明確な谷が現れるような双峰性のヒストグラムが得られる場合において有効な方法である。
一方、可変閾値による処理方法は、照明条件などによって撮影ムラが生じ、背景の濃度値と対象部分の濃度値が画像全体で一定でない場合に有効な方法である。この可変閾値処理法は、注目している画素を中心とする局所領域の平均濃度値を閾値とする方法である。この方法の欠点は、背景領域の微妙な濃淡変化に応じて、例えばひび割れ以外のノイズが多い画像となってしまう点である。
従来の画像処理方法は、撮影された入力画像に対して閾値を決定し、2値化処理をおこないながらひび割れの抽出をおこなうものである。すなわち、この一般的な処理の流れは次のようになる。1)撮影画像をコンピュータに取り込んで入力画像を作成する。2)入力画像の濃度の補正をする。3)2値化処理をおこなってひび割れの抽出をおこなう。
上記する従来の画像処理法は、濃度が一様なコンクリート表面上のひび割れの検出においては比較的高精度のひび割れ検出が可能である。しかし、実際のコンクリート構造物の表面は様々な汚れを含んでおり、さらにはひび割れの濃度も、ひび割れの幅や深度などに応じてばらつきがあるのが一般的である。このようなコンクリート表面に対して従来の画像処理法を用いると、ひび割れの抽出に際しては様々な問題が生じ得る。例えば、固定閾値処理の場合において、コンクリート表面上の汚れ領域とひび割れ領域が同程度の濃度値である場合には、これらを2値化処理することが極めて困難となる。濃度ヒストグラムが双峰性を呈していて、閾値を容易に決定できたとしても、ひび割れ領域と判断される範囲には汚れ領域が含まれる可能性が極めて高くなる。また、逆に、ひび割れ周辺部の汚れ領域を含ませないような閾値をあらたに設定しようとすると、今度は他のひび割れ領域を除外してしまうことになってしまう。
可変閾値処理の場合には、コンクリート表面上の汚れが多くなるにしたがって、ひび割れ抽出画像中にひび割れ以外のノイズが多く含まれることになり、場合によってはひび割れ抽出画像を一見しても、どの部分がひび割れ領域なのか全く判別できないこととなる。
上記する従来手法の問題に対して本発明者等は、撮影されたコンクリート表面の汚れや照明条件などによってひび割れの検出が困難な場合においても、簡易に高精度のひび割れ検出をおこなうことのできるひび割れ検出方法を発案し、特許文献1~7にその開示をおこなっている。
これらのひび割れ検出方法はいずれも、対比される2つの濃度に対応したウェーブレット係数を算定するとともに、2つの濃度をそれぞれ変化させた場合のウェーブレット係数を算定してウェーブレット係数テーブルを作成し、ひび割れ検出対象であるコンクリート表面の撮影画像をコンピュータに入力して入力画像とし、この入力画像をウェーブレット変換することによってウェーブレット画像を作成するステップ、ウェーブレット係数テーブル内において局所領域内の近傍画素の平均濃度と注目画素の濃度に対応するウェーブレット係数を閾値とし、注目画素のウェーブレット係数が閾値よりも大きな場合はこの注目画素をひび割れと判定し、閾値よりも小さな場合は注目画素をひび割れでないと判定し、局所領域および注目画素を変化させながら注目画素のウェーブレット係数と閾値との比較をおこなうことでひび割れ抽出画像を作成するステップ、を少なくとも含む検出方法である。
これらのひび割れ検出方法により、簡易に高精度のひび割れ検出をおこなうことが可能になったものの、ひび割れ検出対象であるコンクリート表面の全面が解析対象であることから、解析に時間を要するといった新たな課題が生じている。
特許第4006007号公報 特許第4870016号公報 特許第4980739号公報 特許第5385593号公報 特許第5421192号公報 特許第5705711号公報 特許第5812705号公報
本発明は上記する問題に鑑みてなされたものであり、ひび割れ検出の解析時間の短縮を図りながら、高精度のひび割れ検出をおこなうことのできるひび割れ検出方法を提供することを目的としている。
前記目的を達成すべく、本発明によるひび割れ検出方法は、コンクリート表面に生じているひび割れの検出をおこなうひび割れ検出方法であって、ひび割れ検出対象であるコンクリート表面の撮影画像をコンピュータに入力して入力画像とし、該入力画像においてひび割れと想定されるひび割れ想定線を含む幅が数ピクセルのパスを作成し、該パスに付したひび割れを覆う面を作成してパス画像とする第1のステップ、対比される2つの濃度に対応したウェーブレット係数を算定するとともに、該2つの濃度をそれぞれ変化させた際のウェーブレット係数を算定してウェーブレット係数テーブルを作成し、前記パス画像におけるひび割れ上でウェーブレット係数を計算する第2のステップ、前記ウェーブレット係数テーブル内において局所領域内の近傍画素の平均濃度と注目画素の濃度に対応するウェーブレット係数をウェーブレット係数に関する閾値とし、注目画素のウェーブレット係数が前記閾値よりも大きな場合は注目画素をひび割れと判定し、小さな場合は注目画素をひび割れでないと判定し、局所領域および注目画素を変化させながら注目画素のウェーブレット係数と前記閾値との比較をおこなうことでひび割れ画像を作成する第3のステップ、前記ひび割れ画像に対して細線化処理をおこなってその中心線で構成されるひび割れ細線化画像を作成する第4のステップ、からなるものである。
ウェーブレット(wavelet)とは、小さな波という意味であり、局在性を持つ波の基本単位を、ウェーブレット関数を用いた式で表現することができる。このウェーブレット関数を拡大または縮小することにより、時間情報や空間情報と周波数情報を同時に解析することが可能となる。このウェーブレット係数を、ひび割れを有するコンクリート表面に適用する場合のこの係数の特徴としては、コンクリート表面の濃度と、ひび割れの濃度と、ひび割れ幅に依存するということである。例えば、ひび割れ幅が大きくなるにつれてウェーブレット係数の値は大きくなる傾向があり、また、ひび割れの濃度が濃くなるにつれて(黒色に近づくにつれて)ウェーブレット係数の値は大きくなる傾向がある。
ウェーブレット変換によって算定されるウェーブレット係数を用いて、ひび割れの検出をおこなうアルゴリズムは以下のようになる。まず、コンクリート表面の撮影画像とウェーブレット関数との内積よりウェーブレット係数を求める。このウェーブレット係数を256階調に変換することで、連続量を持ったウェーブレット画像が作成できる。
ウェーブレット係数は、上記するようにひび割れ幅やひび割れの濃度、コンクリート表面の濃度によって変化することから、擬似的に作成されたデータを用いてひび割れの濃度とコンクリート表面の濃度に関するウェーブレット係数を各階調ごとに算定しておき、ウェーブレット係数テーブルを作成しておく。このウェーブレット係数テーブルにある各階調ごとのウェーブレット係数が、ひび割れ検出の際の閾値となる。例えば、対比される2つの濃度(一方の濃度をコンクリート表面の濃度、他方の濃度をひび割れの濃度と仮定することができる)に対応するウェーブレット係数(閾値)がウェーブレット係数テーブルを参照すれば一義的に決定される。したがって、後述するように、撮影画像において対比される2つの濃度間のウェーブレット係数を算定した際に、このウェーブレット係数がウェーブレット係数テーブルの閾値よりも大きな場合は、ひび割れであると判断できるし、閾値よりも小さな場合はひび割れでないと判断することができる。
このウェーブレット係数テーブルを作成する際の擬似的なデータは特に限定するものではないが、例えば、ひび割れ幅が1画素(1ピクセル)~5画素(5ピクセル)までの中で、各画素幅のひび割れごとに、コンクリート表面の階調とひび割れの階調に対応するウェーブレット係数を算定する。閾値の設定に際しては、例えば、ひび割れ幅が1画素の場合のウェーブレット係数のうち、ひび割れに対応するウェーブレット係数を選定し、ひび割れ幅が5画素の場合のウェーブレット係数のうち、ひび割れ領域でない箇所のウェーブレット係数を選定し、これら2つのウェーブレット係数の平均値をもって任意の階調における閾値とすることができる。
本発明のひび割れ検出方法においては、まず、第1のステップにおいて、ひび割れ検出対象であるコンクリート表面の撮影画像をコンピュータに入力して入力画像とし、該入力画像においてひび割れと想定されるひび割れ想定線を含む幅が数ピクセルのパスを作成し、該パスに付したひび割れを覆う面を作成してパス画像とする。
このように、ひび割れ想定線に沿ってパスを作成し、パスに付したひび割れを覆う面を作成することにより、入力画像においてひび割れ位置を指定することができる。具体的には、市販の画像編集ソフトを使用し、ひび割れと想定されるひび割れ想定線を含む幅が数ピクセルの線や面からなるパスを作成する。
ここで、「ひび割れ想定線を含む幅が数ピクセルのパス」とは、ひび割れ想定線の幅よりも数ピクセル幅の大きなパスのことを意味しており、たとえば1~2ピクセル程度の幅のひび割れ想定線に対して、たとえば3~5ピクセル程度の幅の線や面からなるパスが設定できる。このように、ひび割れと想定されるひび割れ想定線を含む幅が数ピクセルのパスを作成することにより、画像上で大まかにひび割れ位置を指定することができる。仮にパスの幅が1ピクセル程度と狭い場合、パスからひび割れ想定線がはみ出す可能性が高くなり、パスからひび割れ想定線がはみ出している箇所はひび割れと判定されない不具合が生じ得る。
本発明のようにひび割れ想定線に沿ってラフなパス(ひび割れ想定線よりも幅が数ピクセル大きなパス)を作成することにより、このような不具合を効果的に解消することができる。なお、ラフに作成するパスの幅を大きく設定すればよいというものではなく、3~5ピクセル程度の幅に設定する場合はノイズが生じ難い一方で、パスの幅を10ピクセル程度と比較的大きな幅に設定するとノイズが生じ易くなる。したがって、ひび割れ想定線の幅を勘案して、それよりも左右に1~2ピクセル程度ずつ大きな幅に設定するのが好ましい。
また、「パスに付したひび割れを覆う面」は、たとえば3~5ピクセル幅の線からなるパスに対し、その左右両側1ピクセル幅もしくは2ピクセル幅を加えた5~9ピクセル幅程度の面である。
次に、第2のステップにおいて、上記するウェーブレット係数テーブルを作成しておくとともに、適宜に設定された広域領域(例えば30×30画素の領域)に対してウェーブレット係数を算定し、この広域領域から一画素移動した広域領域(同じように例えば30×30画素の領域であって、移動前の30×30画素の領域とほとんどの画素が共通している)で、同じようにウェーブレット係数を算定し、ウェーブレット係数の連続量を求める。
次に、第3のステップにおいて、ウェーブレット係数テーブル内において局所領域内の近傍画素の平均濃度と注目画素の濃度に対応するウェーブレット係数をウェーブレット係数に関する閾値とし、注目画素のウェーブレット係数が閾値よりも大きな場合は注目画素をひび割れと判定し(画面上では例えば白色)、小さな場合は注目画素をひび割れでないと判定し(画面上では例えば黒色)、局所領域および注目画素を変化させながら注目画素のウェーブレット係数と閾値との比較をおこなうことでひび割れ画像を作成する。
最後に、第4のステップにおいて、ひび割れ画像に対して細線化処理をおこなってその中心線で構成されるひび割れ細線化画像を作成する。ここでは、たとえばひび割れ全体が1ピクセル幅を有するひび割れ細線化画像が作成される。
本発明のひび割れ検出方法によれば、まず、入力画像においてひび割れ位置を指定した後に解析を実行することから、ひび割れ画像解析範囲を可及的に縮小することができるため、従来の方法に比して解析領域を大幅に縮小することができ、このことによって解析時間を大幅に短縮することが可能になる。さらに、連続解析が可能となり、データの読み込み時間や書き出し時間の短縮も可能になる。
また、解析者が予めひび割れ位置を指定していることから、ひび割れ以外のものが解析対象になり難い。
また、ひび割れ想定線に沿ってラフなパス(ひび割れ想定線の幅よりも数ピクセル幅の大きなパス)を作成してひび割れ位置を指定していることから、パスからひび割れ想定線がはみ出してひび割れと判定されないといった不具合が生じることも解消され、パスがラフゆえにパスの作成も効率的におこなうことができ、パスの作成(ひび割れ位置の特定)に長い時間を必要としない。
さらに、ひび割れ想定線にラフなパスを作成し、ウェーブレット変換した後に細線化処理をおこなうことで、ひび割れをより一層精度よく特定することが可能になる。
また、本発明によるひび割れ検出方法の好ましい実施の形態は、前記第4のステップで作成されたひび割れ細線化画像に対し、予め設定されたピクセル数以下の長さのひび割れから分岐した線をひび割れでないノイズと判定して除去するヒゲ削除処理をおこない、さらに、予め設定されたピクセル数以下の孤立した長さの線をひび割れでないノイズと判定して除去する輪郭線追跡処理をおこなう第5のステップをさらに含むものである。
この実施の形態のひび割れ検出方法では、第4のステップで作成されたひび割れ細線化画像に対して二種類のノイズ除去処理を実行し、より一層精度の高いひび割れの特定をおこなう方法である。本発明のひび割れ検出方法では、第1のステップにおいて、ひび割れ想定線を含むようにラフなパスを作成することから、実際のひび割れ線の周辺にノイズが生じる可能性がある。たとえば、5ピクセル程度以下の長さの線がひび割れでないことが本発明者等の経験則上分かっていることから、このような線をヒゲ(ノイズ)と特定し、ヒゲ削除処理をおこなう。すなわち、5ピクセル程度の長さを閾値に設定し、この閾値未満の長さの線をヒゲと特定して除去するものである。
また、輪郭線追跡処理は、ある任意の画素(ひび割れと判断されている画素)から出発して、隣接する画素がひび割れ箇所の場合には出発画素と接続し、さらに隣接する画素がひび割れ箇所の場合にはさらに双方を接続し、最終的に出発画素に閉合した場合(例えば、第1画素、第2画素、…、第n-1画素、第n画素、第1画素の順に接続される場合)や、次に繋がるひび割れ箇所が存在しなくなった場合に終了するものである。この輪郭線追跡処理によれば、ループ状に閉合するようなひび割れラインや、複数の屈曲部を備えて線状に伸びるひび割れラインなど、適宜のひび割れラインが作成されることになる。この際、繋げられる画素数の最小数を予め設定しておくことにより、この設定数以下の画素はすべてひび割れでないとして、画面のひび割れ表示から削除することができる。
このように二種類のノイズ除去処理をおこなうことで、第1のステップにおいてひび割れ想定線を含むようにラフなパスを作成することによって生じ得るノイズを効果的に除去することができる。
また、本発明によるひび割れ検出方法の他の実施の形態は、前記第5のステップに続いて、特定されたひび割れをひび割れ幅ごとに分類するとともにひび割れ幅ごとに色分けされたひび割れ展開図を作成し、かつ、ひび割れ幅ごとのひび割れ延長ヒストグラムを作成する第6のステップをさらに備えている。
以上の説明から理解できるように、本発明のひび割れ検出方法によれば、入力画像においてひび割れ位置を指定した後に解析を実行することから、ひび割れ画像解析範囲を可及的に縮小することができるため、従来の方法に比して解析領域を大幅に縮小することができ、解析時間を大幅に短縮することが可能になり、ひび割れ検出の解析時間の短縮を図りながら、高精度のひび割れ検出をおこなうことができる。また、ひび割れ位置の指定に際し、ひび割れ想定線を含むラフなパスを作成することにより、パスからひび割れ想定線がはみ出してひび割れと判定されないといった不具合は生じない。
入力画像と局所領域の関係を示した模式図である。 局所領域と注目画素の関係を示した模式図である。 本発明のひび割れ検出方法の実施の形態1のフロー図である。 パスを付したひび割れを覆う面の作成方法を説明した図である。 擬似画像を示した図である。 図5の擬似画像のウェーブレット係数の鳥瞰図である。 ウェーブレット係数テーブルの一実施の形態を示した図である。 本発明のひび割れ検出方法の実施の形態2のフロー図である。 本発明のひび割れ検出方法において10ピクセル幅のパスを作成した際の各ステップにて作成された画像を示した図であって、(a)はコンクリート表面の撮影画像をコンピュータに入力してなる入力画像を示した図であり、(b)はパス画像を示した図であり、(c)はひび割れ画像を示した図であり、(d)はひび割れ細線化画像がヒゲ除去処理および輪郭線追跡処理された後の画像を示した図である。 図9で示すひび割れ検出方法にて作成されたひび割れの延長ヒストグラムの実施例を示した図である。
以下、図面を参照して本発明のひび割れ検出方法の実施の形態1,2を説明する。
(ひび割れ検出方法の実施の形態1)
図1は、入力画像と局所領域の関係を示した模式図である。本発明のひび割れ検出方法では、入力画像1における広域領域2の中心である局所領域3においてウェーブレット変換をおこない、当該局所領域3の中心でひび割れの検出をおこなうものである。入力画像1内でくまなく広域領域2を上下左右に平行移動して、入力画像1内におけるひび割れの検出をおこなう。この方法により、従来の固定閾値法のように、例えば入力画像1内で一つの閾値を決める方法に比べて、精度のよいひび割れの検出をおこなうことができる。
図2は、局所領域3を拡大した図であり、図示する実施形態では、たとえば3×3の9つの画素(8つの近傍画素31,31,…と、中央に位置する注目画素32)の中心でひび割れ判定をおこなう。なお、ウェーブレット係数の算定は、図1における局所領域3を対象としておこなわれる。
ここで、ウェーブレット関数(マザーウェーブレット関数)を用いたウェーブレット変換をおこなうことでウェーブレット係数を算定する算定式を以下に示す。
Figure 0007008408000001
Figure 0007008408000002
Figure 0007008408000003
ここで、f(x,y)は入力画像(ここで、x,yは2次元入力画像中の任意の座標である)を、Ψはマザーウェーブレット関数(ガボール関数)を、(x,y)はΨの平行移動量を、aはΨの拡大や縮小を(ここで、aは周波数の逆数であって、幾つかの周波数領域について計算するための周波数幅を整数kで示した値)、fは中心周波数を、σはガウス関数の標準偏差を、θは波の進行方向を表す回転角を、(x’,y’)は(x,y)を角度θだけ回転させた座標を、それぞれ示している。
ここで、数式1を用いて計算した複数のθ、kに対して、ウェーブレット係数Ψの累計値C(x,y)を求めたのが数式4となる。
Figure 0007008408000004
上記のパラメータは、任意に設定できるが、例えば、σを0.5~2に、aは0~5に、fは0.1に、回転角θは0~180度に、それぞれ設定できる。
数式4における平行移動量(x,y)は、注目画素の位置に対応するものであり、注目画素の位置を順次移動させることによって、ウェーブレット係数の連続量(C(x,y))が算定できる。
局所領域3を構成する全画素に対して、ウェーブレット係数を上算定式に基づいて算定した後、注目画素を一つ左右または上下に移動させてできる広域領域2の全画素において同様にウェーブレット係数を算定する。
次に、図3に基づいて、ひび割れ検出方法の実施の形態1を説明する。
CCDカメラ等のデジタルカメラで撮影されたコンクリート表面の撮影画像をコンピュータに取り込むことにより、入力画像の作成(ステップS10)がおこなわれる。
次に、入力画像においてひび割れと想定されるひび割れ想定線に沿ってパスを作成し、パスに付したひび割れを覆う面を作成してパス画像とする(ステップS20)。すなわち、解析者が予めひび割れ位置の指定をおこなう。
具体的には、図4で示すように、市販の画像編集ソフトを使用し、ひび割れ想定線に沿って、ひび割れ想定線よりも幅が数ピクセル大きな線や面からなるパスを作成する。たとえば1ピクセル程度のひび割れ想定線に対し、その左右に1~5ピクセル程度の幅を有する、計3~10ピクセル程度の幅のパスを作成する。そして、パスに付したひび割れを覆う面は、たとえば3~10ピクセル幅の線からなるパスに対し、その左右両側1ピクセル幅もしくは2ピクセル幅を加えた5ピクセル幅~14ピクセル幅からなる面である(以上、第1のステップ)。
次に、入力画像とは何らの関係もない、対比する2つの濃度からなる擬似画像に対して、ウェーブレット係数の算定をおこなう。例えば、図5に示すように、コンクリート表面と仮定される背景色a(例えば、背景色のR、G、Bを、255,255,255とする)と、ひび割れと仮定される線分b1~b5からなる擬似画像のウェーブレット係数を求める。ここで、線分b1~b5は、線幅が順に1ピクセル~5ピクセルまで変化しており、さらに、各線分は、3種類の濃度を備えている(例えば、線分b1では、濃度の濃い順に、b11(黒色)、b12(薄い黒色)、b13(灰色)と変化している)。この擬似画像に対してウェーブレット変換をおこなうことで算定されるウェーブレット係数の鳥瞰図を示したのが図6である。図6において、X軸は線分の幅を、Y軸は線分の色の濃度を、Z軸はウェーブレット係数をそれぞれ示している。パス画像におけるひび割れ上でウェーブレット係数が算定される(図3のステップS30)。
同時に、対比する2つの濃度の組み合わせをそれぞれ0~255の256階調でおこなうことで、図7に示すようなウェーブレット係数テーブルの作成(図3のステップS30)がおこなわれる(以上、第2のステップ)。
次に、ウェーブレット係数テーブル内において局所領域内の近傍画素の平均濃度と注目画素の濃度に対応するウェーブレット係数をウェーブレット係数に関する閾値とし、注目画素のウェーブレット係数が閾値よりも大きな場合は注目画素をひび割れと判定し(画面上では例えば白色)、小さな場合は注目画素をひび割れでないと判定し(画面上では例えば黒色)、局所領域および注目画素を変化させながら注目画素のウェーブレット係数と閾値との比較をおこなうことにより、ひび割れ画像が作成される(ステップS40、第3のステップ)。
最後に、ひび割れ画像に対して細線化処理を実行することにより、その中心線で構成され、ひび割れ全体がたとえば1ピクセル幅を有するひび割れ細線化画像を作成する(ステップS50、第4のステップ)。
図示するひび割れ検出方法の実施の形態1によれば、まず、入力画像においてひび割れ位置を指定した後に解析を実行することから、ひび割れ画像解析範囲を可及的に縮小することができるため、従来の方法に比して解析領域を大幅に縮小することができ、このことによって解析時間を大幅に短縮することが可能になる。
また、連続解析が可能となり、データの読み込み時間や書き出し時間の短縮も可能になる。
また、解析者が予めひび割れ位置を指定していることから、ひび割れ以外のものが解析対象になり難い。
また、ひび割れ想定線に沿ってラフなパス(ひび割れ想定線の幅よりも数ピクセル幅の大きなパス)を作成してひび割れ位置を指定していることから、パスからひび割れ想定線がはみ出してひび割れと判定されないといった不具合が生じることも解消され、パスがラフゆえにパスの作成も効率的におこなうことができ、パスの作成(ひび割れ位置の特定)に長い時間を必要としない。
さらに、ひび割れ想定線にラフなパスを作成し、ウェーブレット変換した後に細線化処理をおこなうことで、ひび割れをより一層精度よく特定することが可能になる。
(ひび割れ検出方法の実施の形態2)
図8は本発明のひび割れ検出方法の実施の形態2のフロー図である。
ひび割れ検出方法の実施の形態2は、ひび割れ細線化画像を作成する(ステップS50、第4のステップ)まではひび割れ検出方法の実施の形態1と同様であるが、さらに、二種類のノイズ除去処理を実行する方法である。
まず、第4のステップで作成されたひび割れ細線化画像に対し、予め設定されたピクセル数以下の長さのひび割れから分岐した線をひび割れでないノイズと判定して除去するヒゲ削除処理をおこなう(ステップS60)。
5ピクセル程度以下の長さの線がひび割れでないことが本発明者等の経験則上分かっていることに基づき、たとえば、5ピクセル程度の長さを閾値に設定し、この閾値未満の長さの線をヒゲと特定して除去する。
次に、予め設定されたピクセル数以下の孤立した長さの線をひび割れでないノイズと判定して除去する輪郭線追跡処理をおこなう(ステップS70、以上、第5のステップ)。
輪郭線追跡処理は、ある任意の画素(ひび割れと判断されている画素)から出発して、隣接する画素がひび割れ箇所の場合には出発画素と接続し、さらに隣接する画素がひび割れ箇所の場合にはさらに双方を接続し、最終的に出発画素に閉合した場合(例えば、第1画素、第2画素、…、第n-1画素、第n画素、第1画素の順に接続される場合)や、次に繋がるひび割れ箇所が存在しなくなった場合に終了する。この輪郭線追跡処理では、ループ状に閉合するようなひび割れラインや、複数の屈曲部を備えて線状に伸びるひび割れラインなど、適宜のひび割れラインが作成される。この際、繋げられる画素数の最小数を所定の値に設定しておくことで、この設定数以下の画素はすべてひび割れでないとして、画面のひび割れ表示から削除する。
このように二種類のノイズ除去処理をおこなうことで、第1のステップにおいてひび割れ想定線を含むようにラフなパスを作成することによって生じ得るノイズを効果的に除去することができる。たとえば、パスの幅が3ピクセル程度であればノイズ除去処理を必ずしも必要とせず、パスの幅が10ピクセル程度と広幅になった際に上記二種類のノイズ除去処理を必須とするのが好ましい。
次に、特定されたひび割れをひび割れ幅ごとに分類するとともにひび割れ幅ごとに色分けされたひび割れ展開図を作成し、さらに、ひび割れ幅ごとのひび割れ延長ヒストグラムを作成する(ステップS80、第6のステップ)。
(本発明のひび割れ検出方法において10ピクセル幅のパスを作成した際の各ステップにて作成された画像を説明した実施例)
本発明者等は本発明のひび割れ検出方法を実施し、第1のステップにおいて10ピクセル幅のパスを作成した場合における各ステップで作成される画像を作成するとともに、ひび割れの延長ヒストグラムを作成した。ここで、図9は10ピクセル幅のパスを作成した際の各ステップにて作成された画像を示した図であって、(a)はコンクリート表面の撮影画像をコンピュータに入力してなる入力画像を示した図であり、(b)はパス画像を示した図であり、(c)はひび割れ画像を示した図であり、(d)はひび割れ細線化画像を示した図である。なお、図9(d)は詳細には、ひび割れ細線化画像がヒゲ除去処理および輪郭線追跡処理された後の画像を示した図である。また、図10は図9で示すひび割れ検出方法にて作成されたひび割れの延長ヒストグラムの実施例を示した図である。
図9(a)、(d)より、最終的に得られる図9(d)で示すひび割れ細線化画像およびノイズ処理後の画像はいずれも、コンクリート表面の撮影画像上のひび割れを忠実に再現し、精度よく特定していることが分かる。
すなわち、10ピクセル幅のパスを作成した場合では、ノイズ除去処理をおこなうことでひび割れを精度よく特定できることが分かる。
また、図10より、ひび割れ長さの総計は0.34mであった。
(解析時間(コンピュータのCPU処理時間)の検証とその結果)
本発明者等は、本発明のひび割れ検出方法と、従来のひび割れ検出方法において解析に要する所要時間の検証をおこなった。ここで、解析画像サイズは513×391ピクセル、使用したPCスペックはIntel社製のCore i5-4460CPU 3.20GHzである。検証結果を以下の表1に示す。
Figure 0007008408000005
表1において、「ひび割れ画像からひび割れ以外の成分を除去する時間」とは、従来方法では、解析対象の全画像に対してウェーブレット変換を展開してひび割れ画像を作成することから、作成されたひび割れ画像にはコンクリート表面上の汚れやしみといったひび割れ以外の成分が多分に含まれており、このようなひび割れ以外の成分を除去する必要があることから、この作業に要する時間が必要になる。一方、本発明のひび割れ検出方法では、第1のステップにて解析者がひび割れ想定線に沿ってパスを作成し、ひび割れ位置を特定した後にウェーブレット変換を展開してひび割れ画像を作成することから、このようにコンクリート表面上の汚れやしみといったひび割れ以外の成分を除去する作業は不要となる。なお、上記する従来の方法における、コンクリート表面上の汚れやしみといったひび割れ以外の成分を除去する作業と、本発明の方法の実施の形態2における第5のステップのノイズ除去処理は全く異なるものである。
表1より、従来のひび割れ検出方法に要した解析時間が458秒であったのに対して、本発明のひび割れ検出方法に要した解析時間は52秒と、従来の方法の10%程度にまで短縮できることが分かった。
これは、解析者が予めひび割れ位置を指定した後に解析を実行することから、ひび割れ画像解析範囲が縮小され、解析領域が大幅に縮小されたことによって解析時間が短縮されたことと、連続解析が可能となったことでデータの読み込み時間や書き出し時間が短縮されたこと、さらには、ひび割れ想定線に沿ってラフなパスを作成することでこのパスの作成に多大な時間を要しないことによるものである。
以上、本発明の実施の形態を、図面を用いて詳述してきたが、具体的な構成はこの実施形態に限定されるものではなく、本発明の要旨を逸脱しない範囲における設計変更等があっても、それらは本発明に含まれるものである。
1…入力画像、2…広域領域、3…局所領域、31…近傍画素、32…注目画素

Claims (3)

  1. コンクリート表面に生じているひび割れの検出をおこなうひび割れ検出方法であって、
    ひび割れ検出対象であるコンクリート表面の撮影画像をコンピュータに入力して入力画像とし、該入力画像においてひび割れと想定されるひび割れ想定線に沿って、該ひび割れ想定線の左右に1~2ピクセルずつ幅の大きなラフなパスを作成し、該パスに付したひび割れを覆う面を作成してパス画像とする第1のステップ、
    ひび割れの濃度とコンクリート表面の濃度を擬似的に設定し、対比される2つの濃度に対応したウェーブレット係数を算定するとともに、該2つの濃度をそれぞれ変化させた場合のそれぞれのウェーブレット係数を算定してウェーブレット係数テーブルを作成し、前記パス画像におけるひび割れ上でウェーブレット係数を計算する第2のステップ、
    前記ウェーブレット係数テーブル内において、前記コンクリート表面の濃度と仮定する局所領域内の近傍画素の平均濃度と、前記ひび割れの濃度と仮定する注目画素の濃度に対応するウェーブレット係数をウェーブレット係数に関する閾値とし、任意の近傍画素における任意の注目画素のウェーブレット係数が前記閾値よりも大きな場合は該近傍画素における該注目画素をひび割れと判定し、任意の近傍画素における任意の注目画素のウェーブレット係数が前記閾値よりも小さな場合は注目画素をひび割れでないと判定し、前記局所領域および前記注目画素を変化させながら注目画素のウェーブレット係数と前記閾値との比較をおこなうことでひび割れ画像を作成する第3のステップ、
    前記ひび割れ画像に対して細線化処理をおこなってその中心線で構成されるひび割れ細線化画像を作成する第4のステップ、からなるひび割れ検出方法。
  2. 前記第4のステップで作成された前記ひび割れ細線化画像に対し、予め設定されたピクセル数以下の長さのひび割れから分岐した線をひび割れでないノイズと判定して除去するヒゲ削除処理をおこない、さらに、予め設定されたピクセル数以下の孤立した長さの線をひび割れでないノイズと判定して除去する輪郭線追跡処理をおこなう第5のステップをさらに含む、請求項1に記載のひび割れ検出方法。
  3. 前記第5のステップに続いて、特定されたひび割れをひび割れ幅ごとに分類するとともにひび割れ幅ごとに色分けされたひび割れ展開図を作成し、かつ、ひび割れ幅ごとのひび割れ延長ヒストグラムを作成する第6のステップをさらに備えている、請求項2に記載のひび割れ検出方法。
JP2017020499A 2017-02-07 2017-02-07 ひび割れ検出方法 Active JP7008408B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2017020499A JP7008408B2 (ja) 2017-02-07 2017-02-07 ひび割れ検出方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2017020499A JP7008408B2 (ja) 2017-02-07 2017-02-07 ひび割れ検出方法

Publications (2)

Publication Number Publication Date
JP2018128315A JP2018128315A (ja) 2018-08-16
JP7008408B2 true JP7008408B2 (ja) 2022-01-25

Family

ID=63172404

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2017020499A Active JP7008408B2 (ja) 2017-02-07 2017-02-07 ひび割れ検出方法

Country Status (1)

Country Link
JP (1) JP7008408B2 (ja)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110929757B (zh) * 2019-10-24 2023-04-18 东南大学 一种沥青路面裂缝类型快速分类方法
CN111008956B (zh) * 2019-11-13 2024-06-28 武汉工程大学 基于图像处理的梁底裂纹检测方法、系统、装置及介质
KR102647135B1 (ko) * 2022-10-31 2024-03-15 이관훈 인공지능 기반 객체 탐지 알고리즘을 이용하는 시공 현장을 위한 실시간 균열 탐지 시스템 및 이를 위한 동작 방법

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002310920A (ja) 2001-04-19 2002-10-23 Keisoku Kensa Kk コンクリート壁のひび割れ検出方法およびその装置
JP2006162477A (ja) 2004-12-08 2006-06-22 Ohbayashi Corp 構造物表面のひび割れ検出方法
JP2010121992A (ja) 2008-11-18 2010-06-03 Taisei Corp ひび割れ検出方法
JP2013195074A (ja) 2012-03-15 2013-09-30 Keisoku Res Consultant:Kk クラック検出方法、クラック表示装置
CN103745477A (zh) 2014-01-22 2014-04-23 汕头大学 基于机器视觉的建筑结构裂缝检测及修复方法
JP2016170042A (ja) 2015-03-12 2016-09-23 阪神高速技術株式会社 ひび割れ検出方法及びひび割れ検出プログラム

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2554516B2 (ja) * 1988-01-14 1996-11-13 東京電力株式会社 ひび割れ自動計測装置
JPH087157B2 (ja) * 1988-01-30 1996-01-29 鹿島建設株式会社 コンクリート表面のひび割れ測定方法
JPH04120495A (ja) * 1990-09-12 1992-04-21 Taisei Corp トンネル前方地質判定システム
JP2815263B2 (ja) * 1992-06-04 1998-10-27 松下電器産業株式会社 細線画像整形方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002310920A (ja) 2001-04-19 2002-10-23 Keisoku Kensa Kk コンクリート壁のひび割れ検出方法およびその装置
JP2006162477A (ja) 2004-12-08 2006-06-22 Ohbayashi Corp 構造物表面のひび割れ検出方法
JP2010121992A (ja) 2008-11-18 2010-06-03 Taisei Corp ひび割れ検出方法
JP2013195074A (ja) 2012-03-15 2013-09-30 Keisoku Res Consultant:Kk クラック検出方法、クラック表示装置
CN103745477A (zh) 2014-01-22 2014-04-23 汕头大学 基于机器视觉的建筑结构裂缝检测及修复方法
JP2016170042A (ja) 2015-03-12 2016-09-23 阪神高速技術株式会社 ひび割れ検出方法及びひび割れ検出プログラム

Also Published As

Publication number Publication date
JP2018128315A (ja) 2018-08-16

Similar Documents

Publication Publication Date Title
JP6823486B2 (ja) ひび割れ検出方法
JP5385593B2 (ja) ひび割れ検出方法
JP4006007B2 (ja) ひび割れ検出方法
JP6099479B2 (ja) ひび割れ検出方法
JP7008409B2 (ja) ひび割れ検出方法
JP5812705B2 (ja) ひび割れ検出方法
JP5852919B2 (ja) ひび割れ検出方法
JP5421192B2 (ja) ひび割れ検出方法
JP6620477B2 (ja) コンクリートのひび割れ検出方法及び検出プログラム
JP5705711B2 (ja) ひび割れ検出方法
JP6422198B2 (ja) 画像処理装置、画像処理方法、及び画像処理プログラム
CN113109368A (zh) 玻璃裂纹检测方法、装置、设备及介质
JP4870016B2 (ja) ひび割れ検出方法
JP4908440B2 (ja) 画像処理装置及び方法
CN110706224B (zh) 基于暗场图像的光学元件弱划痕检测方法、系统和装置
JP2005037378A (ja) 奥行計測方法と奥行計測装置
JP7008408B2 (ja) ひび割れ検出方法
CN118351097A (zh) 一种胶路质量检测方法、装置、电子设备及存储介质
JP6637823B2 (ja) ひび割れ検出方法
CN118657740A (zh) 一种五金件的加工质检方法、电子设备及系统
JP5157575B2 (ja) 欠陥検出方法
Nguyen et al. A novel automatic concrete surface crack identification using isotropic undecimated wavelet transform
Prabha et al. Defect detection of industrial products using image segmentation and saliency
JP6909090B2 (ja) コンクリート表面上の損傷部の定量評価方法、定量評価装置、定量評価システム及びプログラム
JP4980739B2 (ja) ひび割れ検出方法

Legal Events

Date Code Title Description
RD03 Notification of appointment of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7423

Effective date: 20170704

RD04 Notification of resignation of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7424

Effective date: 20170824

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20200128

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20201125

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20201201

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20210112

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20210629

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20220104

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20220111

R150 Certificate of patent or registration of utility model

Ref document number: 7008408

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150