买专利,只认龙图腾
首页 专利交易 科技果 科技人才 科技服务 商标交易 会员权益 IP管家助手 需求市场 关于龙图腾
 /  免费注册
到顶部 到底部
清空 搜索

【发明授权】一种桥梁裂缝宽度自动测量方法_西安电子科技大学;陕西康虹交通科技有限公司;西安品码电子科技有限公司_201910590721.5 

申请/专利权人:西安电子科技大学;陕西康虹交通科技有限公司;西安品码电子科技有限公司

申请日:2019-07-02

公开(公告)日:2022-09-20

公开(公告)号:CN110136196B

主分类号:G06T7/62

分类号:G06T7/62;G06T7/13;G06T7/155;G06T5/40;G01B11/02;G01B11/14

优先权:["20190130 CN 2019100888631"]

专利状态码:有效-授权

法律状态:2022.09.20#授权;2019.09.10#实质审查的生效;2019.08.16#公开

摘要:本发明公开了一种桥梁裂缝宽度自动测量方法,主要解决现有的基于图像检测混凝土桥梁裂缝宽度信息中存在误检的问题;其实现方案为:计算机读取原始桥梁图像并进行预处理;基于Sobel算子获取裂缝主干信息,并获取主干的带状区域图;对带状区域进行灰度化并进行直方图均衡化处理;基于Sobel算子获取均衡化后的带状区域图片的梯度,并提取裂缝的共生边缘和裂缝点集合;提取裂缝点集合中位于裂缝内部的点;获取裂缝宽度信息;对裂缝进行分条处理并统计每条裂缝的宽度信息存储到计算机终端。本发明能够以高准确性和实时性完成对混凝土桥梁裂缝的检测,可用于桥梁混凝土裂缝的宽度信息的获取。

主权项:1.一种桥梁裂缝宽度自动测量方法,其特征在于,包括如下:1读取原始混凝土图像,并基于高斯卷积核对其进行平滑处理;2将平滑后的混凝土图像通过Sobel算子获取其裂缝主干图Q;3根据主干图Q,获取裂缝的带状区域图N;4对裂缝带状区域图N做灰度化处理,并进行直方图均衡化,获取图N均衡化后的带状区域图片D;5基于Sobel算子获取均衡化后的裂缝带状区域图片D的梯度值G和梯度方向θ;6根据均衡化后的裂缝带状区域图像D的梯度值G和梯度方向θ,提取均衡化后的裂缝带状区域图片的共生边缘Z和裂缝点集合C;7滤除裂缝点集合C中不在裂缝内部的像素点;8利用共生边缘Z和位于裂缝内部的裂缝点集合C获取裂缝宽度信息:8a按照从上至下,从左至右的顺序遍历共生边缘Z和位于内部的裂缝点C,获取一对共生边缘点Mxm,ym,Kxk,yk和对应的裂缝点Txt,yt;8b依据共生边缘点M,K的坐标计算搜索裂缝边界的方向8c从裂缝点Txt,yt出发在8x8的领域范围内,沿着方向查找裂缝带状区域图片中梯度值G相差最大的像素点U,将该点作为裂缝的第一个边界点,并记录该点的坐标位置xu,yu;8d从裂缝点Txt,yt出发在8x8的领域范围内,沿着的反方向查找裂缝带状区域图片中梯度值G相差最大的像素点H,将该点作为裂缝的第二个边界点,并记录该点的坐标位置xh,yh;8e按照欧式距离计算裂缝宽度w: 其中xh,yh,xu,yu分别代表边界点U点和H点的坐标;9基于主干图Q和裂缝宽度信息w,获取每条裂缝的宽度。

全文数据:一种桥梁裂缝宽度自动测量方法技术领域本发明属于测试技术领域,特别设计一种桥梁裂缝自动测量方法,可用于桥梁混凝土裂缝的宽度信息的获取。技术背景衡量桥梁混凝土病害程度的重要指标包括裂缝的长度、宽度、数量等数据信息,其中裂缝宽度信息为衡量病害程度最重要的指标,现有检测桥梁混凝土裂缝宽度信息的手段包括:人工测量法、红外分析法、和图像处理分析法。其中:人工测量法是人工利用游标卡尺测量裂缝宽度,该方法存在测量精度差,效率较低、还存在一定的危险。红外分析法是利用红外线检测裂缝宽度,该方法具有检测精度高,检测速率快等优势,但是该仪器也存在成本较高,而且需专业的人员操作,维护难度较高,在使用时还存在诸多条件限制等缺点。图像处理分析法是基于图像处理技术检测图像裂缝宽度信息,该方法具有自动化测量,效率高,测量精度高的优势,是目前混凝土检测领域研究最多的技术。程彬在其发表的论文《基于数字图像处理的表面裂缝宽度测量》赤峰学院学报自然科学版,2018,3410:96-98.中提出了两种检测裂缝的算法,手动点取边界法和框选平均值法,这两种方法均是通过手动在图像上标注裂缝的边界点,并计算标记出的边界点之间的欧式距离来获取裂缝的宽度信息。因为该方法近乎是通过人力去测量裂缝宽度,所以实时性极差,且耗费大量的时间与精力。肖利芳在其发表的论文《Sobel算子改进边缘检测算法在混凝土裂缝识别中的应用》软件导刊,2017年01期,112-114提出了一种基于Sobel算子的改进边缘检测算法对图像进行处理,并基于获取到的边缘信息提取图像裂缝的宽度信息。该方法是通过改进Sobel算子提升图像裂缝宽度信息的准确性。但是通过Sobel算子在获取图像边缘信息的过程中存在时间复杂度较高,实时性较差的缺点,而且该方法没有针对实际应用中所存在的诸多干扰加以考虑。所以在实际应用中,检测结果的准确性,实时性均很难保证。综上,目前在国内外提出许多利用图像处理技术检测混凝土裂缝宽度信息的方法,大都存在复杂度较高,抗干扰能力差,实时性差及准确性低的缺点,影响工程应用。发明内容本发明的目的在于针对上述现有技术的不足,提出一种桥梁裂缝宽度自动测量方法,以提高裂缝宽度信息检测的准确性。为实现上述目的,本发明的技术方案包括如下:1读取原始混凝土图像,并基于高斯卷积核对其进行平滑处理;2将平滑后的混凝土图像通过Sobel算子获取其裂缝主干图Q;3根据主干图Q,获取裂缝的带状区域图N;4对裂缝带状区域图N做灰度化处理,并进行直方图均衡化,获取图N均衡化后的带状区域图片D;5基于Sobel算子获取均衡化后的裂缝带状区域图片D的梯度值G和梯度方向θ;6根据均衡化后的裂缝带状区域图像D的梯度值G和梯度方向θ,提取均衡化后的裂缝带状区域图片的共生边缘Z和裂缝点集合C;7滤除裂缝点集合C中不在裂缝内部的像素点;8利用共生边缘Z和位于裂缝内部的裂缝点集合C获取裂缝宽度信息:8a按照从上至下,从左至右的顺序遍历共生边缘Z和位于内部的裂缝点C,获取一对共生边缘点Mxm,ym,Kxk,yk和对应的裂缝点Txt,yt;8b依据共生边缘点M,K的坐标计算搜索裂缝边界的方向8c从裂缝点Txt,yt出发在8x8的领域范围内,沿着方向查找裂缝带状区域图片中梯度值G相差最大的像素点U,将该点作为裂缝的第一个边界点,并记录该点的坐标位置xu,yu;8d从裂缝点Txt,yt出发在8x8的领域范围内,沿着的反方向查找裂缝带状区域图片中梯度值G相差最大的像素点H,将该点作为裂缝的第二个边界点,并记录该点的坐标位置xh,yh;8e按照欧式距离计算裂缝宽度w:其中xh,yh,xu,yu分别代表边界点U点和H点的坐标。9基于主干图Q和裂缝宽度信息w,获取每条裂缝的宽度。进一步,9中基于主干图Q和裂缝宽度信息w,获取每条裂缝的宽度,其实现如下:9a按照从上至下,从左至右的顺序遍历裂缝主干图Q,获取像素点三色通道R、G、B数据均非0的像素点S,且将该像素点S处的B通道的数据赋值为255,R、G通道数据赋值为0,同时判断该点在3x3邻域范围内,是否存在除S点以外另一个非0点P,若存在,则执行9b,若不存在,则继续遍历主干图Q;9b以点P为中心,按照从左至右,从上至下的顺序在其3x3的邻域范围内进行搜索,判断是否存在三色通道R、G、B数据均为非0的像素点L,如果存在,则将像素点L处的B通道数据赋值为255,R、G通道数据赋值为0,继续遍历,如果不存在,则执行9c;9c记录遍历过的裂缝主干点的位置信息,并加以标记裂缝编号,判断是否将裂缝主干图像中所有像素点遍历完毕,若遍历完毕则终止搜索,得到分条裂缝的位置信息,执行9d,如没有遍历完毕,则裂缝编号加1,返回9a;9d遍历分条后每一条裂缝的位置信息和裂缝宽度信息w,并将裂缝宽度信息w按照每条裂缝的位置进行存储,获取每条裂缝的宽度信息。本发明与现有技术相比,具有如下优点:1本发明采用基于灰度直方图均衡化实现了对裂缝边缘的锐化,使得裂缝边缘更加明显,克服了目前检测算法在检测边缘出现的由光照导致的图像边缘虚化的缺点,提高了宽度检测的准确性。2本发明使用基于自适应阈值来判定该裂缝点是位于裂缝内部还是外部点,在检测裂缝宽度时能够更好去除检测错误的裂缝点,克服了现有检测方法中存在的误检率高的问题,提升了检测的准确性。3本发明使用共生边缘的梯度方向作为种子点出发搜索裂缝边缘的方向,克服了目前检测裂缝宽度方向存在极大不确定性、导致图像边缘宽度计算不准确,进一步提高了宽度信息的检索准确率。附图说明图1是本发明的实现流程图;图2是本发明的仿真结果图。具体实施方式下面参照附图1对本发明的实施例及效果进一步的详细描述。参照图1,本发明的实施步骤如下:步骤1,读取原始混凝土图像,并对其进行平滑处理。1a用数字摄相机采集待检测混凝土裂缝图像,并将摄像机中采集到的裂缝图像存储到计算机中;1b按照高斯滤波公式对计算机存储裂缝图像进行平滑处理,得到平滑后的裂缝图像:其中P1代表原始裂缝图像的像素值,P2代表经过高斯卷积核得到平滑后的图像像素值,*代表卷积。[·]表示高斯卷积核。步骤2,基于Sobel算子获取平滑后的裂缝图像的裂缝主干图Q。目前计算图像梯度所使用的算子包括:Roberts、Prewitt、Sobel和Lapacian算子,本实例中采用但不限于Sobel算子,其实现如下:2a获取步骤1b中得到的平滑后裂缝图像;2b获取经过高斯卷积平滑后的裂缝图像的梯度值E和梯度方向α:2b1通过纵向Sobel算子计算平滑后裂缝图像的纵向梯度值Gy:其中P1表示原始灰度图像中像素点的值;*表示卷积操作,[·]表示纵向Sobel算子矩阵;2b2通过横向Sobel算子计算平滑后裂缝图像的横向梯度值Gx:其中[·]表示横向Sobel算子矩阵;2b3通过公式和公式获取到的Gx和Gy,计算裂缝图像的梯度值E和梯度方向α:2b4将步骤2b3获取到的梯度值E和梯度方向α分别作为裂缝图像梯度图的R通道和B通道的值,得到平滑后的裂缝图像的梯度图;2c根据2b获取到的原始裂缝图的梯度值E和梯度方向α,获取裂缝点集合T:2c1本实例中设置阈值L=20,将梯度图像中梯度值G小于阈值L的像素点处的R通道赋值为L-1,大于阈值L的像素点值不变,获取新的梯度图;2c2按照从左至右,从上至下的顺序遍历新的梯度图,判断是否遍历完新的梯度图:如果新的梯度图遍历结束,执行2d;如果未遍历完新的梯度图,再判断新梯度图上的像素点P处R通道值是否为L-1:若是,则遍历下一个像素点,并判断R通道的值,若不是,则先获取该像素点P处梯度值E1和梯度方向α1,再执行2c3;2c3在新的梯度图上,沿着像素点P的梯度方向α1遍历8个像素点,并判断新梯度图上的像素点S处的R通道的值是否为L-1:若是,则执行2c4,若不是,则先获取像素点S处的梯度值E2,再执行2c5;2c4在新的梯度图上,沿着像素点P的梯度方向α1的反方向遍历8个像素点,判断新梯度图上的像素点S处R通道的值是否为L-1:若是,则返回2c2,若不是,则先获取该像素点S处的梯度值E2,再执行2c5;2c5本实例设置梯度差阈值为g=26,若|E1-E2|<g,则像素点S和像素点P处的梯度值保持不变,并取像素点S和像素点P两点之间的中点位置W,并将像素点W存储到裂缝点集合T中;否则,将像素点S和像素点P处的三色通道的值赋值为L-1,并准备运算下一个像素点,返回2c2;2d连接裂缝点集合T中的像素点:2d1设定最小生成树集合为E;2d2从裂缝点集合T中随机取一点u作为路径的开始节点,并从点集合T中选取与点u欧氏距离最小的像素点z;之后,将像素点u和像素点z的连线放入最小生成树集合E中;2d3判断裂缝点集合T是否全部遍历完成,如果未遍历完成,则从像素点z出发,并重复步骤2d2;如果全部遍历结束,则最小生成树集合E获取完成,并执行2e;2e获取裂缝主干图Q:本实例设置边长阈值r=100,先将2d获取到的完整最小生成树集合E中欧式距离大于边长阈值r的连线删除,遍历最小生成树集合E中的剩余的连线,并获取连线上的像素点,这些像素点组成裂缝主干图Q。步骤3,获取裂缝的带状区域图N。3a按照从左至右,从上到下的顺序遍历裂缝主干图片Q,获取其中像素值不为0的像素点集合A;3b遍历像素点集合A中的像素点U,从原始裂缝图中获取与U点坐标位置相同的像素点B及其周围20x20领域内的像素点,将像素点B及其周围20x20领域内的像素点作为带状区域内部的像素点,继续遍历像素点集合A,直至集合A遍历结束,执行3c。3c将3b获取到的所有带状区域内部的像素点组成裂缝带状区域图N。步骤4,获取裂缝带状区域图N均衡化后的带状区域图片D。4a对步骤3获取到的裂缝带状区域图N进行灰度化处理,其实现如下:按照从左至右,从上至下的顺序遍历裂缝带状区域图N,并计算遍历到的像素点U的灰度值gray:gray=0.30×R+0.59×G+0.11×B其中,R表示像素点U的红色通道值,G表示像素点U的绿色通道值,B表示像素点U的蓝色通道值,将计算出的灰度值gray赋值到像素点U的三色通道中,实现对于裂缝区域图N的灰度化;4b对灰度化后的裂缝带状区域图N进行直方图均衡化:4b1按照从左至右,从上至下的顺序遍历裂缝带状区域图N,并获取带状区域裂缝图片中各灰度级的像素数目ni,1≤i<256,其中i代表灰度级;4b利用灰度概率计算公式计算每一个灰度级出现的概率,其公式为:其中i为灰度图像N中的灰度级数,ni为裂缝带状区域中灰度级为i的像素数目,A为裂缝带状区域图像所包含的像素数目,Pni是灰度级数i对应的灰度概率;4c利用概率的累计分布函数计算裂缝带状区域图像的归一化直方图:其中i代表映射前的灰度级数,k代表映射后的灰度级数,Pni为第i灰度级的灰度概率,Fk代表第k个灰度级对应的累计分布概率;4d利用灰度直方图均衡化公式获取不同灰度级对应映射后的灰度级数:r=255*Fbb=1,......,255,其中b代表像素点映射前的灰度级数,r代表像素点映射后的灰度级数,Fb代表第b灰度级对应的累计分布概率;4e将裂缝带状区域图片中的像素点按照从左至右,从上至下的顺序遍历,获取其中像素值不为0的像素点,再通过4d获取这些不为0的像素点经过直方图均衡化后的新的像素值,把新的像素值赋值到这些不为0的像素点,并将这些新像素值的像素点组成均衡化后的裂缝带状区域图D。步骤5,基于Sobel算子获取均衡化后的裂缝带状区域图片D的梯度值G和梯度方向θ。目前计算图像梯度所使用的算子包括:Roberts、Prewitt、Sobel和Lapacian算子,本实例中采用但不限于Sobel算子,其实现如下:5a获取4e得到的均衡化后的裂缝带状区域图片D;5b通过纵向Sobel算子计算图像D的纵向梯度值fy:其中P1表示原始灰度图像中像素点的值;*表示卷积操作,[·]表示纵向Sobel算子矩阵;5c通过横向Sobel算子计算图像D的横向梯度值fx:其中[·]表示横向Sobel算子矩阵;5d通过公式和公式获取到的fx和fy,计算裂缝图像的梯度值G和梯度方向θ:5e将步骤2d获取到的梯度值G和梯度方向θ分别作为裂缝图像梯度图的R通道和B通道的值,得到均衡化后的裂缝带状区域图片D的梯度图。步骤6,提取均衡化后的裂缝带状区域图D的共生边缘Z和提取裂缝点集合C。6a本实例中设置阈值L=20,将均衡化后裂缝带状区域的图像D中梯度值G小于阈值L的像素点处的R通道赋值为L-1,大于阈值L的像素点值不变,获取新的梯度图;6b按照从左至右,从上至下的顺序遍历新的梯度图,获取梯度值不为L-1的像素点,并将获取到的这些像素点组成像素点集合A;6c本实例设置梯度差阈值g=26,在新的梯度图上,遍历像素点集合A中像素点U,并获取点U处的梯度值G1和梯度方向θ1,从点U出发并沿着U的梯度方向θ1和θ1的反方向各遍历8个像素点,获取其中梯度值不为L-1的像素点J,且像素点J处的梯度值G2满足|G1-G2|<g,则将像素点U和像素点J记做一对共生边缘点,并将像素点U和像素点J的中点位置处的像素点V记做裂缝点,继续遍历像素点集合A,直至集合A中的像素点遍历结束;6d将6c获取到的所有共生边缘点组成裂缝带状区域图片的共生边缘集合Z,并将用这些有裂缝点组成裂缝带状区域图片的裂缝点集合C。步骤7,滤除步骤6获取的裂缝点集合C中不在裂缝内部的像素点。7a计算自适应阈值Y:7a1统计经过直方图均衡化后的带状区域图片D中不同灰度级的像素数目,其实现如下:suman=sumbn+1,其中sumbn代表遍历当前像素点之前的第n灰度级的像素点个数的总和,suman代表当前第n灰度级的像素点个数的总和;7a2通过解方程获取自适应阈值Y,其中sumn代表第n灰度级对应的像素点的数目;7b遍历裂缝点集合C中的像素点A,如果集合C中像素点A的像素值小于阈值Y,则保留该像素点A,否则将该像素点A从集合C中删除,直至裂缝点集合C遍历完成。步骤8,获取裂缝宽度信息。8a按照从上至下,从左至右的顺序遍历步骤6得到的共生边缘Z和步骤7得到的位于内部的裂缝点C,获取一对共生边缘点Mxm,ym,Kxk,yk和对应的裂缝点Txt,yt;8b利用角度公式计算搜索裂缝边界方向其中xm,ym,xk,yk代表一对共生边缘点M,K的坐标,arctan代表反正切;8c从裂缝点Txt,yt出发在8x8的领域范围内,沿着方向查找裂缝带状区域图片的梯度图G中梯度值相差最大的像素点U,将该点作为裂缝的一个边界点,并记录该点的坐标位置xu,yu,执行8d;8d从裂缝点Txt,yt出发在8x8的领域范围内,沿着的反方向查找裂缝带状区域图片的梯度图G中梯度值相差最大的像素点H,将该点作为裂缝的另一个边界点,并记录该点的坐标位置xh,yh,执行8e;8e按照欧式距离计算裂缝宽度w:其中xh,yh,xu,yu分别代表边界点U点和H点的坐标。步骤9,分条统计裂缝宽度信息。9a按照从上至下,从左至右的顺序遍历步骤2获取到的裂缝主干图Q,获取像素点三色通道R、G、B数据均非0的像素点S,且将该像素点S处的B通道的数据赋值为255,R、G通道数据赋值为0,同时判断该点在3x3邻域范围内,是否存在除S点以外另一个非0点P,若存在,则执行9b,若不存在,则继续遍历主干图Q;9b以点P为中心,按照从左至右,从上至下的顺序在其3x3的邻域范围内进行搜索,判断是否存在三色通道R、G、B数据均为非0的像素点L,如果存在,则将像素点L处的B通道数据赋值为255,R、G通道数据赋值为0,继续遍历,如果不存在,则执行9c;9c记录遍历过的裂缝主干点的位置信息,并加以标记裂缝编号,判断是否将裂缝主干图像中所有像素点遍历完毕,若遍历完毕则终止搜索,得到分条裂缝的位置信息,执行9d,如没有遍历完毕,则裂缝编号加1,返回9a;9d遍历分条后每一条裂缝的位置信息和裂缝宽度信息w,并将裂缝宽度信息w按照每条裂缝的位置进行存储,获取每条裂缝的宽度信息。下面结合附图1对本发明的仿真效果做进一步描述。1.仿真条件:本发明的仿真实验是在CPU主频2.7GHz、内存7.85GB的硬件环境和VisualStudio2013的软件环境下进行的。首先,在混凝土桥梁表面选取四个典型裂缝场景,并通过游标卡尺测量裂缝的长度、宽度并获取裂缝数量等参数;然后,用图像采集设备CCD相机获取经过游标卡尺测量过的四组混凝土裂缝图像,图像的分辨率为5760*3840,每次采集面积为450mm*300mm。2.仿真实验内容:实验1,使用本发明对采集的第一组混凝土裂缝图像进行检测,结果如图2,其中:图2a为混凝土桥梁裂缝图像,该图像采集于陕西省西安市灞河上的一座桥梁,图2b、2c、2d、2e为本发明对图2a标注的四个区域进行混凝土桥梁裂缝宽度检测得到的结果图。从图2可以看出,本发明在复杂的背景干扰中,仍然能够准确地获取桥梁混凝土裂缝图像中裂缝的宽度信息,说明本发明具有较好的抗干扰能力。实验2,使用本发明对采集的第二组混凝土裂缝图像进行检测。实验3,使用本发明对采集的第三组混凝土裂缝图像进行检测。实验4,使用本发明对采集的第四组混凝土裂缝图像进行检测。对上述实验1,实验2,实验3和实验4的四组实测数据进行统计,结果如表1所示。表1混凝土桥梁检测一览表表1中,“测量相对误差”中的“+”代表测量结果偏大,“-”代表测量结果偏小。通过表1的裂缝平均宽度检测相对误差可以看出,本发明能够准确获取裂缝的宽度信息。

权利要求:1.一种桥梁裂缝宽度自动测量方法,其特征在于,包括如下:1读取原始混凝土图像,并基于高斯卷积核对其进行平滑处理;2将平滑后的混凝土图像通过Sobel算子获取其裂缝主干图Q;3根据主干图Q,获取裂缝的带状区域图N;4对裂缝带状区域图N做灰度化处理,并进行直方图均衡化,获取图N均衡化后的带状区域图片D;5基于Sobel算子获取均衡化后的裂缝带状区域图片D的梯度值G和梯度方向θ;6根据均衡化后的裂缝带状区域图像D的梯度值G和梯度方向θ,提取均衡化后的裂缝带状区域图片的共生边缘Z和裂缝点集合C;7滤除裂缝点集合C中不在裂缝内部的像素点;8利用共生边缘Z和位于裂缝内部的裂缝点集合C获取裂缝宽度信息:8a按照从上至下,从左至右的顺序遍历共生边缘Z和位于内部的裂缝点C,获取一对共生边缘点Mxm,ym,Kxk,yk和对应的裂缝点Txt,yt;8b依据共生边缘点M,K的坐标计算搜索裂缝边界的方向8c从裂缝点Txt,yt出发在8x8的领域范围内,沿着方向查找裂缝带状区域图片中梯度值G相差最大的像素点U,将该点作为裂缝的第一个边界点,并记录该点的坐标位置xu,yu;8d从裂缝点Txt,yt出发在8x8的领域范围内,沿着的反方向查找裂缝带状区域图片中梯度值G相差最大的像素点H,将该点作为裂缝的第二个边界点,并记录该点的坐标位置xh,yh;8e按照欧式距离计算裂缝宽度w:其中xh,yh,xu,yu分别代表边界点U点和H点的坐标。9基于主干图Q和裂缝宽度信息w,获取每条裂缝的宽度。2.根据权利要求1所述的方法,其特征在于,9中基于主干图Q和裂缝宽度信息w,获取每条裂缝的宽度,其实现如下:9a按照从上至下,从左至右的顺序遍历裂缝主干图Q,获取像素点三色通道R、G、B数据均非0的像素点S,且将该像素点S处的B通道的数据赋值为255,R、G通道数据赋值为0,同时判断该点在3x3邻域范围内,是否存在除S点以外另一个非0点P,若存在,则执行9b,若不存在,则继续遍历主干图Q;9b以点P为中心,按照从左至右,从上至下的顺序在其3x3的邻域范围内进行搜索,判断是否存在三色通道R、G、B数据均为非0的像素点L,如果存在,则将像素点L处的B通道数据赋值为255,R、G通道数据赋值为0,继续遍历,如果不存在,则执行9c;9c记录遍历过的裂缝主干点的位置信息,并加以标记裂缝编号,判断是否将裂缝主干图像中所有像素点遍历完毕,若遍历完毕则终止搜索,得到分条裂缝的位置信息,执行9d,如没有遍历完毕,则裂缝编号加1,返回9a;9d遍历分条后每一条裂缝的位置信息和裂缝宽度信息w,并将裂缝宽度信息w按照每条裂缝的位置进行存储,获取每条裂缝的宽度信息。3.根据权利要求1所述的方法,其特征在于,1中基于高斯卷积核对计算机读入的原始混凝土图像进行平滑处理,通过如下公式进行:其中P1代表原始裂缝图像的像素值,P2代表经过高斯卷积核得到平滑后的图像像素值,*代表卷积,[·]表示高斯卷积核。4.根据权利要求1所述的方法,其特征在于,2基于Sobel算子获取裂缝主干图Q,其实现如下:2a获取平滑后原始裂缝图的梯度值E和梯度方向α:其中:表示为是灰度图像像素点在横向的梯度大小;表示为是灰度图像中像素点的纵向的梯度大小;P1表示原始灰度图像中像素点的值;*表示卷积操作,[·]表示Sobel算子矩阵;2b根据原始裂缝图的梯度值E和梯度方向α,获取裂缝点集合T:2b1设置阈值L,将梯度图像中梯度值E小于阈值L的像素点处的R通道赋值为L-1,大于阈值L的像素点值不变,获取新的梯度图;2b2按照从左至右,从上至下的顺序遍历新的梯度图,判断是否遍历完新的梯度图;如果新的梯度图遍历结束,执行2c;如果未遍历完新的梯度图,再判断新梯度图上的像素点处R通道的值是否为L-1:若是,则继续遍历下一个像素点,并判断像点R通道的值,若不是,则先获取该像素点处梯度值E1和梯度方向α1,再执行2b3;2b3在新的梯度图上,沿着方向α1遍历8个像素点,判断新梯度图上的像素点处R通道的值是否为L-1:若是,则执行2b4,若不是,则先获取该像素点处梯度值E2,再执行2b5;2b4在新的梯度图上,沿着α1反方向遍历8个像素点,判断新梯度图上的像素点处R通道的值是否为L-1:若是,则返回2b2,若不是,则先获取该像素点处梯度值E2,再执行2b5;2b5设置梯度差阈值为常数g,若|E1-E2<g则保留E1和E2所在的像素点的梯度值,并取E1和E2两点之间的中点位置W,并将像素点W存储到裂缝点集合T中;否则,将E1和E2所在的像素点处的三色通道的值赋值为L-1,并准备运算下一个像素点,返回2b2;2c连接裂缝点:2c1设定最小生成树集合为E;2c2从裂缝点集合T中随机取一点u作为路径的开始节点,并从点集合T中选取与点u欧氏距离最小的裂缝点z,之后将裂缝点u和裂缝点z的连线放入最小生成树集合E中;2c3从裂缝点z出发,重复步骤2c2,直到裂缝点集合T中的裂缝点全部遍历结束;2d提取裂缝主干:将最小生成树集合E中欧式距离大于阈值r=100的边删除,剩余的边组成裂缝主干图Q。5.根据权利要求1所述的方法,其特征在于,3中获取裂缝带状区域图N,是按照从左至右,从上到下的顺序遍历裂缝主干图片Q,获取像素值不为0的像素点A,并将原始裂缝图中与A点坐标位置相同的像素点B及其周围20x20领域内的像素点的像素值,按照对应位置赋值到裂缝带状区域图N中。6.根据权利要求1所述的方法,其特征在于,4中对裂缝带状区域图N进行直方图均衡化,其实现下:4a统计带状区域裂缝图片中各灰度级的像素数目ni,0≤i<256;4b利用灰度概率计算公式计算每一个灰度级出现的概率,其公式为:其中i为灰度图像N中的灰度级数,ni为裂缝带状区域中灰度级为i的像素数目,A为裂缝带状区域图像所包含的像素数目,Pni是灰度级数i对应的灰度概率;4c利用概率的累计分布函数计算裂缝带状区域图像的归一化直方图:其中i代表映射前的灰度级数,k代表映射后的灰度级数,Pni为第i灰度级的灰度概率,Fk代表第k个灰度级对应的累计分布概率;4d利用灰度直方图均衡化公式获取不同灰度级对应映射后的灰度级数:r=255*Fbb=1,......,255其中b代表像素点映射前的灰度级数,r代表像素点映射后的灰度级数,Fb代表第b灰度级对应的累计分布概率。4e将裂缝带状区域图片中的像素点按照从左至右,从上至下的顺序遍历,获取其中像素值不为0的像素点,再通过4d获取这些不为0的像素点经过直方图均衡化后的新的像素值,把新的像素值赋值到这些不为0的像素点,用这些具有新像素值的像素点组成均衡化后的裂缝带状区域图片D。7.根据权利要求1所述的方法,其特征在于,5中基于Sobel算子获取均衡化后的裂缝带状区域图片D的梯度值G和梯度方向θ,其计算公式如下:其中:表示为是裂缝带状区域图片D像素点在横向的梯度大小;表示为是裂缝带状区域图片D中像素点的纵向的梯度大小;P1表示裂缝带状区域图片D中像素点的值;*表示卷积操作,[·]表示Sobel算子矩阵。8.根据权利要求1所述的方法,其特征在于,6中提取均衡化后的裂缝带状区域图片的共生边缘Z和裂缝点集合C,其实现如下:6a设置阈值L,将均衡化后裂缝带状区域的图像D中梯度值G小于阈值L的像素点处的R通道赋值为L-1,大于阈值L的像素点值不变,获取新的梯度图;6b按照从左至右,从上至下的顺序遍历新的梯度图,获取梯度值不为L-1的像素点,将这些像素点存储到像素点集合A中;6c在新的梯度图上,设置梯度差阈值为常数g,遍历像素点集合A中像素点U,并获取点U处的梯度值G1和梯度方向θ1,从点U出发并沿着U的梯度方向θ1和θ1的反方向各遍历8个像素点,获取其中梯度值不为L-1的像素点J,且像素点J处的梯度值G2满足|G1-G2|<g,则将像素点U和像素点J记做一对共生边缘点,并将像素点U和像素点J的中点位置处的像素点V记做裂缝点,继续遍历像素点集合A,直至集合A中的像素点全部遍历完成;6d将6c获取到的所有共生边缘点组成裂缝带状区域图片的共生边缘集合Z,并将6c获取到的所有裂缝点组成裂缝带状区域图片的裂缝点集合C。9.根据权利要求1所述的方法,其特征在于,7中滤除裂缝点集合C中不在裂缝内部的像素点,其实现方式如下:7a获取自适应阈值Y:7a1统计经过直方图均衡化后的带状区域图片D中不同灰度级的像素数目,其公式如下:suman=sumbn+1,其中sumbn代表遍历当前像素点之前的第n灰度级的像素点个数的总和,suman代表当前第n灰度级的像素点个数的总和;7a2通过解方程获取自适应阈值Y,其中sumn代表第n灰度级对应的像素点的数目;7b遍历裂缝点集合C中的像素点A,如果集合C中像素点A的像素值小于阈值Y,则保留该像素点A;否则将该像素点A从集合C中移除。10.根据权利要求1所述的方法,其特征在于,8b依据共生边缘点计算搜索裂缝边界的方向其计算公式如下:其中xm,ym,xk,yk分别共生边缘点M,K的坐标,arctan代表反正切,代表搜索裂缝边界方向。

百度查询: 西安电子科技大学;陕西康虹交通科技有限公司;西安品码电子科技有限公司 一种桥梁裂缝宽度自动测量方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息,力求客观、公正,但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解,仅供参考使用,不能作为本公司承担任何法律责任的依据或者凭证。