一种基于机器视觉的表面瑕疵检测的纹理特征降维方法
【专利摘要】一种基于机器视觉的表面瑕疵检测的纹理特征降维方法，去除训练样本中的噪声样本和安全样本，以边界样本代替随机选取样本作为特征权值迭代时所用的样本集合，特征权值计算时，根据重要程度的差异赋予三个最近邻样本不同的样本权值，增强了对类别相关性强的特征选择的针对性，降低噪声的干扰程度，具有更强的适应性；求取相关系数矩阵，设定自适应阈值，剔除冗余特征，提取了与类别相关性强的特征，在保证分类识别准确率的同时，降低了特征维数，大大提高了在线特征提取及分类识别的速度，解决了高维特征可能含有冗余特征，甚至是噪声特征，从而导致的在线特征提取时间长、预测准确率降低的问题。
【专利说明】—种基于机器视觉的表面瑕疵检测的纹理特征降维方法
【技术领域】
[0001]本发明涉及基于机器视觉的工业图像识别分类方法，具体是指一种工业产品图像的纹理特征降维方法，属于机器视觉表面瑕疵自动检测研究【技术领域】。
【背景技术】
[0002]基于机器视觉的自动检测技术与人工视觉相比，具有速度快、精度高、永不疲劳的优点，在工业生产流水线上，正逐步代替人工视觉检测，这不仅降低人力成本而且可实现对产品质量的严格控制。
[0003]特征选取是工件瑕疵识别分类的关键步骤之一，所选取特征的好坏在很大程度上影响着最终识别分类的效果。为了对样本进行较为全面的描述，通常采用多维的复合特征而不是单一的特征来表示。高维的特征数据可能存在一定相关性，含有许多冗余特征，甚至是噪声特征，对识别正确率产生负面影响。此外，工业现场对实时性要求较高，特征维数过大，会导致运算量变大，学习时间变长，不利于对工件快速识别分类。如何选取与类别相关性强的特征，同时又能尽可能得降低图像特征的维数，减少在线预测时间，显得至关重要。
[0004]特征降维技术可以降低特征维数，提高分类器分类性能。特征降维可以分为特征提取和特征选择。近年来有不少经典的特征提取方法被提出，如主成分分析(PCA)、独立成分分析(ICA)、核主成分分析(KPCA)等。经过这些方法处理后，虽然特征维数下降了，但是新的特征失去了原有的物理意义，由此带来的一个问题是:在线检测时仍然需要计算降维前的所有特征值，从而使得在线预测计算量大大提高，时间消耗大。有研究者提出用于两类样本特征降维的Re I i ef算法和与之对应的多类样本特征降维的Re I i efF方法，这类算法虽然能找到与类别相关性强的特征子集且不改变特征的原有物理意义，但是仍然存在一些不足之处，如:未考虑不同样本的差异性、无法去除冗余特征、容易受到噪声的干扰适应性能较差等。
【发明内容】
[0005]本发明为了解决高维特征可能含有的冗余特征，甚至是噪声特征，从而导致的在线特征提取时间长、预测准确率降低的问题，提供一种适用于纹理工件表面瑕疵检测的离线特征降维方法。
[0006]为达此目的，本发明通过如下技术方案实现:整个流程可分为离线训练和在线预测两个过程，具体如下:
[0007]离线训练过程:
[0008](I)构造Gabor滤波器组，对每个训练图像利用Gabor滤波器组进行卷积变换获得Gabor子图，并对相同尺度不同方向的子图进行融合。存储Gabor滤波器组G，用于在线特征提取过程。
[0009](2)提取融合后不同尺度的图像的灰度均值、方差和灰度共生矩阵的对比度、相关性、信息熵等统计量作为特征，组成训练样本特征数据矩阵。[0010](3)对提取特征进行降维。
[0011](3.1)找出与类别相关性强的特征。去除噪声点，寻找所有训练样本中的边界样本，根据所得边界样本数据对每一维特征权值进行迭代，权值越大表示其对应的特征与类别相关性越强。比较每一维特征权值与阈值THl的大小:若权值小于THl，则将该维特征所对应的状态标志位置O,否则置I。由此得到所有特征状态标志向量markl。
[0012](3.2)去除冗余特征。将训练样本特征数据矩阵按照权值大小降序排列，计算相关系数矩阵RH0，并根据相关系数矩阵计算阈值TH2。找出与最大权值特征相关系数最大的特征，并比较该相关系数与TH2的大小。若大于TH2，则将该维特征删除，并将其对应状态标志位置O;否则置I。再找出与次最大特征相关系数最大的特征，比较与TH2的大小。以此类推，直到把权值最小的冗余特征去除。由此得到所有特征状态标志向量mark2。
[0013](3.3)将markl和mark2按位与运算得到mark即为最终所有特征的状态标志向量。将mark中状态标志为O对应的特征删除，保留标志位为I的特征。存储特征状态标志向量mark，用于在线特征提取过程。
[0014](4)根据降维后所得的训练样本特征矩阵，训练分类器模型，并存储分类器模型系数model。本发明采用分类精度高、实时性能好、推广能力强的最小二乘支持向量机作为分类器。
[0015]在线预测过程:
[0016](I)下载Gabor滤波器组G，对预测图像进行卷积变换获得Gabor子图，并对相同尺度不同方向的子图进行融合。
[0017](2)下载特征状态标志向量mark并判断每一维特征的标志位是O或1:如果是0，则不进行该维特征值的计算，否则计算该维特征值。
[0018](3)下载离线训练所得分类器模型model，利用所提取的预测图像特征，对表面缺陷进行识别分类。
[0019]本发明的有益效果是:本发明提供了 一种基于机器视觉的表面瑕疵检测的纹理特征降维方法，去除训练样本噪声点和安全样本点，以边界样本代替随机选取样本作为特征权值迭代时所用的样本集合，特征权值计算时，以k(k取3)个最近邻代替一个最近邻且根据三个最近邻的重要程度赋予不同的样本权值，增强样本选择的针对性，减少所受噪声的干扰，适应性更好；此外弥补了 Relief等算法无法剔除冗余特征的不足，提取了与类别相关性强的特征，在保证分类识别准确率的同时，降低了特征维数，大大提高了在线特征提取及分类识别的速度，解决了高维特征可能含有许多冗余特征，甚至是噪声特征，从而导致的在线特征提取时间长、预测准确率降低的问题。
【专利附图】
【附图说明】
[0020]图1是本发明整体流程。
[0021]图2是离线训练纹理特征降维流程图
【具体实施方式】
[0022]为使本发明的目的、技术方案和优点更加清楚明白，下面结合具体实施例，并参照附图，对本发明作进一步详细说明。[0023]本发明一种基于机器视觉的表面瑕疵检测的纹理特征降维方法，整个流程可分为离线训练和在线预测两个过程。如图1所示，离线训练部分主要由Gabor小波变换、图像融合、特征提取和特征降维和分类器学习5个子部分组成；在线预测过程与离线过程基本相同，但只需要下载并利用离线训练时所得的Gabor滤波器组系数G、每维特征状态标志向量mark和分类器模型model。
[0024]进一步的，离线训练过程具体实现步骤为:
[0025]步骤一、对图像进行Gabor变换。
[0026](I)读取工业智能相机拍照所得工件表面图像η张，其中有瑕疵工件和合格工件图像数量大致相等。
[0027](2)构造二维Gabor滤波器组。Gabor小波变换是一种非常规范的图像处理算法。采用具有 4 个方向(Θ = [O。，45°，90° ,135° ])和 5 个尺度(λ = [0,1,2,3,4])的Gabor滤波器组。存储Gabor滤波器组G，用于在线特征提取过程。分别用不同方向不同尺度的Gabor滤波器与图像进行卷积，获得20副Gabor分量子图。
[0028](3)图像融合。以f (v，u)表示原图像经过Gabor小波变换后获得的20副分量子图，其中ve {0,1,2,3,4}表示5个不同尺度，u e {1,2,3,4}表示4个不同方向。将相同尺度不同方向的子图按照下式进行融合:
【权利要求】
1.一种基于机器视觉的表面瑕疵检测的纹理特征降维方法，其特征是，提取出与类别相关性强的特征，去除冗余特征，在保证分类识别准确率的同时，降低了特征维数，大大提高了在线特征提取及分类识别的速度，解决了高维特征可能含有冗余特征，甚至是噪声特征导致的在线特征提取时间长、预测准确率降低的问题；包含如下几个步骤: (1)找出与类别相关性强的特征，得到所有特征状态标志向量markl； (2)去除冗余特征，计算相关系数矩阵RHO，进一步得到所有特征状态标志向量mark2； (3)将markl和mark2按位与运算得到mark即为最终所有特征的状态标志向量mark,存储状态标志向量mark，用于在线特征提取过程。
2.根据权利要求1所述基于机器视觉的表面瑕疵检测的纹理特征降维方法，其特征是:所述步骤(1)中找出与类别相关性强的特征，包含如下步骤: 第一步、剔除安全样本和噪声样本；对于n*m维的训练样本特征矩阵X中的每一个样本Xn，计算出它与其余样本之间的欧式距离；找出与Xn距离最小的5个样本训练，判断这5个训练样本类别标签是否与Xn的类别标签一致，如果5个最近邻样本的类别均与xn的类别相同或者均与Xn不同，则将Xn视为安全样本或者噪声样本剔除；由此获得只含边界样本的特征矩阵Xnew ； 第二步、计算模糊差异度量；对于Xnrat中的每个样本Xn，分别找出其k(k = 3)个最近邻同类、异类样本 :每一维特征对应初始权值
3.根据权利要求1所述基于机器视觉的表面瑕疵检测的纹理特征降维方法，其特征是:所述步骤(2)中去除冗余特征，包含如下步骤: 第一步、将训练样本Xm的所有特征按照特征权值大小降序排列，计算相关系数矩阵L:RHO
4.根据权利要求1所述基于机器视觉的表面瑕疵检测的纹理特征降维方法，其特征是:所述步骤(3),将markl和mark2按位与运算得到mark即为最终所有特征的状态标志向量，将mark中状态标志为O对应的特征删除,保留标志位为I的特征,存储状态标志向量mark,用于在线特征提取过程。
【文档编号】G06K9/46GK103544499SQ201310485150
【公开日】2014年1月29日 申请日期:2013年10月12日 优先权日:2013年10月12日
【发明者】白瑞林, 张振尧, 姜利杰, 李新 申请人:江南大学, 无锡信捷电气股份有限公司