首页 / 专利库 / 视听技术与设备 / 均值漂移过程 / 一种用于水损害检测的探地雷达图谱自适应选取方法

一种用于损害检测的探地雷达图谱自适应选取方法

阅读:564发布:2020-05-08

专利汇可以提供一种用于损害检测的探地雷达图谱自适应选取方法专利检索,专利查询,专利分析的服务。并且本 发明 提供了一种用于 水 损害检测的探地雷达图谱自适应选取方法,该方法根据探地雷达图谱的 对比度 来自适应选取探地雷达图谱,该方法按照以下步骤进行:步骤1,读取预处理后的GPR数据;步骤2,图片 分辨率 调整;步骤3,将数据集输入到识别模型;步骤4,输出水损害结果;步骤5,利用初始随机图谱数据集,进行检测目标的有无判断;步骤6,增量式随机生成图谱,进行合适对比度图谱的选取;本发明的识别 算法 与 深度学习 模型(或者图像分类模型)相结合,可有效、快速、自动寻找出合适的B-scan图谱,实现基于GPR图谱的自动识别和检测,同时提高识别 精度 。,下面是一种用于损害检测的探地雷达图谱自适应选取方法专利的具体信息内容。

1.一种用于损害检测的探地雷达图谱自适应选取方法,其特征在于,该方法根据探地雷达数据来自适应选取合适对比度的探地雷达图谱,该方法按照以下步骤进行:
步骤1,读取预处理后的GPR数据:
对GPR数据进行预处理后,在设定的对比度取值范围内,随机生成不同对比度的雷达图谱,构建初始随机图谱数据集,所述的初始随机图谱数据集的大小为N张图片;
步骤2,图片分辨率调整:
初始随机图谱数据集定义为RID数据集,将RID数据集直接缩放到224×224,定义为RBD数据集;
将水损害初始图谱数据集的分辨率直接缩放到224×224,获得RBD数据集;
步骤3,将数据集输入到识别模型:
将步骤二获得的RBD数据集输入到识别模型,经过识别模型运算后,执行步骤四;
所述的识别模型的图片输入分辨率大小为224×224,图片输出分辨率大小为224×
224;
所述的识别模型为混合式深度学习模型,所述的混合式深度学习模型由两部分构成,特征提取采用ResNet50,目标检测采用YOLO V2框架
步骤4,输出水损害结果:
对步骤三的识别模型的输出结果进行后处理,所述的后处理步骤为:
步骤S41,判断输出结果中的图谱的候选框BBoxes数量,如果候选框BBoxes数量>1,则进行步骤S42,否则直接输出结果;
步骤S42:判断候选框BBoxes是否有重叠,若有重叠,则进行步骤S43,否则直接输出结果;
步骤S43:判断重叠候选框对应的标签名是否相同,若相同,则合并的候选框对应的标签名不变;若不同,则表明同时含有水损害和桥接缝两类标签名,标签名输出为桥接缝Joint;
步骤S44:合并候选框,将相交的候选框x和y方向取最小值,w和h取最大值,合并的候选框坐标为[xmin,ymin,wmax,hmax];
步骤S45:输出结果,在识别模型的输出结果中,将输出图片分辨率调整至等于所述的损害的初始图谱数据集的图片分辨率,输出结果为带目标的标签名及目标对应的候选框位置BBoxes(x,y,w,h)的图谱;
步骤5,利用初始随机图谱数据集,进行检测目标的有无判断:
步骤S51,将步骤四的输出结果转为图片上像素点对应的矩阵Ai,Ai定义为:
where 1≤m≤H0,1≤n≤W0
式中:H0为识别模型输出图谱的图片高度;W0为识别模型输出图谱的图片宽度;
将RID数据集中的N张图片对应的矩阵Ai进行求和,并求其均值,获得均值矩阵A,A定义为:
步骤S52,设定k1=0.8,θ0=0.5,对均值矩阵A进行按照以下公式进行更新,获得更新后的均值矩阵A;
A(A式中:
k1为目标关联系数;
θ0为含有目标时,矩阵A中的最小值,低于此值,则无目标;
max(max(A))为均值矩阵中的最大值;
步骤S53,根据更新后的均值矩阵A,按照以下公式获得目标有无的判断条件T,当T=1,则表示有目标,当T=0,则表示没有目标;
步骤6,增量式随机生成图谱,进行合适对比度图谱的选取:
当图片含有目标后,进行初始化判断:
步骤S61,若Flag=0,则表示为第一次生成随机图谱样本集,即初始化样本集阶段,不进入后续选取判断,设置Flag=1,然后额外增加5%·N张图片作为样随机图谱数据集的本,样本中图片的总数为N=(1+5%)N,返回到步骤2;
步骤S62,若Flag≠0,则表示非初始化阶段,设置一个图片相关系数,选择与平均矩阵A最大相关系数的图片,即为合适对比度的图谱;
相关系数Ri定义为:
式中:Ri为第i张图谱对应的矩阵Ai与均值矩阵A的相关系数;m为高度方向上的坐标值;
n为宽度方向上的坐标值;μA为均值矩阵A的总均值; 为矩阵Ai的总均值;
选取过程的终止条件为:
STOP=abs(F1-F1pre)<0.01&&F1>0.8;
式中:F1为深度学习的评价指标;F1Pre为上一次深度学习的评价指标,初始时为0;TP为真实的目标区域;FP为错误识别的真值,代表未识别的真值,错误判断为负值或背景;FN为错误识别的负值,即把背景当成目标;
当终止条件不满足时,将该次计算的指标F1付给变量F1Pre,然后返回步骤S61,增大样本集,重新进行选取;
当终止条件满足时,系统将输出合适对比度的图谱。
2.如权利要求1所述的用于水损害检测的探地雷达图谱自适应选取方法,其特征在于,所述的GPR数据的获取方法为:采用GPR系统对沥青路面进行现场数据采集,现场数据采集过程中,确定路面出现唧浆或泛白的损害区域,获取损害区域对应的GPR数据。
3.如权利要求1所述的用于水损害检测的探地雷达图谱自适应选取方法,其特征在于,所述的现场数据采集过程中,采样参数要求采样间距<15cm,天线频率>1.6GHz,采样频率为天线主频的10~20倍。
4.如权利要求1所述的用于水损害检测的探地雷达图谱自适应选取方法,其特征在于,所述的预处理过程为采用直流漂移校正算法、地面校正算法、背景扣除算法、带通滤波算法和进行滑动平均算法进行预处理。
5.如权利要求1所述的用于水损害检测的探地雷达图谱自适应选取方法,其特征在于,所述的设定的对比度取值范围为0.5~1.8。
6.如权利要求1所述的用于水损害检测的探地雷达图谱自适应选取方法,其特征在于,所述的N=100。

说明书全文

一种用于损害检测的探地雷达图谱自适应选取方法

技术领域

[0001] 本发明属于自动检测领域,涉及探地雷达图谱,具体涉及一种用于水损害检测的探地雷达图谱自适应选取方法。

背景技术

[0002] 探地雷达(Ground Penetrating Radar,简称GPR)是以雷达脉冲波探测地表以下状况并成像的仪器,其原理是利用天线发射和接收高频电磁波来探测介质内部物质特性。GPR使用高频率且通常被极化的无线电波,将电波发射入地表之下,当电磁波撞击到埋在地表下的物体或到达介电常数变化的边界时,天线接收到的反射波会记录下反射回波的讯号差异。由于探地雷达探测的高精度、高效率以及无损的特点,GPR被广泛用于考古、矿产勘查、灾害地质调查、岩土工程勘察、工程质量检测、建筑结构检测以及军事目标探测等众多领域。
[0003] 由于GPR可进行连续、快速、无破损检测,目前GPR已在公路交通上应用,如在隧道衬底脱空、路面脱空、建筑或者桥面板传杆或筋的识别、地下管线上得到了广泛的应用。然而由于探地雷达成像解释困难,GPR图谱解释往往依赖GPR专家的经验。并且,在实际项目中,路面结构众多,缺陷种类众多,缺陷分布的深度不同,如沥青路面主要是水损害、破碎和裂缝,集中在沥青层(深度0~22cm);水泥路面的缺陷主要是脱空,集中在半刚性(深度>24cm);“白加黑”路面兼备二者缺陷,深度覆盖沥青层、水泥路面、半刚性和路基,深度0~1.5m。为实现不同缺陷的检测,往往需要雷达天线阵列组合。同时公路里程长,路面调查将产生海量数据,依靠人工经验识别的方法将难以满足公路智能化检测要求,迫切需要建立路面缺陷的数据库,引入深度学习等智能识别方法,实现高速公路缺陷(或其他地下目标体)的智能检测。
[0004] 已有相关研究将探地雷达图谱与深度学习方法相结合进行自动识别,目前主要集中在双曲线目标(如钢筋、裂缝等孤立目标体)的识别,采用深度学习模型集中在雷达图谱双曲线特征的识别以及定位,已达到较高的识别精度。数据集依靠专家经验优选对比图强的GPR图片来构建数据集。往往,GPR数据(A-SCAN数据)会经过去除零漂、增益处理、带通滤波、道均衡、二维滤波、偏移归位、反褶积等方法后,将目标体在背景图谱中突出后再截取雷达图谱(B-SCAN),构建合适的数据集。在此基础上,训练对应的深度学习模型(或者图像分类模型),提高了雷达图谱的自动识别效率和精度,为GPR的自动应用提供了一种新的分析方法。
[0005] 然而,现有基于雷达图谱目标体识别的深度学习研究中,所有GPR图谱的对比度(plot scale)选择均依靠人工经验选择,并且集中在特征明显的双曲线目标。好的对比度的图谱可利于专家或者深度模型正确检测出目标体,而不合理对比度的图谱将得到错误的结果。采用手动选择雷达图谱的方式阻碍了GPR的智能检测应用。
[0006] 以常见路面缺陷为例,水损害是造成沥青路面早期损害的主要原因之一。水损害的快速检测一直是一个难点,笔者提出的发明专利(一种基于探地雷达信号时频统计特征的水损害识别方法,201910100046.3)实现了水损害的检测和自动分析,从GPR数据度分析,未从图谱角度分析。利用GPR图谱是一种有效判别水损害和桥接缝的方法,水损害的缺陷并不是双曲线特征,而是呈现更加复杂的图谱特征,未有相关利用图谱的检测研究。
[0007] 由于缺陷含水率不同的原因,同一套GPR设备在同一条路上进行的一次检测中,不同水损害区域对应合适图谱的对比度值均不同,这种方式阻碍了GPR在自动检测中的应用。因此迫切需要一种根据GPR数据自动选择合适对比度雷达图谱的方法,然后将图谱输入给识别模型(如深度模型),实现目标体的自动识别。

发明内容

[0008] 针对上述现有技术的不足与缺陷,本发明的目的在于提供一种用于水损害检测的探地雷达图谱自适应选取方法,解决现有技术中基于GPR图谱依旧依赖人工选择合适对比度GPR图谱的问题。
[0009] 为了解决上述技术问题,本发明采用如下技术方案予以实现:
[0010] 一种用于水损害检测的探地雷达图谱自适应选取方法,该方法根据探地雷达图谱的对比度来自适应选取探地雷达图谱,该方法按照以下步骤进行:
[0011] 步骤1,读取预处理后的GPR数据:
[0012] 对GPR数据进行预处理后,在设定的对比度取值范围内,随机生成不同对比度的雷达图谱,构建初始随机图谱数据集,所述的初始随机图谱数据集的大小为N张图片;
[0013] 步骤2,图片分辨率调整:
[0014] 初始随机图谱数据集定义为RID数据集,将RID数据集直接缩放到224×224,定义为RBD数据集;
[0015] 将水损害初始图谱数据集的分辨率直接缩放到224×224,获得RBD数据集;
[0016] 步骤3,将数据集输入到识别模型:
[0017] 将步骤二获得的RBD数据集输入到识别模型,经过识别模型运算后,执行步骤四;
[0018] 所述的识别模型的图片输入分辨率大小为224×224,图片输出分辨率大小为224×224;
[0019] 所述的识别模型为混合式深度学习模型,所述的混合式深度学习模型由两部分构成,特征提取采用ResNet50,目标检测采用YOLO V2框架
[0020] 步骤4,输出水损害结果:
[0021] 对步骤三的识别模型的输出结果进行后处理,所述的后处理步骤为:
[0022] 步骤S41,判断输出结果中的图谱的候选框BBoxes数量,如果候选框BBoxes数量>1,则进行步骤S42,否则直接输出结果;
[0023] 步骤S42:判断候选框BBoxes是否有重叠,若有重叠,则进行步骤S43,否则直接输出结果;
[0024] 步骤S43:判断重叠候选框对应的标签名是否相同,若相同,则合并的候选框对应的标签名不变;若不同,则表明同时含有水损害和桥接缝两类标签名,标签名输出为桥接缝Joint;
[0025] 步骤S44:合并候选框,将相交的候选框x和y方向取最小值,w和h取最大值,合并的候选框坐标为[xmin,ymin,wmax,hmax];
[0026] 步骤S45:输出结果,在识别模型的输出结果中,将输出图片分辨率调整至等于所述的损害的初始图谱数据集的图片分辨率,输出结果为带目标的标签名及目标对应的候选框位置BBoxes(x,y,w,h)的图谱;
[0027] 步骤5,利用初始随机图谱数据集,进行检测目标的有无判断:
[0028] 步骤S51,将步骤四的输出结果转为图片上像素点对应的矩阵Ai,Ai定义为:
[0029] where1≤m≤H0,1≤n≤W0
[0030] 式中:H0为识别模型输出图谱的图片高度;W0为识别模型输出图谱的图片宽度;
[0031] 将RID数据集中的N张图片对应的矩阵Ai进行求和,并求其均值,获得均值矩阵A,A定义为:
[0032] 步骤S52,设定k1=0.8,θ0=0.5,对均值矩阵A进行按照以下公式进行更新,获得更新后的均值矩阵A;
[0033] A(A
[0034] 式中:
[0035] k1目标关联系数;
[0036] θ0为含有目标时,矩阵A中的最小值,低于此值,则无目标;
[0037] max(max(A))为均值矩阵中的最大值;
[0038] 步骤S53,根据更新后的均值矩阵A,按照以下公式获得目标有无的判断条件T,当T=1,则表示有目标,当T=0,则表示没有目标;
[0039]
[0040] 步骤6,增量式随机生成图谱,进行合适对比度图谱的选取:
[0041] 当图片含有目标后,进行初始化判断:
[0042] 步骤S61,若Flag=0,则表示为第一次生成随机图谱样本集,即初始化样本集阶段,不进入后续选取判断,设置Flag=1,然后额外增加5%·N张图片作为样随机图谱数据集的本,样本中图片的总数为N=(1+5%)N,返回到步骤2;
[0043] 步骤S62,若Flag≠0,则表示非初始化阶段,设置一个图片相关系数,选择与平均矩阵A最大相关系数的图片,即为合适对比度的图谱;
[0044] 相关系数Ri定义为:
[0045] 式中:Ri为第i张图谱对应的矩阵Ai与均值矩阵A的相关系数;m为高度方向上的坐标值;n为宽度方向上的坐标值;μA为均值矩阵A的总均值;μAi为矩阵Ai的总均值;
[0046] 选取过程的终止条件为:
[0047] STOP=abs(F1-F1pre)<0.01&&F1>0.8;
[0048]
[0049]
[0050]
[0051] 式中:F1为深度学习的评价指标;F1Pre为上一次深度学习的评价指标,初始时为0;TP为真实的目标区域;FP为错误识别的真值,代表未识别的真值,错误判断为负值或背景;
FN为错误识别的负值,即把背景当成目标;
[0052] 当终止条件不满足时,将该次计算的指标F1付给变量F1Pre,然后返回步骤S61,增大样本集,重新进行选取;
[0053] 当终止条件满足时,系统将输出合适对比度的图谱。
[0054] 本发明还具有如下技术特征:
[0055] 所述的GPR数据的获取方法为:采用GPR系统对沥青路面进行现场数据采集,现场数据采集过程中,确定路面出现唧浆或泛白的损害区域,获取损害区域对应的GPR数据。
[0056] 所述的现场数据采集过程中,采样参数要求采样间距<15cm,天线频率>1.6GHz,采样频率为天线主频的10~20倍。
[0057] 所述的预处理过程为采用直流漂移校正算法、地面校正算法、背景扣除算法、带通滤波算法和进行滑动平均算法进行预处理。
[0058] 所述的设定的对比度取值范围为0.5~1.8。
[0059] 所述的N=100。
[0060] 本发明与现有技术相比,具有如下技术效果:
[0061] (Ⅰ)基于粒子滤波的思想,读取待检测的GPR数据,随机生成不同对比度的GPR图谱,以此结果构建随机数据集。将生成的图谱输入给识别模型,并基于全局统计结果设定一个值,基于此来寻找合适的图谱,比较前后两次寻找到的合适图谱的一致性作为搜寻判断条件。若不满足停止条件,继续增大随机样本图谱。这种搜寻将包括所有可能的对比度的图谱,可快速有效寻找到的合适B-scan图谱,解决了合适雷达图谱自动选取的问题,为GPR图谱自动识别奠定基础。
[0062] (Ⅱ)本发明的识别算法与深度学习模型(或者图像分类模型)相结合,可有效、快速、自动寻找出合适的B-scan图谱,实现基于GPR图谱的自动识别和检测,同时提高识别精度。
[0063] (Ⅲ)本发明通过实验验证,所提出方法和基于YOLO的水损害检测模型,实现了从GPR原始数据生成合适的GPR图谱,然后完成水损害缺陷的自动检测,提高了沥青路面水损害检测的智能性。附图说明
[0064] 图1是基于GPR的水损害缺陷检测方法流程图
[0065] 图2是水损害缺陷数据集典型图谱特征。
[0066] 图3是识别模型结果的后处理流程图。
[0067] 图4是混合式模型在不同分辨率数据集下的指标对比图。
[0068] 图5是混合式模型下的水损害检测结果图。
[0069] 图6是采用ACF算法的检测结果图。
[0070] 图7是采用Cifar10模型的检测结果图。
[0071] 图8是不同对比度对应图谱图。
[0072] 图9是合适对比度的GPR雷达图谱识别算法。
[0073] 图10是相关系数随采样样本数的变化规律。
[0074] 图11是深度模型相关指标随采样样本数的变化规律。
[0075] 图12优选图谱与随机数据集总体检测的heatmap图。
[0076] 图13是优选结果与均值A的heatmap图。
[0077] 图14是正常路面与随机数据集总体检测的heatmap图。
[0078] 图15是测试数据集的随机样本数分布规律。
[0079] 图16是所提算法(IRS)与随机样本(RS)在水损害测试集的对比结果。
[0080] 图17是所提算法(IRS)与随机样本(RS)增加正常路面的对比结果。
[0081] 图中各个标号的含义为:1-1为合适对比度对应的图谱,1-2.为过小对比度对应的图谱,1-3为过大对比度对应的图谱,1-4为合适对比度对应GPR图谱中真实的水损害缺陷范围。
[0082] 以下结合实施例对本发明的具体内容作进一步详细解释说明。

具体实施方式

[0083] 对比度(plot scale)对GPR图谱的影响原因:采用GPR设置进行沥青路面调查获得雷达数据,对雷达数据进行后处理以增加目标体与背景的差异,然后将处理后的雷达数据转换为GPR图谱。图8是不同对比度对应的GPR图谱,其中1-1对应合适对比度的雷达图谱,1-4是该处GPR数据中检测到的水损害缺陷范围。1-2为过小对比度对应的雷达图谱,但1-4所对应区域的缺陷特征在1-2图谱上未被突出,不易被GPR专家或模型识别;1-3为过大对比度对应雷达图谱,除了1-4对应区域显示为水损害缺陷外,其余正路区域也为高亮区域,会被误判为缺陷区域。从图8分析可知,合适对比度对GPR图谱目标识别非常重要。
[0084] 以下给出本发明的具体实施例,需要说明的是本发明并不局限于以下具体实施例,凡在本申请技术方案基础上做的等同变换均落入本发明的保护范围。
[0085] 实施例1:
[0086] 本实施例给出一种沥青路面水损害检测方法,如图1至图3所示,按照以下步骤进行:
[0087] 步骤一,通过GPR路面调查,获取水损害图谱数据集:
[0088] 步骤S11,GPR路面调查、数据采集:采用GPR系统对沥青路面进行现场数据采集,现场数据采集过程中,确定路面出现唧浆或泛白的损害区域;
[0089] 步骤S11中,所述的现场数据采集过程中,采样参数要求采样间距<15cm,天线频率>1.6GHz,采样频率为天线主频的10~20倍。
[0090] 这些标记将在GPR图谱上方以小正方形出现,如图2中雷达图谱的标记“□”,其下方对应的是水损害缺陷区域,这些标记对应的GPR图谱作为水损害的真值,用于确定水损害缺陷的特征;
[0091] 步骤S12,获取损害的初始图谱数据集:将损害区域对应的GPR数据进行预处理后,设定GPR图谱的对比度并按照长度为5~6m截取GPR图谱,构建具有水损害、桥接缝和正常路面的损害的初始图谱数据集,并对各自特征进行标记;
[0092] 所述的损害的初始图谱数据集的图片分辨率为1090×300;
[0093] 步骤S12中,所述的预处理过程为采用直流漂移校正算法、地面校正算法、背景扣除算法、带通滤波算法和进行滑动平均算法进行预处理。
[0094] 步骤S12中,设定GPR图谱的对比度为1.2~1.6,本实施例优选1.4。
[0095] 图2是水损害缺陷数据集典型图谱,左侧是现场照片,右侧是对应的GPR图谱,标签中Moisture为水损害,Joint为桥接缝。
[0096] 获取水损害的图谱数据集的过程:GPR天线经过水损害区域时,会在数据采集软件中打一个标记,通过大量的实例调查,确定了水损害的主要特征:
[0097] 1)在沥青层中存在连续或者不连续高亮区域;
[0098] 2)图谱区域宽高比Width/Height不定,与水损害的严重程度呈正相关。
[0099] 图2中最下方图谱为桥接缝图谱,其特征是从路面往下呈现连续高亮区域,与水损害高亮区域主要区别在于:
[0100] 1)该特征从路表面往下呈现高亮,两侧会有双曲线特征;
[0101] 2)高亮区域特征连续且从表面到下方的深度Depth≥0.1m;
[0102] 3)图谱区域宽高比Width/Height<4且面积Area>1000pixel2。
[0103] 步骤二,图片分辨率调整:
[0104] 研究发现,不同分辨率图谱在识别模型中的准确度不同,图片分辨率直接影响模型识别效果;
[0105] 损害的初始图谱数据集定义为ID数据集,将ID数据集直接缩放到224×224,定义为BD数据集;
[0106] 将损害的初始图谱数据集的分辨率直接缩放到224×224,获得BD数据集;
[0107] 步骤三,将数据集输入到识别模型:
[0108] 将步骤二获得的BD数据集输入到识别模型,经过识别模型运算后,执行步骤四;
[0109] 所述的识别模型的图片输入分辨率大小为224×224,图片输出分辨率大小为224×224;
[0110] 所述的识别模型为混合式深度学习模型,所述的混合式深度学习模型由两部分构成,特征提取采用ResNet50,目标检测采用YOLO V2框架;
[0111] 所述的ResNet50和YOLO V2框架均为已知的深度学习模型。
[0112] 所述的特征提取由4阶段构成,实现16倍降采样,将输入的224×224转为14×14×1024,为后续YOLO检测提供CNN原始数据;
[0113] 所述的YOLO v2框架中,提供目标检测和候选框,YOLO Class Conv设置栅格Grid=14×14,Anchor boxes=6。YOLO Transform设置的loss function为mse。
[0114] 所述的混合式深度学习模型采用步骤二中获得的图谱分为训练集和测试集,分配比例为70%和30%,具体的模型训练方法为:利用TL方法对设计好的混合式深度学习模型进行训练。模型的loss function采用MSE方法,Anchor boxes数量是按照K-means方法对根据样本集水损害和桥接缝的Height/Width的分类获得。
[0115] 所述的混合式深度学习模型采用F1,Recall和Precision等3个指标来衡量模型的性能。
[0116] 步骤四,输出水损害结果:
[0117] 识别模型给出的结果存在相互重叠的现象,包括:
[0118] 1)一个较长的水损害缺陷,会出现多个预测结果且相互重叠;
[0119] 2)部分桥接缝,除判断为多个结果,还被误判为水损害;
[0120] 为此,图3是针对图谱带有坐标轴的后处理流程图,具体的:
[0121] 对步骤三的识别模型的输出结果进行后处理,所述的后处理步骤为:
[0122] 步骤S41,判断输出结果中的图谱的候选框BBoxes数量,如果候选框BBoxes数量>1,则进行步骤S42,否则直接输出结果;
[0123] 步骤S42:判断候选框BBoxes是否有重叠,若有重叠,则进行步骤S43,否则直接输出结果;
[0124] 步骤S43:判断重叠候选框对应的标签名是否相同,若相同,则合并的候选框对应的标签名不变;若不同,则表明同时含有水损害和桥接缝两类标签名,标签名输出为桥接缝Joint;
[0125] 步骤S44:合并候选框,将相交的候选框x和y方向取最小值,w和h取最大值,合并的候选框坐标为[xmin,ymin,wmax,hmax];
[0126] 步骤S45:输出结果,在识别模型的输出结果中,将输出图片分辨率调整至等于所述的损害的初始图谱数据集的图片分辨率,输出结果为带目标的标签名及目标对应的候选框位置BBoxes(x,y,w,h)的图谱。
[0127] (A)本发明突破了现有GPR领域自动检测中集中在双曲线特征目标的检测,实现了具有复杂目标体特征下沥青路面水损害缺陷的自动检测,可为沥青路面精准预养护和水损害缺陷的自动定位提供依据。
[0128] (B)由于本发明考虑了图片分辨率缩放的影响,并利用混合式模型进行水损害缺陷的自动检测,解决了现有水损害缺陷依靠专家经验识别存在费时费力、人为因素的影响。
[0129] (C)本发明的训练样本集均来自现场试验数据,样本具有广泛的代表性,解决了现有GPR领域研究用FDTD仿真软件合成数据集存在样本不具有代表性的问题,并突破了GPR领域自动识别仅集中在双曲线特征自动检测的局限。
[0130] (D)由于本发明提供的方法可实现水损害缺陷的自动检测和精准定位,因此可以为后期基于无人巡检车的自动检测提供识别模型,实现缺陷区域的定期检测和巡检,进一步实现精准养护和智能养护。
[0131] 对比例1:
[0132] 本对比例给出一种沥青路面水损害检测方法,该方法的其它步骤与实施例1相同,区别仅仅在于步骤二不同以及步骤三中输入的图谱不同。
[0133] 步骤二,图片分辨率调整:
[0134] 损害的初始图谱数据集定义为ID数据集,将ID数据集按照224×224尺寸进行裁剪,将裁剪后的图谱包含水损害和桥接缝的图谱构建的图谱定义为SD数据集;
[0135] 将损害的初始图谱数据集按照224×224尺寸进行裁剪,获得SD数据集。
[0136] 对比例2:
[0137] 本对比例给出一种用于水损害检测的沥青路面水损害检测方法,该方法的其它步骤与实施例1相同,区别仅仅在于步骤二不同以及步骤三中输入的图谱不同。
[0138] 步骤二,图片分辨率调整:
[0139] 损害的初始图谱数据集定义为ID数据集,将ID数据集按照224×224尺寸进行裁剪,将裁剪后的图谱包含水损害和桥接缝的图谱构建的图谱定义为SD数据集,将BD数据集和SD数据集混合构建的图谱定义为MD数据集;
[0140] 将损害的初始图谱数据集的分辨率调整,获得MD数据集。
[0141] 对实施例1、对比例1和对比例2进行对比分析,按照所提算法,构建了原始图谱数据集1431张图谱,按照算法分别构建BD,SD和MD数据集。图4是训练模型的结果,由图可知,所提出的混合式深度模型在各数据集训练的网络在测试集上均有较好的结果,表明所提的混合式深度模型可行。其中,BD数据集上训练的模型最优,模型识别精度为F1=91.97%,Recall=94.53%,以及Precision=91.00%。为此本发明选择训练模型为BD模型,优选分辨率缩放方法为直接将原始图谱进行等比例缩放。
[0142] 对比例3:
[0143] 本对比例给出一种沥青路面水损害检测方法,该方法采用ACF(Aggregate Channel Features)算法对沥青路面水损害进行检测。
[0144] 对比例4:
[0145] 本对比例给出一种沥青路面水损害检测方法,该方法采用Cifar10模型对沥青路面水损害进行检测。
[0146] 对实施例1、对比例3和对比例4进行对比分析,图5至图7是所提深度模型与ACF(Aggregate Channel Features)和Cifar10的对比结果,图中Ground Truth为水损害的真值。通过对比发现,对比的两种方法均存在冗余检测区域或大量漏检区域,对比结果进一步验证了所提方法的准确性。
[0147] 实施例2:
[0148] 本实施例给出一种用于水损害检测的探地雷达图谱自适应选取方法,如图9所示,该方法根据探地雷达图谱的对比度来自适应选取探地雷达图谱,该方法按照以下步骤进行:
[0149] 步骤1,读取预处理后的GPR数据:
[0150] 对GPR数据进行预处理后,在设定的对比度取值范围内,随机生成不同对比度的雷达图谱,构建初始随机图谱数据集,所述的初始随机图谱数据集的大小为N张图片,初始随机图谱数据集用于判断是否含有目标;
[0151] 所述的GPR数据的获取方法为:采用GPR系统对沥青路面进行现场数据采集,现场数据采集过程中,确定路面出现唧浆或泛白的损害区域,获取损害区域对应的GPR数据。
[0152] 所述的现场数据采集过程中,采样参数要求采样间距<15cm,天线频率>1.6GHz,采样频率为天线主频的10~20倍。
[0153] 所述的预处理过程为采用直流漂移校正算法、地面校正算法、背景扣除算法、带通滤波算法和进行滑动平均算法进行预处理。
[0154] 所述的设定的对比度取值范围为0.5~1.8。
[0155] 所述的N=100。
[0156] 步骤2,图片分辨率调整:
[0157] 研究发现,不同分辨率图谱在识别模型中的准确度不同,图片分辨率直接影响模型识别效果;
[0158] 初始随机图谱数据集定义为RID数据集,将RID数据集直接缩放到224×224,定义为RBD数据集;
[0159] 将水损害初始图谱数据集的分辨率直接缩放到224×224,获得RBD数据集;
[0160] 步骤3,将数据集输入到识别模型:
[0161] 将步骤二获得的RBD数据集输入到识别模型,经过识别模型运算后,执行步骤四;
[0162] 所述的识别模型的图片输入分辨率大小为224×224,图片输出分辨率大小为224×224;
[0163] 所述的识别模型为混合式深度学习模型,所述的混合式深度学习模型由两部分构成,特征提取采用ResNet50,目标检测采用YOLO V2框架;
[0164] 所述的ResNet50和YOLO V2框架均为已知的深度学习模型。
[0165] 所述的特征提取由4阶段构成,实现16倍降采样,将输入的224×224转为14×14×1024,为后续YOLO检测提供CNN原始数据;
[0166] 所述的YOLO v2框架中,提供目标检测和候选框,YOLO Class Conv设置栅格Grid=14×14,Anchor boxes=6。YOLO Transform设置的loss function为mse。
[0167] 所述的混合式深度学习模型采用步骤二中获得的图谱分为训练集和测试集,分配比例为70%和30%,具体的模型训练方法为:利用TL方法对设计好的混合式深度学习模型进行训练。模型的loss function采用MSE方法,Anchor boxes数量是按照K-means方法对根据样本集水损害和桥接缝的Height/Width的分类获得。
[0168] 所述的混合式深度学习模型采用F1,Recall和Precision等3个指标来衡量模型的性能。
[0169] 步骤4,输出水损害结果:
[0170] 步骤四,输出水损害结果:
[0171] 对步骤三的识别模型的输出结果进行后处理,所述的后处理步骤为:
[0172] 步骤S41,判断输出结果中的图谱的候选框BBoxes数量,如果候选框BBoxes数量>1,则进行步骤S42,否则直接输出结果;
[0173] 步骤S42:判断候选框BBoxes是否有重叠,若有重叠,则进行步骤S43,否则直接输出结果;
[0174] 步骤S43:判断重叠候选框对应的标签名是否相同,若相同,则合并的候选框对应的标签名不变;若不同,则表明同时含有水损害和桥接缝两类标签名,标签名输出为桥接缝Joint;
[0175] 步骤S44:合并候选框,将相交的候选框x和y方向取最小值,w和h取最大值,合并的候选框坐标为[xmin,ymin,wmax,hmax];
[0176] 步骤S45:输出结果,在识别模型的输出结果中,将输出图片分辨率调整至等于所述的损害的初始图谱数据集的图片分辨率,输出结果为带目标的标签名及目标对应的候选框位置BBoxes(x,y,w,h)的图谱;
[0177] 步骤5,利用初始随机图谱数据集,进行检测目标的有无判断:
[0178] 步骤S51,将步骤四的输出结果转为图片上像素点对应的矩阵Ai,Ai定义为:
[0179] where1≤m≤H0,1≤n≤W0
[0180] 式中:H0为识别模型输出图谱的图片高度;W0为识别模型输出图谱的图片宽度;
[0181] 将RID数据集中的N张图片对应的矩阵Ai进行求和,并求其均值,获得均值矩阵A,A定义为:
[0182] 步骤S52,由于对比度的范围是优选过的,目标与背景的差异比较大,若被测GPR数据中含有目标,那么大部分图谱的输出结果中应该包括了该目标区域,则均值矩阵A在该区域的值大;若被测GPR数据中没有目标,则只有少数不恰当的对比度对应图谱有目标,则均值矩阵A在该区域的值小。
[0183] 设定k1=0.8,θ0=0.5,对均值矩阵A进行按照以下公式进行更新,获得更新后的均值矩阵A;
[0184] A(A
[0185] 式中:
[0186] k1目标关联系数,用于调整均值最大值,以实现不同目标的判断;
[0187] θ0为含有目标时,矩阵A中的最小值,低于此值,则无目标;
[0188] max(max(A))为均值矩阵中的最大值;
[0189] 步骤S53,根据更新后的均值矩阵A,按照以下公式获得目标有无的判断条件T,当T=1,则表示有目标,当T=0,则表示没有目标;
[0190]
[0191] 步骤6,增量式随机生成图谱,进行合适对比度图谱的选取:
[0192] 当图片含有目标后,进行初始化判断:
[0193] 步骤S61,若Flag=0,则表示为第一次生成随机图谱样本集,即初始化样本集阶段,不进入后续选取判断,设置Flag=1,然后额外增加5%·N张图片作为样随机图谱数据集的本,样本中图片的总数为N=(1+5%)N,返回到步骤2;
[0194] 步骤S62,若Flag≠0,则表示非初始化阶段,设置一个图片相关系数,选择与平均矩阵A最大相关系数的图片,即为合适对比度的图谱;
[0195] 相关系数Ri定义为:
[0196] 式中:Ri为第i张图谱对应的矩阵Ai与均值矩阵A的相关系数;m为高度方向上的坐标值;n为宽度方向上的坐标值;μA为均值矩阵A的总均值;μAi为矩阵Ai的总均值;
[0197] 选取过程的终止条件为:
[0198] STOP=abs(F1-F1pre)<0.01&&F1>0.8;
[0199]
[0200]
[0201]
[0202] 式中:F1为深度学习的评价指标;F1Pre为上一次深度学习的评价指标,初始时为0;TP为真实的目标区域;FP为错误识别的真值,代表未识别的真值,错误判断为负值或背景;
FN为错误识别的负值,即把背景当成目标;
[0203] 当终止条件不满足时,将该次计算的指标F1付给变量F1Pre,然后返回步骤S61,增大样本集,重新进行选取;
[0204] 当终止条件满足时,系统将输出合适对比度的图谱。
[0205] 效果测试对比:
[0206] 采用人工构建的水损害图谱数据集,利用YOLO检测框架和transfer learning训练出深度学习模型,结合图9中的算法进行识别。GPR数据与图8中的数据一致。为了观测所提方法各指标随随机样本数的变化规律,图9中的初始样本集N=20,对比度取值范围0.6~1.8,更新矩阵A中的参数k1=0.8,θ0=0.5,实际应用时N=100.
[0207] 图10是相关系数随采样样本数的变化规律,图中Referenced F1为F1指标,其余两条曲线为优选图谱与均值A以及真值(Ground Truth)的相关系数,从图可知,随着样本数的增加,相关系数Ri均在增加,并且稳定在一个固定值,不随样本数增大而增加,表明所提算法已寻找到优选的GPR图谱。
[0208] 图11是深度模型相关指标随采样样本数的变化规律,F1,Precision和Recall均为深度模型评价指标,同样在找到优选结果后,各指标稳定,表明已找到合适的图谱。
[0209] 图12是160个随机样本集下优选图谱与随机数据集总体检测的heatmap,可知在103.5m对应的区域,160张图片均在此处有水损害(目标)结果,与真值(图8中1-4)区域一致,并且3大指标结果大于0.93,表明优选的图谱与真值很接近,验证算法的正确性。
[0210] 图13是优选结果与均值A的heatmap图(即将随机图谱集中所有检测结果进行累积,并通过颜色体现累加结果大小),由于均值A在更新过程中将A矩阵中最大值×0.8以下的结果全部取出,且是平均值,通过对比检测结果与Ground Truth可知,二者吻合度很高,进一步验证所提算法的正确性。
[0211] 图14是正常路面与随机数据集总体检测的heatmap图,由于初始化样本数大小为100,而此图中最大值只有28(均值计算后为0.28),小于θ0=0.5的阈值,该测试数据被分类为正常路面,与实际情况一致。
[0212] 图15是测试数据集的随机样本数分布规律,通过所提算法对31个样本(11个正常路面,20个为水损害样本)进行测试,结果表明所提算法可有效将11个正常路面分类,其所用样本数为100;而水损害区域所需要的随机样本主要集中在200-300之间。由于路面调查的GPR数据中95%为正常路面,因此该算法可有效节约计算开支。
[0213] 为进一步说明所提方法(增量式采样法,记为IRS)的有效性,对比随机选取法(RS)的结果,图16是针对有水损害的数据集(20个)的对比结果,图17增加了11个正常路面的结果。从结果可知所提IRS方法能有效识别缺陷,且精度较高。
[0214] 通过实验表明,所提的增量式采样方法与深度模型联合使用,可有效从GPR原始数据中选取合适对比度的雷达图谱,为GPR的自动应用提供一种有效的方法。
[0215] 所提方法虽然是在水损害缺陷识别进行了验证,但该方法不限定于此案例,利用该方法进行其他雷达图谱中国目标体识别也在本发明保护范围内。
[0216] 实施例3:
[0217] 本实施例给出一种基于图谱灰度自适应选取的沥青路面水损害检测方法,如图1至图17所示,该方法与实施例1的沥青路面水损害检测方法基本相同,区别仅仅在于步骤S12中,将“设定GPR图谱的对比度并按照长度为5~6m截取GPR图谱”替换为“选取合适对比度的GPR图谱并按照长度为5~6m截取GPR图谱”。
[0218] 所述的合适对比度的GPR图谱的选取采用的方法为探地雷达图谱自适应选取方法;
[0219] 所述的探地雷达图谱自适应选取方法采用实施例2中所述的用于水损害检测的沥青路面水损害检测方法相同。
[0220] 实施例1和实施例2中所述的识别模型与相同;实施例1和实施例2中所述的后处理步骤相同。
[0221] 本实施例的方法可以为每一处GPR数据实现图谱优选,然后将优选图谱输入给深度模型得到检测结果。所提方法解决了水损害的图谱优选和图谱识别,真正实现了水损害缺陷检测的自动化和智能化工作。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈