1.本发明属于工业图像质量优化技术领域，具体涉及一种基于优化工业图像质量的缺陷识别方法。
背景技术：2.工业图像经过数字化后，由于光敏材料、射线能量和辐射剂量等众多因素的影响，导致数字化后的图像存在着对比度低、噪声高等缺点。
3.目前，对于工业图像的工件质量评估主要通过人为工作进行，然而，人工检查可能缺乏客观性、一致性、准确性和高效性，不同层次的操作人员的评估结果也会有所不同，而且，评估工作依赖于大量的本领域专业知识，通常是主观的和耗时的。此外，新手操作人员通常缺乏对管道焊缝缺陷进行准确评估所需的经验，受人为因素的影响较大。评价过程主要是通过人眼观察来进行，因此，算子无法获得精确的几何参数，如长度、面积和密度等等。这些因素构成了准确评估被测工件质量的主要障碍。
4.随着工业图像成像技术的发展，基于工业图像的质量增强技术得到了广泛的应用。在工业图像中面临的挑战可以总结为：质量差、亮度不均匀、脏斑点、roi区域与背景对比度低等。一般来说，每种质量增强技术都能很好地应对其中的某个问题，但是对于其他的一些问题可能会处理效果不理想。此外，大多数最近的方法都需要手动设置参数或者影响因素，而不是自适应的计算工业射线图像的相关因子。这也使得应用射线图像质量增强技术来辅助人工检查的成效并不显著。
技术实现要素：5.基于上述问题，本发明提供一种基于优化工业图像质量的缺陷识别方法，包括：
6.步骤1：依据光电技术和数字处理技术将工件信息转换成数字信号，得到原始工件图像；
7.步骤2：将工件图像依据roi区域类型进行分类：非线型区域和线型区域；
8.步骤3：对于非线型roi区域，采用灰度平衡化函数和对比度有限的自适应直方图均衡化方法相结合来提升图像质量；包括：
9.步骤3-1：计算被测工件图像的直方图，然后计算直方图的均值与方差值；
[0010][0011]
式中，e(x)和v(x)表示图像灰度直方图的均值和方差，g(x)表示灰度值为x的像素点个数，n表示总的像素点个数；
[0012]
步骤3-2：使用均值与方差来推导与平衡化系数的关系，公式如下：
[0013][0014]
式中，p代表平衡化系数，并通过p确定要平衡化的直方图的边界值(最小值、最大值)。该值越大，颜色失真就越大，该值越小，对比度增强的程度就越小。在公式(2)中，s是调整均值与方差比值的尺度因子。
[0015]
步骤3-3：利用平衡化系数和像素总数的比值计算期望的最小值和最大值，如公式(3)所示。在图像的累积直方图中，当直方图的像素点累计值小于n
×
p时，将此时直方图所对应的灰度值设置为期望的最小值。当直方图的像素点累计值大于n
×
(1-p)时，则将此时直方图所对应的灰度值设置为期望的最大值。
[0016][0017]
式中，hj表示图像直方图累计的像素点。v
min
和v
max
分别代表期望的最小值和最大值。
[0018]
步骤3-4：使用灰度平衡化公式(4)将原始直方图的值映射到0～255的范围。它可以将大多数像素分布的灰度范围映射到整个像素区间。
[0019][0020]
在公式(4)中，x是像素点的灰度值，v
min
、v
max
分别表示最小的和最大的灰度值。
[0021]
步骤3-5：将灰度平衡化后的图像划分为多个不重叠的区域；
[0022]
步骤3-6：分别计算每个区域的直方图；
[0023]
步骤3-7：计算一个剪切直方图的剪切限制β，公式如下：
[0024][0025]
其中β为剪切限制，m
×
n为每个区域的像素数，l为灰度级，α为0～100的剪切因子，s
max
为允许的最大斜率；
[0026]
步骤3-8：将每个直方图以其高度不超过剪切限制的方式重新进行分配。从公式(5)中可以得出，如果α＝0，则剪切限制
[0027]
步骤4：对于线型roi区域，采用背景衰减的自适应对比度增强算法来提升图像质量；具体表述为：
[0028]
步骤4-1：使用公式(6)对被测工件图像进行灰度值归一化处理用于获取浮点型图像，其灰度级范围在0到1之间；
[0029][0030]
其中，i(x,y)和i
max
分别为整个图像和输入的工件图像的最大值。d(x,y)和r(x,y)分别为细节分量和可变分量，这里，(x,y)是一个空间域索引，细节分量和可变分量，都在0到1之间；
[0031]
步骤4-2：为了寻找归一化图像的局部最大值分量g(y)和局部最小值分量t(y)，使用公式(7)计算局部最大值分量，并使用公式(8)计算背景成分，即局部最小分量。
[0032][0033][0034]
其中，loy表示像素点(x,y)周围的局部区域，(x,y)表示局部区域内的像素。t(x,y)分量可以认为是背景成分，是可以进行衰减的分量，以增强图像的对比度；
[0035]
步骤4-3：利用局部最大值分量和局部最小值分量根据公式(9)生成一个衰减因子λ(x,y)。
[0036][0037]
在公式(9)中，g(x,y)和t(x,y)为局部最大分量和局部最小分量，var(t(x,y))为局部最小分量的方差；
[0038]
步骤4-4：使用公式(10)生成可移动分量r(x,y)，它可以自适应的调节，以找到一个合适的可消除背景成分，以更好地增强图像的对比度。
[0039]
r(x,y)≡λ
·
t(x,y)
ꢀꢀꢀ
(10)
[0040]
衰减因子λ调节去除成分中背景分量的比例。因此，通过调节衰减因子λ来确定的可移动部分r(x,y)和可去除的背景分量t(x,y)，以有效地增强焊缝缺陷图像。
[0041]
步骤4-5：使用公式(11)确定另一个可调亮度参数ψ(x,y)，来控制增强图像的亮度。
[0042][0043]
其中，r(x,y)和g(x,y)为输入射线图像的可变分量和局部最大值；
[0044]
步骤4-6：最后利用公式(12)得到了增强的射线图像e(x,y)。
[0045][0046]
步骤5：将步骤3或步骤4处理过的射线图像e(x,y)进一步执行多尺度细节强化；具体表述为：
[0047]
步骤5-1：通过对射线图像e(x,y)应用三种尺度的高斯卷积核，得到三种不同程度的模糊图像，如公式(13)所示。
[0048][0049]
其中g1(x,y)、g2(x,y)和g3(x,y)分别为不同尺度的高斯卷积核。
[0050]
步骤5-2：我们提取精细细节d1(x,y)、中间细节d2(x,y)和粗糙细节d3(x,y)，如公式(14)所示；
[0051][0052]
步骤5-3：通过设计运算规则来合并这三层细节，以此来生成整体的细节图像，如公式(15)所示。
[0053]d*
(x,y)＝(1-ω1×
sgn(d1(x,y)))
×
d1(x,y)+ω2×
d2(x,y)+ω3×
d3(x,y)
ꢀꢀꢀ
(15)
[0054]
其中，ω1、ω2和ω3分别为精细细节、中间细节和粗糙细节的调节因子；
[0055]
步骤5-4：我们将公式(15)中得到的整体细节d
*
(x,y)添加到对比度增强后的射线图像e(x,y)中，得到多尺度细节强化图像。
[0056]
步骤6：阈值自适应二值化处理，具体表述为：
[0057]
步骤6-1：利用滑动窗口计算强化图像每个窗口的自适应标准差σa，计算方法由公式(16)给出。
[0058][0059]
式中，σw为窗口像素的标准差，σ
min
和σ
max
分别为待处理图像中所有窗口的最小标准差值和最大标准差值；
[0060]
步骤6-2：计算二值化的阈值，计算公式如式(17)所示：
[0061][0062]
式中，t为阈值，μw为窗口像素的平均值，σw为窗口像素的标准差，μg为图像中所有像素的平均值，σa为窗口的自适应标准差。
[0063]
步骤6-3：基于此阈值t，在式(18)中定义了二值化过程。
[0064][0065]
其中，i(x,y)为二值化图像，i(x,y)为原始工件图像的像素值，t为图像中局部窗口的阈值。
[0066]
步骤7：利用背景减法标注被测工件图像的潜在异常区域。
[0067]
步骤7-1：将步骤5所得的细多尺度节强化图像进行大尺度中值滤波，得到背景模型；
[0068]
步骤7-2：应用背景减法得到差分图像；
[0069]
步骤7-3：对差分图像进行步骤6所述阈值自适应二值化处理；
[0070]
步骤7-4：通过面积过滤和形态学操作来获得被测工件图像上的潜在异常区域；
[0071]
步骤7-5：对二值图像进行轮廓跟踪绘制潜在异常区域的最小外接矩形。
[0072]
本发明的有益效果是：
[0073]
本发明提供了一种基于优化工业图像质量的缺陷识别方法，根据被测工件图像的感兴趣区域(即roi区域)是否为线型区域进行不同的对比度增强手段，之后应用多尺度细节强化对roi区域细节信息更加显著，结合阈值自适应的二值化方法进行背景减法获取被测工件图像上的潜在异常区域，并进行标注，从而辅助提高人工检查效率。
附图说明
[0074]
图1为本发明实施例提供的一种工业图像质量提升方法的流程图；
[0075]
图2为本发明实施例提供的对于非线型roi区域图像质量提升的流程图；
[0076]
图3为本发明实施例提供的对于线型roi区域图像质量提升的流程图。
具体实施方式
[0077]
下面结合附图和具体实施实例对发明做进一步说明。以下实施例用于说明本发明，但不用来限制本发明的范围。
[0078]
本实施例中，一种工业图像质量提升方法，如图1所示：包括以下步骤：
[0079]
步骤1：依据光电技术和数字处理技术将工件信息转换成数字信号，得到原始工件图像；本实施例中直接通过工业x射线探伤机采集工件的图像数据。
[0080]
步骤2：将工件图像依据roi区域类型进行分类：非线型区域和线型区域；
[0081]
步骤3：对于非线型roi区域，采用灰度平衡化函数和对比度有限的自适应直方图均衡化方法相结合来提升图像质量，如图2所示。
[0082]
步骤3-1：计算被测工件图像的直方图，然后计算直方图的均值与方差值：
[0083][0084]
式中，e(x)和v(x)表示图像灰度直方图的均值和方差，g(x)表示灰度值为x的像素点个数，n表示总的像素点个数；
[0085]
步骤3-2：我们通过使用均值与方差来推导与平衡化系数的关系，公式如下：
[0086][0087]
式中，p代表平衡化系数，并通过p确定要平衡化的直方图的边界值(最小值、最大值)。该值越大，颜色失真就越大，该值越小，对比度增强的程度就越小。在公式(2)中，s是调整均值与方差比值的尺度因子，在本实施例中将其固定为0.4。
[0088]
步骤3-3：利用平衡化系数和像素总数的比值计算期望的最小值和最大值，如公式(3)所示。在图像的累积直方图中，当直方图的像素点累计值小于n
×
p时，将此时直方图所对应的灰度值设置为期望的最小值。当直方图的像素点累计值大于n
×
(1-p)时，则将此时直方图所对应的灰度值设置为期望的最大值。
[0089][0090]
式中，hj表示图像直方图累计的像素点。v
min
和v
max
分别代表期望的最小值和最大值。
[0091]
步骤3-4：使用灰度平衡化公式(4)将原始直方图的值映射到0～255的范围。它可以将大多数像素分布的灰度范围映射到整个像素区间。
[0092][0093]
在公式(4)中，x是像素点的灰度值，v
min
、v
max
分别表示最小的和最大的灰度值。
[0094]
步骤3-5：将灰度平衡化后的图像划分为多个不重叠的区域；本实施例中将整个图像划分为8*8个不重叠的区域；
[0095]
步骤3-6：分别计算每个区域的直方图；
[0096]
步骤3-7：计算一个剪切直方图的剪切限制β，公式如下：
[0097][0098]
其中，β为剪切限制，m
×
n为每个区域的像素数，l为灰度级，α为0～100的剪切因子，s
max
为允许的最大斜率；
[0099]
步骤3-8：将每个直方图以其高度不超过剪切限制的方式重新进行分配。从公式(5)中可以得出，如果α＝0，则剪切限制
[0100]
步骤4：对于线型roi区域，采用背景衰减的自适应对比度增强算法来提升图像质量，如图3所示。
[0101]
步骤4-1：使用公式(6)对被测工件图像进行灰度值归一化处理用于获取浮点型图像，其灰度级范围在0到1之间；
[0102][0103]
其中，i(x,y)和i
max
分别为整个图像和输入的工件图像的最大值。d(x,y)和r(x,y)分别为细节分量和可变分量，这里，(x,y)是一个空间域索引，细节分量和可变分量，都在0到1之间；
[0104]
步骤4-2：为了寻找归一化图像的局部最大值分量g(y)和局部最小值分量t(y)，我们使用15x15的窗口大小在每个像素点的周围创建局部区域，使用公式(7)计算局部最大值分量，并使用公式(8)计算背景成分，即局部最小分量。
[0105][0106][0107]
其中，loy表示像素点(x,y)周围的局部区域，(x,y)表示局部区域内的像素。t(x,y)分量可以认为是背景成分，是可以进行衰减的分量，以增强图像的对比度；
[0108]
步骤4-3：利用局部最大值分量和局部最小值分量根据公式(9)生成一个衰减因子λ(x,y)。
[0109][0110]
在公式(9)中，g(x,y)和t(x,y)为局部最大分量和局部最小分量，var(t(x,y))为局部最小分量的方差；
[0111]
步骤4-4：使用公式(10)生成可移动分量r(x,y)，它可以自适应的调节，以找到一个合适的可消除背景成分，以更好地增强图像的对比度。
[0112]
r(x,y)≡λ
·
t(x,y)
ꢀꢀꢀ
(10)
[0113]
衰减因子λ调节去除成分中背景分量的比例。因此，通过调节衰减因子λ来确定的可移动部分r(x,y)和可去除的背景分量t(x,y)，以有效地增强焊缝缺陷图像。
[0114]
步骤4-5：使用公式(11)确定另一个可调亮度参数ψ(x,y)，来控制增强图像的亮度。
[0115][0116]
其中，r(x,y)和g(x,y)为输入射线图像的可变分量和局部最大值；
[0117]
步骤4-6：最后利用公式(12)得到增强的射线图像e(x,y)。
[0118][0119]
步骤5：将步骤3或步骤4处理过的射线图像e(x,y)进一步执行多尺度细节强化；具体表述为：
[0120]
步骤5-1：通过对射线图像e(x,y)应用三种尺度的高斯卷积核，得到三种不同程度的模糊图像，如公式(13)所示。
[0121][0122]
其中，g1(x,y)、g2(x,y)和g3(x,y)分别为不同尺度的高斯卷积核，其标准差分别为σ1＝1.0、σ2＝2.0和σ3＝4.0。
[0123]
步骤5-2：我们提取精细细节d1(x,y)、中间细节d2(x,y)和粗糙细节d3(x,y)，如公式(14)所示；
[0124][0125]
步骤5-3：通过设计运算规则来合并这三层细节，以此来生成整体的细节图像，如公式(15)所示。
[0126]d*
(x,y)＝(1-ω1×
sgn(d1(x,y)))
×
d1(x,y)+ω2×
d2(x,y)+ω3×
d3(x,y)
ꢀꢀꢀ
(15)
[0127]
其中，ω1、ω2和ω3分别为精细细节、中间细节和粗糙细节的调节因子，固定为0.5、0.5和0.25；
[0128]
步骤5-4：我们将公式(15)中得到的整体细节d
*
(x,y)添加到对比度增强后的射线图像e(x,y)中，得到多尺度细节强化图像。
[0129]
步骤6：阈值自适应二值化处理，具体表述为：
[0130]
步骤6-1：利用滑动窗口计算强化图像每个窗口的自适应标准差σa，计算方法由公式(16)给出。
[0131]
[0132]
式中，σw为窗口像素的标准差，σ
min
和σ
max
分别为待处理图像中所有窗口的最小标准差值和最大标准差值；
[0133]
步骤6-2：计算二值化的阈值，计算公式如式(17)所示：
[0134][0135]
式中，t为阈值，μw为窗口像素的平均值，σw为窗口像素的标准差，μg为图像中所有像素的平均值，σa为窗口的自适应标准差。
[0136]
步骤6-3：基于此阈值t，在式(18)中定义了二值化过程。
[0137][0138]
其中，i(x,y)为二值化图像，i(x,y)为原始工件图像的像素值。
[0139]
步骤7：利用背景减法标注被测工件图像的潜在异常区域，具体表述为：
[0140]
步骤7-1：将步骤5所得的细多尺度节强化图像进行大尺度中值滤波，得到背景模型；
[0141]
步骤7-2：应用背景减法得到差分图像；
[0142]
步骤7-3：对差分图像进行步骤6所述阈值自适应二值化处理；
[0143]
步骤7-4：通过面积过滤和形态学操作来获得被测工件图像上的潜在异常区域；
[0144]
步骤7-5：对二值图像进行轮廓跟踪绘制潜在异常区域的最小外接矩形。
[0145]
最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明权利要求所限定的范围。
技术特征：1.一种基于优化工业图像质量的缺陷识别方法，其特征在于，包括：步骤1：依据光电技术和数字处理技术将工件信息转换成数字信号，得到原始工件图像；步骤2：将工件图像依据roi区域类型分类为非线型区域和线型区域；步骤3：对于非线型roi区域，采用灰度平衡化函数和对比度有限的自适应直方图均衡化方法相结合来提升图像质量；步骤4：对于线型roi区域，采用背景衰减的自适应对比度增强算法来提升图像质量；步骤5：将步骤3或步骤4处理过的射线图像e(x,y)执行多尺度细节强化；步骤6：进行阈值自适应二值化处理；步骤7：利用背景减法标注被测工件图像的潜在异常区域。2.根据权利要求1所述的一种基于优化工业图像质量的缺陷识别方法，其特征在于，所述步骤3包括：步骤3-1：计算被测工件图像的直方图，然后计算直方图的均值与方差值；式中，e(x)和v(x)表示图像灰度直方图的均值和方差，g(x)表示灰度值为x的像素点个数，n表示总的像素点个数；步骤3-2：使用均值与方差来推导与平衡化系数的关系，公式如下：式中，p代表平衡化系数，并通过p确定要平衡化的直方图的边界值，s是调整均值与方差比值的尺度因子；步骤3-3：利用平衡化系数和像素总数的比值计算期望的最小值和最大值，如公式(3)所示：式中，h
j
表示图像直方图累计的像素点，v
min
和v
max
分别代表期望的最小值和最大值；在图像的累积直方图中，当直方图的像素点累计值小于n
×
p时，将此时直方图所对应的灰度值设置为期望的最小值；当直方图的像素点累计值大于n
×
(1-p)时，则将此时直方图所对应的灰度值设置为期望的最大值；步骤3-4：使用灰度平衡化公式(4)将原始直方图的值映射到0～255的范围；式中，x是像素点的灰度值，v
min
、v
max
分别表示最小的和最大的灰度值；步骤3-5：将灰度平衡化后的图像划分为多个不重叠的区域；
步骤3-6：分别计算每个区域的直方图；步骤3-7：计算一个剪切直方图的剪切限制β，公式如下：其中，β为剪切限制，m
×
n为每个区域的像素数，l为灰度级，α为0～100的剪切因子，s
max
为允许的最大斜率；步骤3-8：将每个直方图以其高度不超过剪切限制的方式重新进行分配；从公式(5)中可以得出，如果α＝0，则剪切限制3.根据权利要求1所述的一种基于优化工业图像质量的缺陷识别方法，其特征在于，所述步骤4包括：步骤4-1：使用公式(6)对被测工件图像进行灰度值归一化处理用于获取浮点型图像，其灰度级范围在0到1之间；其中，i(x,y)和i
max
分别为整个图像和输入的工件图像的最大值；d(x,y)和r(x,y)分别为细节分量和可变分量，这里，(x,y)是一个空间域索引，细节分量和可变分量，都在0到1之间；步骤4-2：为了寻找归一化图像的局部最大值分量g(y)和局部最小值分量t(y)，使用公式(7)计算局部最大值分量，使用公式(8)计算背景成分，即局部最小分量；式(7)计算局部最大值分量，使用公式(8)计算背景成分，即局部最小分量；其中，lo
y
表示像素点(x,y)周围的局部区域，(x,y)表示局部区域内的像素；步骤4-3：利用局部最大值分量和局部最小值分量根据公式(9)生成一个衰减因子λ(x,y)；式中，g(x,y)和t(x,y)为局部最大分量和局部最小分量，var(t(x,y))为局部最小分量的方差；步骤4-4：使用公式(10)生成可移动分量r(x,y)：其中，λ表示衰减因子去除成分中背景分量的比例；步骤4-5：使用公式(11)确定另一个可调亮度参数ψ(x,y)，来控制增强图像的亮度；步骤4-6：利用公式(12)得到增强的射线图像e(x,y)；
4.根据权利要求1所述的一种基于优化工业图像质量的缺陷识别方法，其特征在于，所述步骤5包括：步骤5-1：通过对射线图像e(x,y)应用三种尺度的高斯卷积核，得到三种不同程度的模糊图像，如公式(13)所示：其中，g1(x,y)、g2(x,y)和g3(x,y)分别为不同尺度的高斯卷积核；步骤5-2：提取精细细节d1(x,y)、中间细节d2(x,y)和粗糙细节d3(x,y)，如公式(14)所示；步骤5-3：通过设计运算规则来合并这三层细节，以此来生成整体的细节图像，如公式(15)所示：d
*
(x,y)＝(1-ω1×
sgn(d1(x,y)))
×
d1(x,y)+ω2×
d2(x,y)+ω3×
d3(x,y)
ꢀꢀꢀꢀ
(15)其中，d
*
(x,y)为整体细节，ω1、ω2和ω3分别为精细细节、中间细节和粗糙细节的调节因子；步骤5-4：将公式(15)中得到的整体细节d
*
(x,y)添加到对比度增强后的射线图像e(x,y)中，得到多尺度细节强化图像。5.根据权利要求1所述的一种基于优化工业图像质量的缺陷识别方法，其特征在于，所述步骤6包括：步骤6-1：利用滑动窗口计算强化图像每个窗口的自适应标准差σ
a
，如公式(16)所示：式中，σ
w
为窗口像素的标准差，σ
min
和σ
max
分别为待处理图像中所有窗口的最小标准差值和最大标准差值；步骤6-2：计算二值化的阈值，如公式(17)所示：式中，t为阈值，μ
w
为窗口像素的平均值，σ
w
为窗口像素的标准差，μ
g
为图像中所有像素的平均值，σ
a
为窗口的自适应标准差；步骤6-3：设计二值化过程生成二值化图像，如公式(18)所示：其中，i(x,y)为二值化图像，i(x,y)为原始工件图像的像素值，t为图像中局部窗口的
阈值。6.根据权利要求1所述的一种基于优化工业图像质量的缺陷识别方法，其特征在于，所述步骤7包括：步骤7-1：将步骤5所得的细多尺度节强化图像进行大尺度中值滤波，得到背景模型；步骤7-2：应用背景减法得到差分图像；步骤7-3：对差分图像进行步骤6所述阈值自适应二值化处理；步骤7-4：通过面积过滤和形态学操作来获得被测工件图像上的潜在异常区域；步骤7-5：对二值图像进行轮廓跟踪绘制潜在异常区域的最小外接矩形。
技术总结本发明提供一种基于优化工业图像质量的缺陷识别方法，涉及工业图像质量优化技术领域；本发明针对工业图像上的ROI区域类型是否为线型分别进行对比度增强处理，再应用多尺度细节强化方法将对比度增强后的射线图像的ROI区域细节更加显著化，最后结合阈值自适应的二值化方法进行背景减法获取被测工件图像上的潜在异常区域，并进行标注，从而辅助提高人工检查效率。检查效率。检查效率。
技术研发人员：王赢博 鲁丹宇 那新 贾迪
受保护的技术使用者：沈阳智谷科技有限公司
技术研发日：2022.06.24
技术公布日：2022/11/1
