CN106770323A - 基于层次聚类和Gabor滤波的纺织品瑕疵检测方法 - Google Patents
基于层次聚类和Gabor滤波的纺织品瑕疵检测方法 Download PDFInfo
- Publication number
- CN106770323A CN106770323A CN201611156224.7A CN201611156224A CN106770323A CN 106770323 A CN106770323 A CN 106770323A CN 201611156224 A CN201611156224 A CN 201611156224A CN 106770323 A CN106770323 A CN 106770323A
- Authority
- CN
- China
- Prior art keywords
- grid pattern
- class grid
- row
- calculate
- class
- 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
- 239000004753 textile Substances 0.000 title claims abstract description 75
- 238000001514 detection method Methods 0.000 title claims abstract description 18
- 238000000034 method Methods 0.000 claims abstract description 53
- 238000004458 analytical method Methods 0.000 claims abstract description 12
- 238000000605 extraction Methods 0.000 claims abstract description 12
- 230000011218 segmentation Effects 0.000 claims abstract description 11
- 239000011159 matrix material Substances 0.000 claims description 22
- 230000007547 defect Effects 0.000 claims description 14
- 230000002950 deficient Effects 0.000 claims description 14
- 230000006641 stabilisation Effects 0.000 claims description 8
- 238000011105 stabilization Methods 0.000 claims description 8
- 238000006073 displacement reaction Methods 0.000 claims description 6
- 238000010276 construction Methods 0.000 claims description 4
- 238000001914 filtration Methods 0.000 claims description 2
- 102000002274 Matrix Metalloproteinases Human genes 0.000 claims 1
- 108010000684 Matrix Metalloproteinases Proteins 0.000 claims 1
- 238000005286 illumination Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 22
- 210000000988 bone and bone Anatomy 0.000 description 8
- 238000004364 calculation method Methods 0.000 description 8
- 230000001174 ascending effect Effects 0.000 description 7
- 239000004744 fabric Substances 0.000 description 5
- 230000006870 function Effects 0.000 description 5
- 238000012360 testing method Methods 0.000 description 4
- 238000000205 computational method Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000009394 selective breeding Methods 0.000 description 2
- 241000894007 species Species 0.000 description 2
- HUTDUHSNJYTCAR-UHFFFAOYSA-N ancymidol Chemical compound C1=CC(OC)=CC=C1C(O)(C=1C=NC=NC=1)C1CC1 HUTDUHSNJYTCAR-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000012217 deletion Methods 0.000 description 1
- 230000037430 deletion Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000012854 evaluation process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 210000003746 feather Anatomy 0.000 description 1
- 244000144992 flock Species 0.000 description 1
- 230000015654 memory Effects 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 230000007115 recruitment Effects 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 235000015170 shellfish Nutrition 0.000 description 1
- 238000009987 spinning Methods 0.000 description 1
- 230000017105 transposition Effects 0.000 description 1
- 238000005303 weighing Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/84—Systems specially adapted for particular applications
- G01N21/88—Investigating the presence of flaws or contamination
- G01N21/8851—Scan or image signal processing specially adapted therefor, e.g. for scan signal adjustment, for detecting different kinds of defects, for compensating for structures, markings, edges
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/84—Systems specially adapted for particular applications
- G01N21/88—Investigating the presence of flaws or contamination
- G01N21/8851—Scan or image signal processing specially adapted therefor, e.g. for scan signal adjustment, for detecting different kinds of defects, for compensating for structures, markings, edges
- G01N2021/8887—Scan or image signal processing specially adapted therefor, e.g. for scan signal adjustment, for detecting different kinds of defects, for compensating for structures, markings, edges based on image processing techniques
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Landscapes
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Image Analysis (AREA)
Abstract
本发明提供一种基于层次聚类和Gabor滤波的纺织品瑕疵检测方法，通过结合层次聚类算法和Gabor滤波器组，分析基于照明光源下平坦纺织品表面的数字图像像素灰度信息，自动定位纺织品表面瑕疵，主要步骤包括类格图案分割，特征提取和特征比较三部分。本发明特别适用于自动识别在稳定照明光源下采集的纺织品平坦表面的数字图像中的纺织品表面瑕疵，是一种基于层次聚类算法的从纺织品图像中自动分割类格图案，并基于Gabor滤波器组对类格图案进行特征提取与瑕疵识别的方法。
Description
技术领域
本发明涉及一种基于层次聚类和Gabor滤波的纺织品瑕疵检测方法。
背景技术
传统的纺织品瑕疵人工识别准确率只有60-75％(K.Srinivasan,P.H.Dastoor,P.Radhakrishnaiah,et al..FDAS:a knowledge-based framework for analysis ofdefects in woven textiles tructures,J.Text.Inst.83(1992)431–448.)，机器自动识别纺织品瑕疵的方法具有实际应用需求。平坦纺织品表面的数字图像采样(以下简称纺织品图像)属于二维纹理，二维纹理已被证明可根据17种壁纸群(wallpaper group)定义的图案排列方法生成(K.Srinivasan,P.H.Dastoor,P.Radhakrishnaiah,et al..FDAS:aknowledge-based framework for analysis of defects in woven textilestructures,J.Text.Inst.83(1992)431–448.)，用于生成二维纹理的图案称为格(lattice)，格的内部图案称为motif。多数纺织品瑕疵自动检测方法只能处理壁纸群中p1类型的纺织品图像(H.Y.T.Ngan,G.K.H.Pang,N.H.C.Yung.Automated fabric defectdetection—A review,Image and Vision Computing 29(7)(2011)442–458.)，仅有少数方法能处理p1类型以外的纺织品图像，例如基于小波预处理的基准图像差分方法(wavelet-pre-processed golden image subtraction，以下简称WGIS，出自文献：H.Y.T.Ngan,G.K.H.Pang,N.H.C.Yung,et al.,Wavelet based methods on patternedfabric defect detection,PatternRecognit.38(4)(2005)559–576)，共生矩阵方法，布林带方法(Bollinger bands，以下简称BB，出自文献：H.Y.T.Ngan,G.K.H.Pang,Novelmethodforpatternedfabricinspectionusingbollinger bands,Opt.Eng.45(8)(2006)087202-1–087202-15)，规则带方法(regular bands，以下简称RB，出自文献：H.Y.T.Ngan,G.K.H.Pang,Regularityanalysisforpatternedtextureinspection,IEEETrans.Autom.Sci.Eng.6(1)(2009)131–144)，Elo评估方法(Elo rating method，以下简称ER，出自文献C.S.C.Tsang,H.Y.T.Ngan,G.K.H.Pang,Fabric inspection based on theElo rating method,Pattern Recognit.51(2016)378–394.)等。尽管这些方法可以处理p1以外的纺织品图像，但它们的计算方法多是建立在基于人工选择的类似格的图案(以下简称类格图案)之上。例如WGIS要求人工选择类格图案的尺寸和纹理，BB，RB和ER要求人工定义类格图案的尺寸。这些先验知识在一定程度上降低了机器识别纺织品瑕疵的自动化程度。
发明内容
本发明要解决的技术问题是：为了克服现有纺织品瑕疵自动检测方法仍基于人工选择或人工定义，自动化程度不高的不足，本发明提供一种基于层次聚类和Gabor滤波的纺织品瑕疵检测方法，通过结合基于层次聚类(hierarchical clustering,HC)算法和Gabor滤波器组分析基于照明光源下平坦纺织品表面的数字图像像素灰度信息，自动定位纺织品表面瑕疵，本发明特别适用于自动识别在稳定照明光源下采集的纺织品平坦表面的数字图像中的纺织品表面瑕疵。
本发明解决其技术问题所采用的技术方案是：一种基于层次聚类和Gabor滤波的纺织品瑕疵检测方法，包括以下步骤：
输入灰度化的纺织品图像；
类格图案分割：分割纺织品图像产生类格图案，类格图案满足：相对于纺织品图像的行和列，类格图案按图像行的方向横向排列，并按照列的方向纵向排列；在形态成分分析方法的纺织品图案卡通成分Ic中，类格图案具有几何形状并与背景像素在灰度上有显著差异；
特征提取：将类格图案与Gabor滤波器组的滤波进行卷积，对卷积结果降维得到一维投影并计算其能量和振幅，根据能量和振幅计算特征向量；计算一维投影之间的车贝雪夫距离，挑选每行类格图案的典型无瑕疵类格图案，根据典型无瑕疵类格图案的特征向量计算理想特征向量；每个类格图案特征向量与理想特征向量间的车贝雪夫距离为所提取的特征；和
特征比较：分析特征提取步骤中得到的基于特征向量的车贝雪夫距离直方图，获取表示有瑕疵的像素索引集合。
基于层次聚类和Gabor滤波的纺织品瑕疵检测方法，具体步骤为：
步骤1：输入灰度化的纺织品图像，通过形态学成分分析方法计算输入图像的卡通成分Ic；
步骤2：使用阈值fc·max(Ic)二值化卡通成分Ic得到二值图像Itc；
步骤3：计算二值图像Itc中各行背景像素数的多重集各列背景像素数的多重集计算各行背景像素数的多重集和各列背景像素数的多重集的峰值和
步骤4：分别计算基于步骤(3)中的峰值和的HC算法轮廓系数，以最大轮廓系数对应的聚类个数初始化HC算法，并对峰值和分别进行聚类，得到的聚类中心构成多重集和
步骤5：根据多重集和计算阈值和分别计算中不小于的峰值所对应行索引的间距和中不小于的峰值所对应列索引的间距，即用阈值对中的峰值进行筛选，只计算不小于的峰值所对应的行索引间距，用阈值对中的峰值进行筛选，只计算不小于的峰值所对应的列索引间距；计算构成稳定行间距的最长连续行索引的集合Sh以及构成稳定列间距的最长连续列索引的集合Sv；
其中，将中元素按降序排列，从开始，计算其与下一个元素的差的绝对值，若该值不大于DCT尺寸的高，则继续计算当前元素与下一个元素的差的绝对值并与DCT尺寸的高比较，若不大于则继续，若大于则下一个元素为并终止；如果一直没有出现大于DCT尺寸的高的情况，则即为
将中元素按降序排列，从开始，计算其与下一个元素的差的绝对值，若该值不大于DCT尺寸的高，则继续计算当前元素与下一个元素的差的绝对值并与DCT尺寸的宽比较，若不大于则继续，若大于则下一个元素为并终止；如果一直没有出现大于DCT尺寸的宽的情况，则即为
DCT尺寸：形态学成分分析方法在图像局部应用离散余弦变换，即首先将图像划分为不重叠且具有固定大小的矩形区域，然后对每个区域应用离散余弦变换，矩形区域的大小称为DCT尺寸，单位为像素，区域内一行的像素数称为DCT尺寸的宽，一列的像素数称为DCT尺寸的高；
步骤6：扩展集合Sh和集合Sv以覆盖图像的卡通成分Ic的大部分区域，并计算类格图案的理想行数和理想列数
步骤7：根据集合Sh和集合Sv将图像的卡通成分Ic分割为类格图案(其中，)，使用Gabor滤波器组与类格图案卷积并计算一维卷积投影Gi，j|s，θ，根据一维卷积投影Gi，j|s，θ的能量和振幅构建特征向量
其中，i和j分别是类格图案以类格图案为单位的行与列索引；
表示尺度参数，表示角度参数；
步骤8：计算类格图案与第i行类格图案之间基于一维卷积投影Gi，j|s，θ的车贝雪夫距离，结果保存为矩阵的第j行，根据矩阵计算第i行类格图案的典型无瑕疵类格图案
的意思是：对于固定i，s，θ，然后对于特定的j，计算Gi，j|s，θ与第i行所有类格图案一维投影的距离，这些距离保存为的第j行，这一行表示特定的j所对应的类格图案与第i行中其它类格图案基于一维投影的距离，而整个表示第i行所有类格图案一维投影之间的距离。
步骤9：计算第i行类格图案的典型无瑕疵类格图案与第i行类格图案之间基于特征向量的车贝雪夫距离，结果保存为矩阵的第i行，根据矩阵计算第i行的典型无瑕疵类格图案的距离之和di；
步骤8根据得到每行类格图案的步骤9计算所对应特征向量的距离，找出整幅图的理想特征向量步骤9和步骤8的距离计算区别在于步骤8是计算一维投影的距离，而步骤9则是特征向量之间的距离。
步骤10：根据距离之和di筛选第i行的典型无瑕疵类格图案基于筛选结果计算理想特征向量
其中，理想特征向量为S*的均值；
为每行的典型类格图案的特征向量，i为类格图案行索引；
其中表示按i的索引顺序进行连接；di为距离向量d中的第i个元素，表示与所有典型图案的车贝雪夫距离之和，典型图案是指每行的典型无瑕疵类格图案；
步骤11：计算特征向量和理想特征向量的车贝雪夫距离，结果保存为矩阵中索引为(i，j)的元素其中i为类格图案行索引，j为类格图案列索引；
步骤12：计算矩阵的直方图
步骤13：基于距离阈值d*和直方图计算缺口值t′和断崖值t″，根据缺口值t′和断崖值t″近似表示真实阈值t*；
步骤14：所有对应车贝雪夫距离>真实阈值t*的类格图案波标记为有瑕疵的类格图案输出其中i为类格图案行索引，j为类格图案列索引。
本发明的有益效果是，本发明的基于层次聚类和Gabor滤波的纺织品瑕疵检测方法，通过结合基于层次聚类算法和Gabor滤波器组分析基于照明光源下平坦纺织品表面的数字图像像素灰度信息，自动定位纺织品表面瑕疵，本发明特别适用于自动识别在稳定照明光源下采集的纺织品平坦表面的数字图像中的纺织品表面瑕疵，是一种基于层次聚类算法的从纺织品图像中自动分割类格图案的方法，并且针对类格图案的基于Gabor滤波器组的特征提取与瑕疵识别。
附图说明
下面结合附图和实施例对本发明进一步说明。
图1是本发明的基于层次聚类和Gabor滤波的纺织品瑕疵检测方法的算法流程图。
图2是算法1的假设条件展示图。
图3是算法1的概念图。
图4是算法4中投影的概念图。
图5是有瑕疵的类格图案与无瑕疵类格图案在特征向量各元素取值上存在差异的展示图。
图6是典型类格图案的计算过程示意图。
图7是距离矩阵的计算过程和瑕疵识别过程示意图。
图8是箱形图像TPR-FPR散点图：(a)为箱形图像断端散点图；(b)为箱形图像孔洞散点图；(c)为箱形图像网纹散点图；(d)为箱形图像粗条纹散点图；(e)为箱形图像细条纹散点图。
图9是星形图像TPR-FPR散点图：(a)为星形图像断端散点图；(b)为星形图像孔洞散点图；(c)为星形图像网纹散点图；(d)为星形图像粗条纹散点图；(e)为星形图像细条纹散点图。
具体实施方式
现在结合附图对本发明作进一步详细的说明。
为使陈述清楚明了，现集中定义本发明所涉及的部分符号和概念。
1.表示正整数集合。
2.表示包括零的整数集合。
3.表示包括零的正实数集合。
4.表示包括零的实数集合。
5.T表示矩阵或向量转置。
6.表示比小的最大整数，例如
7.表示按操作数顺序连接产生向量，例如标量v1＝1和向量 对于标量s1＝8，s2＝1，s3＝5，对于向量
8.其中
9.Cb(v1，v2)表示维数相同的向量v1与v2的车贝雪夫距离(Chebychev distance)。
10.{ai}表示由索引i确定的由元素ai组成的集合或多重集。
11.|S|表示集合S中的元素个数。
12.avg(S)：计算集合或多重集S的均值，S的元素均为实数。
13.std(S)：计算集合或多重集S的标准差，S的元素均为实数。
14.med(S)：计算集合或多重集S的中位值，S的元素均为实数。
15.mod(S)：计算多重集S的众数，S的元素均为实数。
16.max(S)表示找出集合或多重集S的元素最大值，例如max(Ic)代表Ic中像素的最大灰度值。
17.max(s「条件)表示找出符合条件的的最大值。
18.min(S)表示找出集合或多重集S的元素最小值，例如min(Ic)代表Ic中像素的最小灰度值。
19.arg maxs f(s)表示在函数f的定义域内变量s的取值范围中，使得函数f(s)取最大值的s。
20.arg mins f(s)表示在函数f的定义域内变量s的取值范围中，使得函数f(s)取最小值的s。
21.arg modi({ai})表示对应多重集{ai}众数mod({ai})的索引。
22.dimx(I)表示二维图像I的总行数，dimy(I)表示I的总列数。
23.I(x，y)表示在二维图像I中具有行列索引(x，y)的像素值。行索引1≤x≤dimx(I)；列索引1≤y≤dimy(I)。
24.纺织品图像卡通成分Ic：对一幅灰度化的纺织品图像，应用基于曲波(curvelet)和离散余弦变换(local discrete cosine transform，以下简称DCT)的形态成分分析方法(morphological component analysis，以下简称MCA)所计算的具有光滑边缘图案的图像称为卡通成分Ic，Ic是一幅灰度图像。
25.阈值系数fc：用于二值化Ic的参数，该参数由算法3计算得到。
26.二值化卡通成分Itc：使用fc·max(Ic)作为阈值二值化Ic所得到的二值图像，其中1表示前景像素，即Ic中灰度值不小于阈值的像素，0表示背景像素。Itc与Ic的行数和列数相同。
27.横向投影是的多重集，其中1≤k≤dimy(I)，即表示行索引为x的背景像素个数。
28.纵向投影是的多重集，其中1≤l≤dimx(I)，即表示列索引为y的背景像素个数。
29.表示的峰值多重集，中的元素称为峰值，峰值是指中满足和的元素其中x表示行索引。
30.表示的峰值多重集，中的元素称为峰值，峰值是指中满足和的元素其中y表示列索引。
31.表示对中元素使用HC算法所得到的聚类中心的多重集。
32.表示对中元素使用HC算法所得到的聚类中心的多重集。
33.DCT尺寸：MCA在图像局部应用DCT，即首先将图像划分为不重叠且具有固定大小的矩形区域，然后对每个区域应用DCT，矩形区域的大小称为DCT尺寸，单位为像素，区域内一行的像素数称为DCT尺寸的宽，一列的像素数称为DCT尺寸的高。
34.将中元素按降序排列，从开始，计算其与下一个元素的差的绝对值，若该值不大于DCT尺寸的高，则继续计算当前元素与下一个元素的差的绝对值并与DCT尺寸的高比较，若不大于则继续，若大于则下一个元素为并终止；如果一直没有出现大于DCT尺寸的高的情况，则将即为
35.将中元素按降序排列，从开始，计算其与下一个元素的差的绝对值，若该值不大于DCT尺寸的高，则继续计算当前元素与下一个元素的差的绝对值并与DCT尺寸的宽比较，若不大于则继续，若大于则下一个元素为并终止；如果一直没有出现大于DCT尺寸的宽的情况，则将即为
36.K：HC算法的聚类个数，用于指定HC算法的聚类个数，
37.S′h：中不小于的元素的行索引x的集合，S′h中x按从小到大排列。
38.S′v：中不小于的元素的列索引y的集合，S′v中y按从小到大排列。
39.集合S′h中相邻两元素之间后一元素与前一元素差的多重集，即行索引间距的多重集，中行索引间距按产生该间距的x∈S′h中的较大值从小到大排列。
40.集合S′v中相邻两元素之间后一元素与前一元素差的多重集，即列索引间距的多重集，中列索引间距按产生该间距的y∈S′v中的较大值从小到大排列。
41.稳定行间距：中与之差的绝对值小于的行索引间距。即里的元素是行索引间距，这些行索引间距都与故差，差的绝对值小于的行索引间距称为稳定行间距。
42.稳定列间距：中与之差的绝对值小于的列索引间距。即里的元素是列索引间距，这些列索引间距都与做差，差的绝对值小于的列索引间距称为稳定列间距。
43.满足行索引连续的稳定行间距组成的多重集，中的稳定行间距按产生该间距的两个行索引x∈S′h中的较大值升序排列，行索引连续是指中排列顺序相邻的稳定行间距所对应的行索引xi∈S′h中的i值连续，例如i可以取2，3，4，但不能只取2和4。
44.满足列索引连续的稳定列间距组成的多重集，中的稳定列间距按产生该间距的两个列索引y∈S′v中的较大值升序排列，列索引连续是指中排列顺序相邻的稳定列间距所对应的yj∈S′v的的j值连续，例如j可以取2，3，4，但不能只取2和4。
45.具有最多元素个数的
46.具有最多元素个数的
47.Sh：Sh的初始值为产生中稳定行间距所对应行索引x∈S′h的集合，后经算法1扩展，Sh表示类格图案分界处像素行索引的集合。
48.Sv：Sv的初始值为产生中稳定列间距所对应列索引y∈S′v的集合，后经算法1扩展，Sv表示类格图案分界处像素列索引的集合。
49.类格图案分界：Sh中行索引所对应的Ic中的行与Sv中列索引所对应的Ic中的列。
50.理想行数
51.理想列数
52.类格图案：在Ic中，根据Sh中行索引对应的行与Sv中列索引对应的列，将Ic分割为矩形区域，矩形区域称为类格图案，其中Sh中行索引对应的行与Sv中列索引对应的列不包含在类格图案里。
在以上定义的基础上，现介绍本发明的技术方案。正如附图1所示，本发明方法由三部分构成：类格图案分割，特征提取和特征比较，为了方便表述，这里分别定义这三部分的大体步骤为算法1，算法4和算法6，本发明涉及6个自定义算法(算法1-6)，其他涉及到的方法都是调用现成的程序，比如MCA，HC之类。以下按图1中从上至下的顺序依次介绍三部分的内容。
类格图案分割的功能是自动分割纺织品图像产生类格图案。由于纺织品图案丰富多样，相应的类格图案种类繁多。本发明的类格图案分割方法建立在对类格图案的假设上，即：相对于纺织品图像的行和列，类格图案按图像行的方向横向排列，并按列的方向纵向排列；在MCA的卡通成分Ic中，类格图案具有几何形状并与背景像素在灰度上有显著差异。例如附图2所示的三种情况，图2中每一行显示了一种情况，每行第一列是纺织品图像，第二列是Ic，第三列是Ic的三维Mesh图，第四列是二值化卡通成分每行背景像素的个数分布，第四列图的横坐标是行索引，纵坐标是背景像素个数。图2中第一行纺织品图像的类格图案没有几何形状，这导致了背景像素分布缺乏明显的周期性；第二行纺织品图像的类格图案虽然有几何形状，但类格图案中的形状与背景在Ic中的差异小，即相应的Mesh图大部分区域几乎是平坦的，这导致背景像素数量过多，背景像素分布缺乏明显的周期性；第三行纺织品图像类格图案具有几何形状并在Ic中与背景的差异大，其背景像素分布具有周期性。
如附图3所示，对于一幅给定的纺织品图像，类格图案分割(算法1)使用MCA计算Ic和纹理成分，根据算法3计算得到的阈值系数fc，使用阈值fc·max(Ic)二值化Ic得到Itc。图3中显示了Ic的Mesh图，Ic中的二维图案在Mesh图中显示为三维“山峰”，二值化Ic相当于用一个灰色平面截断山峰，山峰平面上方的部分所对应的像素保存为1，山峰平面下方的部分所对应的像素保存为0，这个二值化结果为Itc，即图3右下角箭头“使用阈值二值化Ic得到Itc”所指的图案。
假设纺织品图像至少由4个类格图案构成，那么Itc中对应类格图案的对象尺寸应小于图像尺寸的一半，因此如果出现了尺寸过大的情况，那么这个对象则不是类格图案，应从Itc中删除尺寸过大的对象，即：由Moore-Neighbor跟踪算法(Moore-Neighbortracingalgorithm)获取Itc中对象的闭合边缘。对于每个具有闭合边缘的对象，找出该对象行列索引的极值，如果该对象行索引极值之差的绝对值超过0.75·dimx(Itc)，或列索引极值之差的绝对值超过0.75·dimy(Itc)，则从Itc中删除该对象，即将尺寸过大的对象的像素置为0。
纺织品图像类格图案的几何形状被Itc中的二值对象所描述，类格图案的丰富多样导致了二值对象几何形状的多样性，但二值对象之间背景像素的分布受其形状影响小，即：不同形状的二值对象，如果其沿相同方向的分布相同，那么背景像素在该方向上的分布相似。如附图3所示，统计二值化卡通成分每行和每列的背景像素个数，按行和列的次序分别排列背景像素个数即构成了背景像素的横向投影和纵向投影和的峰值分别记为多重集和即图3中横向投影和纵向投影图示中靠近标签“横向投影”和“纵向投影”的深色小点，这些峰值反映了Itc中类格图案的分界。为获取这些峰值而过滤其他峰值，对和分别使用HC算法进行聚类，本发明方法中的HC算法采用常见的基于单关联(Singlelinkage)的合并层次聚类(Agglomerative Hierarchical Clustering)方式处理数据，即：对于给定的一维离散实数，HC算法计算每两个数之间的欧氏距离，距离最小的两个数聚为一类；对于没有归类的数，每个数视为单独一类，类与类之间的距离定义为两类之间欧氏距离最小的两数间距离，将距离最小的两类聚为一类，如此反复直至所有数都聚为一类。这个过程可以由树状图(dendrogram)描述，例如附图3中标有“的HC聚类”和“的HC聚类”的图例显示了对峰值和的聚类过程树状图，两幅树状图的横轴表示峰值，纵轴表示类间距离。在树状图中，HC算法采取自底向上的方式两两聚类，被聚为一类的两类由倒置的“U”型线条连接，最上方的U型线条表示所有峰值最终被聚为一类。对于给定的聚类个数K，HC算法返回类数为K时各类的聚类中心(聚类中心定义为该类中所有元素的平均值)，即在图3的树状图中，用平行于横轴的直线与U型线条中平行于纵轴的部分相交，交点个数为K时，自底向上连接到同一交点的峰值视为一类，从而产生K个类。
对和的聚类首先根据基于欧氏距离的轮廓系数估计最优聚类个数K，以具有最大轮廓系数的K初始化HC算法并对扣聚类，即：计算K为的HC算法对分类的轮廓系数，用对应最大轮廓系数的K初始化HC算法并对进行分类，分类所得的聚类中心保存为多重集计算K为的HC算法对分类的轮廓系数，用对应最大轮廓系数的K初始化HC算法并对进行分类，分类所得的聚类中心保存为多重集由于数据的随机性，和有可能包含多个相近的聚类中心，这些相近聚类中心中的最小值被分别选为阈值和即：将降序排列，从开始，计算前后两元素之差的绝对值，差大于DCT尺寸高的的第一个元素为类似地，将降序排列，从开始，计算前后两元素之差的绝对值，差大于DCT尺寸宽的的第一个元素为中不小于的峰值按所对应的行索引，记为S′h；中不小于的峰值按所对应的列索引，记为S′v。
由于瑕疵等因素的干扰，S′h和S′v中的行列索引不一定准确反映类格图案的分界。因此，需要评估S′h中是否存在具有稳定行间距的行索引，以及S′v中是否存在具有稳定列间距的列索引，这些行列索引作为类格图案的分界以分割类格图案。对于S′h，将S′h中元素做升序排列，行索引间距的多重集定义为具有行索引连续的稳定行间距的多重集定义如下。
其中xi∈S′h，i为数值连续的正整数，例如i可以取2，3，4，但不能只取2和4。由于可能存在多个符合定义的连续间距，所以可能存在多个其中具有最多元素个数的(记作)所对应的行索引作为类格图案分界集合Sh的初始值，其定义如下。
类似地，可以计算和Sv的初始值，定义如下。
因为Sh与Sv初始值所对应的和分别包含稳定行间距和稳定列间距，因此与中位值定义为理想行数和理想列数Sh与Sv的初始值示例见附图3，如图3所示，纺织品图像仅部分区域被Sh和Sv的初始值同时覆盖，即图3标有“Sh和Sv的初始值”图示中的网格，只有竖直直线的区域仅有Sv覆盖。Sh和Sv的扩展基于和进行。由于Sh中的行索引按升序排列，从min(Sh)开始，以步长向Itc行索引最小值1移动，即计算检查是否存在 满足如果存在x′，则将x′添加进Sh，否则将x添加进Sh并保持Sh中元素升序排列，再次计算并重复上述步骤；类似地，从max(Sh)开始，以步长向Itc行索引最大值dimx(Ic)移动，即计算检查是否存在满足如果存在x′，则将x′添加进Sh，否则将x添加进Sh，再次计算并重复上述步骤。对Sv以步长做类似扩展，详见算法1。
扩展得到的Sh与Sv基本覆盖了纺织品图像大部分区域，如附图3所示。根据Sh与Sv分别包含的行与列索引，可以将Ic按这些索引所在的行与列进行分割，分割所得的区域定义为类格图案，其定义如下。
其中1≤i≤|Sh|-1且1≤j≤|Sv|-1，因此Ic分割出(|Sh|-1)·(|Sv|-1)个而i和j分别是以类格图案为单位的行与列索引。类格图案分割的一个重要参数是fc，如图3所示，Ic的二值化是基于阈值fc·max(Ic)完成，而该阈值取决于fc。如果存在多个无瑕疵纺织品图像，对于每幅图像应用基于不同fc取值的算法1能得到多个和计算所得的和的直方图，其中出现次数最多的和及其对应的fc取值对确定fc的最终取值具有参考意义。假设存在n个无瑕疵纺织品图像I1，I2...In和m个fc的可选值c1，c2…cm，令和分别表示输入为Ik，k＝1，2...n且fc＝cl，l＝1，2...m的算法1计算得到的和对于每个Ik都存在m个和m个在关于Ik的多重集和多重集中出现次数最多元素所对应的fc取值分别记为ch(k，l)和定义如下。
对于Ik，如果ch(k，l)与相同，则相应索引(k，l)保存在集合中，的定义如下。
对于每个即对于Ik，如果至少存在一个l′使得成立，则与med({cl})距离最近的ch(k，l′)参与fc的计算，fc的定义如下。
其中ch(k，l)的索引表示取符合定义中的k值。至此，类格图案分割及相关算法(算法1和2)的描述结束。下面介绍特征提取部分(算法3)的内容，如图1所示，特征提取采用了Gabor滤波器组对类格图案分割生成的1≤i≤|Sh|-1且1≤j≤|Sv|-1进行卷积，卷积结果通过投影方法降维，计算降维结果的能量和振幅，组合能量和振幅构造特征向量，最终产生特征向量距离矩阵，矩阵中的元素即为算法3提取的特征。
Gabor滤波器包含实部和虚部，本发明仅采用了Gabor滤波器的虚数部分，其定义如下。
其中表示尺度参数，表示角度参数，如果s取值的个数记作θ的取值个数记作可构成ns·nθ个Gabor滤波器，这些滤波器的集合称为Gabor滤波器组。每个分别与ns·nθ个Gabor滤波器进行卷积，可产生ns·nθ个卷积结果，将这些卷积结果旋转-θ，然后按行累加，即可将二维卷积结果变为一维，达到降维目的。令Gi，j(x|s，θ)表示行索引为x的行投影，其定义如下。
其中和为Gabor滤波器系数的索引，表示Gabor滤波器系数索引的取值范围，Gabor滤波器系数可排列为行数和列数均为(2·SG+1)的矩阵。对于将其卷积结果的各行投影Gi，j(x|s，θ)按x升序排列则得到一维的卷积投影Gi，j|s，θ，其定义如下。
其中表示按x升序排列的连接，连接方式与定义7相同。附图4以单个类格图案为例显示了本发明计算Gi，j|s，θ的阶段性结果，其中Gabor滤波器的s＝1，θ取值为0°，45°，90°和135°，因此ns＝1，nθ＝4。对于纺织品图像中的每个类格图案都会产生ns·nθ＝4个Gi，j|s，θ，计算每个Gi，j|s，θ的能量和振幅，并按s和θ的取值顺序连接构造的特征向量(其中1≤i≤|Sh|-1，1≤j≤|Sv|-1)，其定义如下。
其中Ei，j和Ai，j分别表示的能量和振幅。若Ei，j|s，θ表示Gi，j|s，θ的能量，则Ei，j|s，θ＝∑x Gi，j(x|s，θ)2为Gi，j|s，θ中所有元素的平方和，其中x为Gi，j|s，θ的元素索引，那么 其中表示按先s后θ的顺序连接。本发明中 类似地，定义Ai，j|s，θ表示Gi，j|s，θ的振幅，则Ai，j|s，θ＝∑x|Gi，j(x|s，θ)|为Gi，j|s，θ中所有元素的绝对值之和，其中x为Gi，j|s，θ的元素索引，那么本发明中对于含瑕疵的纺织品图像，如果瑕疵出现在部分类格图案中，那么有瑕疵的类格图案与无瑕疵类格图案在特征向量各元素取值上存在差异。例如附图5所示，图5显示了一幅有瑕疵纺织品图像的部分特征向量。图5中，纺织品图像由算法1分割产生14×11的类格图案，即中1≤i≤14，1≤j≤11，其中加粗边框显示的和包含瑕疵。图5显示了第7行类格图案的Gi，j|s，θ，标有“G7，j|s，θ”的图示中，从左向右所示的4列分别是G7，j|1，0°，G7，j|1，45°，G7，j|1，90°和G7，j|1，135°，从上到下所示的11行分别对应可以观察到加粗显示的G7，2|s，θ和G7，3|s，θ(其中s＝1，θ为0°，45°，90°和135°)与其他G7，j|s，θ(其中j≠2或3)存在差异。图5中标有的图示，从上到下分别是可以观察到加粗显示的和与其他(其中j≠2或3)存在差异。
由于有瑕疵的类格图案与无瑕疵类格图案在Gi，j|s，θ与(其中1≤i≤|Sh|-1，1≤j≤|Sv|-1)两者上差异大，而无瑕疵类格图案之间在Gi，j|s，θ与两者上差异小，本发明通过分析类格图案在Gi，j|s，θ上的差异，找出差异小的类格图案，在这些类格图案的基础上构建理想特征向量并比较与的差距，在一定程度上区别有瑕疵和无瑕疵的类格图案。本发明将构建的问题转换为找出每行的典型无瑕疵类格图案，基于每行的典型类格图案的特征向量构建对于每行类格图案，计算每个类格图案与其所在行中所有类格图案关于卷积投影的车贝雪夫距离(Chebychevdistance)，即对于给定的i，s和θ，对于每个(其中1≤j≤|Sv|-1)的取值，计算Cb(Gi，j|s，θ，Gi，1|s，θ)，Cb(Gi，j|s，θ，Gi，2|s，θ)...Cb(Gi，j|s，θ，Gi，|Sv|-1|s，θ)，并按此计算顺序排列成一行，一共有行，再将这些行按j的取值大小从上到下排列，构成如下关于第i行的距离矩阵
对于第i行，存在ns·nθ个对于本发明的每行类格图案，则需计算 和附图6显示了典型类格图案的计算结果，图中标有“G7，j|s，θ”的图例是图5中标有“G7，j|s，θ”图例的简化，表示的都是图5所示的第7行的卷积投影G7，j|s，θ。根据G7，j|s，θ计算得到的和以重叠方式显示为图6中标有的图例，(其中s＝1，θ为0°，45°，90°和135°)中的每个元素以一个色块的形式显示，色块越亮则元素值越大，色块越暗则元素值越小，所有关于车贝雪夫距离的矩阵在附图中都以这种方式图形化。的第j行表示与第i行中所有类格图案卷积投影的车贝雪夫距离，由于对角线元素表示与它自己卷积投影的车贝雪夫距离，所以对角线上的元素恒为0。移除该的对角线得到矩阵的定义如下。
的第j行表示与第i行中除外的其他类格图案卷积投影的车贝雪夫距离。计算每行的标准差，具有最小标准差的行的行索引j′对应的类格图案定义为第i行类格图案关于s和θ的典型无瑕疵类格图案，记作j′的定义如下。
其中表示中的第j行。对于第i行类格图案，如果存在多个对应相同s和θ的则任取其中一个作为唯一的关于s和θ的典型无瑕疵类格图案。对于不同的s和θ，的索引j′有可能相同，即不同的s和θ的典型无瑕疵类格图案是同一个类格图案，重合次数最多的类格图案定义为第i行的典型无瑕疵类格图案，记作其定义如下。
对于第i行类格图案，如果没有出现不同的s和θ的典型无瑕疵类格图案是同一个类格图案的情况，即对于给定的i和不同的s和θ，或立，则任取中的一个。在某些情况下不一定能代表第i行的无瑕疵类格图案，例如一行中瑕疵类格图案较多或瑕疵不明显时，有可能是有瑕疵类格图案中的一个。附图6显示了这种情况，图中标有“各行的典型图案”的图例以白色边框的形式显示1≤i≤14，其中第7行(即i＝7)整行以加粗边框形式突出显示，可以观察到是第7行中有瑕疵的类格图案。为了排除有瑕疵的(其中1≤i≤|Sh|)，可计算特征向量之间的车贝雪夫距离，即对于每个计算并将距离排列为如下所示的矩阵
附图6显示了图中标有的图例显示了根据标有“各行的典型图案”图例中(其中1≤i≤14)所计算的其中有瑕疵的在中所对应的第7行色块比其他行要亮，同时第7列色块比其他列要亮。中的第i行和第i列都表示与纺织品图像中所有典型图案特征向量的车贝雪夫距离，的对角线元素为因此，如果对应的是有瑕疵的相比较没有瑕疵的典型图案在中所对应的行和列，对应的第i行和第i列中除对角线以外元素的值较大，例如图6中明亮的第7行和第7列，按行累加中的元素，则可以量化中关于(其中1≤i≤|Sh|-1)的车贝雪夫距离的大小情况，即计算如下定义的距离向量d。
其中1≤i′≤|Sh|-1，表示按i的索引顺序进行连接，连接方式与定义7相同。距离向量d中的第i个元素di表示与所有典型图案的车贝雪夫距离之和。只有满足条件avg(d)-std(d)＜di＜avg(d)+std(d)的可以参与理想特征向量的计算，满足该条件的组成的多重集S*定义如下。
定义为S*的均值，即以作为标准，通过计算纺织品图像中类格图案特征向量与的车贝雪夫距离并分析的分布，识别有瑕疵的类格图案。将作为一个矩阵第i行和第j列的元素，得到距离矩阵其定义如下。
附图7显示了以为标准计算的简略过程，标有“纺织品图像中所有类格图案的特征向量的图例以图形的形式象征性地显示以特征向量表示的纺织品图像的类格图案，标有“理想特征向量的图例图形化表示了这两幅图例之间用标有“比较”的箭头连接，该箭头简略表示的计算方法。图7中标有“距离矩阵”的图例图形化显示了该图例被标有“等价”的箭头与标有“的Mesh图”的图例相连，即以三维图形的形式显示了的值，这两个图例中，有瑕疵的类格图案对应的和与其他相比具有较大的值，即Mesh图中的最高峰。
至此，特征提取及相关算法(算法3)的描述结束。下面介绍特征比较(算法5和6)的内容，如附图7所示，中有瑕疵类格图案对应的和的值很大，在的直方图，即标有“中元素分布的直方图”的图例中，和所对应的分布用加粗边框突出显示，其分布位于直方图的最右边。图7中标有“有瑕疵的类格图案”的图例显示了对应和的和因此，有瑕疵的类格图案所对应的往往靠近直方图中最大值分布的最右端区域，分离直方图中有瑕疵和无瑕疵类格图案所对应的可以转化为基于阈值的判断问题，即选择距离阈值d*，比较与d*，对应且符合一定条件的类格图案被标记为有瑕疵的类格图案。本发明中d*取一组无瑕疵纺织品图像最大值的平均值，计算过程由算法5描述。
结合距离阈值d*与的直方图可以分析直方图中的分布情况，从而确定有瑕疵的类格图案。由于d*是基于一组无瑕疵纺织品图像的阈值，对于一幅特定的纺织品图像，d*与准确分离该图像中有瑕疵和无瑕疵类格图案的真实阈值t*可能存在差异。本发明方法选择缺口值t′和断崖值t″为t*的近似值。令t表示的横轴刻度，即的取值范围，h(t)表示的纵轴刻度，即取值为t的的个数，t′表示d*＜t时h(t)的“第一个缺口”，t″表示d*＜t时h(t)的“第一个断崖”，其定义如下。
实际情况中，t′和t″并不一定存在，因此当t′存在时，则用t′近似t*，如果t′不存在而t″存在，则使用t″近似t*，如果t′和t″都不存在，则定义t*为d*。所有对应的类格图案被标记为有瑕疵的类格图案计算过程由算法6描述。
本发明方法的瑕疵检测效果评估中使用了香港大学电气和电子工程系工业自动化实验室提供的106幅像素大小为256×256的24位彩色纺织品图像，在实验中这些图像被转换为8位的灰度图像。106幅图像包括两种图案：箱形图像和星形图像，其中箱形图像包括26幅无瑕疵和30幅有瑕疵图像；星形图像包括25幅无瑕疵和25幅有瑕疵图像。两种图案的有瑕疵图像包括5种瑕疵类型：断端(broken end)，孔洞(hole)，网纹(netting multiple)，粗条纹(thick bar)和细条纹(thin bar)，每种瑕疵类型的具体数量详见表1与表2的第一列。所有瑕疵图像都有相同大小的瑕疵基准图(ground-truth image)，瑕疵基准图为2值图像，其中1表示瑕疵，0表示背景。用于比较的算法包括WGIS，BB，RB和ER。用于评估的指标包括真阳性(true positive，以下简称TP)，假阳性(positiverate，以下简称FPR)，真阳性率(truepositiverate，以下简称TPR)，假阳性率(positiverate，以下简称FPR)，阳性预测值(positivepredictivevalue，以下简称PPV)和阴性预测值negativepredictivevalue，以下简称NPV)。
TPR衡量瑕疵基准图中表示瑕疵的像素被算法正确标定为瑕疵的比例，FPR衡量瑕疵基准图中表示背景的像素被算法错误标定为瑕疵的比例，PPV衡量算法输出的瑕疵中瑕疵基准图中的瑕疵所占比例，NPV衡量算法输出的背景中瑕疵基准图中的背景所占比例。对于TPR，PPV和NPV，指标值越大越好，对于FPR则越小越好。
相关数学定义可以在文献(M.K.Ng，H.Y.T.Ngan，X.Yuan，，et al.，Patternedfabricinspectionandvisualizationbythemethodofimagedecomposition，IEEETrans.Autom.Sci.Eng.11(3)(2014)943-947)中找到。实验硬件平台为含处理器IntelCoreTMi7-3610QM 230-GHz和8.00GB内存的笔记本电脑，软件为Windows 8.1和Maltab8.4。
BB和RB输出二值图像(其中1表示瑕疵，0表示背景)，虽然大小与输入图像一致，但输出图像中靠近边缘的狭长区域未被处理，这些区域的像素值被设为0。ER和WGIS虽然也输出二值图像(其中1表示瑕疵，0表示背景)，但尺寸比输入图像小，因此其处理结果用最近邻插值将尺寸转换为输入图像尺寸。经过处理后的BB，RB，ER和WGIS的输出图像可以通过逻辑运算直接与瑕疵基准图进行比较。这四种算法都要求人工输入参数，其中ER和WGIS要求输入一幅图案模版。对于WGIS，图案模版从无瑕疵图像按名称排序的第一幅图像的左上角截取，箱形图像的图案模版的像素大小为27×25，星形图像的图案模版的像素大小为22×18。对于ER，图案模版像素大小恒为28×26，表示比赛次数的参数设为15。对于RB，规则带像素大小恒为25。对于BB，行带(row band)和列带(column band)的像素大小分别为15和25，标准差数量为2。
由于本发明方法并不输出二值图像，而是输出以行列索引形式表示的瑕疵区域，因此输出结果不能直接与瑕疵基准图比较。为计算评价指标，评估过程中采取生成二值图像的方法将瑕疵区域集合转换为二值图像(其中1表示瑕疵，0表示背景)，具体方法是将瑕疵基准图中被瑕疵区域覆盖的部分拷贝到一幅像素全为0且大小与瑕疵图相同的二值图像的相同位置，所合成的二值图像作为本发明方法参与评估的输出。对于箱形图像，本发明方法的参数fc为0.4722，d*为3415；对于星形图像，本发明方法的参数fc为0.5，d*为2184。
表1和表2显示了参与评估的几种算法的检测结果，其中每行(除标有“概况”的最后五行)都表示一种算法在特定瑕疵种类上的指标平均值，第一列表示瑕疵类型(除标有“概况”的最后五行)，其中括号中的数字表示该类型的纺织品图像数量，第二列表示算法标识为瑕疵的像素平均数量，第三列至第六列表示TPR平均值，FPR平均值，PPV平均值和NPV平均值(平均值或平均数量是算法对特定瑕疵类型的所有图像检测结果的平均值)，最后一列显示了算法名称，各列中的最优值以加粗字体形式显示。表1和表2标有“概况”的最后五行分别表示对箱形图像和星形图像所有纺织品图像的指标平均值。表1显示了箱形图像的检测结果，本发明方法的总体表现与RB相近，但TPR比RB和WGIS低。本发明方法的TPR在断端类型达到最大值。
表1
表2罗列了星形图像的检测结果，本发明方法的TPR平均值在孔洞，网纹和细条纹均达到最优，FPR平均值在除细条纹和断端外的所有瑕疵类型中均为最优，PPV平均值和NPV平均值在所有瑕疵类型中均为最优。在标有“概况”的最后五行中，本发明方法的总体TPR平均值，FPR平均值，PPV平均值和NPV平均值均为最优。
表2
图8和图9是评估涉及的各算法TPR和FPR的散点图(scatter plot)，两幅图都显示了根据不同瑕疵类型每幅纺织品图像检测结果计算的TPR和FPR组成的坐标点，例如箱形图像的瑕疵类型断端有5幅图像，则图8和图9中标有“箱形图像断端”的图例中，每个算法在TPR-FPR坐标系中都有5个以TPR和FPR为坐标值的点，每个点的TPR和FPR坐标分别表示算法对1幅瑕疵类型为断端的箱形图像检测结果的TPR和FPR指标值。有些散点图中存在个别算法的坐标点少于图像个数的情况，这是由于有些图像未被该算法识别为含瑕疵的图像，因此不存在相应的指标值。TPR和FPR的理想值分别是1.00和0.00，对应TPR-FPR坐标系的左上角。图8显示了箱形图像各算法检测结果的TPR和FPR值，其中本发明方法的TPR较低，WGIS虽然有较高的TPR但FPR很高。对于星形图像，即图9所示，本发明方的TPR-FPR点非常接近坐标系左上角，说明本发明方法具有较好的星形图像检测效果。
以上述依据本发明的理想实施例为启示，通过上述的说明内容，相关工作人员完全可以在不偏离本项发明技术思想的范围内，进行多样的变更以及修改。本项发明的技术性范围并不局限于说明书上的内容，必须要根据权利要求范围来确定其技术性范围。
Claims (2)
1.一种基于层次聚类和Gabor滤波的纺织品瑕疵检测方法，其特征在于，包括以下步骤：
输入灰度化的纺织品图像；
类格图案分割：分割纺织品图像产生类格图案，类格图案满足：相对于纺织品图像的行和列，类格图案按图像行的方向横向排列，并按列的方向纵向排列；在形态成分分析方法的纺织品图案卡通成分Ic中，类格图案具有几何形状并与背景像素在灰度上有显著差异；
特征提取：将类格图案与Gabor滤波器组的滤波进行卷积，对卷积结果降维得到一维投影并计算其能量和振幅，根据能量和振幅计算特征向量；计算一维投影之间的车贝雪夫距离，挑选每行类格图案的典型无瑕疵类格图案，根据典型无瑕疵类格图案的特征向量计算理想特征向量；每个类格图案特征向量与理想特征向量间的车贝雪夫距离为所提取的特征；和
特征比较：分析特征提取步骤中得到的基于特征向量的车贝雪夫距离直方图，获取表示有瑕疵的像素索引集合。
2.如权利要求1所述的基于层次聚类和Gabor滤波的纺织品瑕疵检测方法，其特征在于，具体步骤为：
步骤1：输入灰度化的纺织品图像，通过形态学成分分析方法计算输入图像的卡通成分Ic；
步骤2：使用阈值fc·max(Ic)二值化卡通成分Ic得到二值图像Itc；
步骤3：计算二值图像Itc中各行背景像素数的多重集各列背景像素数的多重集计算各行背景像素数的多重集和各列背景像素数的多重集的峰值和
步骤4：分别计算基于步骤(3)中的峰值和的HC算法轮廓系数，以最大轮廓系数对应的聚类个数初始化HC算法，并对峰值和分别进行聚类，得到的聚类中心构成多重集和
步骤5：根据多重集和计算阈值和分别计算中不小于的峰值所对应行索引的间距和中不小于的峰值所对应列索引的间距，计算构成稳定行间距的最长连续行索引的集合Sh以及构成稳定列间距的最长连续列索引的集合Sv；
其中，将中元素按降序排列，从开始，计算其与下一个元素的差的绝对值，若该值不大于DCT尺寸的高，则继续计算当前元素与下一个元素的差的绝对值并与DCT尺寸的高比较，若不大于则继续，若大于则下一个元素为并终止；如果一直没有出现大于DCT尺寸的高的情况，则即为
将中元素按降序排列，从开始，计算其与下一个元素的差的绝对值，若该值不大于DCT尺寸的高，则继续计算当前元素与下一个元素的差的绝对值并与DCT尺寸的宽比较，若不大于则继续，若大于则下一个元素为并终止；如果一直没有出现大于DCT尺寸的宽的情况，则即为
DCT尺寸：形态学成分分析方法在图像局部应用离散余弦变换，即首先将图像划分为不重叠且具有固定大小的矩形区域，然后对每个区域应用离散余弦变换，矩形区域的大小称为DCT尺寸，单位为像素，区域内一行的像素数称为DCT尺寸的宽，一列的像素数称为DCT尺寸的高；
步骤6：扩展集合Sh和集合Sv以覆盖图像的卡通成分Ic的大部分区域，并计算类格图案的理想行数和理想列数
步骤7：根据集合Sh和集合Sv将图像的卡通成分Ic分割为类格图案其中1≤i≤|Sh|-1，1≤j≤|Sv|-1，使用Gabor滤波器组与类格图案卷积并计算一维卷积投影Gi,j|s,θ，根据一维卷积投影Gi,j|s,θ的能量和振幅构建特征向量
其中i和j分别是类格图案以类格图案为单位的行与列索引；
表示尺度参数，表示角度参数；
步骤8：计算类格图案与第i行类格图案之间基于一维卷积投影Gi,j|s,θ的车贝雪夫距离，结果保存为矩阵的第j行，根据矩阵计算第i行类格图案的典型无瑕疵类格图案
步骤9：计算第i行类格图案的典型无瑕疵类格图案与第i行类格图案之间基于特征向量的车贝雪夫距离，结果保存为矩阵的第i行，根据矩阵计算第i行的典型无瑕疵类格图案的距离之和di；
步骤10：根据距离之和di筛选第i行的典型无瑕疵类格图案基于筛选结果计算理想特征向量
其中，理想特征向量为S*的均值；
为每行的典型类格图案的特征向量，i为类格图案行索引；
其中1≤i′≤|Sh|-1，表示按i的索引顺序进行连接；di为距离向量d中的第i个元素，表示与所有典型图案的车贝雪夫距离之和，典型图案是指每行的典型无瑕疵类格图案；
步骤11：计算特征向量和理想特征向量的车贝雪夫距离，结果保存为矩阵中索引为(i,j)的元素其中i为类格图案行索引，j为类格图案列索引；
步骤12：计算矩阵的直方图
步骤13：基于距离阈值d*和直方图计算缺口值t′和断崖值t″，根据缺口值t′和断崖值t″近似表示真实阈值t*；
步骤14：所有对应车贝雪夫距离＞真实阈值t*的类格图案被标记为有瑕疵的类格图案输出其中i为类格图案行索引，j为类格图案列索引。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611156224.7A CN106770323B (zh) | 2016-12-15 | 2016-12-15 | 基于层次聚类和Gabor滤波的纺织品瑕疵检测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611156224.7A CN106770323B (zh) | 2016-12-15 | 2016-12-15 | 基于层次聚类和Gabor滤波的纺织品瑕疵检测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106770323A true CN106770323A (zh) | 2017-05-31 |
| CN106770323B CN106770323B (zh) | 2019-05-28 |
Family
ID=58888181
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611156224.7A Active CN106770323B (zh) | 2016-12-15 | 2016-12-15 | 基于层次聚类和Gabor滤波的纺织品瑕疵检测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106770323B (zh) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107895371A (zh) * | 2017-11-24 | 2018-04-10 | 常州大学 | 基于峰值覆盖值和Gabor特征的纺织品瑕疵检测方法 |
| CN107945165A (zh) * | 2017-11-24 | 2018-04-20 | 常州大学 | 基于峰值覆盖值和面积计算的纺织品瑕疵检测方法 |
| CN107945164A (zh) * | 2017-11-24 | 2018-04-20 | 常州大学 | 基于峰值阈值、旋转校准和混合特征的纺织品瑕疵检测方法 |
| CN107967680A (zh) * | 2017-11-24 | 2018-04-27 | 常州大学 | 基于峰值阈值和混合特征的纺织品瑕疵检测方法 |
| CN107977961A (zh) * | 2017-11-24 | 2018-05-01 | 常州大学 | 基于峰值覆盖值和混合特征的纺织品瑕疵检测方法 |
| CN108010029A (zh) * | 2017-12-27 | 2018-05-08 | 江南大学 | 基于深度学习和支持向量数据描述的织物疵点检测方法 |
| CN109816631A (zh) * | 2018-12-25 | 2019-05-28 | 河海大学 | 一种基于新代价函数的图像分割方法 |
| CN111680750A (zh) * | 2020-06-09 | 2020-09-18 | 创新奇智(合肥)科技有限公司 | 图像识别方法、装置和设备 |
| CN115082482A (zh) * | 2022-08-23 | 2022-09-20 | 山东优奭趸泵业科技有限公司 | 一种金属表面缺陷检测方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102706881A (zh) * | 2012-03-19 | 2012-10-03 | 天津工业大学 | 基于机器视觉的布匹瑕疵检测方法 |
| CN102879401A (zh) * | 2012-09-07 | 2013-01-16 | 西安工程大学 | 基于模式识别和图像处理的纺织品瑕疵自动检测及分类方法 |
| CN103234976A (zh) * | 2013-04-03 | 2013-08-07 | 江南大学 | 基于Gabor变换的经编机布匹瑕疵在线视觉检测方法 |
-
2016
- 2016-12-15 CN CN201611156224.7A patent/CN106770323B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102706881A (zh) * | 2012-03-19 | 2012-10-03 | 天津工业大学 | 基于机器视觉的布匹瑕疵检测方法 |
| CN102879401A (zh) * | 2012-09-07 | 2013-01-16 | 西安工程大学 | 基于模式识别和图像处理的纺织品瑕疵自动检测及分类方法 |
| CN103234976A (zh) * | 2013-04-03 | 2013-08-07 | 江南大学 | 基于Gabor变换的经编机布匹瑕疵在线视觉检测方法 |
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107977961B (zh) * | 2017-11-24 | 2019-10-11 | 常州大学 | 基于峰值覆盖值和混合特征的纺织品瑕疵检测方法 |
| CN107945164B (zh) * | 2017-11-24 | 2019-07-26 | 常州大学 | 基于峰值阈值、旋转校准和混合特征的纺织品瑕疵检测方法 |
| CN107945164A (zh) * | 2017-11-24 | 2018-04-20 | 常州大学 | 基于峰值阈值、旋转校准和混合特征的纺织品瑕疵检测方法 |
| CN107967680A (zh) * | 2017-11-24 | 2018-04-27 | 常州大学 | 基于峰值阈值和混合特征的纺织品瑕疵检测方法 |
| CN107977961A (zh) * | 2017-11-24 | 2018-05-01 | 常州大学 | 基于峰值覆盖值和混合特征的纺织品瑕疵检测方法 |
| CN107895371B (zh) * | 2017-11-24 | 2021-10-01 | 常州大学 | 基于峰值覆盖值和Gabor特征的纺织品瑕疵检测方法 |
| CN107945165A (zh) * | 2017-11-24 | 2018-04-20 | 常州大学 | 基于峰值覆盖值和面积计算的纺织品瑕疵检测方法 |
| CN107967680B (zh) * | 2017-11-24 | 2019-07-09 | 常州大学 | 基于峰值阈值和混合特征的纺织品瑕疵检测方法 |
| CN107895371A (zh) * | 2017-11-24 | 2018-04-10 | 常州大学 | 基于峰值覆盖值和Gabor特征的纺织品瑕疵检测方法 |
| CN108010029A (zh) * | 2017-12-27 | 2018-05-08 | 江南大学 | 基于深度学习和支持向量数据描述的织物疵点检测方法 |
| CN109816631A (zh) * | 2018-12-25 | 2019-05-28 | 河海大学 | 一种基于新代价函数的图像分割方法 |
| CN111680750A (zh) * | 2020-06-09 | 2020-09-18 | 创新奇智(合肥)科技有限公司 | 图像识别方法、装置和设备 |
| CN115082482A (zh) * | 2022-08-23 | 2022-09-20 | 山东优奭趸泵业科技有限公司 | 一种金属表面缺陷检测方法 |
| CN115082482B (zh) * | 2022-08-23 | 2022-11-22 | 山东优奭趸泵业科技有限公司 | 一种金属表面缺陷检测方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106770323B (zh) | 2019-05-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106770323B (zh) | 基于层次聚类和Gabor滤波的纺织品瑕疵检测方法 | |
| CN107274385A (zh) | 基于类格图案及其Gabor特征的纺织品瑕疵检测方法 | |
| CN107248152A (zh) | 基于类格图案及其面积特征的纺织品瑕疵检测方法 | |
| Ng et al. | Patterned fabric inspection and visualization by the method of image decomposition | |
| CN105550660B (zh) | 一种机织物组织结构种类的识别方法 | |
| CN108154519A (zh) | 眼底图像中血管的分割方法、装置及存储介质 | |
| CN108288271A (zh) | 基于三维残差网络的图像检测系统及方法 | |
| Liu et al. | Fabric defect detection based on information entropy and frequency domain saliency | |
| CN107895371A (zh) | 基于峰值覆盖值和Gabor特征的纺织品瑕疵检测方法 | |
| Chang et al. | Fabric defect detection based on pattern template correction | |
| CN110245593A (zh) | 一种基于图像相似度的手势图像关键帧提取方法 | |
| CN111862027B (zh) | 基于低秩稀疏矩阵分解的纺织品瑕疵检测方法 | |
| CN107945164B (zh) | 基于峰值阈值、旋转校准和混合特征的纺织品瑕疵检测方法 | |
| Rezaei et al. | Texture classification approach based on combination of random threshold vector technique and co-occurrence matrixes | |
| Zhang et al. | A method for defect detection of yarn-dyed fabric based on frequency domain filtering and similarity measurement | |
| Techniková et al. | Weaving density evaluation with the aid of image analysis | |
| Jing et al. | Automatic recognition of weave pattern and repeat for yarn-dyed fabric based on KFCM and IDMF | |
| CN107977961B (zh) | 基于峰值覆盖值和混合特征的纺织品瑕疵检测方法 | |
| CN109507193A (zh) | 一种基于局部对比度增强与二值模式的织物瑕疵检测方法 | |
| CN107967680B (zh) | 基于峰值阈值和混合特征的纺织品瑕疵检测方法 | |
| CN112802074A (zh) | 一种基于光照校正与视觉显著特征的纺织品瑕疵检测方法 | |
| Bao et al. | Low-rank decomposition fabric defect detection based on prior and total variation regularization | |
| Mohanty et al. | Detection and classification of fabric defects in textile using image mining and association rule miner | |
| Jiang et al. | Automatic identification of jacquard warp-knitted fabric patterns based on the wavelet transform | |
| CN107945165B (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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20170531 Assignee: Changzhou Bai Jia Textile Technology Co.,Ltd. Assignor: CHANGZHOU University Contract record no.: X2023980049372 Denomination of invention: Textile defect detection method based on hierarchical clustering and Gabor filtering Granted publication date: 20190528 License type: Common License Record date: 20231203 |
|
| EE01 | Entry into force of recordation of patent licensing contract | ||
| EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20170531 Assignee: Shandong Hongde Yuheng Information Technology Co.,Ltd. Assignor: CHANGZHOU University Contract record no.: X2023980051060 Denomination of invention: Textile defect detection method based on hierarchical clustering and Gabor filtering Granted publication date: 20190528 License type: Common License Record date: 20231209 |
|
| EE01 | Entry into force of recordation of patent licensing contract |