CN103745234A - 一种带钢表面缺陷的特征提取与分类方法 - Google Patents
一种带钢表面缺陷的特征提取与分类方法 Download PDFInfo
- Publication number
- CN103745234A CN103745234A CN201410038260.8A CN201410038260A CN103745234A CN 103745234 A CN103745234 A CN 103745234A CN 201410038260 A CN201410038260 A CN 201410038260A CN 103745234 A CN103745234 A CN 103745234A
- Authority
- CN
- China
- Prior art keywords
- sample
- defect
- sigma
- size
- pixel
- 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.)
- Granted
Links
Images
Landscapes
- Image Analysis (AREA)
Abstract
一种带钢表面缺陷的特征提取与分类方法,属于模式识别与图像处理领域。提取带钢表面缺陷样本数据库的基准采样尺寸表;获得基准采样图像,构造梯度大小-方向共生矩阵;针对基准采样图像的缺陷内区域,构造灰度大小-方向共生矩阵;生成特征向量样本的训练库;用K-最近邻和R-最近邻相结合的方法,修剪训练样本集并提取倍数因子;并利用修剪后样本的倍数因子改进分类器;获得多类别分类器模型;根据基准采样尺寸表,将缺陷测试样本转为基准采样图像,然后提取25维的特征量,并输入到多类别分类器模型,完成缺陷的自动识别。本发明能够实现尺度和旋转不变,抑制其它不利因素的影响,提高识别的效率与精度。
Description
技术领域
本发明属于模式识别与图像处理领域,具体涉及带钢表面缺陷识别方向的缺陷特征提取与缺陷分类的方法。
背景技术
近年来,随着不断增长的对高品质带钢产品的需求和日益白热化的市场竞争,表面缺陷的检测成为钢铁企业实施带钢质量监测与控制的重要技术支撑。在带钢的生产和加工过程中会出现各种类型的缺陷,比如:裂纹、结疤、孔洞、氧化铁皮、卷边、擦伤等。表面缺陷监测的一个重要过程就是对这些缺陷进行识别与区分,以便对产品质量问题做出快速的处理和控制。带钢表面缺陷识别的过程主要包括四个环节:缺陷预处理、缺陷分割、缺陷特征提取和缺陷分类。其中,缺陷的特征提取和分类是两个关键性环节。针对它们的研究也越来越广泛与深入。
针对带钢表面缺陷,同一个样本可能存在着不同的尺度、旋转、平移、亮度、视角上的图像版本。在实际的带钢生产线,具有固定环境和位置的表面缺陷采集装置避免了亮度与视角上的明显差异,同时,采用图像分割的方法将缺陷区域从图像中提取出来,避免了平移上的差异。但是,同一个缺陷样本的不同图像在尺度、旋转上的差异不可避免。这就要求缺陷特征提取算法具备尺度、旋转的不变性。在传统的表面缺陷特征提取技术中,大部分算法都不具备尺度、旋转不变性,仅有少量的算法考虑了尺度、旋转不变性的要求,如:学者HU提出的不变矩算法。但是,HU不变矩算法仅提取7个具有尺度、旋转不变性的特征量,无法满足后续分类方法依靠多特征量来提高精度的需求。因此,研究其它的具有尺度、旋转不变性的缺陷特征提取算法以提取更多的特征量势在必行。
另外,针对带钢表面缺陷的分类,支持向量机被证实是最为有效的方法之一。但是,在实际的分类问题中,传统的支持向量机存在着分类精度与效率的冲突。幸运的是,随着技术的发展,学者们提出了性能得到提升的新型支持向量机。如:近似支持向量机、孪生支持向量机、最小二乘孪生支持向量机等。因此,选择有效的新型支持向量机,并加以改进,以适合带钢表面缺陷的高精度与高效率的多类别分类要求。
发明内容
针对现有技术存在的缺陷,本发明的目的是提供一种带钢表面缺陷的特征提取与分类方法,采用具有尺度、旋转不变性的特征提取技术,以达到多类别缺陷分类的高效率与高精度要求的目的。
本发明的技术方案是这样实现的:一种带钢表面缺陷的特征提取与分类方法,包括以下步骤:
步骤1:提取带钢表面缺陷样本数据库的基准采样尺寸表;
步骤2:计算缺陷样本数据库中样本图像的高度-宽度比,并检索基准采样尺寸表,然后按照检索到的尺寸对缺陷图像进行采样,获得基准采样图像,基准采样图像保证了后续算法提取的特征具有尺度不变性;
步骤3:针对基准采样图像的缺陷边缘,构造梯度大小-方向共生矩阵,并根据此矩阵提取9个梯度大小-方向共生矩阵的特征,即大小均值、方向均值、大小方差、方向方差、大小熵、方向熵、混合熵、惯性矩、相关性,使提取到的特征具有尺度不变性和旋转不变性;
步骤4:针对基准采样图像的缺陷内区域,构造灰度大小-方向共生矩阵,并根据此矩阵提取9个灰度大小-方向共生矩阵的特征,使这些特征具有尺度不变性和旋转不变性;
步骤5:将缺陷边缘的9个梯度大小-方向共生矩阵的特征量、缺陷内区域的9个灰度大小-方向共生矩阵的特征量和缺陷图像的7个HU不变矩算法的特征量结合,生成具有25维特征向量样本的训练库;
步骤6:针对具有25维特征向量样本的训练库,使用K-最近邻和R-最近邻相结合的方法,修剪训练样本集,同时,提取修剪后样本的倍数因子;
步骤7:采用最小二乘孪生支持向量机作为分类器,并利用修剪后样本的倍数因子改进分类器;
步骤8:将改进的分类器与二叉树结合,建立多类别分类器模型,并利用修剪后的训练样本集对多类别分类器模型进行训练;
步骤9:根据基准采样尺寸表,将缺陷测试样本转为基准采样图像,然后提取25维的特征量,并输入到多类别分类器模型,完成缺陷的自动识别。
步骤1所述的提取带钢表面缺陷样本数据库的基准采样尺寸表,步骤如下:
步骤1-1:对带钢表面缺陷样本数据库的所有样本进行[0°,360°)角度范围内的旋转,并计算所有样本对应的每一个旋转角度下的图像尺寸;
步骤1-2:计算所有样本对应的图像尺寸的高度-宽度比rhw,底标hw是高度-宽度的英文缩写;
步骤1-3:分析每一个rhw对应的所有图像尺寸,确定一个基准采样尺寸Shw;
步骤1-4:建立rhw与Shw一一对应的基准采样尺寸表;
步骤3所述的提取9个梯度大小-方向共生矩阵特征的步骤如下:
步骤3-1:假设基准采样图像中缺陷边缘有nL个像素点,每一个像素点的灰度值为fi L,以fi L像素点为中心,提取3×3窗口中属于缺陷边缘的像素点,组成集合Pi L,且该集合中像素点的个数为这里L符号代表缺陷边缘,i=1,2,…,nL;
步骤3-2:计算Pi L中每一个像素点梯度的大小(GM)i(k0,l0)和方向(GO)i(k0,l0),并定义Pi L中所有像素点相对于3×3窗口中心点的坐标偏移量集合为Wi L,这里k0和l0表示Wi L中的坐标偏移量,G符号代表梯度,M符号代表大小,O符号代表方向;
步骤3-4:为保证旋转不变性,将fi L像素点的梯度方向旋转到0方向,对应的Pi L中所有像素点的梯度方向也随之旋转,最终,缺陷边缘的每一个像素点的梯度方向被定义为Pi L中所有像素点旋转后的梯度方向均值,假设MOD符号表示取余运算,则的公式如下:
步骤3-5:将和的值域进行32等级离散化,且离散化后缺陷边缘像素点的梯度大小和方向在离散域[1,2,…,32]中取值,这里i1,i2=1,2,…,nL,M32符号代表32等级离散化后的大小,O32符号代表32等级离散化后的方向;
步骤3-6:针对缺陷边缘,统计梯度大小为和梯度方向为的缺陷边缘像素点的个数,并将梯度大小-方向共生矩阵HL中k1行l1列的元素HL(k1,l1)设定为此值,构造出32×32梯度大小-方向共生矩阵HL;
步骤3-7:将梯度大小-方向共生矩阵HL的元素转为概率形式,即HL在(k1,l1)元素的概率为p(k1,l1),其公式如下:
步骤3-8:根据缺陷边缘的梯度大小-方向共生矩阵HL,计算对应的具有尺度、旋转不变性的9个梯度大小-方向共生矩阵的特征,即:大小均值μM、方向均值μO、大小方差方向方差大小熵EM、方向熵EO、混合熵EMO、惯性矩IMO、相关性CMO,具体公式如下:
步骤4所述的提取9个灰度大小-方向共生矩阵特征的步骤如下:
步骤4-1:假设基准采样图像中缺陷内区域有nΩ个像素点,每一个像素点的灰度值为ft Ω,以ft Ω像素点为中心,提取3×3窗口中所有的像素点,组成集合Pt Ω,这里Ω符号代表缺陷内区域,t=1,2,…,nΩ;
步骤4-2:计算Pt Ω中像素点梯度的方向(GO)t(k2,l2),这里k2和l2表示Pt Ω中像素点相对于3×3窗口中心点的坐标偏移量,且k2,l2∈{-1,0,1};
步骤4-3:为保证旋转不变性,将ft Ω像素点的梯度方向旋转到0方向,对应的Pt Ω中所有像素点的梯度方向也随之旋转,最终,缺陷内区域的每一个像素点的灰度方向被定义为Pt Ω中所有像素点旋转后的梯度方向均值,其公式如下:
步骤4-5:针对缺陷内区域,统计灰度大小为和灰度方向为的缺陷内区域像素点的个数,并将灰度大小-方向共生矩阵(k3,l3)位置的元素HΩ(k3,l3)设定为此值,构造出32×32灰度大小-方向共生矩阵HΩ;
步骤4-6:将灰度大小-方向共生矩阵HΩ的元素转为概率形式,即HΩ在(k3,l3)的概率为p(k3,l3),其公式如下:
步骤4-7:根据缺陷内区域的灰度大小-方向共生矩阵HΩ,计算对应的具有尺度、旋转不变性的9个灰度大小-方向共生矩阵的特征,即:大小均值、方向均值、大小方差、方向方差、大小熵、方向熵、混合熵、惯性矩、相关性,具体公式与步骤3-8公式的形式一致。
步骤6所述的修剪训练样本集和提取修剪后样本的倍数因子的步骤如下:
步骤6-1:假设训练样本集X=[X1,X2,…,Xl],样本个数为l,每一个样本为25维的特征向量,计算样本距离dab=||Xa-Xb||(a,b=1,2,…,l);
步骤6-2:根据dab和由用户选择的K值,获得每一个样本的K个最近邻样本(K-最近邻法),并确定该样本对应的K-最近邻矩阵的元素kab,其公式如下:
步骤6-7:针对修剪后的样本集XS,重复使用步骤6-5和6-6的方法,直至XS样本集被修剪为空,此时获得所有修剪后的样本Aj和对应的倍数因子矩阵的对角线元素(ρ1)jj,并最终确定针对训练样本集X的被修剪后的样本矩阵A=[A1,A2,…,An1]和具有对角线形式的倍数因子矩阵ρ1,这里j=1,2,…,n1,n1为样本集X被修剪后的样本数;
步骤6-8:使用相同的方法,可获得针对训练样本集Y的被修剪后的样本矩阵B=[B1,B2,…,Bn2]和倍数因子矩阵ρ2,这里n2为样本集Y被修剪后的样本数。
步骤7所述的利用修剪后样本的倍数因子改进最小二乘孪生支持向量机分类器的过程如下:
针对被修剪后的训练样本矩阵A=[A1,A2,…,An1]和B=[B1,B2,…,Bn2],利用相应的倍数因子ρ1和ρ2修改最小二乘孪生支持向量机的模型,其修改后的模型如下:
式中,C=[A′ B′]′,C1、C2为惩罚因子,表示核函数,e1,e2表示全1的向量,ξ1、ξ2、η1、η2为误差量,w1、w2、b1、b2为两个向量机模型的最优化解,倍数因子ρ1,ρ2具有对角阵形式,其中,C1、C2、由用户确定。
对上述两个模型进行推导,可求得对应的解:
步骤1所述的带钢表面缺陷样本数据库,是从现场的带钢表面缺陷检测***中提取的样本,且包含六种类型缺陷样本集,它们是:裂纹、结疤、孔洞、氧化铁皮、卷边、擦伤。
步骤6-1和6-8所述的训练样本集X和Y是六种类型缺陷样本集的任意两个。
本发明的有益效果:本发明带钢表面缺陷的特征提取与分类方法,具有如下优点:
1)利用训练样本库的基准采集尺寸将样本转换为基准采样图像,可以保证特征提取的尺度不变性;
2)利用缺陷边缘的梯度大小-方向共生矩阵模型和缺陷内区域的灰度大小-方向共生矩阵模型进行特征的提取,可以实现尺度和旋转不变性的特征提取,并抑制其它不利因素的影响,从而提高后续识别的精度;
3)利用K-最近邻和R-最近邻相结合的方法修剪训练样本,可以提高识别执行的效率;
4)使用新型的最小二乘孪生支持向量机模型进行分类,可以提高识别的效率与精度;
5)利用修剪后样本的倍数因子改进最小二乘孪生支持向量机模型,可以保证修剪后的样本集训练分类器模型的精度;
6)改进的分类器模型与二叉树结合,实现了多类别缺陷的分类。
附图说明
图1为本发明一种实施方式带钢表面缺陷的特征提取与分类方法流程图;
图2为本发明一种实施方式一个缺陷样本对应16种旋转角度图像的示意图;
图3为本发明一种实施方式一个缺陷图像转换为基准采样图像的流程图;
图4为本发明一种实施方式一个缺陷图像、缺陷边缘和缺陷内区域的示意图;
图5为本发明一种实施方式多类别分类器模型示意图。
具体实施方式
下面结合附图对本发明的实施例作进一步详细的说明。
本发明实施例采用一种带钢表面缺陷的特征提取与分类方法对带钢的缺陷进行处理,流程如图1所示,包括以下步骤:
步骤1:提取带钢表面缺陷样本数据库的基准采样尺寸表;
基准采样尺寸表是在分析带钢表面缺陷训练样本库的基础上获取的,它可以避免尺度对缺陷样本特征提取的影响。本发明实施例中的训练库是从现场的带钢表面缺陷检测***中提取的样本,且包含六种类型缺陷样本集,它们是:裂纹、结疤、孔洞、氧化铁皮、卷边和擦伤。
步骤1-1:对带钢表面缺陷样本数据库的所有样本进行[0°,360°)角度范围内的旋转,并计算所有样本对应的每一个旋转角度下的图像尺寸;
样本库中的一个缺陷样本,会存在不同旋转角度的版本。本发明实施例通过人工旋转[0°,360°)角度的方式来获取不同旋转角度的图像。如果人工旋转的角度为θ,那么这个缺陷样本就具有360/θ个旋转角度的图像。可以看出,θ值决定了人工旋转样本角度的准确度和效率。为了兼顾旋转的准确度和效率,本发明实施例选取θ=22.5°,即一个缺陷样本对应16种旋转角度图像,如图2所示。可以看出,16个旋转角度的图像一共具有8个不同的图像尺寸(高度和宽度),且4个图像的高度和宽度比与另外4个图像的高度和宽度比刚好倒置。因此,对于本实施例中的一个缺陷样本,只需人工旋转出0°、22.5°、45°、67.5°四个旋转角度的图像,就可确定完整的16个旋转角度的图像。
步骤1-2:计算所有样本对应的图像尺寸的高度-宽度比rhw,这里符号r表示比值,符号hw代表高度-宽度;
步骤1-3:分析每一个rhw对应的所有图像尺寸,确定一个基准采样尺寸Shw;
对于相同rhw所对应的图像,会存在不同的图像尺寸(高度和宽度)。分析图像的高度和宽度,确定一个适中的尺寸作为该rhw对应的基准采样尺寸Shw。Shw不应过大,否则转换为基准采样图像时耗费的时间过多;Shw也不应过小,否则转换后的基准采样图像分辨率过低。确定Shw的第一种方法是,将rhw对应的所有图像尺寸进行累加,然后求尺寸均值S,将Shw设定为确定Shw的第二种方法是,将rhw对应的所有图像尺寸进行排序,获得最大尺寸Sm,将Shw定义为3Sm/4。本发明实施例采用的是第二种方法。
步骤1-4:建立rhw与Shw一一对应的基准采样尺寸表。
基准采样尺寸表包含了带钢表面缺陷训练库中每一个样本不同旋转角度下的图像尺寸。对于一个缺陷图像,将其转换为对应的基准采样尺寸下的图像,可有效的保证后续特征提取的尺度不变性。
步骤2:计算缺陷样本数据库中样本图像的高度-宽度比,并检索基准采样尺寸表,然后按照检索到的尺寸对缺陷图像进行采样,获得基准采样图像,而基准采样图像保证了后续算法提取的特征具有尺度不变性;
对于一个缺陷图像,首先通过计算rhw和检索基准采样尺寸表来确定Shw。先计算该图像的高度-宽度比,根据这个比值查基准采样尺寸表获得该比值下图像应该具有的尺寸,这个尺寸比缺陷图像的尺寸大或者小,这就需要对缺陷图像进行上采样或下采样(如,对缺陷图像尺寸进行放大或缩小),获得基准采样尺寸下的图像。基准采样图像保证了后续算法提取的特征具有尺度不变性,因此本发明实施例的特征提取都将基于基准采样图像进行。图3展示了一个缺陷图像转换为基准采样图像的流程。
步骤3:针对基准采样图像的缺陷边缘,构造梯度大小-方向共生矩阵,并根据此矩阵提取9个梯度大小-方向共生矩阵的特征,即大小均值μM、方向均值μO、大小方差方向方差大小熵EM、方向熵EO、混合熵EMO、惯性矩IMO、相关性CMO,使提取到的特征具有尺度不变性和旋转不变性;
对于一个缺陷图像,经过缺陷预处理和缺陷分割方法,可以确定缺陷的边缘和缺陷的内区域。图4展示了一个缺陷图像、缺陷边缘和缺陷内区域。缺陷图像的边缘具有剧烈的灰度变化特征,常采用梯度来描述这种变化。基于梯度的矢量信息和灰度共生矩阵的思想,本发明提出了梯度大小-方向共生矩阵。梯度大小-方向共生矩阵能够完整的反映边缘点的灰度的变化大小统计、方向统计和联合统计特征。
步骤3-1:假设基准采样图像中缺陷边缘有nL个像素点,每一个像素点的灰度值为fi L,以fi L像素点为中心,提取3×3窗口中属于缺陷边缘的像素点,组成集合Pi L,且该集合中像素点的个数为这里L符号代表缺陷边缘,i=1,2,…,nL;
步骤3-2:计算Pi L中每一个像素点梯度的大小(GM)i(k0,l0)和方向(GO)i(k0,l0),并定义Pi L中所有像素点相对于3×3窗口中心点的坐标偏移量集合为Wi L,这里k0和l0表示Wi L中的坐标偏移量,G符号代表梯度,M符号代表大小,O符号代表方向;
本发明实施例定义fi L=f(xi,yi),即缺陷边缘灰度值为fi L的像素点在原缺陷图像中坐标为(xi,yi),且灰度值为f(xi,yi),Pi L中所有像素点坐标偏移量集合Wi L可表示为:
Wi L={k0,l0|k0,l0∈{-1,0,1};f(xi+k0,yi+l0)∈Pi L}
Pi L中每一个像素点梯度的计算采用差分运算替代。梯度的差分运算方法有多种,如:水平垂直差分法、Roberts算子、Sobel算子等。本发明实施例的梯度大小(GM)i(k0,l0)和方向(GO)i(k0,l0)的计算过程如下:
dxi(k0,l0)=(f(xi+k0+1,yi+l0)-f(xi+k0-1,yi+l0))/2
dyi(k0,l0)=(f(xi+k0,yi+l0+1)-f(xi+k0,yi+l0-1))/2
(GM)i(k0,l0)=[dxi(k0,l0)2+dxi(k0,l0)2]0.5 k0,l0∈Wi L
这里,MOD符号表示取余运算。x方向差分dxi(k0,l0)和y方向差分dyi(k0,l0)的计算采用间隔为2的水平垂直差分处理方式,可更好的捕捉边缘的突变量。梯度的大小(GM)i(k0,l0)被确定为梯度矢量的模。梯度的方向(GO)i(k0,l0)被确定为灰度变化最快的方向,且值被调整到[0,2π)范围,以符合实际的角度变化范围。
步骤3-3:缺陷边缘的每一个像素点的梯度大小被定义为Pi L中所有像素点梯度大小的均值,其公式如下:
由于缺陷边缘是在基准采样图像上进行提取的,因此缺陷边缘像素点的梯度大小已经具备了尺度不变性。梯度的大小本身具有旋转不变性,而且使用局部区域的梯度统计均值作为边缘像素点的梯度大小,可以有效的降低其它因素(如:噪声)对梯度大小的影响。因此,基于基准采样图像的边缘点梯度大小,具有尺度和旋转不变性,且具有一定的抑制其它因素影响的能力。
步骤3-4:为保证旋转不变性,将fi L像素点的梯度方向旋转到0方向,对应的Pi L中所有像素点的梯度方向也随之旋转。最终,缺陷边缘的每一个像素点的梯度方向被定义为Pi L中所有像素点旋转后的梯度方向均值,假设MOD符号表示取余运算,则的公式如下:
基于基准采样图像的缺陷边缘像素点,其梯度方向已经具备了尺度不变性。但是,缺陷边缘像素点的绝对梯度方向受旋转因素影响。为此,将Pi L中所有像素点的梯度方向旋转灰度值为fi L像素点的梯度方向角度,这样做的原因是边缘点局部区域内的所有像素点对于边缘点的相对梯度方向具备了旋转不变性。因此,旋转后的Pi L中所有像素点的相对梯度方向保证了尺度和旋转不变性。
步骤3-5:将和的值域进行32等级离散化,且离散化后缺陷边缘像素点的梯度大小和方向在离散域[1,2,…,32]中取值,这里i1,i2=1,2,…,nL,M32符号代表32等级离散化后的大小,O32符号代表32等级离散化后的方向;
和值域的32等级离散化,压缩了值域的区间,虽然精度上有了一定的损失,但是提高了运算效率、加强了抵御其它因素(如:噪声)的不利影响。这里,32等级离散化的过程,就是将值域进行32等分,且每一个等分区间对应离散域[1,2,…,32]中的一个值。本发明实施例假定为梯度大小的最大值,为梯度方向的最大值,INT为取整运算,则离散化后的缺陷边缘梯度大小和方向的公式如下:
步骤3-6:针对缺陷边缘,统计梯度大小为和梯度方向为的缺陷边缘像素点的个数,并将梯度大小-方向共生矩阵HL中k1行l1列的元素HL(k1,l1)设定为此值,构造出32×32梯度大小-方向共生矩阵HL;
32等级离散化后的梯度大小和方向具有32个等级,因此k1和l1的范围为[1,2,…,32]。那么,构造的梯度大小-方向共生矩阵HL的大小为32×32,且HL的元素是对具有相同梯度大小和梯度方向的像素点个数的统计。梯度大小-方向共生矩阵采用了灰度共生矩阵的思想,它既能反映出独立的梯度大小和梯度方向的统计信息,又能反映出梯度大小和方向的联合矢量统计信息。
步骤3-7:将梯度大小-方向共生矩阵HL的元素转为概率形式,即HL在(k1,l1)元素的概率为p(k1,l1),其公式如下:
步骤3-8:根据缺陷边缘的梯度大小-方向共生矩阵HL,计算对应的具有尺度、旋转不变性的9个梯度大小-方向共生矩阵的特征,即:大小均值μM、方向均值μO、大小方差方向方差大小熵EM、方向熵EO、混合熵EMO、惯性矩IMO、相关性CMO,具体公式如下:
大小均值μM和方向均值μO反映了梯度大小和梯度方向的集中趋势或总体水平;大小方差和方向方差反映了梯度大小和梯度方向的离散程度;大小熵EM和方向熵EO是对梯度大小和梯度方向有序化程度的度量,大小熵的值越大表示梯度大小的量越混乱,大小熵越小表示梯度大小的量越有序;混合熵EMO是对梯度大小和梯度方向联合统计量有序程度的度量;惯性矩IMO是梯度大小和梯度方向对图像清晰程度的联合反映;相关性CMO是梯度大小和梯度方向对图像相关程度的联合反映。
步骤4:针对基准采样图像的缺陷内区域,构造灰度大小-方向共生矩阵,并根据此矩阵提取9个灰度大小-方向共生矩阵的特征,使这些特征具有尺度不变性和旋转不变性。
缺陷的内区域具有缺陷的灰度变化特征与纹理分布特征,如图4所示。基于灰度信息和梯度方向信息,本发明提出了灰度大小-方向共生矩阵。灰度大小-方向共生矩阵能够反映缺陷内区域的灰度分布规律,梯度方向特征规律以及纹理特征规律。
步骤4-1:假设基准采样图像中缺陷内区域有nΩ个像素点,每一个像素点的灰度值为ft Ω,以ft Ω像素点为中心,提取3×3窗口中所有的像素点,组成集合Pt Ω,这里Ω符号代表缺陷内区域,t=1,2,…,nΩ;
步骤4-2:计算Pt Ω中像素点梯度的方向(GO)t(k2,l2),这里k2和l2表示Pt Ω中像素点相对于3×3窗口中心点的坐标偏移量,且k2,l2∈{-1,0,1};
本发明的实施例定义ft Ω=f(xt,yt),即缺陷内区域灰度值为ft Ω的像素点在原缺陷图像中坐标为(xt,yt),且灰度值为f(xt,yt)。Pt Ω中每一个像素点梯度的计算采用差分运算替代。本发明实施例采用水平垂直差分法,其梯度方向(GO)t(k2,l2)计算过程如下:
dxt(k2,l2)=f(xt+k2+1,yt+l2)-f(xt+k2,yt+l2)
dyt(k2,l2)=f(xt+k2,yt+l2+1)-f(xt+k2,yt+l2)
步骤4-3:为保证旋转不变性,将ft Ω像素点的梯度方向旋转到0方向,对应的Pt Ω中所有像素点的梯度方向也随之旋转。最终,缺陷内区域的每一个像素点的灰度方向被定义为Pt Ω中所有像素点旋转后的梯度方向均值,其公式如下:
基于基准采样图像的缺陷内区域像素点,其梯度方向已经具备了尺度不变性。但是,缺陷内区域像素点的绝对梯度方向受旋转因素影响。为此,将Pt Ω中所有像素点的梯度方向旋转灰度值为ft Ω像素点的梯度方向角度,这样做的原因是以缺陷内区域点为中心的局部区域像素点对于中心点的相对梯度方向具备了旋转不变性。而且使用局部区域的相对梯度方向统计均值作为缺陷内区域像素点的灰度方向,可以有效的降低其它因素(如:噪声)对灰度方向的影响。最终,Pt Ω中所有像素点的灰度方向保证了尺度和旋转不变性,且具有一定的抑制其它因素影响的能力。
对于灰度ft Ω而言,它所在的缺陷内区域是在基准采样图像上进行提取的,因此缺陷内区域像素点的灰度大小已经具备了尺度不变性。同时,灰度的大小本身具有旋转不变性。所以,中所有像素点的灰度大小保证了尺度和旋转不变性。
这里,32等级离散化的过程,就是将值域进行32等分,且每一个等分区间对应离散域[1,2,…,32]中的一个值。本发明实施例假定为灰度大小的最大值,为灰度方向的最大值,则离散化后的缺陷内区域灰度大小和方向的公式如下:
步骤4-5:针对缺陷内区域,统计灰度大小为和灰度方向为的缺陷内区域像素点的个数,并将灰度大小-方向共生矩阵(k3,l3)位置的元素HΩ(k3,l3)设定为此值,构造出32×32灰度大小-方向共生矩阵HΩ;
32等级离散化后的梯度大小和方向具有32个等级,因此k3和l3的范围为[1,2,…,32]。那么,构造的灰度大小-方向共生矩阵HΩ的大小为32×32,且HΩ的元素是对具有相同灰度大小和灰度方向的像素点个数的统计。灰度大小-方向共生矩阵同样采用了灰度共生矩阵的思想,它既能反映出灰度的纹理信息和灰度的统计信息,又能反映出灰度大小和方向的联合统计信息。
步骤4-6:将灰度大小-方向共生矩阵HΩ的元素转为概率形式,即HΩ在(k3,l3)的概率为p(k3,l3),其公式如下:
步骤4-7:根据缺陷内区域的灰度大小-方向共生矩阵HΩ,计算对应的具有尺度、旋转不变性的9个灰度大小-方向共生矩阵的特征,即:大小均值、方向均值、大小方差、方向方差、大小熵、方向熵、混合熵、惯性矩、相关性,具体公式与步骤3-8公式的形式一致。
步骤5:将缺陷边缘的9个梯度大小-方向共生矩阵的特征量、缺陷内区域的9个灰度大小-方向共生矩阵的特征量和缺陷图像的7个HU不变矩算法的特征量结合,生成具有25维特征向量样本的训练库;
经过上述的特征提取步骤,本发明实施例有六种类型训练样本集(裂纹、结疤、孔洞、氧化铁皮、卷边、擦伤),其样本都被转换为25维特征向量样本,这25维的特征量包括缺陷边缘的9个梯度大小-方向共生矩阵的特征量、缺陷内区域的9个灰度大小-方向共生矩阵的特征量和缺陷图像的7个HU不变矩算法的特征量。转换后的具有25维特征向量样本的六个训练库分别为:
步骤6:针对具有25维特征向量样本的训练库,使用K-最近邻和R-最近邻相结合的方法,修剪训练样本集,同时,提取修剪后样本的倍数因子;
步骤6-1:假设训练样本集X=[X1,X2,…,Xl],样本个数为l,每一个样本为25维的特征向量,计算样本距离dab=||Xa-Xb||(a,b=1,2,…,l);
步骤6-2:根据dab和由用户选择的K值,获得每一个样本的K个最近邻样本(K-最近邻法),并确定K-最近邻矩阵对应该样本的元素kab,其公式如下:
K值的大小影响本发明分类算法的效率与精度。K值越小,分类的效率越低,但分类的精度会越高;K值越大,分类的效率越高,但分类的精度会降低。因此,用户选择K值要考虑到分类精度与效率的平衡,可以采用经验值法或穷举验证法。经验值法就是根据实际的经验给出一个K值;穷举验证法是对不同的K值进行一一实验,根据实验的结果,确定一个最优的K值。本发明实施例中K值被设定为7,就是通过穷举验证法来判断分类精度与效率是否平衡而获取的。
关键样本的R-最近邻的所有样本被修剪,可以有效的降低训练输入样本的数量,从而提高训练的效率。但是,样本的修剪造成样本数量的减少,必然会降低训练的精度。为此,提取关键样本的倍数因子值m1,作为样本的权重,以提高精度。
步骤6-7:针对修剪后的样本集XS,重复使用步骤6-5和6-6的方法,直至XS样本集被修剪为空,此时获得所有修剪后的样本Aj和对应的倍数因子矩阵的对角线元素(ρ1)jj,并最终确定针对训练样本集X的被修剪后的样本矩阵A=[A1,A2,…,An1]和具有对角线形式的倍数因子矩阵ρ1,这里j=1,2,…,n1,n1为样本集X被修剪后的样本数。
步骤6-8:使用相同的方法,可获得针对训练样本集Y的被修剪后的样本矩阵B=[B1,B2,…,Bn2]和倍数因子矩阵ρ2,这里n2为样本集Y被修剪后的样本数。
本发明实施例的修剪样本方法使用K-最近邻算法捕获训练样本集的关键样本,而使用R-最近邻算法修剪关键样本R邻域的其它样本。本发明实施例的修剪样本方法能够降低训练输入样本的数量和提高训练的效率,而且能有效的剔除噪声样本。本发明实施例针对修剪后样本提取的倍数因子,可以用于抑制样本数量的减少而引起分类精度降低的不利影响。
步骤7:采用新型的最小二乘孪生支持向量机作为分类器,并利用修剪后样本的倍数因子改进分类器;
最小二乘孪生支持向量机是一款新型的分类器,它结合了孪生支持向量机和近似支持向量机的思想。最小二乘孪生支持向量机既与孪生支持向量机一样获取的是两个最优分类面,
又与近似支持向量机一样求解线性方程来代替求解二次规划问题。因此,最小二乘孪生支持向量机既适合交叉分类面数据集的分类,又适合大规模数据集的分类。
针对被修剪后的训练样本矩阵A=[A1,A2,…,An1]和B=[B1,B2,…,Bn2],利用相应的倍数因子ρ1和ρ2修改最小二乘孪生支持向量机的模型,其修改后的模型如下:
式中,C=[A′ B′]′,C1、C2为惩罚因子,表示核函数,e1,e2表示全1的向量,ξ1、ξ2、η1、η2为误差量,w1、w2、b1、b2为两个模型的最优化解,倍数因子ρ1,ρ2具有对角阵形式。其中,C1、C2、由用户确定。
参数C1、C2确定的方法有多种,如:经验值法、十折交叉验证法、智能搜索法等。经验值法是根据实际的经验给出C1、C2的值;十折交叉验证法是将训练集的数据划分为10分,轮流将其中9份作为训练数据,1份作为测试数据,依次进行试验,根据实验的结果,确定最优的C1、C2值。智能搜索法则利用智能算法对参数进行最优化搜索。本发明实施例采用十折交叉验证法,且与分类器模型训练过程相结合来确定最优的C1、C2值。核函数是支持向量机中关键的一个参数。一股采用线性函数、多项式非线性函数、径向基非线性函数等。本发明实施例采用径向基非线性函数,其形式为:
这里,||·||表示范数,δ为径向基函数的核半径。参数δ值的确定同样采用十折交叉验证法,且与分类器模型训练过程相结合。
使用修剪后的样本集A和B对分类器模型进行训练,提高了训练的效率。倍数因子ρ1和ρ2可提高修剪后样本的权重,从而保证算法的执行精度。
对上述两个模型进行推导,可求得对应的解:
这里,w1、w2、b1、b2为分类器模型的解,它们按照上面的公式对训练集进行训练获得。结合参数C1、C2、δ获取方法,本发明实施例将分类器训练的过程与十折交叉验证法结合,完成所有分类器模型解和参数的确定。
步骤8:将改进的分类器与二叉树结合,建立多类别分类器模型,并利用修剪后的训练样本集对多类别分类器模型进行训练;
改进的最小二乘孪生支持向量机模型适合于二分类问题。但是,带钢表面缺陷的分类属于多类别分类。本发明的实施例包含了六类缺陷类型的分类。二叉树是实现多类分类的一种有效模型,与改进的最小二乘孪生支持向量机结合后的多类别分类器模型如图5所示。可以看出,实施例的多类别分类器模型由5个改进的最小二乘孪生支持向量机模型组成。通过六种缺陷类型训练样本集可训练获得这5个二类分类器的模型。
步骤9:根据基准采样尺寸表,将缺陷测试样本转为基准采样图像,然后提取25维的特征量,并输入到多类别分类器模型,完成缺陷的自动识别。
本发明实施例实现了六种带钢表面缺陷(裂纹、结疤、孔洞、氧化铁皮、卷边、擦伤)的自动分类识别。根据自动识别的分类缺陷信息,带钢产品工作人员对生产问题做出快速的反应,及时改进带钢生产的工艺和调整带钢生产的过程,从而形成有效的带钢产品质量监测与控制。
虽然以上描述了本发明的具体实施方式,但是本领域内的熟练的技术人员应当理解,这些仅是举例说明,可以对这些实施方式做出多种变更或修改,而不背离本发明的原理和实质。本发明的范围仅由所附权利要求书限定。
Claims (8)
1.一种带钢表面缺陷的特征提取与分类方法,其特征在于:包括以下步骤:
步骤1:提取带钢表面缺陷样本数据库的基准采样尺寸表;
步骤2:计算缺陷样本数据库中样本图像的高度-宽度比,并检索基准采样尺寸表,然后按照检索到的尺寸对缺陷图像进行采样,获得基准采样图像;
步骤3:针对基准采样图像的缺陷边缘,构造梯度大小-方向共生矩阵,并根据此矩阵提取9个梯度大小-方向共生矩阵的特征,即大小均值、方向均值、大小方差、方向方差、大小熵、方向熵、混合熵、惯性矩、相关性,使提取到的特征具有尺度不变性和旋转不变性;
步骤4:针对基准采样图像的缺陷内区域,构造灰度大小-方向共生矩阵,并根据此矩阵提取9个灰度大小-方向共生矩阵的特征,使这些特征具有尺度不变性和旋转不变性;
步骤5:将缺陷边缘的9个梯度大小-方向共生矩阵的特征量、缺陷内区域的9个灰度大小-方向共生矩阵的特征量和缺陷图像的7个HU不变矩算法的特征量结合,生成具有25维特征向量样本的训练库;
步骤6:针对具有25维特征向量样本的训练库,使用K-最近邻和R-最近邻相结合的方法,修剪训练样本集,同时,提取修剪后样本的倍数因子;
步骤7:采用最小二乘孪生支持向量机作为分类器,并利用修剪后样本的倍数因子改进分类器;
步骤8:将改进的分类器与二叉树结合,建立多类别分类器模型,并利用修剪后的训练样本集对多类别分类器模型进行训练;
步骤9:根据基准采样尺寸表,将缺陷测试样本转为基准采样图像,然后提取25维的特征量,并输入到多类别分类器模型,完成缺陷的自动识别。
2.根据权利要求1所述的带钢表面缺陷的特征提取与分类方法,其特征在于:步骤1所述的提取带钢表面缺陷样本数据库的基准采样尺寸表,步骤如下:
步骤1-1:对带钢表面缺陷样本数据库的所有样本进行[0°,360°)角度范围内的旋转,并计算所有样本对应的每一个旋转角度下的图像尺寸;
步骤1-2:计算所有样本对应的图像尺寸的高度-宽度比rhw,底标hw是高度-宽度的英文缩写;
步骤1-3:分析每一个rhw对应的所有图像尺寸,确定一个基准采样尺寸Shw;
步骤1-4:建立rhw与Shw一一对应的基准采样尺寸表。
3.根据权利要求1所述的带钢表面缺陷的特征提取与分类方法,其特征在于:步骤3所述的提取9个梯度大小-方向共生矩阵特征的步骤如下:
步骤3-1:假设基准采样图像中缺陷边缘有nL个像素点,每一个像素点的灰度值为fi L,以fi L像素点为中心,提取3×3窗口中属于缺陷边缘的像素点,组成集合Pi L,该集合中像素点的个数为且有这里L符号代表缺陷边缘,i=1,2,…,nL;
步骤3-2:计算Pi L中每一个像素点梯度的大小(GM)i(k0,l0)和方向(GO)i(k0,l0),并定义Pi L中所有像素点相对于3×3窗口中心点的坐标偏移量集合为Wi L,这里k0和l0表示Wi L中的坐标偏移量,G符号代表梯度,M符号代表大小,O符号代表方向;
步骤3-3:缺陷边缘的每一个像素点的梯度大小被定义为Pi L中所有像素点梯度大小的均值,其公式如下:
步骤3-4:为保证旋转不变性,将fi L像素点的梯度方向旋转到0方向,对应的Pi L中所有像素点的梯度方向也随之旋转,最终,缺陷边缘的每一个像素点的梯度方向被定义为Pi L中所有像素点旋转后的梯度方向均值,假设MOD符号表示取余运算,则的公式如下:
步骤3-5:将和的值域进行32等级离散化,且离散化后缺陷边缘像素点的梯度大小和方向在离散域[1,2,…,32]中取值,这里i1,i2=1,2,…,nL,M32符号代表32等级离散化后的大小,O32符号代表32等级离散化后的方向;
步骤3-6:针对缺陷边缘,统计梯度大小为和梯度方向为的缺陷边缘像素点的个数,并将梯度大小-方向共生矩阵HL中k1行l1列的元素HL(k1,l1)设定为此值,构造出32×32梯度大小-方向共生矩阵HL;
步骤3-7:将梯度大小-方向共生矩阵HL的元素转为概率形式,即HL在(k1,l1)元素的概率为p(k1,l1),其公式如下:
步骤3-8:根据缺陷边缘的梯度大小-方向共生矩阵HL,计算对应的具有尺度、旋转不变性的9个梯度大小-方向共生矩阵的特征,即:大小均值μM、方向均值μO、大小方差方向方差大小熵EM、方向熵EO、混合熵EMO、惯性矩IMO、相关性CMO,具体公式如下:
4.根据权利要求1所述的带钢表面缺陷的特征提取与分类方法,其特征在于:步骤4所述的提取9个灰度大小-方向共生矩阵特征的步骤如下:
步骤4-1:假设基准采样图像中缺陷内区域有nΩ个像素点,每一个像素点灰度值为ft Ω,以ft Ω像素点为中心,提取3×3窗口中所有的像素点,组成集合Pt Ω,这里Ω符号代表缺陷内区域,t=1,2,…,nΩ;
步骤4-2:计算Pt Ω中像素点梯度的方向(GO)t(k2,l2),这里k2和l2表示Pt Ω中像素点相对于3×3窗口中心点的坐标偏移量,且k2,l2∈{-1,0,1};
步骤4-3:为保证旋转不变性,将ft Ω像素点的梯度方向旋转到0方向,对应的Pt Ω中所有像素点的梯度方向也随之旋转,最终,缺陷内区域的每一个像素点的灰度方向被定义为Pt Ω中所有像素点旋转后的梯度方向均值,其公式如下:
步骤4-5:针对缺陷内区域,统计灰度大小为和灰度方向为的缺陷内区域像素点的个数,并将灰度大小-方向共生矩阵(k3,l3)位置的元素HΩ(k3,l3)设定为此值,构造出32×32灰度大小-方向共生矩阵HΩ;
步骤4-6:将灰度大小-方向共生矩阵HΩ的元素转为概率形式,即HΩ在(k3,l3)的概率为p(k3,l3),其公式如下:
步骤4-7:根据缺陷内区域的灰度大小-方向共生矩阵HΩ,计算对应的具有尺度、旋转不变性的9个灰度大小-方向共生矩阵的特征,即:大小均值、方向均值、大小方差、方向方差、大小熵、方向熵、混合熵、惯性矩、相关性。
5.根据权利要求1所述的带钢表面缺陷的特征提取与分类方法,其特征在于:步骤6所述的修剪训练样本集和提取修剪后样本的倍数因子的步骤如下:
步骤6-1:假设训练样本集X=[X1,X2,…,Xl],X1、X2、……、Xl为训练样本集X中的样本,l为样本个数,每一个样本为25维的特征向量,计算样本距离dab=||Xa-Xb||,且有a,b=1,2,…,l;
步骤6-2:根据dab和由用户选择的K值,获得每一个样本的K个最近邻样本,并确定该样本对应的K-最近邻矩阵的元素kab,其公式如下:
步骤6-3:根据kab,确定每一个样本Xb成为K-最近邻的频次pb和等效距离假定ONE符号表示对非零元素置1运算,则pb和的公式如下:
步骤6-7:针对修剪后的样本集XS,重复使用步骤6-5和6-6的方法,直至XS样本集被修剪为空,此时获得所有修剪后的样本Aj和对应的倍数因子矩阵的对角线元素(ρ1)jj,并最终确定针对训练样本集X的被修剪后的样本矩阵A=[A1,A2,…,An1]和具有对角线形式的倍数因子矩阵ρ1,这里j=1,2,…,n1,n1为样本集X被修剪后的样本数;
步骤6-8:使用相同的方法,可获得针对训练样本集Y的被修剪后的样本矩阵B=[B1,B2,…,Bn2]和倍数因子矩阵ρ2,这里n2为样本集Y被修剪后的样本数。
6.根据权利要求1所述的带钢表面缺陷的特征提取与分类方法,其特征在于:步骤7所述的利用修剪后样本的倍数因子改进最小二乘孪生支持向量机分类器的过程如下:
针对被修剪后的训练样本矩阵A=[A1,A2,…,An1]和B=[B1,B2,…,Bn2],利用倍数因子ρ1和ρ2修改最小二乘孪生支持向量机的模型,其修改后的模型如下:
式中,C=[A′B′]′,C′是C的转置,C1、C2为惩罚因子,表示核函数,e1,e2表示全1的向量,ξ1、ξ2、η1、η2为误差量,w1、w2、b1、b2为两个模型的最优化解,倍数因子ρ1,ρ2具有对角阵形式,其中,C1、C2、由用户确定;
对上述两个模型进行推导,求得对应的解:
式中,
7.根据权利要求1所述的带钢表面缺陷的特征提取与分类方法,其特征在于:步骤1所述的带钢表面缺陷样本数据库,是从现场的带钢表面缺陷检测***中提取的样本,且包含六种类型缺陷样本集,它们是:裂纹、结疤、孔洞、氧化铁皮、卷边、擦伤。
8.根据权利要求5所述的带钢表面缺陷的特征提取与分类方法,其特征在于:步骤6-1所述的训练样本集X和6-8所述的训练样本集Y取自缺陷样本集中的任意两个。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410038260.8A CN103745234B (zh) | 2014-01-23 | 2014-01-23 | 一种带钢表面缺陷的特征提取与分类方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410038260.8A CN103745234B (zh) | 2014-01-23 | 2014-01-23 | 一种带钢表面缺陷的特征提取与分类方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103745234A true CN103745234A (zh) | 2014-04-23 |
CN103745234B CN103745234B (zh) | 2017-01-25 |
Family
ID=50502251
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410038260.8A Expired - Fee Related CN103745234B (zh) | 2014-01-23 | 2014-01-23 | 一种带钢表面缺陷的特征提取与分类方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103745234B (zh) |
Cited By (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104866862A (zh) * | 2015-04-27 | 2015-08-26 | 中南大学 | 一种带钢表面面积型缺陷识别分类的方法 |
CN105447524A (zh) * | 2015-11-27 | 2016-03-30 | 优渊商贸(上海)有限公司 | 图像识别方法及装置 |
CN105654121A (zh) * | 2016-04-11 | 2016-06-08 | 李云栋 | 一种基于深度学习的复杂提花织物缺陷检测方法 |
CN105843212A (zh) * | 2016-03-29 | 2016-08-10 | 东北大学 | 一种高炉故障诊断***及方法 |
CN106097360A (zh) * | 2016-06-17 | 2016-11-09 | 中南大学 | 一种带钢表面缺陷识别方法及装置 |
CN106127756A (zh) * | 2016-06-21 | 2016-11-16 | 西安工程大学 | 一种基于多特征信息融合技术的绝缘子识别检测方法 |
CN106770675A (zh) * | 2016-12-06 | 2017-05-31 | 郑州轻工业学院 | 基于声发射信号的金刚石压机顶锤裂纹在线检测方法 |
CN107784323A (zh) * | 2017-10-16 | 2018-03-09 | 武汉科技大学 | 基于局部不相关约束的带钢表面缺陷图像分类方法 |
CN107992894A (zh) * | 2017-12-12 | 2018-05-04 | 北京小米移动软件有限公司 | 图像识别方法、装置及计算机可读存储介质 |
CN108510977A (zh) * | 2018-03-21 | 2018-09-07 | 清华大学 | 语种识别方法及计算机设备 |
CN108846831A (zh) * | 2018-05-28 | 2018-11-20 | 中冶南方工程技术有限公司 | 基于统计特征和图像特征相结合的带钢表面缺陷分类方法 |
CN109767545A (zh) * | 2017-01-10 | 2019-05-17 | 中国人民银行印制科学技术研究所 | 有价票据的缺陷分类方法及缺陷分类*** |
CN109816645A (zh) * | 2019-01-18 | 2019-05-28 | 创新奇智(广州)科技有限公司 | 一种钢卷松卷的自动检测方法 |
CN109827971A (zh) * | 2019-03-19 | 2019-05-31 | 湖州灵粮生态农业有限公司 | 一种无损检测水果表面缺陷的方法 |
CN110516733A (zh) * | 2019-08-23 | 2019-11-29 | 西南石油大学 | 一种基于改进多分类孪生支持向量机的测井岩性识别方法 |
CN110779928A (zh) * | 2019-11-19 | 2020-02-11 | 汪科道 | 缺陷检测装置及方法 |
CN111079832A (zh) * | 2019-12-13 | 2020-04-28 | 辽宁科技大学 | 具有抗特征噪声性能的钢板表面缺陷分类方法 |
CN111402236A (zh) * | 2020-03-17 | 2020-07-10 | 北京科技大学 | 一种基于图像灰度值的热轧带钢表面缺陷分级方法 |
CN111618265A (zh) * | 2020-05-06 | 2020-09-04 | 大连理工大学 | 一种基于k近邻分类的连铸坯纵裂纹在线检测方法 |
CN112116587A (zh) * | 2020-09-29 | 2020-12-22 | 西安热工研究院有限公司 | 基于孪生支持向量机的水轮机转轮叶片裂纹识别方法、***、设备及存储介质 |
CN112200806A (zh) * | 2020-12-03 | 2021-01-08 | 晶芯成(北京)科技有限公司 | 一种晶圆的缺陷分析方法和*** |
CN113962997A (zh) * | 2021-12-22 | 2022-01-21 | 武汉同顺源钢构钢模有限公司 | 基于图像处理的带钢边裂缺陷检测方法及*** |
CN114049354A (zh) * | 2022-01-12 | 2022-02-15 | 山东仲良格环保技术有限公司 | 一种基于金属锈蚀程度的除锈剂优化配比方法及*** |
CN114240942A (zh) * | 2022-02-25 | 2022-03-25 | 武汉市雷成龙液压机械有限公司 | 一种基于图像处理的金属零件表面质量控制方法及*** |
CN115018833A (zh) * | 2022-08-05 | 2022-09-06 | 山东鲁芯之光半导体制造有限公司 | 一种半导体器件的加工缺陷检测方法 |
CN115049668A (zh) * | 2022-08-16 | 2022-09-13 | 江苏众联管业有限公司 | 基于特征抽取的钢带辊印识别方法 |
CN115375676A (zh) * | 2022-10-24 | 2022-11-22 | 微山三利特不锈钢有限公司 | 一种基于图像识别的不锈钢制品质量检测方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101670533A (zh) * | 2009-09-25 | 2010-03-17 | 南京信息工程大学 | 基于工件加工表面图像分析的刀具磨损状态评价方法 |
CN101852768B (zh) * | 2010-05-05 | 2012-02-22 | 电子科技大学 | 磁粉探伤环境下基于复合特征的工件伤痕识别方法 |
CN102253049B (zh) * | 2011-06-30 | 2013-07-10 | 东北大学 | 带钢生产过程表面质量在线精准检测方法 |
-
2014
- 2014-01-23 CN CN201410038260.8A patent/CN103745234B/zh not_active Expired - Fee Related
Cited By (42)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104866862B (zh) * | 2015-04-27 | 2017-12-22 | 中南大学 | 一种带钢表面面积型缺陷识别分类的方法 |
CN104866862A (zh) * | 2015-04-27 | 2015-08-26 | 中南大学 | 一种带钢表面面积型缺陷识别分类的方法 |
CN105447524A (zh) * | 2015-11-27 | 2016-03-30 | 优渊商贸(上海)有限公司 | 图像识别方法及装置 |
CN105843212A (zh) * | 2016-03-29 | 2016-08-10 | 东北大学 | 一种高炉故障诊断***及方法 |
CN105654121A (zh) * | 2016-04-11 | 2016-06-08 | 李云栋 | 一种基于深度学习的复杂提花织物缺陷检测方法 |
CN105654121B (zh) * | 2016-04-11 | 2019-11-19 | 李云栋 | 一种基于深度学习的复杂提花织物缺陷检测方法 |
CN106097360A (zh) * | 2016-06-17 | 2016-11-09 | 中南大学 | 一种带钢表面缺陷识别方法及装置 |
CN106127756B (zh) * | 2016-06-21 | 2019-03-26 | 西安工程大学 | 一种基于多特征信息融合技术的绝缘子识别检测方法 |
CN106127756A (zh) * | 2016-06-21 | 2016-11-16 | 西安工程大学 | 一种基于多特征信息融合技术的绝缘子识别检测方法 |
CN106770675B (zh) * | 2016-12-06 | 2019-05-21 | 郑州轻工业学院 | 基于声发射信号的金刚石压机顶锤裂纹在线检测方法 |
CN106770675A (zh) * | 2016-12-06 | 2017-05-31 | 郑州轻工业学院 | 基于声发射信号的金刚石压机顶锤裂纹在线检测方法 |
CN109767545B (zh) * | 2017-01-10 | 2021-06-08 | 中钞印制技术研究院有限公司 | 有价票据的缺陷分类方法及缺陷分类*** |
CN109767545A (zh) * | 2017-01-10 | 2019-05-17 | 中国人民银行印制科学技术研究所 | 有价票据的缺陷分类方法及缺陷分类*** |
CN107784323A (zh) * | 2017-10-16 | 2018-03-09 | 武汉科技大学 | 基于局部不相关约束的带钢表面缺陷图像分类方法 |
CN107784323B (zh) * | 2017-10-16 | 2019-11-08 | 武汉科技大学 | 基于局部不相关约束的带钢表面缺陷图像分类方法 |
CN107992894A (zh) * | 2017-12-12 | 2018-05-04 | 北京小米移动软件有限公司 | 图像识别方法、装置及计算机可读存储介质 |
CN107992894B (zh) * | 2017-12-12 | 2022-02-08 | 北京小米移动软件有限公司 | 图像识别方法、装置及计算机可读存储介质 |
CN108510977A (zh) * | 2018-03-21 | 2018-09-07 | 清华大学 | 语种识别方法及计算机设备 |
CN108510977B (zh) * | 2018-03-21 | 2020-05-22 | 清华大学 | 语种识别方法及计算机设备 |
CN108846831A (zh) * | 2018-05-28 | 2018-11-20 | 中冶南方工程技术有限公司 | 基于统计特征和图像特征相结合的带钢表面缺陷分类方法 |
CN108846831B (zh) * | 2018-05-28 | 2021-09-28 | 中冶南方工程技术有限公司 | 基于统计特征和图像特征相结合的带钢表面缺陷分类方法 |
CN109816645A (zh) * | 2019-01-18 | 2019-05-28 | 创新奇智(广州)科技有限公司 | 一种钢卷松卷的自动检测方法 |
CN109827971A (zh) * | 2019-03-19 | 2019-05-31 | 湖州灵粮生态农业有限公司 | 一种无损检测水果表面缺陷的方法 |
CN110516733A (zh) * | 2019-08-23 | 2019-11-29 | 西南石油大学 | 一种基于改进多分类孪生支持向量机的测井岩性识别方法 |
CN110779928A (zh) * | 2019-11-19 | 2020-02-11 | 汪科道 | 缺陷检测装置及方法 |
CN110779928B (zh) * | 2019-11-19 | 2022-07-26 | 汪科道 | 缺陷检测装置及方法 |
CN111079832B (zh) * | 2019-12-13 | 2023-04-18 | 辽宁科技大学 | 具有抗特征噪声性能的钢板表面缺陷分类方法 |
CN111079832A (zh) * | 2019-12-13 | 2020-04-28 | 辽宁科技大学 | 具有抗特征噪声性能的钢板表面缺陷分类方法 |
CN111402236A (zh) * | 2020-03-17 | 2020-07-10 | 北京科技大学 | 一种基于图像灰度值的热轧带钢表面缺陷分级方法 |
CN111402236B (zh) * | 2020-03-17 | 2023-04-28 | 北京科技大学 | 一种基于图像灰度值的热轧带钢表面缺陷分级方法 |
CN111618265A (zh) * | 2020-05-06 | 2020-09-04 | 大连理工大学 | 一种基于k近邻分类的连铸坯纵裂纹在线检测方法 |
CN111618265B (zh) * | 2020-05-06 | 2021-07-16 | 大连理工大学 | 一种基于k近邻分类的连铸坯纵裂纹在线检测方法 |
CN112116587A (zh) * | 2020-09-29 | 2020-12-22 | 西安热工研究院有限公司 | 基于孪生支持向量机的水轮机转轮叶片裂纹识别方法、***、设备及存储介质 |
CN112200806A (zh) * | 2020-12-03 | 2021-01-08 | 晶芯成(北京)科技有限公司 | 一种晶圆的缺陷分析方法和*** |
CN113962997A (zh) * | 2021-12-22 | 2022-01-21 | 武汉同顺源钢构钢模有限公司 | 基于图像处理的带钢边裂缺陷检测方法及*** |
CN114049354A (zh) * | 2022-01-12 | 2022-02-15 | 山东仲良格环保技术有限公司 | 一种基于金属锈蚀程度的除锈剂优化配比方法及*** |
CN114240942A (zh) * | 2022-02-25 | 2022-03-25 | 武汉市雷成龙液压机械有限公司 | 一种基于图像处理的金属零件表面质量控制方法及*** |
CN114240942B (zh) * | 2022-02-25 | 2022-04-26 | 武汉市雷成龙液压机械有限公司 | 一种基于图像处理的金属零件表面质量控制方法及*** |
CN115018833B (zh) * | 2022-08-05 | 2022-11-04 | 山东鲁芯之光半导体制造有限公司 | 一种半导体器件的加工缺陷检测方法 |
CN115018833A (zh) * | 2022-08-05 | 2022-09-06 | 山东鲁芯之光半导体制造有限公司 | 一种半导体器件的加工缺陷检测方法 |
CN115049668A (zh) * | 2022-08-16 | 2022-09-13 | 江苏众联管业有限公司 | 基于特征抽取的钢带辊印识别方法 |
CN115375676A (zh) * | 2022-10-24 | 2022-11-22 | 微山三利特不锈钢有限公司 | 一种基于图像识别的不锈钢制品质量检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103745234B (zh) | 2017-01-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103745234A (zh) | 一种带钢表面缺陷的特征提取与分类方法 | |
CN103996018B (zh) | 基于4dlbp的人脸识别方法 | |
CN103530590A (zh) | Dpm二维码识别*** | |
CN101807258B (zh) | 基于核标度切维数约简的合成孔径雷达图像目标识别方法 | |
CN104156734A (zh) | 一种基于随机蕨分类器的全自主在线学习方法 | |
CN104680144A (zh) | 基于投影极速学习机的唇语识别方法和装置 | |
CN103810704A (zh) | 基于支持向量机和判别随机场的sar图像变化检测方法 | |
CN103246894B (zh) | 一种解决光照不敏感问题的地基云图识别方法 | |
CN103544499A (zh) | 一种基于机器视觉的表面瑕疵检测的纹理特征降维方法 | |
CN104063713A (zh) | 一种基于随机蕨分类器的半自主在线学习方法 | |
CN113160139A (zh) | 基于注意力的Faster R-CNN网络的钢板表面缺陷检测方法 | |
CN101561865A (zh) | 基于多参数谱特征的合成孔径雷达图像目标识别方法 | |
CN114723994B (zh) | 一种基于双分类器对抗增强网络的高光谱图像分类方法 | |
CN106780503A (zh) | 基于后验概率信息熵的遥感图像最优分割尺度确定方法 | |
CN112766194A (zh) | 一种海洋中尺度涡检测方法 | |
CN104616319A (zh) | 基于支持向量机的多特征选择目标跟踪方法 | |
CN116883393B (zh) | 一种基于无锚框目标检测算法的金属表面缺陷检测方法 | |
CN103500343A (zh) | 基于mnf变换结合扩展属性滤波的高光谱影像分类方法 | |
CN103426004A (zh) | 基于纠错输出编码的车型识别方法 | |
CN115661072A (zh) | 一种基于改进Faster RCNN算法的圆盘耙片表面缺陷检测方法 | |
CN114758132B (zh) | 一种基于卷积神经网络的果树病虫害识别方法及*** | |
CN104809471A (zh) | 一种基于空间光谱信息的高光谱图像残差融合分类方法 | |
CN107909073A (zh) | 多维局部二值模式和机器学习手写乐谱谱线删除方法 | |
CN105894035B (zh) | 基于sar-sift和dbn的sar图像分类方法 | |
CN104504391A (zh) | 一种基于稀疏特征和马尔科夫随机场的高光谱图像分类方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170125 Termination date: 20220123 |
|
CF01 | Termination of patent right due to non-payment of annual fee |