基于图像距离变换的混凝土裂缝检测方法技术

一种基于图像距离变换的混凝土裂缝检测方法，其步骤为：(1)计算机读取原始混凝土裂缝图像并预处理得到平滑图像(2)计算平滑图像的梯度值(3)标记平滑图像极大值像素点(4)标记平滑图像中裂缝种子点(5)连接裂缝种子点(6)裁剪部分种子点的连线(7)获取裂缝参数信息。本发明专利技术可以应用于混凝土裂缝检测领域，通过图像处理技术实现高速连续地从混凝土裂缝图像中获取裂缝的数量、长度、宽度信息。由于本发明专利技术采用了像素梯度值对称性分析法，并充分利用图像距离变换减少裂缝连接的计算量，因此具有抗背景干扰能力强，检测精度高、检测速度快的特点。

【技术实现步骤摘要】
基于图像距离变换的混凝土裂缝检测方法
本专利技术属于物理
，更进一步涉及图像处理
中的一种基于图像距离变换的混凝土裂缝检测方法。本专利技术可应用于混凝土裂缝检测领域，可以通过图像处理技术实现快速地从混凝土裂缝图像中获取裂缝的长度、宽度、数量等数据，通过这些数据对混凝土的病害程度做出判定。
技术介绍
混凝土裂缝的长度、宽度和数量等信息是衡量混凝土病害程度的重要指标。常见的混凝土裂缝信息获取方法中：游标卡尺测量法是手工操作，效率低、精度较差，在混凝土桥梁场景中还可能对工人造成一定的危险；红外分析法测量速度较快，精度较高，但是相应的设备昂贵，使用条件高，操作复杂，维护困难，并且需要专业的技术人员操作，人力成本很高；图像处理分析法是自动化测量，效率高，测量精度高，是目前混凝土裂缝检测领域研究最多的技术，该方法通过CCD摄像机和图像采集卡完成图像数据采集，计算机对采集的混凝土裂缝图像进行处理和计算得到裂缝的相关信息，从而高效率、低成本地对评定混凝土目标的灾害等级。目前国内外提出了多种利用图像处理技术检测混凝土裂缝的方法，逐渐形成了完整、实用、科学的检测体系、方法及标准，但是目前这些检测方法还存在误检率高、算法复杂度高、计算速度慢等缺点，对实际应用造成了困难。南京理工大学在其申请的专利文献“基于多种低层显著特征的路面裂缝检测方法”(专利申请号：201410141476.7，申请公开号：103903268A)中提出了一种基于多种低层显著特征的路面裂缝检测方法。该方法首先采集路面灰度图像，然后将图像划分为局部分块，计算每个像素的稀少性概率，生成裂缝特征图，通过贝叶斯理论的双邻域扩张提取裂缝，再通过自适应阈值分割进行区域增长增强进一步提取裂缝，最后计算裂缝的面积、位置等各项具体参数。该方法虽然能够检测出裂缝，但是该方法的不足之处在于，通过计算像素的稀少性概率生成裂缝特征图无法排除与裂缝像素类似的背景像素的干扰，通过阈值分割多次提取裂缝导致计算量过大，耗费时间较长。整个算法的检测基于多种低层显著性特征，因此检测的精度受裂缝图像成像质量的制约。许薛军在其发表的论文“基于数字图形的混凝土裂缝检测技术”(湖南大学学报，JHU(NS)，2013年07期)中提出了一种基于改进的平滑滤波和边缘检测的方法。该方法首先对混凝土裂缝进行灰度化处理，然后再进行平滑滤波，再提取图像边缘，获取裂缝骨架，最后计算裂缝的面积、数量的参数。该方法虽然步骤简单，计算复杂度较低，但是该方法仍然存在的不足之处是，平滑滤波对于消除噪声的作用非常有限，在后续处理过程中，图像里包含了大量的背景噪声，对于边缘检测造成了很大的干扰，导致检测的准确率较低，进一步导致裂缝的长度、面积等信息误差很大，对实际应用造成很大的困难。
技术实现思路
本专利技术的目的在于克服上述已有技术的不足，提出一种基于图像距离变换的混凝土裂缝检测方法。本专利技术既可以很好地保留混凝土裂缝的位置信息，又可以较好的去除背景噪声的干扰，具有较高的检测精度，同时由于本算法的复杂度较低，因而可以大幅减少混凝土裂缝检测在图像处理上耗费的时间。实现本专利技术的基本思路是：计算机读取混凝土裂缝图像，采用高斯滤波和梯度计算相结合的方法对图像进行处理得到二值化图像，然后对图像中的像素点计算局部极大值获取种子点、再通过距离变换方法连接种子点、提取裂缝骨架，最后获取裂缝的长度、宽度、数量信息，并输出到计算机终端。为实现上述目的，本专利技术的方法包括如下步骤：(1)计算机对读取的图像进行预处理(1a)计算机读取图像采集卡从数字摄相机采集的待检测的混凝土裂缝图像，(1b)使用灰度转化公式，将待检测的混凝土裂缝图像作灰度化处理，得到灰度图像；(1c)使用高斯内核卷积处理公式，将灰度图像中的像素点，按从左至右、从上到下的顺序依次与高斯内核卷积，得到平滑图像；(2)计算平滑图像的梯度值(2a)计算平滑图像中每一个像素点与其右边相邻及下边相邻的像素点的横向差值和纵向差值；(2b)使用梯度公式，根据横向差值和纵向差值，计算每一个像素点的梯度值，将所有梯度值组成梯度图像；(2c)使用梯度角度公式，根据横向差值和纵向差值，计算每一个像素点的梯度方向角度值，将所有梯度方向角度值组成梯度方向图像；(3)标记平滑图像局部极大值像素点(3a)按照四象限划分法，将平面直角坐标系划分为四个象限，给梯度方向图像中每个像素点赋予一个象限序号，给梯度图像中像素点赋予与梯度方向图像中像素点相同的象限序号；(3b)按照最大值标记法，对梯度图像中的每一个像素点进行标记；(4)标记平滑图像中裂缝种子点(4a)按照八象限划分法，将平面直角坐标系划分为八个象限，给梯度方向图像中每个像素点赋予一个象限序号，给梯度图像中像素点赋予与梯度方向图像中像素点相同的象限序号；(4b)按照对称像素点标记法，对梯度图像中的每一个像素点进行标记；(5)连接裂缝种子点(5a)将梯度图像中的值不为0的像素点作为种子点，按照从左向右，从上到下的顺序依次对每个种子点进行编号；(5b)建立一个点集合、一个长度集合、一个边集合；(5c)将1号种子点放入点集合中；(5d)将点集合中的所有种子点作为计算变换距离的初始种子点，使用距离变换公式，计算初始种子点与其所有相邻种子点的变换距离；(5e)将与初始种子点变换距离最小的相邻种子点放入点集合中；(5f)将最小变换距离放入长度集合中；(5g)将初始种子点与变换距离最小的种子点连成的直线放入边集合中；(5h)判断是否所有的种子点都放入到点集合中，若是，则执行步骤(6)，否则,执行步骤(5d)；(6)裁剪部分种子点的连线(6a)从长度集合中选取所有长度大于30个像素点的长度元素，将所选取的长度元素的值置为平滑图像的宽度值；(6b)从边集合中去除所有长度值与所选取的长度元素值相等的边元素；(7)获取裂缝参数信息(7a)从边集合中选取所有长度值与平滑图像的宽度值不相等的边元素，将梯度图像中所选取的边元素的两个端点的像素点值置为1后形成多个连通域；(7b)找出每个连通域的任意两个端点之间的最长路径，将连通域中不在最长路径上的像素点的值置为0；(7c)将梯度图像中连通域的数量作为混凝土裂缝图像中裂缝的数量，连通域的位置作为混凝土裂缝图像中裂缝的位置；(7d)将梯度图像中连通域包含的像素点总数作为混凝土裂缝图像中裂缝的长度。本专利技术与现有技术相比较，具有如下优点：第一，本专利技术在混凝土裂缝图像预处理部分采用了高斯滤波法和灰度转化法相结合的方法，克服了现有技术中图像预处理对于混凝土裂缝图像细微边缘显示不清的缺点，同时减少了待处理图像数据量，使得本专利技术提高了图像处理速度。第二，本专利技术通过对后处理图像采用梯度计算法和种子点标记法相结合的方法，有效地去除了背景噪声，克服了现有技术噪声干扰大，裂缝位置误差大、误检率高的缺点，使得本专利技术具有检测准确率高的优点。第三，本专利技术采用了距离变换算法，能够以较低的复杂度获取完整的裂缝，克服了现有技术连接裂缝复杂度高、检测结果中裂缝不连续的缺点，使得本专利技术具有处理速度快、误检率低的优点。附图说明图1为本专利技术的流程图。图2为本专利技术的仿真效果图。具体实施方式下面参照附图1对本专利技术的步骤作进一步的详细描述。步骤1，计算机读取原始图像并进行预处理。第一本文档来自技高网...

【技术保护点】
一种基于图像距离变换的混凝土裂缝检测方法，其特征在于，包括如下步骤：(1)计算机对读取的图像进行预处理：(1a)计算机读取图像采集卡从数字摄相机采集的待检测的混凝土裂缝图像，(1b)使用灰度转化公式，将待检测的混凝土裂缝图像作灰度化处理，得到灰度图像；(1c)使用高斯内核卷积处理公式，将灰度图像中的像素点，按从左至右、从上到下的顺序依次与高斯内核卷积，得到平滑图像；(2)计算平滑图像的梯度值：(2a)计算平滑图像中每一个像素点与其右边相邻及下边相邻的像素点的横向差值和纵向差值；(2b)使用梯度公式，根据横向差值和纵向差值，计算每一个像素点的梯度值，将所有梯度值组成梯度图像；(2c)使用梯度角度公式，根据横向差值和纵向差值，计算每一个像素点的梯度方向角度值，将所有梯度方向角度值组成梯度方向图像；(3)标记平滑图像局部极大值像素点：(3a)按照四象限划分法，将平面直角坐标系划分为四个象限，给梯度方向图像中每个像素点赋予一个象限序号，给梯度图像中像素点赋予与梯度方向图像中像素点相同的象限序号；(3b)按照最大值标记法，对梯度图像中的每一个像素点进行标记；(4)标记平滑图像中裂缝种子点：(4a)按照八象限划分法，将平面直角坐标系划分为八个象限，给梯度方向图像中每个像素点赋予一个象限序号，给梯度图像中像素点赋予与梯度方向图像中像素点相同的象限序号；(4b)按照对称像素点标记法，对梯度图像中的每一个像素点进行标记；(5)连接裂缝种子点：(5a)将梯度图像中的值不为0的像素点作为种子点，按照从左向右，从上到下的顺序依次对每个种子点进行编号；(5b)建立一个点集合、一个长度集合、一个边集合；(5c)将1号种子点放入点集合中；(5d)将点集合中的所有种子点作为计算变换距离的初始种子点，使用距离变换公式，计算初始种子点与其所有相邻种子点的变换距离；(5e)将与初始种子点变换距离最小的相邻种子点放入点集合中；(5f)将最小变换距离放入长度集合中；(5g)将初始种子点与变换距离最小的种子点连成的直线放入边集合中；(5h)判断是否所有的种子点都放入到点集合中，若是，则执行步骤(6)，否则,执行步骤(5d)；(6)裁剪部分种子点的连线：(6a)从长度集合中选取所有长度大于30个像素点的长度元素，将所选取的长度元素的值置为平滑图像的宽度值；(6b)从边集合中去除所有长度值与所选取的长度元素值相等的边元素；(7)获取裂缝参数信息：(7a)从边集合中选取所有长度值与平滑图像的宽度值不相等的边元素，将梯度图像中所选取的边元素的两个端点的像素点值置为1后形成多个连通域；(7b)找出每个连通域的任意两个端点之间的最长路径，将连通域中不在最长路径上的像素点的值置为0；(7c)将梯度图像中连通域的数量作为混凝土裂缝图像中裂缝的数量，连通域的位置作为混凝土裂缝图像中裂缝的位置；(7d)将梯度图像中连通域包含的像素点总数作为混凝土裂缝图像中裂缝的长度。...

【技术特征摘要】
1.一种基于图像距离变换的混凝土裂缝检测方法，其特征在于，包括如下步骤：(1)计算机对读取的图像进行预处理：(1a)计算机读取图像采集卡从数字摄相机采集的待检测的混凝土裂缝图像，(1b)使用灰度转化公式，将待检测的混凝土裂缝图像作灰度化处理，得到灰度图像；(1c)使用高斯内核卷积处理公式，将灰度图像中的像素点，按从左至右、从上到下的顺序依次与高斯内核卷积，得到平滑图像；(2)计算平滑图像的梯度值：(2a)计算平滑图像中每一个像素点与其右边相邻及下边相邻的像素点的横向差值和纵向差值；(2b)使用梯度公式，根据横向差值和纵向差值，计算每一个像素点的梯度值，将所有梯度值组成梯度图像；(2c)使用梯度角度公式，根据横向差值和纵向差值，计算每一个像素点的梯度方向角度值，将所有梯度方向角度值组成梯度方向图像；(3)标记平滑图像局部极大值像素点：(3a)按照四象限划分法，将平面直角坐标系划分为四个象限，给梯度方向图像中每个像素点赋予一个象限序号，给梯度图像中像素点赋予与梯度方向图像中像素点相同的象限序号；(3b)按照最大值标记法，对梯度图像中的每一个像素点进行标记；(4)标记平滑图像中裂缝种子点：(4a)按照八象限划分法，将平面直角坐标系划分为八个象限，给梯度方向图像中每个像素点赋予一个象限序号，给梯度图像中像素点赋予与梯度方向图像中像素点相同的象限序号；(4b)按照对称像素点标记法，对梯度图像中的每一个像素点进行标记；(5)连接裂缝种子点：(5a)将梯度图像中的值不为0的像素点作为种子点，按照从左向右，从上到下的顺序依次对每个种子点进行编号；(5b)建立一个点集合、一个长度集合、一个边集合；(5c)将1号种子点放入点集合中；(5d)将点集合中的所有种子点作为计算变换距离的初始种子点，使用距离变换公式，计算初始种子点与其所有相邻种子点的变换距离；(5e)将与初始种子点变换距离最小的相邻种子点放入点集合中；(5f)将最小变换距离放入长度集合中；(5g)将初始种子点与变换距离最小的种子点连成的直线放入边集合中；(5h)判断是否所有的种子点都放入到点集合中，若是，则执行步骤(6)，否则,执行步骤(5d)；(6)裁剪部分种子点的连线：(6a)从长度集合中选取所有长度大于30个像素点的长度元素，将所选取的长度元素的值置为平滑图像的宽度值；(6b)从边集合中去除所有长度值与所选取的长度元素值相等的边元素；(7)获取裂缝参数信息：(7a)从边集合中选取所有长度值与平滑图像的宽度值不相等的边元素，将梯度图像中所选取的边元素的两个端点的像素点值置为1后形成多个连通域；(7b)找出每个连通域的任意两个端点之间的最长路径，将连通域中不在最长路径上的像素点的值置为0；(7c)将梯度图像中连通域的数量作为混凝土裂缝图像中裂缝的数量，连通域的位置作为混凝土裂缝图像中裂缝的位置；(7d)将梯度图像中连通域包含的像素点总数作为混凝土裂缝图像中裂缝的长度。2.根据权利要求1中所述的基于图像距离变换的混凝土裂缝检测方法，其特征在于，步骤(1b)中所述的灰度转化公式如下：gray(x,y)＝0.30R(x,y)+0.59G(x,y)+0.11B(x,y)其中,gray(x,y)表示坐标为(x,y)的像素点的灰度值，R(x,y)表示坐标为(x,y)的像素点的红色分量值，G(x,y)表示坐标为(x,y)的像素点的绿色分量值，B(x,y)表示坐标为(x,y)的像素点的蓝色分量值。3.根据权利要求1中所述的基于图像距离变换的混凝土裂缝检测方...

【专利技术属性】
技术研发人员：李云松，郭展扩，杜建超，胡银林，
申请(专利权)人：西安电子科技大学，
类型：发明
国别省市：陕西,61