一种基于单演小波分析的织物疵点检测方法技术

本发明专利技术涉及一种基于单演小波分析的织物疵点检测方法，首先对织物图像进行单演小波变换得到多分辨率三元小波系数向量，计算向量方向与向量模得到多分辨率方向子带和振幅子带；然后对各分辨率下的方向子带和振幅子带进行预处理后，计算其响应优秀程度，选择最优响应的方向子带和振幅子带；最后对最优响应的方向子带和最优响应的振幅子带分别进行二值化操作，得到两个二值图像；将二者融合得到织物疵点检测结果，值为255的像素表示无疵区域，值为0的像素表示疵点区域。本发明专利技术使用单演小波分析提取的织物图像振动特征避免使用函数响应，具有更好的自适应性并且不受光照不匀影响；对多数检测方法难以识别的微弱纹理变化类疵点具有优秀的检测效果。

【技术实现步骤摘要】

本专利技术属图像分析处理
，涉及一种基于单演小波分析的织物疵点检测方法，应用于纺织品表面质量自动检测与控制领域。
技术介绍
在基于机器视觉和图像处理技术的织物疵点检测方法中，难度最大的问题是对种类繁多形式各异的织物组织和疵点形式难以实现通用性。现有方法通常是使用固定基函数、滤波器、特征值、纹理单元等对象的响应，因此理论上对不同织物组织和疵点形式不具有通用性；使用人工神经网络对如小波变换、Gabor变换中滤波器和参数进行优化也无法从根本解决问题。解析信号分析是一种与上述方法中所用的图像分析技术有着本质不同的方法，具有不依赖基函数响应和高度时频同步优势，但其使用的Hilbert变换的二维扩展存在许多问题，产生了不同的解析信号形式。Krieger等人(1996)提出将原信号沿x轴或y轴进行Hilbert变换后作为虚部与原信号组成单向解析信号。Stark(1971)提出将原信号沿x轴和y轴各进行一次Hilbert变换和作为虚部与原信号组成全向解析信号。Hahn(1992)提出基于单向和全向Hilbert变换结果的单象限解析信号。Bulow等人(2001)提出包含两个单向Hilbert变换结果和全向Hilbert变换结果及原信号的四元解析形式。但上述方法非各向同性，无法避免真实振动方向与坐标轴夹角引起的系统性误差。Felsberg等人(2001)提出一种三元解析信号形式，包含两个单向Hilbert变换结果作为虚部，称为单演信号，并提出计算局部最大Hilbert能量方向分析和该方向上的振幅与频率分析方法，解决了方向引起的系统性误差问题，但该方法仍然缺乏多分辨率分析能力，对图像尺度无适应性。
技术实现思路
本专利技术所要解决的技术问题是克服现有算法的不足，提高算法的对不同纹理和疵点适应性，提供一种基于单演小波分析的织物疵点检测方法。本专利技术的一种基于单演小波分析的织物疵点检测方法：(1)对图像进行多分辨率单演分析，即对图像依次进行下采样，在每个分辨率都进行单演分析，得到各分辨率下的振动方向和振幅计算结果。具体实现方式为使用三元向量形式的单演小波对织物图像进行单演小波变换，得到多分辨率的三元小波系数向量(wi[k]，r1，i[k]，r2，i[k])；然后根据多分辨率的三元小波系数向量方向和向量模，计算各分辨率下最大Hilbert能量 方向，得到多分辨率方向子带；计算各分辨率下瞬时Hilbert域振幅，得到多分辨率振幅子带；(2)对多分辨率方向子带进行预处理后，计算每个分辨率方向子带的响应优秀程度，选择最优响应方向子带；所述的方向信号最优子带判断标准为内部均方差σ最低，其依据为该子带分辨率尺度最接近织物纹理真实尺度，因此方向信号反应信息最为准确，在织物正常纹理处的振动方向一致性最高，疵点区域与无疵区域间差异最明显；对多分辨率振幅子带进行预处理后，计算每个分辨率振幅子带的响应优秀程度，选择最优响应振幅子带；所述的一种基于单演小波分析的织物疵点检测方法，所述的振幅信号最优子带判断标准为最大类间均值差MICMD最高，其依据为该子带分辨率振幅反映疵点纹理模式变化最充分，疵点区域与无疵区域间差异最明显；(3)对最优响应方向子带进行二值化操作，得到最优响应方向子带的二值化矩阵，对应为最优响应方向子带的分割图像SO(x，y)；所述二值化操作使用的阈值化方法可使用图像处理领域中多种常见的方法，因最优响应子带前景背景间反差已较为强烈，对阈值化方法的选择无要求。对最优响应振幅子带进行二值化操作，得到最优响应振幅子带的二值化矩阵，对应为最优响应振幅子带的分割图像SA(x，y)；所述二值化操作使用的阈值化方法可使用图像处理领域中多种常见的方法，因最优响应子带前景背景间反差已较为强烈，对阈值化方法的选择无要求。融合SO(x，y)与SA(x，y)，得到融合图像SF(x，y)＝SO(x，y)*SA(x，y)/255；所述融合图像中，值为255的像素表示无疵区域，值为0的像素表示疵点区域，即为织物图像的疵点检测结果，也即织物疵点检测结果。作为优选的技术方案：如上所述的一种基于单演小波分析的织物疵点检测方法，所述三元向量形式的单演小波为复Riesz-Laplacian小波；所述三元向量形式的单演小波的具体构建如下：结合各向同性的分数阶Laplacian算子(-Δ)α与多重调和样条φλ，定义各向同性分数阶Laplacian样条小波ψ(x)＝(-Δ)λ/2φ2λ(Dx)，对Laplacian样条小波进行Riesz变换得到三元向量形式的复Riesz-Laplacian小波即三元向量形式的单演小波。如上所述的一种基于单演小波分析的织物疵点检测方法，计算各分辨率下最大Hilbert能量方向，得到多分辨率方向子带，具体为：分辨率i下位置k处的最大Hilbert能量方向为 其中Jmn＝[Ji(k)]mn，m，n∈1，2为分辨率i下位置k处的加权结构矩阵；计算各分辨率下瞬时Hilbert域振幅，得到振幅子带，具体为：分辨率i下位置k处的瞬 时Hilbert域振幅为如上所述的一种基于单演小波分析的织物疵点检测方法，对多分辨率方向子带进行预处理是指：对每个分辨率方向子带的数据分布进行统计，当数值位于中的像素点数超过图像面积的75％时，将负值数据+π，正值数据-π。如上所述的一种基于单演小波分析的织物疵点检测方法，对多分辨率振幅子带进行预处理是指将每个分辨率的振幅子带标准化，具体为：将某分辨率振幅子带记为A，A的最大值记为AMAX，A的最小值记为AMIN，标准化后的振幅子带为ASTD＝(A-AMIN)/(AMAX-AMIN)。如上所述的一种基于单演小波分析的织物疵点检测方法，计算每个分辨率方向子带的响应优秀程度，选择最优响应方向子带，具体步骤为：计算预处理后的每个分辨率方向子带标准差σ，选择σ最小并满足σ<0.1的子带作为最优响应方向子带。如上所述的一种基于单演小波分析的织物疵点检测方法，计算每个分辨率振幅子带的响应优秀程度，选择最优响应振幅子带，具体步骤为：对每个分辨率下预处理后的振幅子带ASTD，分别进行OTSU阈值分割，计算分割结果前景与背景间均值差的绝对值，定义为最大类间均值差MICMD，选择MICMD最大的子带作为最优响应振幅子带。如上所述的一种基于单演小波分析的织物疵点检测方法，所述二值化操作是指计算最优响应的方向子带或最优响应的振幅子带中所有元素的均值μ与标准差σ，当残差图像中某一元素的数值介于μ±3σ之间时，赋值为0；反之赋值为255。如上所述的一种基于单演小波分析的织物疵点检测方法，所述织物图像为位深度为8位以上的灰度图像；所述的多分辨率范围为：[1，log2(N)]。有益效果1、本专利技术通过使用单演小波分析提取的织物图像振动方向和振幅，是振动本身特征而非函数响应，具有更好的自适应性；2、本专利技术的方法分析对象为振动模式，不受光照不匀影响；3、本方法具有对振动方向的分析，对多数检测方法难以识别的微弱纹理变化类疵点具有优秀的检测效果。附图说明图1为吊经疵点图像图2为吊经疵点图像的最优响应方向子带图3为吊经疵点图像的疵点检测结果图4为双纬疵点图像图5为双纬疵点图像的最优响应振幅子带图6为双纬疵点图像的疵点检测结果具体实施方式下面结合具体实施方式本文档来自技高网...

【技术保护点】
一种基于单演小波分析的织物疵点检测方法，其特征是：(1)使用三元向量形式的单演小波对织物图像进行单演小波变换，得到多分辨率的三元小波系数向量(wi[k],r1,i[k],r2,i[k])；然后根据多分辨率的三元小波系数向量方向和向量模，计算各分辨率下最大Hilbert能量方向，得到多分辨率方向子带；计算各分辨率下瞬时Hilbert域振幅，得到多分辨率振幅子带；(2)对多分辨率方向子带进行预处理后，计算每个分辨率方向子带的响应优秀程度，选择最优响应方向子带；对多分辨率振幅子带进行预处理后，计算每个分辨率振幅子带的响应优秀程度，选择最优响应振幅子带；(3)对最优响应方向子带进行二值化操作，得到最优响应方向子带的二值化矩阵，对应为最优响应方向子带的分割图像SO(x,y)；对最优响应振幅子带进行二值化操作，得到最优响应振幅子带的二值化矩阵，对应为最优响应振幅子带的分割图像SA(x,y)；融合SO(x,y)与SA(x,y)，得到融合图像SF(x,y)＝SO(x,y)*SA(x,y)/255；所述融合图像中，值为255的像素表示无疵区域，值为0的像素表示疵点区域，即为织物图像的疵点检测结果，也即织物疵点检测结果。...

【技术特征摘要】
1.一种基于单演小波分析的织物疵点检测方法，其特征是：(1)使用三元向量形式的单演小波对织物图像进行单演小波变换，得到多分辨率的三元小波系数向量(wi[k],r1,i[k],r2,i[k])；然后根据多分辨率的三元小波系数向量方向和向量模，计算各分辨率下最大Hilbert能量方向，得到多分辨率方向子带；计算各分辨率下瞬时Hilbert域振幅，得到多分辨率振幅子带；(2)对多分辨率方向子带进行预处理后，计算每个分辨率方向子带的响应优秀程度，选择最优响应方向子带；对多分辨率振幅子带进行预处理后，计算每个分辨率振幅子带的响应优秀程度，选择最优响应振幅子带；(3)对最优响应方向子带进行二值化操作，得到最优响应方向子带的二值化矩阵，对应为最优响应方向子带的分割图像SO(x,y)；对最优响应振幅子带进行二值化操作，得到最优响应振幅子带的二值化矩阵，对应为最优响应振幅子带的分割图像SA(x,y)；融合SO(x,y)与SA(x,y)，得到融合图像SF(x,y)＝SO(x,y)*SA(x,y)/255；所述融合图像中，值为255的像素表示无疵区域，值为0的像素表示疵点区域，即为织物图像的疵点检测结果，也即织物疵点检测结果。2.根据权利要求1所述的一种基于单演小波分析的织物疵点检测方法，其特征在于，所述三元向量形式的单演小波为复Riesz-Laplacian小波；所述三元向量形式的单演小波的具体构建如下：结合各向同性的分数阶Laplacian算子(-Δ)α与多重调和样条φλ，定义各向同性分数阶Laplacian样条小波ψ(x)＝(-Δ)λ/2φ2λ(Dx)，对Laplacian样条小波进行Riesz变换得到三元向量形式的复Riesz-Laplacian小波即三元向量形式的单演小波。3.根据权利要求1所述的一种基于单演小波分析的织物疵点检测方法，其特征在于，计算各分辨率下最大Hilbert能量方向，得到多分辨率方向子带，具体为：分辨率i下位置k处的最大Hilbert能量方向为其中Jmn...

【专利技术属性】
技术研发人员：厉征鑫，高卫东，周建，潘如如，
申请(专利权)人：江南大学，
类型：发明
国别省市：江苏;32