CN115690001A - 一种用于检测钢管焊接数字射线图像中缺陷的方法 - Google Patents

一种用于检测钢管焊接数字射线图像中缺陷的方法 Download PDF

Info

Publication number
CN115690001A
CN115690001A CN202211228352.3A CN202211228352A CN115690001A CN 115690001 A CN115690001 A CN 115690001A CN 202211228352 A CN202211228352 A CN 202211228352A CN 115690001 A CN115690001 A CN 115690001A
Authority
CN
China
Prior art keywords
image
images
steps
feature
training
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.)
Pending
Application number
CN202211228352.3A
Other languages
English (en)
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.)
Sinopec Oilfield Equipment Corp
Sinopec Petroleum Engineering Machinery Co Ltd Shashi Steel Pipe Works Branch
Original Assignee
Sinopec Oilfield Equipment Corp
Sinopec Petroleum Engineering Machinery Co Ltd Shashi Steel Pipe Works Branch
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 Sinopec Oilfield Equipment Corp, Sinopec Petroleum Engineering Machinery Co Ltd Shashi Steel Pipe Works Branch filed Critical Sinopec Oilfield Equipment Corp
Priority to CN202211228352.3A priority Critical patent/CN115690001A/zh
Publication of CN115690001A publication Critical patent/CN115690001A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/30Computing systems specially adapted for manufacturing

Landscapes

  • Image Analysis (AREA)

Abstract

本发明提供一种用于检测钢管焊接数字射线图像中缺陷的方法,包括以下步骤:S1、采集多个DR图像,扩大样本数据集,对多个DR图像进行预处理,并构造训练数据集合;S2、构建基于图像重建的深度卷积网络模型,该网络模型至少包含一个基于卷积网络的编码器和一个对应的解码器;S3、进行模型参数优化训练,利用训练数据集以优化网络参数,并得到优化后的网络模型;S4、将采集的DR图像输入优化后的网络模型进行推理,得到重建的图像;将重建图像与输入图像进行差分运算,并将差分结果进行二值化以获得最终缺陷位置;通过采用以上的方案,获得最终缺陷位置。

Description

一种用于检测钢管焊接数字射线图像中缺陷的方法
技术领域
本发明涉及钢管焊接缺陷检测领域,特别是一种用于检测钢管焊接数字射线图像中缺陷的方法。
背景技术
在钢管焊接过程中,受生产环境和制造工艺的影响,不可避免地出现各种焊接缺陷,如气孔、裂纹、不熔合和溅射等。这些缺陷在很大程度上影响了产品的性能。检 测焊接缺陷并及时修复是非常重要的。目前,行业内常用的焊接检测方法有人工检测 和非破坏性检测,包括射线检测、涡流检测、渗透检测、超声波检测和磁粉检测。虽 然无损检测可以解决大部分重要的表面缺陷,但它需要专门的技术人员进行现场检测。 对于企业来说,在检测效率和经济成本上很难达标。在实际焊接中,应尽量避免剩余 材料的产生。对焊接区域的检测和对缺陷位置的再加工,有利于提高钢管的使用质量 和利用率。因此,提高焊接缺陷检测的质量和效率就显得尤为重要。
目前,钢管焊缝通常采用无损检测方法,其中X射线是最常用的方法。X射线焊 缝检测根据图像的不同分为射线胶片检测和数字检测。前者主观性强,容易产生误判, 而后者利用计算机算法进行图像识别和检测。为了减轻人工分析的困难,目前已经开 发了一些自动化或半自动化技术。特别地,依靠传统计算机视觉或浅层机器学习的自 动分析已经在该领域得到了一些应用。遗憾的是,由于存在各种挑战,如焊缝的几何 形状,较低的对比度和不一致的亮度,以及外观差异较大的缺陷,此类方法在实际焊 接应用中并没有被广泛采用。最近,人们对机器学习和深度学习的兴趣激增,特别是 在图像分析方面。深度学习模型是通过从包含数千或数百万例子的大量数据集中学习 而开发的。对于许多无损检测应用来说,有注释的训练数据资源是有限的,因此,深 度学习的潜力可能没有得到充分的利用。当前该领域常采用目标检测或者语义分割模 型以发现和定位缺陷的位置,但此类方法均属于监督学习方法,自然训练实例的不足 导致模型无法有效地的被优化,在实际中性能较差。
综上所述,大多数相关企业仍然采用传统的人工观测的方法,这种方法主观因素多,检测效率低,容易造成焊接质量问题。而自动方法因为种种原因,例如,特征表 能力不足或人工标记样本数量不足,难以取得满意的性能。因此,有必要提供一种用 于检测钢管焊接数字射线图像中缺陷的方法,以解决上述问题。
中国专利号CN111738991A涉及一种焊缝缺陷的数字射线检测模型的创建方法。其工作原理是:运用数字图像处理技术提取出图像中的焊缝;通过运用图像增强技术 对焊缝图像进行样本数量的扩充;利用搭建的深度学习模型和设计的损失函数来训练 模型,能够从现有的缺陷图像中学习缺陷规律,使用深度神经网络能够更加灵活的对 所有待测图像进行批处理,更灵活且性能更加稳定,保证模型较高的召回率和一定的 准确率;通过对模型的封装保存方便对模型的调用和使用。减少人为工作量,提高工 作效率。但该专利在使用时,特征表能力不足或人工标记样本数量不足,难以取得满 意的性能。使用存在缺陷。
发明内容
本发明所要解决的技术问题是提供一种用于检测钢管焊接数字射线图像中缺陷的方法,以解决人工观测的方法主观因素多,检测效率低,容易造成焊接质量问题; 而自动方法特征表能力不足或人工标记样本数量不足,难以取得满意的性能的问题。
为解决上述技术问题,本发明所采用的技术方案是:一种用于检测钢管焊接数字射线图像中缺陷的方法,包括以下步骤:
S1、采集多个DR图像,扩大样本数据集,对多个DR图像进行预处理,并构造训 练数据集合;
S2、构建基于图像重建的深度卷积网络模型,该网络模型至少包含一个基于卷积网络的编码器和一个对应的解码器;
S3、进行模型参数优化训练,利用训练数据集以优化网络参数,并得到优化后的网络模型;
S4、将采集的DR图像输入优化后的网络模型进行推理,得到重建的图像;将重 建图像与输入图像进行差分运算,并将差分结果进行二值化以获得最终缺陷位置;
通过以上方法,获得最终缺陷位置。
优选的方案中,在步骤S1中,用于构造训练数据集合的DR图像中应不包含缺陷,多个DR图像进行图像增强处理;
对多个DR图像进行预处理的步骤为:
A1:将多通道图像转换为单通道灰度图,转换函数为:
f0(r,c)=0.299×R(r,c)+0.587×G(r,c)+0.114×B(r,c)
其中,f0(r,c)为转换后的灰度图像,R(r,c),G(r,c)和B(r,c)分别为输入图像 的红色、绿色和蓝色通道图像。r和c表示图像的第r行和第c列;
A2:采用直方图均衡化以强调图像的对比度;
A3:采用中值滤波以过滤图像中的脉冲型噪声;
A4:采用其他非线性变化函数对图像对比度进行调整,函数为:
Figure BDA0003880924340000031
其中,f0(x,y)为输入单通道图 像,f(x,y)为输出图像,a和b分别为f0(x,y)的最小灰度值和最大灰度值。
优选的方案中,图像增强处理包括直接和间接扩大样本数量,直接扩大样本数量的方法包括额外进行现场图像采集,生成对抗模型、变分自动编码器模型深度模型生 成额外的图像,间接扩大样本数量的方法包括图像进行平移、尺度缩放、旋转、剪切、 镜像、灰度抖动、模板消除处理。
优选的方案中,卷积网络模型设计网络模块为:ReLU”块为激活层;“MaxPool” 块为最大池化层;“Resize”块为2倍分辨率上采样层;“Concat”块为特征拼接层, “Cluster_N”和“Memory_N”分别为特征聚类层和特征存储层;
编码器用于将输入的图像编码为多尺度的描述特征,解码器用于将描述特征翻译为二维图像;
优选的方案中,步骤S2中基于图像重建的深度卷积网络模型的步骤是:
B1:特征聚类层仅工作在训练过程中,通过损失设计以增强特征关于图像中主要部分的描述;
B2:采用了多个特征存储模块,用于存储训练时选取的有代表性的特征,在测试阶段中利用存储的特征来重建输入图像。
优选的方案中,B1中增强特征关于图像中主要部分的描述的步骤是:
C1:通过变形模块将其变为(H×W)×C,通过聚类算法将H×W个维度为C的一维 向量聚类成K类,得到K个类心,每个类心向量维度为C维;
其中C,H和W分别为特征的通道数;
C2:通过变形模块将其变为(H×W)×C,然后按照下式来计算特征图上每个位置f_i与聚类中心向量μ_j之间的相似度S_ij:
Figure BDA0003880924340000041
其中i∈(1,2,…,W*H),j∈(1,2,…,K),j^'∈(1,2,…,K)。α控制相似度分布 的形状,在本申请的方案中,α=1是一个常规的选择;
C3:相似度归一化:
Figure BDA0003880924340000042
C4:计算损失L_cluster以增强特征对图像中主要部分的描述能力:
Figure BDA0003880924340000043
其中KL为Kullback-Leibler散度。
优选的方案中,B2中用存储的特征来重建输入图像的步骤是:
D1:计算特征图上每个位置q_i与每个存储单元m_j之间的距离:
Figure BDA0003880924340000044
其中i∈(1,2,…,W*H),j∈(1,2,…,N),j^'∈(1,2,…,N);
D2:求出输入向量q_i的重建结果q^^i:
Figure BDA0003880924340000045
D3:对损失函数进行约束:
Figure BDA0003880924340000046
D4:最初输入图像计算重建损失:
Figure BDA0003880924340000051
r和c表示图像的第r行和第c列;
D5:整个系统的损失函数如下式:
Figure BDA0003880924340000052
优选的方案中,S3中模型参数优化训练的步骤是:
E1:计算重建图像和输入图像之间的差值图像:
Figure BDA0003880924340000053
E2:二值化处理:
Figure BDA0003880924340000054
r和c表示图像的第r行和第c列。β为一个常数,β=0.7是一个常规选择。TH_val 也是一个常数;
计算平均灰度值作为TH_val:
Figure BDA0003880924340000055
S为集合中图像张数,H和W分别为图像的高度和宽度。dio(s,r,c)表示第s张 差值图像中的位于第r行、第c列的像素灰度值。
优选的方案中,缺陷位置确定方法为:通过二值化结果图像中非零像素位置将被判断为缺陷所在位置,如果没有非零像素,则该图像中没有缺陷。
本发明提供了一种用于检测钢管焊接数字射线图像中缺陷的方法,通过采用以上的方案,通过定制的预处理算法,提高了DR图像清晰度,增加了缺陷与背景的区分 度;构建新颖的图像重建深度卷积网络模型,提高网络的特征提取能力和图像重建 的精确性,稳定性,使重建图像与输入图像之间的差分结果更稳健,实现了缺陷检测 的稳定性和精确性的统一。本发明能够实时检测钢管焊缝中的裂纹﹑未焊透﹑未熔合 ﹑咬边﹑条形﹑圆形﹑飞溅等各类缺陷,优选的方案中,在保证缺陷误检率的基础上, 获得了较高的焊缝缺陷检测率。在提高检测速度的同时,极大的降低了人工缺陷检测 的工作量。
附图说明
下面结合附图和实施例对本发明作进一步说明:
图1是一张没有缺陷的钢管焊缝DR图像;
图2是一种符合本申请要求的网络模型的实施例;
图3是聚类模块中一种生成聚类中心的示意图;
图4是聚类模块中一种计算聚类损失的示意图;
图5是一种训练流程的示意图;
图6是图5所示训练流程中,一种验证子流程的示意图;
图7是图6中一种计算二值化分割阈值的示意图;
图8是当训练完成后,一种实际工作中缺陷检测的流程示意图。
具体实施方式
实施例1:
一种用于检测钢管焊接数字射线图像中缺陷的方法,包括以下步骤:
S1、采集多个DR图像,扩大样本数据集,对多个DR图像进行预处理,并构造训 练数据集合;
S2、构建基于图像重建的深度卷积网络模型,该网络模型至少包含一个基于卷积网络的编码器和一个对应的解码器;
S3、进行模型参数优化训练,利用训练数据集以优化网络参数,并得到优化后的网络模型;
S4、将采集的DR图像输入优化后的网络模型进行推理,得到重建的图像;将重 建图像与输入图像进行差分运算,并将差分结果进行二值化以获得最终缺陷位置;
通过以上方法,获得最终缺陷位置。
优选的方案中,在步骤S1中,用于构造训练数据集合的DR图像中应不包含缺陷,多个DR图像进行图像增强处理;
对多个DR图像进行预处理的步骤为:
A1:将多通道图像转换为单通道灰度图,转换函数为:
f0(r,c)=0.299×R(r,c)+0.587×G(r,c)+0.114×B(r,c)
其中,f0(r,c)为转换后的灰度图像,R(r,c),G(r,c)和B(r,c)分别为输入图像 的红色、绿色和蓝色通道图像。r和c表示图像的第r行和第c列;
A2:采用直方图均衡化以强调图像的对比度;
A3:采用中值滤波以过滤图像中的脉冲型噪声;
A4:采用其他非线性变化函数对图像对比度进行调整,函数为:
Figure BDA0003880924340000071
其中,f0(x,y)为输入单通道图 像,f(x,y)为输出图像,a和b分别为f0(x,y)的最小灰度值和最大灰度值。
优选的方案中,图像增强处理包括直接和间接扩大样本数量,直接扩大样本数量的方法包括额外进行现场图像采集,生成对抗模型、变分自动编码器模型深度模型生 成额外的图像,间接扩大样本数量的方法包括图像进行平移、尺度缩放、旋转、剪切、 镜像、灰度抖动、模板消除处理。
优选的方案中,卷积网络模型设计网络模块为:ReLU”块为激活层;“MaxPool” 块为最大池化层;“Resize”块为2倍分辨率上采样层;“Concat”块为特征拼接层, “Cluster_N”和“Memory_N”分别为特征聚类层和特征存储层;
编码器用于将输入的图像编码为多尺度的描述特征,解码器用于将描述特征翻译为二维图像;
优选的方案中,步骤S2中基于图像重建的深度卷积网络模型的步骤是:
B1:特征聚类层仅工作在训练过程中,通过损失设计以增强特征关于图像中主要部分的描述;
B2:采用了多个特征存储模块,用于存储训练时选取的有代表性的特征,在测试阶段中利用存储的特征来重建输入图像。
优选的方案中,B1中增强特征关于图像中主要部分的描述的步骤是:
C1:通过变形模块将其变为(H×W)×C,通过聚类算法将H×W个维度为C的一维 向量聚类成K类,得到K个类心,每个类心向量维度为C维;
其中C,H和W分别为特征的通道数;
C2:通过变形模块将其变为(H×W)×C,然后按照下式来计算特征图上每个位置f_i与聚类中心向量μ_j之间的相似度S_ij:
Figure BDA0003880924340000081
其中i∈(1,2,…,W*H),j∈(1,2,…,K),j^'∈(1,2,…,K)。α控制相似度分布 的形状,在本申请的方案中,α=1是一个常规的选择;
C3:相似度归一化:
Figure BDA0003880924340000082
C4:计算损失L_cluster以增强特征对图像中主要部分的描述能力:
Figure BDA0003880924340000083
其中KL为Kullback-Leibler散度。
优选的方案中,B2中用存储的特征来重建输入图像的步骤是:
D1:计算特征图上每个位置q_i与每个存储单元m_j之间的距离:
Figure BDA0003880924340000084
其中i∈(1,2,…,W*H),j∈(1,2,…,N),j^'∈(1,2,…,N);
D2:求出输入向量q_i的重建结果q^^i:
Figure BDA0003880924340000085
D3:对损失函数进行约束:
Figure BDA0003880924340000086
D4:最初输入图像计算重建损失:
Figure BDA0003880924340000091
r和c表示图像的第r行和第c列;
D5:整个系统的损失函数如下式:
Figure BDA0003880924340000092
优选的方案中,S3中模型参数优化训练的步骤是:
E1:计算重建图像和输入图像之间的差值图像:
Figure BDA0003880924340000093
E2:二值化处理:
Figure BDA0003880924340000094
r和c表示图像的第r行和第c列。β为一个常数,β=0.7是一个常规选择。TH_val 也是一个常数;
计算平均灰度值作为TH_val:
Figure BDA0003880924340000095
S为集合中图像张数,H和W分别为图像的高度和宽度。dio(s,r,c)表示第s张 差值图像中的位于第r行、第c列的像素灰度值。
优选的方案中,缺陷位置确定方法为:通过二值化结果图像中非零像素位置将被判断为缺陷所在位置,如果没有非零像素,则该图像中没有缺陷。
实施例2:
结合图1~8进一步说明:
针对钢管焊接的DR图像进行处理。图X给出了本申请实施例提供的钢管焊接缺 陷检测方法的一种示意性流程图,可以包括以下几个步骤:
构造训练数据集合。构成训练的数据集合是现场采集的DR图像,本申请的特征 在于,图像中无需缺陷存在。在现代化生产过程中,以当前的焊接工艺水平,无缺陷 图像的收集难度要远低于有缺陷的图像。因此,在实际应用中可以较为方便地收集到 足够数量的用于训练的图像。注意,足够数量并未指定固定的张数。考虑到模型参数 的训练过程可以在长时间的调试过程中不断迭代,这里可以是多张。图1给出了一张 训练图像的示例,注意,该图像中焊缝的方向,宽度,灰度以及对比度等不应作为唯 一标准。针对采集的图像,需要对图像进行预处理。作为一种实现方式,预处理步骤 可以为:
如果图像为多通道图像,转换为单通道灰度图。以常见3通道RGB图像转换为单 通道灰度图为例,一种转换函数为:
f0(r,c)=0.299×R(r,c)+0.587×G(r,c)+0.114×B(r,c) (1)
其中,f0(r,c)为转换后的灰度图像,R(r,c),G(r,c)和B(r,c)分别为输入图像的红色、 绿色和蓝色通道图像。r和c表示图像的第r行和第c列。注意,该转换函数符合国际电信联盟推荐的BT.601-7建议书,但不应作为唯一转换标准。
采用直方图均衡化以强调图像的对比度;
采用中值滤波以过滤图像中的脉冲型噪声,注意,中值滤波模板的尺寸,形状可以根据不同设备采集的图像来针对设计,不应作为唯一限制。
采用其他非线性变化函数对图像对比度进行调整,例如以下函数:
Figure BDA0003880924340000101
其中,f0(x,y)为输入单通道图像,f(x,y)为输出图像,a和b分别为f0(x,y)的最小灰度值和最大灰度值。
由于在训练过程中,样本数据集中的样本数量越大,其参数泛化能力更高,因此,还需要对图像进行数据增强处理。其中,本申请所述的数据增强,指直接或间接扩大 样本数据集中的样本数量。其中直接扩大样本数量的方法包括额外进行现场图像采集, 或通过生成对抗模型、变分自动编码器模型等深度模型生成额外的图像。间接扩大样 本数量的方法包括但不限于将批注后的图像进行平移、尺度缩放、旋转、剪切、镜像、 灰度抖动、模板消除等处理。
构建基于图像重建的深度卷积网络模型。图2给出了一种可选的实现方式。图2中,Cout×Cin×H×W块为卷积层,其中Cout,Cin,H,W分别为卷积核的输出通道数 (卷积核个数)、输入通道数、高度和宽度;“ReLU”块为激活层;“MaxPool”块为 最大池化层;“Resize”块为2倍分辨率上采样层;“Concat”块为特征拼接层。上 述这些层均为卷积网络模型设计中常见的网络模块。“Cluster_N”和“Memory_N”分 别为特征聚类层和特征存储层,为本申请所示例模型中特有的模块。整个网络由编码 器和解码器两部分组成。编码器用于将输入的图像编码为多尺度的描述特征,解码器 用于将描述特征翻译为二维图像。注意的是,图2中采用了3个尺度的编码器和解码 器,但该层数不应作为唯一配置。
图2中采用了三个特征聚类模块:“Cluster_1”、“Cluster_2”和“Cluster_3”。 聚类模块仅工作在训练过程中,通过损失设计以增强特征关于图像中主要部分的描述。 其工作模式分为两步:
在第一阶段训练中,其工作模式如图3所示。设编码器侧面输出的特征维度为 C×H×W,其中C,H和W分别为特征的通道数,高度和宽度。首先通过变形模块将其 变为(H×W)×C,然后通过聚类算法将H×W个维度为C的一维向量聚类成K类,即得到 K个类心,每个类心向量维度为C维。
在第二阶段训练中,其工作模式如图4所示。设编码器侧面输出的特征维度为 C×H×W,其中C,H和W分别为特征的通道数,高度和宽度。首先通过变形模块将其 变为(H×W)×C,然后按照下式来计算特征图上每个位置fi与聚类中心向量μj之间的相 似度Sij
Figure BDA0003880924340000111
其中i∈(1,2,…,W*H),j∈(1,2,…,K),j′∈(1,2,…,K)。α控制相似度分布的形状,在本申 请的方案中,α=1是一个常规的选择。
然后对于相似度分布进行如下所示的归一化:
Figure BDA0003880924340000112
最后通过计算损失Lcluster以增强特征对图像中主要部分的描述能力。其定义为:
Figure BDA0003880924340000121
其中KL为Kullback-Leibler散度。
图2中采用了三个特征存储模块:“Memory_1”、“Memory_2”和“Memory_3”。 特征存储单元用于存储训练时选取的有代表性的特征。然后在测试阶段中利用存储的 特征来重建输入图像。考虑到本申请方案中在训练阶段只采用了没有缺陷的样本,因 此在测试阶段中若输入有缺陷的图像,则可以认为是将无缺陷的样本组合起来以重建 缺陷图像。在这种情况下,能够有效地的去除图像中的缺陷。具体来说,设M= {m1,m2,…,mN}∈RN×C是一个实值的矩阵,其中N表示存储向量的个数,C表示存储向量 的维度。类似聚类模块,设存储模块的输入特征的维度为C×H×W,其中C,H和W分 别为特征的通道数,高度和宽度。首先通过变形模块将其变为(H×W)×C,然后按照 下式计算特征图上每个位置qi与每个存储单元mj之间的距离:
Figure BDA0003880924340000122
其中i∈(1,2,…,W*H),j∈(1,2,…,N),j′∈(1,2,…,N)。得到所有wi后,输入向量qi的重 建结果
Figure BDA0003880924340000123
可以由下式确定:
Figure BDA0003880924340000124
考虑到存储单元应保留输入特征向量的代表性部分,其数值应该是较为稀疏的,因此,每个存储单元引入下列损失函数进行约束:
Figure BDA0003880924340000125
最后,重建的结果图像
Figure BDA0003880924340000126
需要和最初输入图像x之间计算如下的重建损失:
Figure BDA0003880924340000127
其中,r和c表示图像的第r行和第c列。综上所述,整个系统的损失函数如下式所示:
Figure BDA0003880924340000131
其中,a0+a1+a2+a3+b1+b2+b3=1。在本申请的方案中,a0=0.8,
Figure BDA0003880924340000132
Figure BDA0003880924340000133
b1=b2=b3=0.05是一个常规的选择。
进行模型参数优化训练。一种模型训练的流程如图5所示。该训练流程中需要构建验证样本集合。本申请中构建验证样本集合的过程和训练样本集合保持一致。其中, 一种验证流程如图6所示。下面先对图5中的部分步骤进行详细说明。网络前向推理 后得到的特征参与计算上述损失
Figure BDA0003880924340000134
并通过该损失计算S02,在S02中,一种优化 策略是采用随机梯度下降法,并设置动量值为0.9,权重衰减系数为0.001。S03中 判断是否进入验证阶段的标准可以设置为网络前向传播-反向优化固定次数后进行, 比如500次迭代后开始执行验证步骤S04。当S04的输出满足一定条件时,S05输出 训练过程结束信号。下面对S04对应的一种实施例,即图6中的部分步骤进行详细说 明。S06获取网络前向推理结果,即针对输入图像x的一张重建图像
Figure BDA0003880924340000138
后,首先根据下 式计算重建图像和输入图像之间的差值图像:
Figure BDA0003880924340000135
可以看到,dio
Figure BDA0003880924340000136
具有类似的计算形式。然后对dio进行二值化得到二值图像dbio,一种二值化的方法如下式所示:
Figure BDA0003880924340000137
其中,r和c表示图像的第r行和第c列。β为一个常数,β=0.7是一个常规选择。THval也是一个常数,该常数的计算方法如图7所示。首先,选取验证样本中所有带有 缺陷的图像作为一个图像集合。逐一输入集合中图像至S10中并获取该图像的重建图 像。S11用于计算两张图像之间的差值图像,并保存在计算机中。S12中进行对输入 图像的计数,若达到固定数值则可开始计算所有保存图像的平均灰度值。该固定数值 选择集合中图像总数是一个常规选择。S13按照下式计算所有保存的差值图像的平均 灰度值,并将该值作为THval
Figure BDA0003880924340000141
其中,S为集合中图像张数,H和W分别为图像的高度和宽度。dio(s,r,c)表示第s 张差值图像中的位于第r行、第c列的像素灰度值。
实际工作中,将采集的DR图像输入优化后的网络模型进行推理,得到重建的图像;将得到的重建图像于输入图像进行差分运算,并将差分结果进行二值化以获得最 终缺陷位置。图2给出了一种实现该流程的示意图。实际工作时,前端设备会实时采 集焊缝的DR图像并输入到本申请设计的处理系统。首先,利用训练好的网络模型对 单张图像进行推理,并得到重建图像。计算重建图像和输入图像之间的差值图像,计 算方法如公式(11)所示。对差值图像进行二值化,计算方法如公式(12)所示。最后, 因为此时已经得到二值化结果图像,因此图像中非零像素位置将被判断为缺陷所在位 置。如果没有非零像素,则该图像中没有缺陷。
上述的实施例仅为本发明的优选技术方案,而不应视为对于本发明的限制,本申请中的实施例及实施例中的特征在不冲突的情况下,可以相互任意组合。本发明的保 护范围应以权利要求记载的技术方案,包括权利要求记载的技术方案中技术特征的等 同替换方案为保护范围。即在此范围内的等同替换改进,也在本发明的保护范围之内。

Claims (9)

1.一种用于检测钢管焊接数字射线图像中缺陷的方法,其特征是包括以下步骤:
S1、采集多个DR图像,扩大样本数据集,对多个DR图像进行预处理,并构造训练数据集合;
S2、构建基于图像重建的深度卷积网络模型,该网络模型至少包含一个基于卷积网络的编码器和一个对应的解码器;
S3、进行模型参数优化训练,利用训练数据集以优化网络参数,并得到优化后的网络模型;
S4、将采集的DR图像输入优化后的网络模型进行推理,得到重建的图像;将重建图像与输入图像进行差分运算,并将差分结果进行二值化以获得最终缺陷位置;
通过以上方法,获得最终缺陷位置。
2.根据权利要求1所述的一种用于检测钢管焊接数字射线图像中缺陷的方法,其特征是:在步骤S1中,用于构造训练数据集合的DR图像中应不包含缺陷,多个DR图像进行图像增强处理;
对多个DR图像进行预处理的步骤为:
A1:将多通道图像转换为单通道灰度图,转换函数为:
f0(r,c)=0.299×R(r,c)+0.587×G(r,c)+0.114×B(r,c)
其中,f_0(r,c)为转换后的灰度图像,R(r,c),G(r,c)和B(r,c)分别为输入图像的红色、绿色和蓝色通道图像。r和c表示图像的第r行和第c列;
A2:采用直方图均衡化以强调图像的对比度;
A3:采用中值滤波以过滤图像中的脉冲型噪声;
A4:采用其他非线性变化函数对图像对比度进行调整,函数为:
Figure DEST_PATH_BDA0003880924340000031
其中,f_0(x,y)为输入单通道图像,f(x,y)为输出图像,a和b分别为f_0(x,y)的最小灰度值和最大灰度值。
3.根据权利要求2所述的一种用于检测钢管焊接数字射线图像中缺陷的方法,其特征是:图像增强处理包括直接和间接扩大样本数量,直接扩大样本数量的方法包括额外进行现场图像采集,生成对抗模型、变分自动编码器模型深度模型生成额外的图像,间接扩大样本数量的方法包括图像进行平移、尺度缩放、旋转、剪切、镜像、灰度抖动、模板消除处理。
4.根据权利要求1所述的一种用于检测钢管焊接数字射线图像中缺陷的方法,其特征是:
卷积网络模型设计网络模块为:ReLU”块为激活层;“MaxPool”块为最大池化层;“Resize”块为2倍分辨率上采样层;“Concat”块为特征拼接层,“Cluster_N”和“Memory_N”分别为特征聚类层和特征存储层;
编码器用于将输入的图像编码为多尺度的描述特征,解码器用于将描述特征翻译为二维图像。
5.根据权利要求1所述的一种用于检测钢管焊接数字射线图像中缺陷的方法,其特征是:步骤S2中基于图像重建的深度卷积网络模型的步骤是:
B1:特征聚类层仅工作在训练过程中,通过损失设计以增强特征关于图像中主要部分的描述;
B2:采用了多个特征存储模块,用于存储训练时选取的有代表性的特征,在测试阶段中利用存储的特征来重建输入图像。
6.根据权利要求5所述的一种用于检测钢管焊接数字射线图像中缺陷的方法,其特征是:B1中增强特征关于图像中主要部分的描述的步骤是:
C1:通过变形模块将其变为(H×W)×C,通过聚类算法将H×W个维度为C的一维向量聚类成K类,得到K个类心,每个类心向量维度为C维;
其中C,H和W分别为特征的通道数;
C2:通过变形模块将其变为(H×W)×C,然后按照下式来计算特征图上每个位置f_i与聚类中心向量μ_j之间的相似度S_ij:
Figure FDA0003880924330000021
其中i∈(1,2,…,W*H),j∈(1,2,…,K),j^'∈(1,2,…,K)。α控制相似度分布的形状,在本申请的方案中,α=1是一个常规的选择;
C3:相似度归一化:
Figure FDA0003880924330000031
C4:计算损失L_cluster以增强特征对图像中主要部分的描述能力:
Figure FDA0003880924330000032
其中KL为Kullback-Leibler散度。
7.根据权利要求5所述的一种用于检测钢管焊接数字射线图像中缺陷的方法,其特征是:B2中用存储的特征来重建输入图像的步骤是:
D1:计算特征图上每个位置q_i与每个存储单元m_j之间的距离:
Figure FDA0003880924330000033
其中i∈(1,2,…,W*H),j∈(1,2,…,N),j^’∈(1,2,…,N);
D2:求出输入向量q_i的重建结果q^^i:
Figure FDA0003880924330000034
D3:对损失函数进行约束:
Figure FDA0003880924330000035
D4:最初输入图像计算重建损失:
Figure FDA0003880924330000036
r和c表示图像的第r行和第c列;
D5:整个系统的损失函数如下式:
ltotal=α0×lrec1×lmem1+b1×lcluster1+a2×lmem2+b2×lcluster2+a3×lmem3+b3×lcluster3
8.根据权利要求1所述的一种用于检测钢管焊接数字射线图像中缺陷的方法,其特征是:S3中模型参数优化训练的步骤是:
E1:计算重建图像和输入图像之间的差值图像:
Figure FDA0003880924330000041
E2:二值化处理:
Figure FDA0003880924330000042
r和c表示图像的第r行和第c列。β为一个常数,β=0.7是一个常规选择。TH_val也是一个常数;
计算平均灰度值作为TH_val:
Figure FDA0003880924330000043
S为集合中图像张数,H和W分别为图像的高度和宽度。d_io(s,r,c)表示第s张差值图像中的位于第r行、第c列的像素灰度值。
9.根据权利要求1所述的一种用于检测钢管焊接数字射线图像中缺陷的方法,其特征是:缺陷位置确定方法为:通过二值化结果图像中非零像素位置将被判断为缺陷所在位置,如果没有非零像素,则该图像中没有缺陷。
CN202211228352.3A 2022-10-09 2022-10-09 一种用于检测钢管焊接数字射线图像中缺陷的方法 Pending CN115690001A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211228352.3A CN115690001A (zh) 2022-10-09 2022-10-09 一种用于检测钢管焊接数字射线图像中缺陷的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211228352.3A CN115690001A (zh) 2022-10-09 2022-10-09 一种用于检测钢管焊接数字射线图像中缺陷的方法

Publications (1)

Publication Number Publication Date
CN115690001A true CN115690001A (zh) 2023-02-03

Family

ID=85064815

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211228352.3A Pending CN115690001A (zh) 2022-10-09 2022-10-09 一种用于检测钢管焊接数字射线图像中缺陷的方法

Country Status (1)

Country Link
CN (1) CN115690001A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117788396A (zh) * 2023-12-18 2024-03-29 江苏省特种设备安全监督检验研究院 一种基于深度学习的dr图像缺陷智能识别算法

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117788396A (zh) * 2023-12-18 2024-03-29 江苏省特种设备安全监督检验研究院 一种基于深度学习的dr图像缺陷智能识别算法

Similar Documents

Publication Publication Date Title
CN111027547B (zh) 一种针对二维图像中的多尺度多形态目标的自动检测方法
CN111815601B (zh) 一种基于深度卷积自编码器的纹理图像表面缺陷检测方法
JP2024509411A (ja) 欠陥検出方法、装置及びシステム
CN112967243A (zh) 一种基于yolo的深度学习芯片封装裂纹缺陷检测方法
CN112102229A (zh) 一种基于深度学习的工业ct检测缺陷的智能识别方法
CN111667455A (zh) 一种刷具多种缺陷的ai检测方法
JP2020085546A (ja) 構造物の点検・補修支援システム
CN113469951B (zh) 一种基于级联区域卷积神经网络的轮毂缺陷检测方法
CN111382785A (zh) 一种gan网络模型及实现样本自动清洗、辅助标记的方法
CN114972759A (zh) 基于分级轮廓代价函数的遥感图像语义分割方法
CN115690001A (zh) 一种用于检测钢管焊接数字射线图像中缺陷的方法
Provencal et al. Identification of weld geometry from ultrasound scan data using deep learning
CN110334775B (zh) 一种基于宽度学习的无人机线路故障识别方法和装置
CN114639102A (zh) 基于关键点与尺寸回归的细胞分割方法和装置
CN111709936B (zh) 一种基于多级特征比对的令纸缺陷检测方法
CN110889418A (zh) 一种气体轮廓识别方法
CN114782822A (zh) 电力设备的状态检测方法、装置、电子设备及存储介质
CN116977239A (zh) 缺陷检测方法、装置、计算机设备及存储介质
García et al. A configuration approach for convolutional neural networks used for defect detection on surfaces
CN113269251A (zh) 基于机器视觉和深度学习融合的水果瑕疵分类方法、装置、存储介质及计算机设备
CN113034432A (zh) 一种产品缺陷检测方法、系统、装置及存储介质
CN116109627B (zh) 基于迁移学习和小样本学习缺陷检测方法、装置及介质
CN116977334B (zh) 一种光缆表面瑕疵检测方法及装置
CN117372720B (zh) 一种基于多特征交叉掩码修复的无监督异常检测方法
CN116778398B (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