一种基于平稳小波变换视觉显著性的织物疵点检测方法
【技术领域】
[0001] 本发明涉及纺织品图像处理的技术领域，具体涉及一种基于平稳小波变换（SWT) 视觉显著性的织物疵点检测方法，利用图像的平稳小波变换和显著性分析方法对织物疵点 图像进行疵点的检测和定位。
【背景技术】
[0002] 织物疵点检测是纺织品制造过程中质量控制的重要环节，通过对织物疵点的检测 有利于发现生产过程中的工艺问题，从而及时改进生产工艺、提高织物的质量等级。随着集 成电路和图像处理技术的飞速发展，机器视觉已经在工业表面检测领域中得到了越来越广 泛的应用，以计算机视觉来代替人工操作不仅可以提高检测速度，降低劳动成本，而且通过 布匹疵点自动检测系统可以为布匹质量等级的评定提供双方可信的参考标准，有利于国际 贸易的往来。织物疵点检测与判别方法是该类系统的核心环节，直接影响着系统的性能。
[0003] 目前，关于织物疵点检测算法大致分为以下三类:基于统计的方法、基于模型的方 法和基于频谱分析的方法。（1)基于统计的方法通过研究背景与疵点在纹理、灰度等方面的 不同统计特性实现织物疵点的有效检测，简单易行，得到了广泛的研究，但不同纹理、灰度 统计策略对检测效果影响较大。（2)基于模型的方法通过建立模型、参数估计提取图像纹理 特征，计算量大，实现复杂，在线学习尤为困难，疵点检测效果差目前已很少研究。（3)基于 谱分析的方法主要针对具有均匀纹理的织物图像，将其变换到频域时，其纹理由一些基本 的纹理基元按照周期性规则构成，纹理基元的周期性有利于进行疵点检测。常用的处理方 法有快速傅里叶变换(FFT)、基于Gabor滤波器的方法、Gabor滤波和形态学结合的方法、小 波变换等。基于谱分析的方法检测效果往往受滤波器组选择的影响较大。
[0004] 目前，提出的疵点检测方法在一定程度上达到了疵点检测目的，但仍有许多共性 的问题和新问题仍然未得到解决或仍待进一步研究：1)布匹种类较多，造成表面纹理多样 化(譬如:斜纹、花纹等），大部分方法对纹理比较简单的布匹检测效果好，而对复杂纹理的 织物检测效果较差，不能有效地把疵点与背景分离开来。2)织物疵点种类较多，有横裆疵、 斜纹疵、弓弧、断炜疵、斑点疵、扭结纱疵等90多种，目前的检测方法一般只能检测特定的几 种疯点类型，且检测精度有待提尚。
[0005] 根据大量神经学研究表明，人类视觉具有快速对视觉场景中感兴趣区域快速选择 与定位的能力。织物疵点图像虽然纹理复杂，但疵点在织物图像中较为显著。因此，基于视 觉显著性模型的织物疵点检测具有很好的研究价值。目前已有的视觉显著性模型主要分 为:基于特征融合的视觉显著性模型，基于图论的视觉显著性模型，基于频谱分析的视觉显 著性模型。近些年越来越多的学者倾向于利用频谱分析模型构建显著图。从上个世纪开始 小波理论就已经发展到很重要的地位。小波变换(Wavelet Transform，WT)在1980年代被首 次应用到信号处理上，在过去的几十年里，WT在图像处理领域中一直被认为具有很大的应 用潜力。WT具有多分辨率、空间和频率定位的特性，是进行图像特征提取的理想工具。基于 小波理论的优点，提出了许多传统的基于小波变换的图像检测方法如加权平均法和极大极 小值法。随后Tang等提出了基于多尺度小波变换的融合显著性检测算法，高频系数是通过 加权不规则的维数或者加权梯度特征获得，而低频系数则是通过加权局部能量值得到。 Jian提出了一个新的方案，测量原图像每个小波分解系数的显著值。它反映了小波分解系 数的视觉注意，符合人类的视觉感知。但是基于小波变换的视觉显著性模型还存在一些缺 陷：（1)小波变换存在多种小波基，一旦选定后，其特征就被固定下来。（2)小波系数是通过 尺度和平移变换获得的，每分解一次，逼近图像和细节图像的长度减小一半，重构时会丢失 原来的特征。
【发明内容】
[0006] 为了解决上述技术问题，本发明提供了一种基于平稳小波变换视觉显著性的织物 疵点检测方法，实现了对织物图像疵点的有效检测与定位，且具有较高的检测精度。
[0007] 为了达到上述目的，本发明的技术方案是：一种基于平稳小波变换视觉显著性的 织物疵点检测方法，包括平稳小波变换、背景估计、高斯估计模型及分割四部分，其步骤如 下：
[0008] 步骤一:将输入的原始RGB图像转化到灰度图像，对灰度图像进行图像去噪；
[0009] 步骤二:对去噪后的图像进行平稳小波变换得到一系列小波系数；
[0010] 步骤三:对除去平稳小波变换高频成分的小波系数进行逆平稳小波变换，重构得 到特征图；
[0011] 步骤四：对重构的特征图进行背景估计，得到局部显著图；
[0012] 步骤五:采用高斯估计模型对局部显示图进行全局估计，得到全局显著图；
[0013] 步骤六:利用自适应阈值分割算法对全局显著图进行分割，实现疵点与背景分离。
[0014] 所述对灰度图像进行图像去噪的方法是：利用高斯低通滤波器对灰度图像I进行 滤波其中，1为mXm的2D高斯滤波器，Γ为去除噪声后的灰度图像，*表示卷积 操作。
[0015]所述平稳小波变换选择'coiflet'滤波器作为小波基进行分解：[An，Hs，V s，Ds] = SWT0- )，其中，N表示平稳小波变换分解的层数，s e {1，…，N}; An是低频分解的系数，Hs、Vs、 Ds分别是在s层水平、垂直、对角方向的系数。
[0016]所述重构特征图的方法是:除去平稳小波变换分解产生的高频成分小波系数后进 行逆平稳小波变换，得到特征图：F=ISWT(A，H，V)，其中，A、H、V分别表示图像的低频系数、 水平方向系数和竖直方向系数。
[0017]所述对重构的特征图进行背景估计，得到局部显著图的方法是:将特征图F不重叠 的分成sXs大小相等的图像块Fa，a=l，2,…，K，K为原始图像分成图像块的块数;估计每个 图像块F a的均值:BG = {BG!，BG2，…，BGK}，其中，BGa为图像块?3的背景均值;计算原始图像和 每个图像块F a的欧式距离D(i，j) = {| Il'ajhBGil |，| |l'（i，j)-BG2| |，···| |l'（i，j)-BGK I }，其中，Γ (i，j)是灰度图像I经过mXm的2D高斯滤波器模糊得到；由欧式距离得到局部 显著图：S(i, j)=min{D(i, j)}。
[0018]所述得到全局显著图的方法是:采用正态分布的高斯密度函数对局部显示图S(i， j)进行全局估计：
[0019]
其中，n=l 表 示局部显示图S(i，j)的空间维数，2 = ￡[(€(1，」）4)(汽1，」）4，]是局部显著图5(1」） 的方差，y = E[S]是局部显著图S(i，j)的均值；由局部显著图S(i，j)得到显著图SG(i，j):SG (i，j) = (l〇g(P(f(i，J·))-")1/2。
[0020] 本发明的有益效果为：
[0021] (1)本发明综合考虑织物纹理特征的随机性和疵点种类的多样性，基于频谱分析 的视觉显著性模型能有效利用图像的全局与局部信息，取得了较好的效果。
[0022] (2)本发明不需要提取任何疵点信息，不需要参考样本，自适应性和鲁棒性较强。 [0023] (3)本发明不仅适用于检测表面纹理简单的织物，也适用于检测纹理较为复杂的 织物，扩大了算法的使用范围;且计算速度较快，适合在线检测。
【附图说明】
[0024]为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本 发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以 根据这些附图获得其他的附图。
[0025]图1为本发明的流程框图。
[0026] 图2(a)_(f)是本发明的原始待测疵点图像。
[0027] 图3(a)_(f)是本发明对原始待测疵点图像背景估计产生的局部显著图。
[0028] 图4(a)_(f)是本发明对原始待测疵点图像高斯估计模型估计产生的显著图。
[0029]图5(a)_(f)是本发明对原始待测疵点图像采用文献[1]的检测结果。
[0030]图6(a)_(f)是本发明对原始待测疵点图像采用文献[2]的检测结果。
[0031]图7(a)_(f)是本发明对原始待测疵点图像采用文献[3]的检测结果。
[0032]图8(a)_(f)是本发明对原始待测疵点图像基于本发明阈值分割的检测结果。
【具体实施方式】
[0033]下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完 整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于 本发明中的实施例，本领域普通技术人员在没有付出创造性劳动前提下所获得的所有其他 实施例，都属于本发明保护的范围。
[0034]如图1所示，一种基于平稳小波变换视觉显著性的织物疵点检测方法，包括平稳小 波变换、背景估计、高斯估计模型及分割四部分，其步骤如下：
[0035]步骤一:将输入的原始RGB图像转化到灰度图像，对灰度图像进行图像去噪。
[0036] 首先将输入的RGB图像转化到灰度图像，然后对灰度图像进行图像去噪，避免噪声 对图像分割的干扰，具体方法是:利用高斯低通滤波器对灰度图像I进行滤波：
[0037] l'=l*lmXm (1)；
[0038] 其中，1为mXm的2D高斯滤波器，Γ为去除噪声后的灰度图像，*表示卷积操作。此 处，取m = 5〇
[0039] 步骤二:对去噪后的图像进行平稳小波变换得到一系列小波系数。
[0040] 对一幅织物图像进行平稳小波分解，这样图像就相应的被分成多层次多尺度的具 有原始图像大小的子图像。平稳小波变换选择'coiflet'滤波器作为小波基进行分解：
[0041] [AN,Hs,Vs,Ds] = SffT(r) (2)，
[0042] 其中，N表示平稳小波变换分解的层数，se {1，···，N} ;An是低频分解的系数，HS、VS、 Ds分别是在s层水平、垂直、对角方向的系数。根据织物图像纹理特性，选择N=6。
[0043] 步骤三:对除去平稳小波变换高频成分的小波系数进行逆平稳小波变换，重构得 到特征图。
[0044] 小波系数表示了图像在不同尺度下的图像细节信息，能够用来创建频率带宽增大 的特征图谱，特征图是由逆平稳小波变换计算得到。由于已经对原始图像进行了高斯滤波 处理，因此用除去SWT分解产生的高频成分的小波系数重构即可得到特征图。
[0045] 所述重构特征图的方法是:除去平稳小波变换分解产生的高频成分小波系数后进 行逆平稳小波变换，得到特征图：F