CN101691994B - 一种隧道裂缝最大宽度自动定位和检测方法 - Google Patents

一种隧道裂缝最大宽度自动定位和检测方法 Download PDF

Info

Publication number
CN101691994B
CN101691994B CN2009101535277A CN200910153527A CN101691994B CN 101691994 B CN101691994 B CN 101691994B CN 2009101535277 A CN2009101535277 A CN 2009101535277A CN 200910153527 A CN200910153527 A CN 200910153527A CN 101691994 B CN101691994 B CN 101691994B
Authority
CN
China
Prior art keywords
picture
value
pixel
edge
crack
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.)
Expired - Fee Related
Application number
CN2009101535277A
Other languages
English (en)
Other versions
CN101691994A (zh
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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN2009101535277A priority Critical patent/CN101691994B/zh
Publication of CN101691994A publication Critical patent/CN101691994A/zh
Application granted granted Critical
Publication of CN101691994B publication Critical patent/CN101691994B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

本发明公开了一种隧道裂缝最大宽度自动定位和检测方法。包括如下步骤:1)对摄像头获取到的隧道裂缝图片进行像素级别的放大处理,并通过边缘检测算法来获得裂缝的矢量化边缘,对得到的图片进行适应性阈值的二值化操作,并对结果进行腐蚀和膨胀的去噪处理;2)在步骤1)所得到的低噪音二值化图片基础上,进行欧氏距离变换操作,所得到的结果和步骤1)所得到的矢量化边缘结果进行匹配寻找距离最大的点作为定位点,并以它的值作为裂缝的最大宽度值;本发明使用了一系列计算机图像处理知识来自动定位和检测隧道裂缝的最大宽度,获得了良好的效果,具有较高的准确度和精度。

Description

一种隧道裂缝最大宽度自动定位和检测方法
技术领域
[0001] 本发明涉及一种隧道裂缝最大宽度自动定位和检测方法,尤其涉及基于摄像头获 得近距离隧道裂缝图片条件下,运用计算机图形图像学的知识和技术来获得对于裂缝最大 宽度点的自动定位和检测,属于图形图像应用技术的领域。
背景技术
[0002] 混凝土裂缝的存在和影响是导致目前高速公路施工建设特别是隧道安全的最大 隐患。正因为裂缝的存在,使空气中的二氧化碳极易渗透到混凝土内部,在潮湿的环境作用 下,二氧化碳与水泥中的氢氧化钙、硅酸三钙、硅酸二钙相互作用并转化成碳酸盐,中和水 泥的基本碱性,从而使混凝土的碱度降低,导致钢筋的纯化膜遭受破坏而引起锈蚀;同时由 于混凝土的碳化会加剧混凝土收缩开裂,从而导致桥梁结构破坏。
[0003] 随着高速公路中隧道施工和保养技术手段的翻新,在不影响影响桥梁结构外观的 正常使用和耐久性的基础上,发展出了如下几种隧道裂缝及其宽度的无损检测和监测方 法:
[0004] 1.超声波检测:超声波法用于非破损性检测,就是以超声波为媒介,获得物体内 部信息的一种方法,目前超声法已应用于医疗诊断、钢材探伤、鱼群探测等许多领域。在这 些领域里,由于组成颗粒小密度大,密度分部也很均勻,所以声波能很好地传播,对其内部 缺陷及其位置等都能准确地检测出来。掌握混凝土表面产生的裂缝深度,对耐久性诊断和 研究修补加固对策有重要意义。测定裂缝深度,基本上都是将发射探头和接收探头,布置在 混凝土同一面上的裂缝附近,但由于所选用的波形种类纵波、横波及表面波和声学参数声 速、频率、相位等的不同,已有许多种具体方法。
[0005] 2.声发射检测法:声发射检测法也是利用弹性波进行声学检测的具体检测方法 检测裂缝,和其他方法最大的不同是只能检测正在发生的裂缝,不能检测已发生的旧裂缝, 对正在发生的裂缝可检测裂缝发生的位置声发射源定位,裂缝的大小,扩展情况和种类,以 及裂缝的深度等。
[0006] 3.摄影检测法:摄影检测法主要用作调查混凝土表面的裂缝摄影法包括普通照 相机、摄像机、放射线、红外线摄影等进行检测。由于所依赖的硬件设备大量普及带来的成 本低廉和部署简便,已经成为目前商用部署的主流方法。
[0007] 4.传感仪器监测:利用埋设在混凝土中的仪器进行裂缝监测,常规技术是利用卡 尔逊式或弦式测缝计,其控制范围仅0. 12毫米,属点式检测,由于裂缝出现的空间随机性, 因此往往漏检,为了及时无遗漏地监测裂缝,必须实施大范围的、连续、分布式监测,即所谓 全分布监测。
[0008] 5.光纤传感网络监测:在目前竞相开发的结构监测高科技领域里,光纤传感以其 独特优势居于科研开发的中心地位,它灵巧、精度高、抗电磁干扰,且可靠耐久,易于光纤传 输组成自动化遥测系统。光纤传感应用于结构工程监测始于上世纪九十年代初,如航空航 天器、桥梁等的温度、振动、应变检测等裂缝的发生可以用埋设在混凝土中光纤的光强变化监测,而裂缝的定位可用多模光纤在裂缝处的光强突然下降或诊断完成,通过衰减曲线上 的裂缝损耗突变点,可以准确地确定裂缝的位置,针对混凝土裂缝检测的特点,研制出基于 光时域反射技术的光纤裂缝传感网络,可实现桥梁混凝土结构的分布检测,凡裂缝与光纤 传感网络相交,均可感知,并可定宽、定位、定向。
发明内容
[0009] 本发明的目的是克服现有技术的不足,提供一种隧道裂缝最大宽度自动定位和检 测方法。
[0010] 隧道裂缝最大宽度自动定位和检测方法包括如下步骤:
[0011] 1)对摄像头获取到的隧道裂缝图片进行像素级别的放大处理,并通过边缘检测算 法来获得裂缝的矢量化边缘,对得到的图片进行适应性阈值的二值化操作,并对结果进行 腐蚀和膨胀的去噪处理;
[0012] 2)在步骤1)所得到的低噪音二值化图片基础上,进行欧氏距离变换操作,所得到 的结果和步骤1)所得到的矢量化边缘结果进行匹配寻找距离最大的点作为定位点,并以 它的值作为裂缝的最大宽度值;
[0013] 所述的对摄像头获取到的隧道裂缝图片进行像素级别的放大处理,并通过边缘检 测算法来获得裂缝的矢量化边缘,对得到的图片进行适应性阈值的二值化操作,并对结果 进行腐蚀和膨胀的去噪处理步骤包括:
[0014] a)对于所获取的隧道裂缝图片,获得它对应的像素点信息集合为{Pn,m},其中η为 图片的横坐标像素点数量,m为图片的纵坐标像素点数量,这里的Py= (Ri,,., GiijjBi,^是 用来描述像素点信息的一组向量,包括了红色通道的值氏,」,绿色通道的值Gm,蓝色通道的 值Biij,这些值可以通过对裂缝图片的读取获得;
[0015] 按照输入的图片尺寸大小,按照图片的格式新建一个新的相同格式的空白图片, 其尺寸为原始图片尺寸大小的两倍,其对应的像素点信息集合为{p2n,2m’ },这里的P1,/ = (R1,/,Gi,/,B1,/ }是用来描述新的图片像素点信息的一组向量,定义方式如Pi,」。其中Pu 和Plk’之间存在着一一对应的关系,公式如下:
[0016] P2*i,2*j(R,G,B),= Pij j(R, G,B);
[0017] P2»(i+1),州(R,G,B),= Pi+1, j (R,G,B);
[0018] Ρ2„,2*(」+1) (R,G,B),= Pi, J+1 (R,G,B);
[0019] P2»(i+1),2*(J+1) (R,G,B),= Pi+1, J+1 (R,G,B);
[0020] P2*i+1,2*j (R,G,B),= (Pi, j (R,G,B) +Pi+1,』(R,G,B)) /2 ; 1
[0021] P2»ij2»J+1 (R, G,B),=們,」0?,G,B)+PijJ+1 (R, G,B))/2 ;
[0022] P2»(i+1) ,2*j+1 (R,G,B),= (Pi+1,」(R,G,B) +Pi+1, J+1 (R,G,B)) /2 ;
[0023] P2*i+1,2*(j+1) (R,G,B),= (Pi, J+1 (R,G,B) +Pi+1, J+1 (R,G,B)) /2 ;
[0024] P2»i+lj2»J+1 (R,G, B)' = (Pi, j (R,G,B) +Pi+1, j (R,G,B) +Pi, J+1 (R,G,B) +Pi+1, J+1 (R,G, B))/4;
[0025] 其中R,G和B三个通道的值都是按照公式1所示进行获得和赋值操作;
[0026] b)对放大得到的图片,调用以动态链接库形式存在的Matlab2007b中的edge函 数,按照 “sobel” “prewitt” “roberts” “log” “zerocross” “canny” 方法及其所对应的参数进行边缘检测;对于获得的已经矢量化边缘的图片,按照整个图片的{P2n,2m’ }像素统 计信息,对于P1,/的欧式距离按照阈值Ther进行过滤来完成二值化,公式如下:
[0027]
Figure CN101691994BD00061
[0028]
[0029] 其中IlP1Y Il为计算R,G和B三个通道分量的欧氏距离,变量Hpct为经验性的 数据,预先指定来控制Ther的取值;接下来对于得到的二值化图片,进行腐蚀和膨胀操作, 消除由于二值化带来的噪声,获得完整的二值化后的矢量化边缘图片。
[0030] 所述的在步骤1)所得到的低噪音二值化图片基础上,进行欧氏距离变换操作,所 得到的结果和步骤1)所得到的矢量化边缘结果进行匹配寻找距离最大的点作为定位点, 并以它的值作为裂缝的最大宽度值步骤包括:
[0031] c)对在步骤1)得到的低噪音二值化图片基础上,进行欧氏距离变换操作,欧式距 离变换的公式如下:
[0032] D(ρ) := min{d(p, q) |q e 0C} = min{d(p, q) 11 (q) = 0} 3
[0033] 其中D(p)表示在新生成的图D中每一点像素ρ的值,(T表示在图像矩阵中的背 景,I(q)表示图像I中像素点q的值,公式表明了欧氏距离变换是一种变换,它生成一个图 D,其中的每一点像素ρ的值为它距离背景集合的最小值;
[0034] 将所获得的欧氏距离变换结果图D和在步骤1)中所得的矢量化边缘图片进行匹 配处理,只有当像素AP1,/处于矢量化边缘之间,并且D^k')取到最大值的时候,将此 P1,/取为定位点,并且以该点P1,/的欧氏距离的值为该点的像素宽度,经过像素-实际尺 寸之间的比例换算关系,获得真实隧道裂缝的宽度。
[0035] 本发明使用了一系列计算机图像处理知识来自动定位和检测隧道裂缝的最大宽 度,获得了良好的效果,具有较高的准确度和精度。
附图说明
[0036] 下面结合附图和具体实施例对本发明作进一步的说明。
[0037] 图1为本发明所述的结果图;
具体实施方式
[0038] 隧道裂缝最大宽度自动定位和检测方法包括如下步骤:
[0039] 1)对摄像头获取到的隧道裂缝图片进行像素级别的放大处理,并通过边缘检测算 法来获得裂缝的矢量化边缘,对得到的图片进行适应性阈值的二值化操作,并对结果进行 腐蚀和膨胀的去噪处理;
[0040] 2)在步骤1)所得到的低噪音二值化图片基础上,进行欧氏距离变换操作,所得到 的结果和步骤1)所得到的矢量化边缘结果进行匹配寻找距离最大的点作为定位点,并以 它的值作为裂缝的最大宽度值;
[0041] 所述的对摄像头获取到的隧道裂缝图片进行像素级别的放大处理,并通过边缘检测算法来获得裂缝的矢量化边缘,对得到的图片进行适应性阈值的二值化操作,并对结果 进行腐蚀和膨胀的去噪处理步骤包括:
[0042] a)对于所获取的隧道裂缝图片,获得它对应的像素点信息集合为{Pn,m},其中η为 图片的横坐标像素点数量,m为图片的纵坐标像素点数量,这里的Py= (Ri,,., GiijjBi,^是 用来描述像素点信息的一组向量,包括了红色通道的值氏,」,绿色通道的值Gm,蓝色通道的 值Biij,这些值可以通过对裂缝图片的读取获得;
[0043] 按照输入的图片尺寸大小,按照图片的格式新建一个新的相同格式的空白图片, 其尺寸为原始图片尺寸大小的两倍,其对应的像素点信息集合为{p2n,2m’ },这里的P1,/ = (R1,/,G1,/,B1,/ }是用来描述新的图片像素点信息的一组向量,定义方式如Pi,」。其中Pu 和Plk’之间存在着一一对应的关系,公式如下:
[0044] P2„,2<y (R, G, B)' = Pi,」(R,G,B);
[0045] P2»(i+1) j2<iJ (R,G, B) ' = Pi+1, j (R,G,B);
[0046] Ρ2„,2*(」+1) (R,G, B) ' = Pi, J+1 (R,G,B);
[0047] P2»(i+1),2*(J+1) (R,G,B),= Pi+1, J+1 (R,G,B);
[0048] P2*i+1,2*j (R,G,B),= (Pi,』(R,G,B) +Pi+1,』(R,G,B)) /2 ; 1
[0049] P測,2<y+1 (R,G, B) ' = (Pi, j (R,G,B) +Pi, J+1 (R,G,B)) /2 ;
[0050] P2»(i+1) ,2*j+1 (R,G,B),= (Pi+1,」(R,G,B) +Pi+1, J+1 (R,G,B)) /2 ; [0051 ] P2»i+1,2»(J+1) (R,G,B),= (Pi, J+1 (R,G,B) +Pi+1, J+1 (R,G,B)) /2 ;
[0052] P2»i+lj2»J+1 (R,G, B)' = (Pi, j (R,G,B) +Pi+1, j (R,G,B) +Pi, J+1 (R,G,B) +Pi+1, J+1 (R,G, B))/4;
[0053] 其中R,G和B三个通道的值都是按照公式1所示进行获得和赋值操作;
[0054] b)对放大得到的图片,调用以动态链接库形式存在的Matlab2007b中的edge函 数,按照 “sobel” “prewitt” “roberts” “log” “zerocross” “canny” 方法及其所对应的 参数进行边缘检测;对于获得的已经矢量化边缘的图片,按照整个图片的{P2n,2m’ }像素统 计信息,对于P1,/的欧式距离按照阈值Ther进行过滤来完成二值化,公式如下:
Figure CN101691994BD00071
[0057] 其中IlP1Y Il为计算R,G和B三个通道分量的欧氏距离,变量Hpct为经验性的 数据,预先指定来控制Ther的取值;接下来对于得到的二值化图片,进行腐蚀和膨胀操作, 消除由于二值化带来的噪声,获得完整的二值化后的矢量化边缘图片。
[0058] 所述的在步骤1)所得到的低噪音二值化图片基础上,进行欧氏距离变换操作,所 得到的结果和步骤1)所得到的矢量化边缘结果进行匹配寻找距离最大的点作为定位点, 并以它的值作为裂缝的最大宽度值步骤包括:
[0059] c)对在步骤1)得到的低噪音二值化图片基础上,进行欧氏距离变换操作,欧式距 离变换的公式如下:
[0060] D (ρ) : = min {d (ρ, q) | q e 0C} = min {d (ρ, q) 11 (q) = 0} 3[0061] 其中D(p)表示在新生成的图D中每一点像素ρ的值,(T表示在图像矩阵中的背 景,I(q)表示图像I中像素点q的值,公式表明了欧氏距离变换是一种变换,它生成一个图 D,其中的每一点像素ρ的值为它距离背景集合的最小值;
[0062] 将所获得的欧氏距离变换结果图D和在步骤1)中所得的矢量化边缘图片进行匹 配处理,只有当像素AP1,/处于矢量化边缘之间,并且D^k')取到最大值的时候,将此 P1,/取为定位点,并且以该点P1,/的欧氏距离的值为该点的像素宽度,经过像素-实际尺 寸之间的比例换算关系,获得真实隧道裂缝的宽度。
[0063] 实施例
[0064] (1)对于通过紧附在隧道裂缝上方高精度摄像头所获取的隧道裂缝图片,获得它 对应的像素点信息集合为{pn, J,其中η为图片的横坐标像素点数量,m为图片的纵坐标像 素点数量,在本实施例中η = 800,m= 700。这里的Piij= (RijjiGijjjBijjI是用来描述像 素点信息的一组向量,包括了红色通道的值I^j,绿色通道的值Gi,」,蓝色通道的值4,」,这些 值可以通过对裂缝图片的读取获得;
[0065] 按照输入的图片尺寸大小,按照图片的格式新建一个新的相同格式的空白图片, 其尺寸为原始图片尺寸大小的两倍,其对应的像素点信息集合为{P2n,2m’}。这里的P1,/ = (R1,/,G1,/,B1,/ }是用来描述新的图片像素点信息的一组向量,定义方式如Pi,」。其中Pu 和Plk’之间存在着一一对应的关系,公式如下:
[0066] P2„,2<y (R, G, B)' = Pi,」(R,G,B);
[0067] P2»(i+1) j2<iJ (R,G, B) ' = Pi+1, j (R,G,B);
[0068] Ρ2„,2*(」+1) (R,G, B) ' = Pi, J+1 (R,G,B);
[0069] P2»(i+1),2*(J+1) (R,G,B),= Pi+1, J+1 (R,G,B);
[0070] P2*i+1,2*j (R,G,B),= (Pi, j (R,G,B) +Pi+1,』(R,G,B)) /2 ; 1
[0071] P2„,2*j+1(R,G,B),= (Ρ^^Κ,Ο, B)+PijJ+1(R,G, B))/2 ;
[0072] P2»(i+1) ,2*j+1 (R,G,B),= (Pi+1,」(R,G,B) +Pi+1, J+1 (R,G,B)) /2 ;
[0073] P2*i+1,2*(j+1) (R,G,B),= (Pi, J+1 (R,G,B) +Pi+1, J+1 (R,G,B)) /2 ;
[0074] P2»i+lj2»J+1 (R,G, B)' = (Pi, j (R,G,B) +Pi+1, j (R,G,B) +Pi, J+1 (R,G,B) +Pi+1, J+1 (R,G, B))/4;
[0075] 其中R,G和B三个通道的值都是按照公式1所示进行获得和赋值操作;
[0076] (2)对放大得到的图片,调用以动态链接库形式存在的Matlab2007b中的edge函 数,按照 “sobel” “prewitt” “roberts” “log” “zerocross” “canny” 方法及其所对应的 参数进行边缘检测;
[0077] 其中Sobel方法对应的数学表达式为:
[0078] g = [G2X+G2J
[0079] = {[(z7 + 2z8 + Z9) - (ζ, + Iz2 + z3)f
[0080] +[(z3+2Z6+Z9)-(Z1+2Z4+Z7)]2}^ 2
[0081] 其中(^和Gy为一阶导数值,Z1〜%为中心点Z5的领域点,Sobel公式用掩模来数 字化的近似了一阶倒数和Gy,若中心点的g值超过某个指定的阈值,则认为该位置的一 个像素是一个边缘像素。其中I^ewitt方法对应的数学表达式为
Figure CN101691994BD00091
[0086] 其中和Gy为一阶导数值,Z1〜%为中心点Z5的领域点,Prewitt公式用掩模来 数字化的近似了一阶倒数和Gy,同Sobel方法相比,在计算上要简单一些,但是由于没有 平滑系数的作用,容易产生一些噪声。
Figure CN101691994BD00092
其中Roberts方法对应的数学表达式为 g = [G2x+G2J
Figure CN101691994BD00093
[0090] 其中和Gy为一阶导数值,Z1〜%为中心点Z5的领域点,Roberts公式用掩模来 数字化的近似了一阶倒数和Gy,同Sobel方法相比,它是非对称的,而且不能检测诸如45 度倍数的边缘,但是由于它的简单快速,还是经常用于硬件实现中。
[0091] 其中LoG方法对应的数学表达式为:
[0092]
Figure CN101691994BD00094
[0093] 其中r2 = x2+y2,σ是标准偏差,LoG对图像卷积会有两个效果:使图像变得更加 平滑;计算拉普拉斯算子,以便产生双边缘图像。然后,定位边缘就是找到两个边缘之间的 零交叉。
[0094] 对于获得的已经矢量化边缘的图片,按照整个图片的{&„》’}像素统计信息,对于 P1,/的欧式距离按照阈值Ther进行过滤来完成二值化,公式如下:
‘P1 k{R = 255,G = 0,忍=0>/化 t 'J > Ther
[0095]
Plk{R = O5G = 0,5 = J ^ Ther 6
2*n,2*m
[0096]
[0097]
Σ Il^
Ther = Hn
f^r 4*rt*m
其中IlP1Y Il为计算R,G和B三个通道分量的欧氏距离,变量Hpct为经验性的 数据,预先指定来控制Ther的取值,在本实施例中,该变量设定为0. 67〜0. 75之间;接下 来对于得到的二值化图片,进行腐蚀和膨胀操作,消除由于二值化带来的噪声,获得完整的 二值化后的矢量化边缘图片。
[0098] (3)对在步骤1)得到的低噪音二值化图片基础上,进行欧氏距离变换操作,欧式 距离变换的公式如下:
[0099] D (ρ) : = min {d (ρ, q) | q e 0C} = min {d (ρ, q) 11 (q) = 0} 7
[0100] 其中D(p)表示在新生成的图D中每一点像素ρ的值,(T表示在图像矩阵中的背 景,I(q)表示图像I中像素点q的值。公式表明了欧氏距离变换是一种变换,它生成一个
9图D,其中的每一点像素ρ的值为它距离背景集合的最小值。
[0101] 将所获得的欧氏距离变换结果图D和在步骤1)中所得的矢量化边缘图片进行匹 配处理,只有当像素AP1,/处于矢量化边缘之间,并且D^k')取到最大值的时候,将此 Pu’取为定位点,并且以该点Pu’的欧氏距离的值为该点的像素宽度,经过像素-实际尺 寸之间的比例换算关系,获得真实隧道裂缝的宽度,本方法的模块图和最终结果图见图1。

Claims (2)

1. 一种隧道裂缝最大宽度自动定位和检测方法,其特征在于包括如下步骤:1)对摄像头获取到的隧道裂缝图片进行像素级别的放大处理,并通过边缘检测算法来 获得裂缝的矢量化边缘,对得到的图片进行适应性阈值的二值化操作,并对结果进行腐蚀 和膨胀的去噪处理;2)在步骤1)所得到的低噪音二值化图片基础上,进行欧氏距离变换操作,所得到的结 果和步骤1)所得到的矢量化边缘结果进行匹配寻找距离最大的点作为定位点,并以它的 值作为裂缝的最大宽度值;所述的对摄像头获取到的隧道裂缝图片进行像素级别的放大处理,并通过边缘检测算 法来获得裂缝的矢量化边缘,对得到的图片进行适应性阈值的二值化操作,并对结果进行 腐蚀和膨胀的去噪处理步骤包括:a)对于所获取的隧道裂缝图片,获得它对应的像素点信息集合为{Pn,m},其中η为图片 的横坐标像素点数量,m为图片的纵坐标像素点数量,这里的Piij= {Ry,Giij,BiJ是用来 描述像素点信息的一组向量,包括了红色通道的值氏,」,绿色通道的值Gm,蓝色通道的值Bi, 」,这些值可以通过对裂缝图片的读取获得;按照输入的图片尺寸大小,按照图片的格式新建一个新的相同格式的空白图片,其尺 寸为原始图片尺寸大小的两倍,其对应的像素点信息集合为{P2n,2m’ },这里的P1,/ = IR1, k’,G1,/, B1,/ }是用来描述新的图片像素点信息的一组向量,定义方式SPi,」,其中Pij和 Plk’之间存在着一一对应的关系,公式如下: P2*i,2*j(R,G,B),= Pi, j (Rj G,B); P2*(i+i),2*j (R,G,B),= Pi+1,j (R,G,B); P2*i,2*(j+i) (R,G,B),= Pi, j+i (R,G,B); P2*(i+l),2*(j+l) (R,G,B),= Pi+1, j+1 (R? G,B); P2*i+1,2*j (R,G,B),= (Pi, j (R,G,B) +Pi+1, j (R,G,B)) /2 ; 1P2*ij2*J+1 (R,G, B)' = (Pijj(R,G,B) +Pi,J+1 (R,G,B))/2 ; P2*(i+i),2*j+i (R,G,B),= (Pi+i, j (R' G, B) +Pi+lj J+1 (R7 G, B)) /2 ; P2*i+i,2*(j+i) (R,G, B)' = (Pi, j+i (R» G, B) +Pi+lj J+1 (R7 G, B)) /2 ;P2*i+1,2*J+1 (R,G,B)' = (Pi,」(R,G,B) +Pi+ljJ (R,G,B) +PijJ+1 (R,G,B) +Pi+1, j+1 (R,G,B))/4 ; 其中R,G和B三个通道的值都是按照公式1所示进行获得和赋值操作;b)对放大得到的图片,调用以动态链接库形式存在的Matlab2007b中的edge函数,按 照“sobel” “prewitt” “roberts” “log” “zerocross” “canny”方法及其所对应的参数进 行边缘检测;对于获得的已经矢量化边缘的图片,按照整个图片的{P2n,2m’ }像素统计信息, 对于P1,/的欧式距离按照阈值Ther进行过滤来完成二值化,公式如下:P1AR = 255,G = = 0)if\\PIJc 'J > Ther ^ ‘P,k(R = 0,G = 0,B = 0>/||^ ’I -< Ther2*n,2*mΣ MlTher 二 Hper * -其中I Ip1^ 11为计算r,G和Β三个通道分量的欧氏距离,变量Hpct为经验性的数据,预先指定来控制Ther的取值;接下来对于得到的二值化图片,进行腐蚀和膨胀操作,消除 由于二值化带来的噪声,获得完整的二值化后的矢量化边缘图片。
2.根据权利要求1所述的一种隧道裂缝最大宽度自动定位和检测方法,其特征在于: 所述的在步骤1)所得到的低噪音二值化图片基础上,进行欧氏距离变换操作,所得到的结 果和步骤1)所得到的矢量化边缘结果进行匹配寻找距离最大的点作为定位点,并以它的 值作为裂缝的最大宽度值步骤包括:c)对在步骤1)得到的低噪音二值化图片基础上,进行欧氏距离变换操作,欧式距离变 换的公式如下:D(p) := min{d(p, q) q e 0C} = min{d(p, q) I (q) = 0} 3其中D(p)表示在新生成的图D中每一点像素ρ的值,(f表示在图像矩阵中的背景,I (q) 表示图像I中像素点q的值,公式表明了欧氏距离变换是一种变换,它生成一个图D,其中的 每一点像素P的值为它距离背景集合的最小值;将所获得的欧氏距离变换结果图D和在步骤1)中所得的矢量化边缘图片进行匹配处 理,只有当像素点P1,/处于矢量化边缘之间,并且D^k')取到最大值的时候,将此P1,/ 取为定位点,并且以该点P1,/的欧氏距离的值为该点的像素宽度,经过像素-实际尺寸之 间的比例换算关系,获得真实隧道裂缝的宽度。
CN2009101535277A 2009-09-30 2009-09-30 一种隧道裂缝最大宽度自动定位和检测方法 Expired - Fee Related CN101691994B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2009101535277A CN101691994B (zh) 2009-09-30 2009-09-30 一种隧道裂缝最大宽度自动定位和检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2009101535277A CN101691994B (zh) 2009-09-30 2009-09-30 一种隧道裂缝最大宽度自动定位和检测方法

Publications (2)

Publication Number Publication Date
CN101691994A CN101691994A (zh) 2010-04-07
CN101691994B true CN101691994B (zh) 2011-06-01

Family

ID=42080702

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2009101535277A Expired - Fee Related CN101691994B (zh) 2009-09-30 2009-09-30 一种隧道裂缝最大宽度自动定位和检测方法

Country Status (1)

Country Link
CN (1) CN101691994B (zh)

Families Citing this family (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102735169B (zh) * 2011-03-31 2015-06-03 上海西门子医疗器械有限公司 一种计算剂量曲线带宽度的系统和方法
CN102521802A (zh) * 2011-11-28 2012-06-27 广东省科学院自动化工程研制中心 一种数学形态学和LoG算子结合的边缘检测算法
CN102636125A (zh) * 2012-04-18 2012-08-15 长沙金码高科技实业有限公司 表面裂缝宽度检测仪
CN103061237B (zh) * 2013-01-31 2014-12-24 长安大学 一种水泥混凝土路面错台量的检测算法
CN103486971B (zh) * 2013-08-14 2016-07-13 北京交通大学 一种地铁隧道裂缝宽度检测及校正算法
CN103440657B (zh) * 2013-08-27 2015-11-04 武汉大学 一种在线道路裂缝筛查方法
CN103942783A (zh) * 2014-04-01 2014-07-23 同济大学 一种岩体裂隙信息的计算方法
CN103955923B (zh) * 2014-04-18 2017-03-01 南京理工大学 一种快速的基于图像的路面病害检测方法
CN104330023B (zh) * 2014-10-15 2017-02-08 浙江大学 混凝土表面初始开裂信息的采集系统与识别方法
CN105741285B (zh) * 2016-01-28 2018-09-04 北京工业大学 用于隧道病害的图像处理装置
CN106446905B (zh) * 2016-07-25 2019-05-28 广东省智能制造研究所 渗透算法和自适应Canny算法相融合的表面裂纹纹理的提取方法
CN106225702A (zh) * 2016-08-03 2016-12-14 武汉科技大学 裂缝宽度检测装置和方法
CN106504246B (zh) * 2016-11-08 2019-04-30 太原科技大学 隧道裂缝检测的图像处理方法
CN106651893A (zh) * 2016-12-23 2017-05-10 贵州电网有限责任公司电力科学研究院 一种基于边缘检测的墙体裂缝识别方法
CN108007355B (zh) * 2017-10-20 2019-06-18 西安电子科技大学 基于图像距离变换的混凝土裂缝检测方法
CN108596869B (zh) * 2018-01-08 2020-06-16 东北大学 一种隧道衬砌裂缝快速检测方法
CN109820537A (zh) * 2018-12-05 2019-05-31 余姚市华耀工具科技有限公司 结石分布探测平台
CN110057308A (zh) * 2019-04-30 2019-07-26 北京理工大学 一种形变监测方法及系统

Also Published As

Publication number Publication date
CN101691994A (zh) 2010-04-07

Similar Documents

Publication Publication Date Title
CN101691994B (zh) 一种隧道裂缝最大宽度自动定位和检测方法
Kudela et al. Identification of cracks in thin-walled structures by means of wavenumber filtering
Dobson et al. Independent component analysis for improved defect detection in guided wave monitoring
CN105549016A (zh) 超声矩阵检验
Castellini et al. The laser doppler vibrometer as an instrument for nonintrusive diagnostic of works of art: application to fresco paintings
CA2460901C (en) System and method for inspecting an industrial furnace or the like
CN110231409B (zh) 一种地下管道损伤的检测方法及系统
US6205859B1 (en) Method for improving defect detectability with magnetostrictive sensors for piping inspection
US20130088724A1 (en) Method and apparatus for the inspection of sandwich structures using laser-induced resonant frequencies
Liu et al. Locating and imaging contact delamination based on chaotic detection of nonlinear Lamb waves
Shi et al. Non-destructive testing of full-length bonded rock bolts based on HHT signal analysis
CN101479597B (zh) 检测和表征介质的声学方法和装置
Chong et al. Statistical threshold determination method through noise map generation for two dimensional amplitude and time-of-flight mapping of guided waves
Castellini et al. On field validation of non-invasive laser scanning vibrometer measurement of damaged frescoes: experiments on large walls artificially aged
Liu et al. Fuzzy pattern recognition of impact acoustic signals for nondestructive evaluation
US10883969B2 (en) System and method for inspecting a structure with coda acoustic waves
Dackermann et al. A new non-destructive testing system based on narrow-band frequency excitation for the condition assessment of pole structures using frequency response functions and principle component analysis
JP2008032474A (ja) 巨視的探査理論に基づいたコンクリート電柱地中部の広帯域超音波探査方法
JP6371575B2 (ja) 超音波探傷検査方法及び超音波探傷検査装置
Ambrozinski et al. Damage detection in plate-like structures based on mode-conversion sensing with 3D laser vibrometer
Kudela et al. Damage visualization enhancement by the wave field filtering and processing
Giri Integrated approaches for non-destructive testing of construction materials and structures
Golato et al. Extended defect localization in sparsity-based guided wave structural health monitoring
Gao et al. Enhanced wavefield imaging method for impact damage detection in composite laminates via laser-generated Lamb waves
McKenna et al. Infrasound measurements of a railroad bridge

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20110601

Termination date: 20120930