CN114998096B - 一种路面图像中裂缝比例检测方法及系统 - Google Patents
一种路面图像中裂缝比例检测方法及系统 Download PDFInfo
- Publication number
- CN114998096B CN114998096B CN202210796996.6A CN202210796996A CN114998096B CN 114998096 B CN114998096 B CN 114998096B CN 202210796996 A CN202210796996 A CN 202210796996A CN 114998096 B CN114998096 B CN 114998096B
- Authority
- CN
- China
- Prior art keywords
- image
- crack
- pavement
- pavement crack
- center
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000001914 filtration Methods 0.000 claims abstract description 15
- 230000004927 fusion Effects 0.000 claims abstract description 14
- 238000012545 processing Methods 0.000 claims abstract description 14
- 238000004458 analytical method Methods 0.000 claims abstract description 12
- 238000005498 polishing Methods 0.000 claims abstract description 9
- 238000005070 sampling Methods 0.000 claims description 75
- 238000013528 artificial neural network Methods 0.000 claims description 9
- 238000001514 detection method Methods 0.000 claims description 8
- 238000012163 sequencing technique Methods 0.000 claims description 8
- 238000005286 illumination Methods 0.000 claims description 5
- 238000004590 computer program Methods 0.000 claims description 3
- 238000000605 extraction Methods 0.000 abstract description 10
- 238000004364 calculation method Methods 0.000 abstract description 4
- 238000004088 simulation Methods 0.000 abstract 1
- 238000005516 engineering process Methods 0.000 description 5
- 230000000877 morphologic effect Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 238000004422 calculation algorithm Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000002372 labelling Methods 0.000 description 2
- 230000001795 light effect Effects 0.000 description 2
- 238000012549 training Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000010339 dilation Effects 0.000 description 1
- 230000008676 import Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/04—Context-preserving transformations, e.g. by using an importance map
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/10—Image enhancement or restoration using non-spatial domain filtering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/20—Image enhancement or restoration using local operators
- G06T5/30—Erosion or dilatation, e.g. thinning
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/187—Segmentation; Edge detection involving region growing; involving region merging; involving connected component labelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/20—Image preprocessing
- G06V10/26—Segmentation of patterns in the image field; Cutting or merging of image elements to establish the pattern region, e.g. clustering-based techniques; Detection of occlusion
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10004—Still image; Photographic image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20081—Training; Learning
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20084—Artificial neural networks [ANN]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20092—Interactive image processing based on input by user
- G06T2207/20104—Interactive definition of region of interest [ROI]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20212—Image combination
- G06T2207/20216—Image averaging
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20212—Image combination
- G06T2207/20221—Image fusion; Image merging
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Multimedia (AREA)
- Geometry (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种路面图像中裂缝比例检测方法及系统,包括:采集路面裂缝图像,对二值化路面裂缝图像进行闭运算;基于三点布光原理在模拟器中设置三点光源,对路面裂缝图像进行打光;利用闭运算结果以及打光模拟得到的图像进行引导滤波;基于暗通道先验与图像融合,对滤波后图像进行处理,得到路面裂缝灰度图像;依次进行二值化、图像细化以及连通域分析、顶帽运算,去裂缝根节点,以各连通域的端点为中心构建ROI,连接断裂的裂缝,得到完整裂缝图像,进行膨胀操作,计算裂缝像素所占图像的比例。本发明能够提高裂缝提取精度和裂缝比例计算精度。
Description
技术领域
本发明涉及计算机视觉、裂缝提取技术领域,具体涉及一种路面图像中裂缝比例检测方法及系统。
背景技术
路面的情况复杂多样,一般的人力检测不仅消耗严重而且时效性和可靠性都比较差,因此采用计算机自动检测公路路面裂缝的技术成为该技术领域的研究热点之一。
现有裂缝检测容易受到路面光照、阴影等因素的影响,且裂缝的提取通常为断续的多条曲线,这对于路面裂缝比例的检测存在很大的影响。
发明内容
本发明的目的在于针对上述现有技术存在的缺陷,提出一种路面图像中裂缝比例检测方法及系统。
一种路面图像中裂缝比例检测方法,该方法包括:
步骤1,采集路面裂缝图像,提取路面裂缝图像中的裂缝,得到二值化路面裂缝图像;
步骤2,对二值化路面裂缝图像进行闭运算,得到精细二值化路面裂缝图像;
步骤3,基于三点布光原理在模拟器中设置三点光源,对路面裂缝图像进行打光,得到三幅突显纹理的路面裂缝图像I1、I2、I3;
步骤4,精细二值化路面裂缝图像作为引导图,分别对I1、I2、I3进行引导滤波,得到滤波后路面裂缝图像J1、J2、J3;
步骤5,基于暗通道先验与图像融合,对J1、J2、J3进行处理,得到路面裂缝灰度图像M;
步骤6,对路面裂缝灰度图像M依次进行二值化、图像细化以及连通域分析,得到图像N;
步骤7,对图像N进行顶帽运算,得到去裂缝根节点图像P;
步骤8,以去裂缝根节点图像P中各连通域的端点为中心构建ROI,连接图像P中断裂的裂缝,得到完整裂缝图像Q;
步骤9,对完整裂缝图像Q进行膨胀操作,计算裂缝像素所占图像的比例。
进一步地,所述步骤5包括:
步骤5a,基于暗通道先验分别对J1、J2、J3进行处理,得到三张单通道图像L1、L2、L3;
步骤5b,对L1、L2、L3进行逐像素融合得到路面裂缝灰度图像M:
M(i,j)=Average(L1(i,j),L2(i,j),L3(i,j))
其中,Average表示均值运算,(i,j)表示像素位置。
进一步地,所述以去裂缝根节点图像P中各连通域的端点为中心构建ROI,连接图像P中断裂的裂缝包括:
将各连通域按照面积大小进行排序,从排序序列中依次取连通域作为待处理连通域,以待处理连通域的端点为中心,构建半径为R的ROI;
搜索ROI内是否存在其他连通域的端点,如果存在,则将距离中心端点最近的端点与中心端点连接。
进一步地,将距离中心端点最近的端点与中心端点连接包括:
步骤8a,分别对待处理连通域、距离中心端点最近的端点所在连通域内点进行采样,得到第一采样点坐标序列、第二采样点坐标序列;
步骤8b,将第一采样点坐标序列与第二采样点坐标序列接合,输入全连接神经网络进行分析,输出预测点坐标序列;
步骤8c,将中心端点、预测点坐标序列中的预测点、距离中心端点最近的端点连接。
进一步地,第一采样点坐标序列与第二采样点坐标序列中采样点的个数相同。
进一步地,所述第一采样点坐标序列、第二采样点坐标序列中采样点个数按照以下方式计算:
步骤a1,分别计算待处理连通域、距离中心端点最近的端点所在连通域的最小外接矩形的长宽比;
步骤a2,计算两个长宽比的比值ratio,第一采样点坐标序列与第二采样点坐标序列中采样点的个数比值与两个长宽比的比值相同;
步骤a3,将全连接神经网络的输入层数据的尺寸乘以ratio,作为第一采样点坐标序列中采样点个数,将输入层数据的尺寸减去第一采样点坐标序列中采样点个数,作为第二采样点坐标序列中采样点个数。
进一步地,所述计算裂缝像素所占图像的比例具体为:
其中,图像为膨胀后的完整裂缝图像Q。
一种路面图像中裂缝比例检测系统,该系统包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现以下步骤:
步骤1,采集路面裂缝图像,提取路面裂缝图像中的裂缝,得到二值化路面裂缝图像;
步骤2,对二值化路面裂缝图像进行闭运算,得到精细二值化路面裂缝图像;
步骤3,基于三点布光原理在模拟器中设置三点光源,对路面裂缝图像进行打光,得到三幅不同光照角度的路面裂缝图像I1、I2、I3;
步骤4,精细二值化路面裂缝图像作为引导图,分别对I1、I2、I3进行引导滤波,得到滤波后路面裂缝图像J1、J2、J3;
步骤5,基于暗通道先验与图像融合,对J1、J2、J3进行处理,得到路面裂缝灰度图像M;
步骤6,对路面裂缝灰度图像M依次进行二值化、图像细化以及连通域分析,得到图像N;
步骤7,对图像N进行顶帽运算,得到去裂缝根节点图像P;
步骤8,以去裂缝根节点图像P中各连通域的端点为中心构建ROI,连接图像P中断裂的裂缝,得到完整裂缝图像Q;
步骤9,对完整裂缝图像Q进行膨胀操作,计算裂缝像素所占图像的比例。
本发明与现有技术相比,具有以下有益效果:
1.本发明通过三点光源打光,使得路面裂缝图像光效果主次分明,互相补充,使得后续图像处理结果更具层次感,纹理更明显;进一步结合引导滤波技术,对三幅突显纹理的图像进行滤波,使图像本身突显的纹理部分更突出,并强化了其缺少以及不明显的纹理信息;更进一步结合暗通道先验技术、图像融合技术进行处理,融合三幅图像中的关键纹理信息,提高了路面裂缝的提取精度。
2.本发明对路面裂缝灰度图像进行二值化、图像细化、连通域分析去除图像中的噪声,进一步结合顶帽运算,去除裂缝根节点,无需判断区分断裂位置、非断裂位置,降低了裂缝连接过程的计算量。
3.本发明通过裂缝端点ROI启发式连接,对待处理连通域、距离中心端点最近的连通域进行采样点提取,结合神经网络技术输出预测点坐标,使得图像中的路面裂缝更加完整、更加顺滑、更加符合原始裂缝图像中裂缝的走向。根据待处理连通域、距离中心端点最近的连通域的长宽比设置采样点个数,能够进一步提高预测点的准确率,使裂缝连接更加符合原始裂缝走向。而且,本发明通过计算机视觉算法,使得路面裂缝的提取更加鲁棒。通过连通域面积筛选掉小的连通域,避免了噪声对提取结果的影响。
附图说明
图1为本发明的流程图;
图2为连通域端点示意图,附图标记1、2表示连通域端点;
图3为预测点示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图以及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
本发明提供一种路面图像中裂缝比例检测方法及系统。图1为本发明的流程图。下面通过具体实施例来进行说明。
实施例1:
路面图像中裂缝比例检测方法:
步骤1,采集路面裂缝图像,提取路面裂缝图像中的裂缝,得到二值化路面裂缝图像。
首先利用数据采集设备,对路面图像进行采集,得到路面裂缝图像。对于路面裂缝图像进行路面裂缝提取,具体的方法有很多,如基于Beamlet变换的路面图像裂缝检测方法、基于加权融合纹理的路面裂缝检测算法等,同样实施者也可基于DNN语义分割来提取路面裂缝,提取效果更鲁棒,但通常存在裂缝中断的状况。经过路面裂缝提取后,得到二值化路面裂缝图像。
步骤2,对二值化路面裂缝图像进行闭运算,得到精细二值化路面裂缝图像。
对于二值化路面裂缝图像,利用形态学操作闭运算填平小孔,弥合小裂缝,得到精细二值化路面裂缝图像。
步骤3,基于三点布光原理在模拟器中设置三点光源,对路面裂缝图像进行打光,得到三幅突显纹理的路面裂缝图像I1、I2、I3。
基于三点布光原理,对原图进行打光,得到三幅突显纹理的路面裂缝图像I1、I2、I3。具体的步骤可将图像导入相关软件(如Maya等)或引擎中,进行模拟三点布光,根据虚拟相机位置在路面裂缝不同位置部署主体光、辅助光、轮廓光,优选地,分别开启主体光、辅助光、轮廓光,当开启某一位置布光时,其他位置的布光关闭,得到三幅突显纹理的光照化路面裂缝图像。实施者也可以将三点布光全部打开,依次关闭主体光、辅助光、轮廓光,当关闭某一位置布光时,其他位置布光开启,得到三幅突显纹理的光照化路面裂缝图像。实施者也可基于光照分布变化来对图像进行滤波处理,常见于照片的滤镜。但是,第一种优选的实施方式,不仅纹理更明显,而且三幅图像的纹理更具层次,有助于提高后续裂缝的提取精度。
通过三点布光使得路面图像光效果主次分明,互相补充,使得图像中的裂缝具有透视感、空间感。
进一步的利用引导滤波,对三幅突显纹理的路面裂缝图像进行处理。
步骤4,精细二值化路面裂缝图像作为引导图,分别对I1、I2、I3进行引导滤波,得到滤波后路面裂缝图像J1、J2、J3。
引导滤波是一种图像滤波技术,通过一张引导图,对初始输入图像进行滤波处理,使得最后的输出图像大体上与初始图像相似,但是纹理部分与引导图相似。
此处,对于每一幅突显纹理的路面裂缝图像进行引导滤波处理,引导图为精细二值化路面裂缝图像。通过上述原理,得到三张引导滤波后的图像J1、J2、J3。
步骤5,基于暗通道先验与图像融合,对J1、J2、J3进行处理,得到路面裂缝灰度图像M。
对于J1、J2、J3里的每一张图像,都为三通道图像。然后此处基于暗通道先验对图像进行处理。暗通道原则的表达式如下:
P(i,j)=min{R(i,j),G(i,j),B(i,j)}
其中(i,j)是指输入图像的像素位于输入图像某通道的第i行第j列,R(i,j),G(i,j),B(i,j)分别代表输入图像在位置(i,j)处的像素中R(红),G(绿),B(蓝)通道分量值,P(i,j)是指位置(i,j)处的输入图像的R,G,B通道中最小亮度值。
至此,得到三张引导滤波后的图像J1、J2、J3每张图像的暗通道图像L1、L2、L3。
进一步的对于三张暗通道图像(单通道灰度图),进行图像融合,求平均操作。
对L1、L2、L3进行逐像素融合得到路面裂缝灰度图像M:
M(i,j)=Average(L1(i,j),L2(i,j),L3(i,j))
其中,Average表示均值运算,(i,j)表示像素位置。
对于每一个相同像素位置的三张暗通道图像像素值进行求平均,最终得到一张灰度图M,该图像具有和线扫描相机得到的光照均一结果一样,能够明显地得到缺陷信息,且可以直接阈值化。
步骤6,对路面裂缝灰度图像M依次进行二值化、图像细化以及连通域分析,得到图像N。
进一步的,对路面裂缝灰度图像M进行二值化处理,此处可采用OTSU、Kittle等阈值化算法。最终得到二值化图像M’。
对二值化图像M’进行图像细化操作(Image Thinning),使裂缝变得比较细,得到图像M”。
M”图像中包含了许多噪点,因此需要对图像进行进一步处理,以得到更加清晰的路面裂缝结果N。本实施例采用连通域面积筛选处理,将小于阈值的噪点连通域过滤。对二值化图像M”进行连通域分析,得到图像中的每个连通域,同时将小于阈值的噪点连通域过滤,以去除很小的连通域。其中连通区域的邻接方式(4邻接、8邻接),实施者可自由选择,效果较佳即可。连通域面积可直接通过相关接口实现,当面积小于阈值时,进行过滤,该阈值为人为经验,与图像大小也有关。
步骤7,对图像N进行顶帽运算,得到去裂缝根节点图像P。
对图像N进行顶帽运算,去除图像中的裂缝根节点,以减少后续计算量,并防止裂缝的误连接。裂缝的根节点即裂缝交叉点,对于顶帽运算,建议使用椭圆结构的模板来去除裂缝根节点,使得图像更加平滑。
步骤8,以去裂缝根节点图像P中各连通域的端点为中心构建ROI,连接图像P中断裂的裂缝,得到完整裂缝图像Q。
对步骤7结果中的连通域的端点进行启发式连接,以连接间断的路面裂缝。由于经过步骤7的处理,祛除了裂缝根节点,利用算法例如最小外接凸多边形容易得到图像P中的裂缝连通域的端点。将各连通域按照面积大小进行排序,依次从排序序列中取连通域作为待处理连通域,以待处理连通域的端点为中心,构建半径为R的ROI;搜索ROI内是否存在其他连通域的端点,如果存在,则将距离中心端点最近的端点与中心端点连接。具体的方法如下:
1)首先求得连通域端点坐标,可通过在连通域遍历每个像素并分析其邻接的像素的值来判断,连通域端点示意图如图2所示,图中附图标记1、2为裂缝端点,是连通域的起点或终点;
2)选取面积最大的连通域作为初始操作连通域;
3)初始操作连通域中,将所有端点作为圆心形成一个半径为R的圆形,当其他连通域的端点位于该圆形区域时连接端点,该圆形区域称为ROI区域;
4)当该圆形区域存在连通域的多个端点时,选取距离圆心最近的端点进行连接;
5)然后选取第二面积的连通域重复上述操作,直至所有连通域判断结束。
需要说明的时,上述所述圆形半径R,对于512*512的图像来说,半径R的经验值为12,即12个像素长。
众所周知,通过拟合方法得到的函数通常每个x坐标只对应1个y坐标,因此当两个连通域存在一个x坐标对应多个y坐标或一个y坐标对应多个x坐标时,拟合方法将会非常不准确,极大影响后续裂缝比例的计算。
对于上述端点连接,分别对待处理连通域、距离中心端点最近的端点所在连通域内点进行采样,得到第一采样点坐标序列、第二采样点坐标序列;将第一采样点坐标序列与第二采样点坐标序列接合,输入全连接神经网络进行分析,输出预测点坐标序列;将中心端点、预测点坐标序列中的预测点、距离中心端点最近的端点连接。具体方法如下:当上述3)、4)判断成功时,即可得知哪两个连通域进行端点连接,由于直接连接为直线连接,不平滑且不符合逻辑,因此此处使用全连接网络进行两端点间点坐标的预测。对于全连接方法,由于全连接网络的输入为定长,因此首先对两个连通域的点进行采样。然后采样连通域内的N个点坐标,N为超参数,本实施例选择为10个。采样的方法有很多种,如随机采样、基于特定规则的采样等,实施者可自由选择。基于特定规则的采样包括,利用最小外接矩形将连通域分块,每个分块内采样一个点。至此通过上述方法,即可得到两个连通域总共20个点坐标。然后训练全连接神经网络,该网络输入为[B,20],输出为[B,M],M为超参数,同样为人为设置,M决定了两个要连接的连通域端点间预测的点个数。该网络损失函数为回归型损失函数,如均方差等,标签数据同样基于人为标注,根据原始图像中断裂位置的走向标注预测点,预测点标注示意如图3所示,点的个数与M应保持一致。最终网络输出得到M个坐标,然后进行连接即可,连接的规则实施者可自由设置,这里不做约束。进一步地,为了隔离工况,降低全连接网络训练复杂度,以待处理连通域中心端点为原点,坐标轴方向仍为原始图像的坐标轴方向,构建新的坐标系,将采样点的坐标变换到新的坐标系。
为了进一步提高预测点精度,第一采样点坐标序列、第二采样点坐标序列中采样点个数按照以下方式计算:对于每个连通域寻找最小外接矩形端点,(Xmin,y0)、(Xmax,y1)、(x0,Ymin)、(x1,Ymax),基于四个端点构建最小外接矩形,分别计算待处理连通域、距离中心端点最近的端点所在连通域的最小外接矩形的长宽比;计算两个长宽比的比值ratio,第一采样点坐标序列与第二采样点坐标序列中采样点的个数比值与两个长宽比的比值相同;将全连接神经网络的输入层数据的尺寸乘以ratio,作为第一采样点坐标序列中采样点个数,将输入层数据的尺寸减去第一采样点坐标序列中采样点个数,作为第二采样点坐标序列中采样点个数。此外,还有一种计算方式:分别计算待处理连通域、距离中心端点最近的端点所在连通域拟合得到的曲线的极值个数,极值包括极大值、极小值;分别计算两个连通域端点的最大x轴距离,即将坐标中最大x值减去最小x值,然后计算极值个数与最大x轴距离的比值得到r1、r2,该比值表征极值在x轴方向上的密度;第一采样点坐标序列与第二采样点坐标序列中采样点的个数比值与r1/r2相同;将全连接神经网络的输入层数据的尺寸乘以r1/r2,作为第一采样点坐标序列中采样点个数,将输入层数据的尺寸减去第一采样点坐标序列中采样点个数,作为第二采样点坐标序列中采样点个数。
通过步骤8,即可得到完整的、不间断的路面裂缝二值化图像—完整裂缝图像Q,但该图像的裂缝比较细。
步骤9,对完整裂缝图像Q进行膨胀操作,计算裂缝像素所占图像的比例。
对于图像Q构建形态学结构元素,对图像Q进行形态学膨胀操作,使得裂缝变得宽大,得到图像Q’。
构建形态学结构元素。一种实施例是设计形态学结构元素K,形态元结构元素的尺寸是3×3,具体可以采用十字结构进行膨胀操作,使得路面裂缝连接线变得较宽。
然后计算图像Q’中的裂缝比例。对于路面裂缝二值化图像Q’,将代表路面裂缝的像素数量与图像的总像素数相除,其结果作为路面图像的裂缝比例。公式如下:
至此,即可得到每张路面图像的裂缝比例。
实施例2:
一种路面图像中裂缝比例检测系统,该系统包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现以下步骤:
步骤1,采集路面裂缝图像,提取路面裂缝图像中的裂缝,得到二值化路面裂缝图像;
步骤2,对二值化路面裂缝图像进行闭运算,得到精细二值化路面裂缝图像;
步骤3,基于三点布光原理在模拟器中设置三点光源,对路面裂缝图像进行打光,得到三幅突显纹理的路面裂缝图像I1、I2、I3;
步骤4,精细二值化路面裂缝图像作为引导图,分别对I1、I2、I3进行引导滤波,得到滤波后路面裂缝图像J1、J2、J3;
步骤5,基于暗通道先验与图像融合,对J1、J2、J3进行处理,得到路面裂缝灰度图像M;
步骤6,对路面裂缝灰度图像M依次进行二值化、图像细化以及连通域分析,得到图像N;
步骤7,对图像N进行顶帽运算,得到去裂缝根节点图像P;
步骤8,以去裂缝根节点图像P中各连通域的端点为中心构建ROI,连接图像P中断裂的裂缝,得到完整裂缝图像Q;
步骤9,对完整裂缝图像Q进行膨胀操作,计算裂缝像素所占图像的比例。
以上实施例仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种路面图像中裂缝比例检测方法,其特征在于,该方法包括:
步骤1,采集路面裂缝图像,提取路面裂缝图像中的裂缝,得到二值化路面裂缝图像;
步骤2,对二值化路面裂缝图像进行闭运算,得到精细二值化路面裂缝图像;
步骤3,基于三点布光原理在模拟器中设置三点光源,对路面裂缝图像进行打光,得到三幅突显纹理的路面裂缝图像I1、I2、I3;
步骤4,精细二值化路面裂缝图像作为引导图,分别对I1、I2、I3进行引导滤波,得到滤波后路面裂缝图像J1、J2、J3;
步骤5,基于暗通道先验与图像融合,对J1、J2、J3进行处理,得到路面裂缝灰度图像M;
步骤6,对路面裂缝灰度图像M依次进行二值化、图像细化以及连通域分析,得到图像N;
步骤7,对图像N进行顶帽运算,得到去裂缝根节点图像P;
步骤8,以去裂缝根节点图像P中各连通域的端点为中心构建ROI,连接图像P中断裂的裂缝,得到完整裂缝图像Q;
步骤9,对完整裂缝图像Q进行膨胀操作,计算裂缝像素所占图像的比例。
2.如权利要求1所述的方法,其特征在于,所述步骤5包括:
步骤5a,基于暗通道先验分别对J1、J2、J3进行处理,得到三张单通道图像L1、L2、L3;
步骤5b,对L1、L2、L3进行逐像素融合得到路面裂缝灰度图像M:
M(i,j)=Average(L1(i,j),L2(i,j),L3(i,j))
其中,Average表示均值运算,(i,j)表示像素位置。
3.如权利要求1所述的方法,其特征在于,所述以去裂缝根节点图像P中各连通域的端点为中心构建ROI,连接图像P中断裂的裂缝包括:
将各连通域按照面积大小进行排序,从排序序列中依次取连通域作为待处理连通域,以待处理连通域的端点为中心,构建半径为R的ROI;
搜索ROI内是否存在其他连通域的端点,如果存在,则将距离中心端点最近的端点与中心端点连接。
4.如权利要求3所述的方法,其特征在于,将距离中心端点最近的端点与中心端点连接包括:
步骤8a,分别对待处理连通域、距离中心端点最近的端点所在连通域内的点进行采样,得到第一采样点坐标序列、第二采样点坐标序列;
步骤8b,将第一采样点坐标序列与第二采样点坐标序列接合,输入全连接神经网络进行分析,输出预测点坐标序列;
步骤8c,将中心端点、预测点坐标序列中的预测点、距离中心端点最近的端点连接。
5.如权利要求4所述的方法,其特征在于,第一采样点坐标序列与第二采样点坐标序列中采样点的个数相同。
6.如权利要求4所述的方法,其特征在于,所述第一采样点坐标序列、第二采样点坐标序列中采样点个数按照以下方式计算:
步骤a1,分别计算待处理连通域、距离中心端点最近的端点所在连通域的最小外接矩形的长宽比;
步骤a2,计算两个长宽比的比值ratio,第一采样点坐标序列与第二采样点坐标序列中采样点的个数比值与两个长宽比的比值相同;
步骤a3,将全连接神经网络的输入层数据的尺寸乘以ratio,作为第一采样点坐标序列中采样点个数,将输入层数据的尺寸减去第一采样点坐标序列中采样点个数,作为第二采样点坐标序列中采样点个数。
7.如权利要求1所述的方法,其特征在于,所述计算裂缝像素所占图像的比例具体为:
其中,图像为膨胀后的完整裂缝图像Q。
8.一种路面图像中裂缝比例检测系统,该系统包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述程序时实现以下步骤:
步骤1,采集路面裂缝图像,提取路面裂缝图像中的裂缝,得到二值化路面裂缝图像;
步骤2,对二值化路面裂缝图像进行闭运算,得到精细二值化路面裂缝图像;
步骤3,基于三点布光原理在模拟器中设置三点光源,对路面裂缝图像进行打光,得到三幅不同光照角度的路面裂缝图像I1、I2、I3;
步骤4,精细二值化路面裂缝图像作为引导图,分别对I1、I2、I3进行引导滤波,得到滤波后路面裂缝图像J1、J2、J3;
步骤5,基于暗通道先验与图像融合,对J1、J2、J3进行处理,得到路面裂缝灰度图像M;
步骤6,对路面裂缝灰度图像M依次进行二值化、图像细化以及连通域分析,得到图像N;
步骤7,对图像N进行顶帽运算,得到去裂缝根节点图像P;
步骤8,以去裂缝根节点图像P中各连通域的端点为中心构建ROI,连接图像P中断裂的裂缝,得到完整裂缝图像Q;
步骤9,对完整裂缝图像Q进行膨胀操作,计算裂缝像素所占图像的比例。
9.如权利要求8所述的系统,其特征在于,所述步骤5包括:
步骤5a,基于暗通道先验分别对J1、J2、J3进行处理,得到三张单通道图像L1、L2、L3;
步骤5b,对L1、L2、L3进行逐像素融合得到路面裂缝灰度图像M:
M(i,j)=Mean(L1(i,j),L2(i,j),L3(i,j))
其中,Mean表示均值运算,(i,j)表示像素位置。
10.如权利要求8所述的系统,其特征在于,所述以去裂缝根节点图像P中各连通域的端点为中心构建ROI,连接图像P中断裂的裂缝包括:
将各连通域按照面积大小进行排序,依次从排序序列中取最大面积连通域作为待处理连通域,以待处理连通域的端点为中心,构建半径为R的ROI;
搜索ROI内是否存在其他连通域的端点,如果存在,则将距离中心端点最近的端点与中心端点连接。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210796996.6A CN114998096B (zh) | 2022-07-06 | 2022-07-06 | 一种路面图像中裂缝比例检测方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210796996.6A CN114998096B (zh) | 2022-07-06 | 2022-07-06 | 一种路面图像中裂缝比例检测方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114998096A CN114998096A (zh) | 2022-09-02 |
CN114998096B true CN114998096B (zh) | 2023-11-03 |
Family
ID=83019106
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210796996.6A Active CN114998096B (zh) | 2022-07-06 | 2022-07-06 | 一种路面图像中裂缝比例检测方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114998096B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2019066263A (ja) * | 2017-09-29 | 2019-04-25 | 清水建設株式会社 | ひび割れ検出装置、ひび割れ検出方法、および、コンピュータプログラム |
CN111489352A (zh) * | 2020-04-29 | 2020-08-04 | 安徽国钜工程机械科技有限公司 | 一种基于数字图像处理的隧道缝隙检测与测量方法及装置 |
CN112419250A (zh) * | 2020-11-13 | 2021-02-26 | 湖北工业大学 | 路面裂缝数字图像提取、裂缝修补及裂缝参数计算方法 |
-
2022
- 2022-07-06 CN CN202210796996.6A patent/CN114998096B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2019066263A (ja) * | 2017-09-29 | 2019-04-25 | 清水建設株式会社 | ひび割れ検出装置、ひび割れ検出方法、および、コンピュータプログラム |
CN111489352A (zh) * | 2020-04-29 | 2020-08-04 | 安徽国钜工程机械科技有限公司 | 一种基于数字图像处理的隧道缝隙检测与测量方法及装置 |
CN112419250A (zh) * | 2020-11-13 | 2021-02-26 | 湖北工业大学 | 路面裂缝数字图像提取、裂缝修补及裂缝参数计算方法 |
Non-Patent Citations (2)
Title |
---|
Effect of stress ratio or mean stress on fatigue delamination growth in composites: Critical review;Rafiullah Khan;《Composite Structures》;215-227 * |
线阵相机拍摄路面图像的灰度校正及应用;孟祥成;《计算机应用研究》;355-357 * |
Also Published As
Publication number | Publication date |
---|---|
CN114998096A (zh) | 2022-09-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111145209B (zh) | 一种医学图像分割方法、装置、设备及存储介质 | |
WO2018107939A1 (zh) | 一种基于边缘完备度的图像分割最优识别方法 | |
CN111626146B (zh) | 一种基于模板匹配的合并单元格表格分割识别方法 | |
CN108986122B (zh) | 室内泊车引导地图智能重构方法 | |
CN112837290B (zh) | 一种基于种子填充算法的裂缝图像自动识别方法 | |
CN103295199B (zh) | 古壁画的裂缝智能修复辅助系统 | |
CN113240623B (zh) | 一种路面病害检测方法及装置 | |
CN111582294A (zh) | 一种构建用于表面缺陷检测的卷积神经网络模型的方法及其利用 | |
CN113096121B (zh) | 基于交叉断裂力学和图像处理的路面裂缝检测方法和系统 | |
CN113012124B (zh) | 一种鞋印孔洞和嵌入物特征检测及描述方法 | |
CN101106716A (zh) | 一种分水岭图象分割处理方法 | |
CN110969620A (zh) | 一种磁瓦波纹类缺陷的检测方法及装置 | |
CN114005081A (zh) | 一种烟丝异物智能检测装置及方法 | |
CN116539619A (zh) | 产品缺陷检测方法、系统、装置及存储介质 | |
CN110246139B (zh) | 基于双阈值的浮游生物原位图像roi快速提取方法 | |
CN114998096B (zh) | 一种路面图像中裂缝比例检测方法及系统 | |
WO2020164042A1 (zh) | 基于边界提取的区域合并图像分割算法 | |
CN113935961A (zh) | 一种鲁棒性乳腺钼靶mlo视角图像胸肌分割方法 | |
JP3889361B2 (ja) | 核領域認識法および細胞系譜作成法 | |
KR19990074775A (ko) | 지문 특징 추출방법 | |
CN105894501B (zh) | 一种高分辨率遥感图像单木探测和树冠描绘方法 | |
CN105069766A (zh) | 一种基于汉字图像轮廓特征描述的碑文修复方法 | |
CN116863134A (zh) | 一种隧道衬砌裂缝长度与宽度的检测分割方法及系统 | |
CN116258660A (zh) | 口腔鳞癌侵袭深度智能测量系统 | |
CN110544260B (zh) | 融合自学习语义特征与设计特征的遥感影像目标提取方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |