王傳桐， 胡 峰， 徐啟永， 吳雨川， 余聯(lián)慶

(武漢紡織大學(xué) 機(jī)械工程與自動(dòng)化學(xué)院， 湖北 武漢 430074)

采用Gabor濾波簇和等距映射算法的織物疵點(diǎn)檢測方法

王傳桐， 胡 峰， 徐啟永， 吳雨川， 余聯(lián)慶

(武漢紡織大學(xué) 機(jī)械工程與自動(dòng)化學(xué)院， 湖北 武漢 430074)

為提高織物疵點(diǎn)檢測率，將Gabor濾波法與等距映射方法進(jìn)行融合，克服疵點(diǎn)檢測過程中存在的問題。首先用由3個(gè)尺度和5個(gè)方向組成的15個(gè)Gabor濾波器簇對織物疵點(diǎn)圖像進(jìn)行濾波，減少疵點(diǎn)圖像光照不均和對比度低的影響；然后將濾波圖像劃分成面積相等且互不重合的鄰域，并從鄰域中提取高維特征向量。采用等距映射方法對高維特征向量進(jìn)行降維，剔除高維特征中冗余信息，強(qiáng)化分類器擬合能力；再用低維嵌入模型提取新增樣本低維特征向量，用于概率神經(jīng)網(wǎng)絡(luò)分類器分類，檢測是否存在疵點(diǎn)；最后用2種不同紋理的織物進(jìn)行檢測實(shí)驗(yàn)。結(jié)果表明，本文方法能有效提高疵點(diǎn)的檢測精度。

疵點(diǎn)檢測； Gabor濾波簇； 等距映射； 圖像處理

中國是紡織品生產(chǎn)大國，提高紡織品疵點(diǎn)檢測準(zhǔn)確率對提升紡織品質(zhì)量和競爭力具有重要意義。小波分析法是常用的疵點(diǎn)檢測手段之一[1]。

離散正交小波具有快速M(fèi)allat算法支持，適合紋理的實(shí)時(shí)分析[2]。它從3個(gè)方向(0°、45°、90°)進(jìn)行多尺度織物圖像分解。在小波基和尺度選擇合適的情況下，織物背景紋理將被抑制，疵點(diǎn)邊緣能顯著增強(qiáng)，但是，正交離散小波進(jìn)行疵點(diǎn)圖像分解時(shí)會(huì)出現(xiàn)明顯相位畸變，不利于疵點(diǎn)檢測[3]。復(fù)輪廓變換法比離散正交小波更具優(yōu)勢，該方法核心是雙樹復(fù)小波變換和方向?yàn)V波器[4]，可在6個(gè)方向(±15°、±45°、±75°)下對疵點(diǎn)圖像進(jìn)行不同尺度分解，并且雙樹小波是復(fù)小波，回避了正交離散小波相位畸變的缺點(diǎn)。與前2種小波方法相比，Gabor小波的方向和尺度不受限制，并且和雙樹小波一樣也屬于復(fù)小波。用多尺度、多方向的Gabor濾波簇對疵點(diǎn)圖像進(jìn)行濾波，能夠提高疵點(diǎn)圖像對比度，抑制正常紋理信息和噪聲[5-6]。通常情況下，Gabor小波實(shí)部用于平滑疵點(diǎn)圖像，虛部用于檢測疵點(diǎn)邊緣[7]。Gabor濾波函數(shù)通過剔除直流分量，使其具有對光照變化不敏感的特點(diǎn)，但是，最優(yōu)Gabor小波濾波器設(shè)計(jì)困難，基于Gabor濾波簇的疵點(diǎn)特征向量又具有維數(shù)高和冗余信息大的缺點(diǎn)[7]。等距映射算法(Isomap)是非線性降維方法，它能保持?jǐn)?shù)據(jù)在空間中的幾何結(jié)構(gòu)，可獲得高維特征向量本質(zhì)維數(shù)[8]，有助于克服基于Gabor濾波簇特征向量的維數(shù)高和冗余信息大的缺點(diǎn)。

本文將2種方法結(jié)合用于織物疵點(diǎn)檢測，并提出Isomap結(jié)構(gòu)參數(shù)的優(yōu)化方法和新增樣本的低維嵌入模型，以用于提高疵點(diǎn)的檢測準(zhǔn)確率。

1 織物疵點(diǎn)檢測算法

圖1示出疵點(diǎn)檢測流程圖。

圖1 疵點(diǎn)檢測流程圖Fig.1 Flow chart of fabric detection

步驟1：利用Gabor濾波簇進(jìn)行高維特征向量提取。

2)將疵點(diǎn)灰度圖I和Gabor濾波器進(jìn)行卷積，生成15張濾波圖像。濾波后的疵點(diǎn)圖像可表示為

(1)

式中?表示卷積算子。

3)將每張濾波圖像劃分成N個(gè)尺寸為M×M互不重疊的子窗口。

4)針對灰度圖子窗口I(a,b)(第b行a列子窗口)，將15張濾波圖像中相同位置子窗口Fλθ(a,b)的灰度值依次連接組成高維特征向量，向量維數(shù)為M×M×15。

步驟2：訓(xùn)練樣本的降維。

使用Isomap算法[10]對高維特征向量進(jìn)行降維，消除冗余信息，有助于降低分類器泛化能力和提高分類器擬合能力，達(dá)到提高分類器分類精度的目的。降維過程中，嵌入維數(shù)d和鄰域大小k會(huì)影響降維精度。因此，采用離散粒子群算法[11-12]進(jìn)行參數(shù)優(yōu)化。適應(yīng)度函數(shù)采用Ncut準(zhǔn)則[13-14]：

(2)

(3)

式中，織物樣本分為2類(正常樣本和疵點(diǎn)樣本)，Ci是屬于類i的樣本集合，S(xi,xj)是樣本xi和xj間的邊權(quán)值，d(xi,xj)是xi和xj間的歐式距離，σi(σj)是xi(xj)與鄰域點(diǎn)的平均距離。如果樣本xi和xj之間沒有邊，則邊權(quán)值為零。由式(3)可知，S∈[0,1]，它使不同維度空間計(jì)算得到的邊權(quán)值具有可比性，從而使用于衡量類可分性的Ncut值在不同維度空間具有可比性。當(dāng)類Ci邊界點(diǎn)與其他類邊界點(diǎn)相距越遠(yuǎn)，邊權(quán)值S越小，式(2)的分子越??；當(dāng)類內(nèi)Ci數(shù)據(jù)點(diǎn)相互越靠近時(shí)，邊權(quán)值S越大，式(1)的分母越大。Ncut值越小時(shí)，類可分性越好。

1)初始化離散粒子群參數(shù)。設(shè)置種群大小為22，粒子維數(shù)為30，最大迭代次數(shù)為100。

2)利用Isomap算法計(jì)算訓(xùn)練樣本的低維輸出向量。計(jì)算每個(gè)粒子適應(yīng)度值，獲取最佳嵌入維數(shù)和鄰域大小。運(yùn)用Isomap算法，使用得到的嵌入維數(shù)和鄰域值對訓(xùn)練樣本進(jìn)行降維，得到訓(xùn)練樣本的低維嵌入結(jié)果。

步驟3：驗(yàn)證樣本的降維。

輸入新樣本xnew，它在訓(xùn)練樣本集中k個(gè)最近鄰點(diǎn)為xj，j=1,2,…,k。yj為xj的低維輸出向量。X=[x1,x2,…,xk]為高維訓(xùn)練樣本集合，Y=[y1,y2,…,yk]為訓(xùn)練樣本低維特征向量。首先找出新樣本在訓(xùn)練樣本集中的k個(gè)鄰近點(diǎn)。然后計(jì)算xnew和xj間的權(quán)值wij[15]。最后利用式(4)計(jì)算新樣本的低維特征向量ynew。

(4)

步驟4：疵點(diǎn)檢測。

1)使用訓(xùn)練樣本對概率神經(jīng)網(wǎng)絡(luò)(PNN)[16]進(jìn)行訓(xùn)練，將訓(xùn)練樣本中正常樣本的低維特征向量、疵點(diǎn)樣本的低維特征向量輸入概率神經(jīng)網(wǎng)絡(luò)進(jìn)行訓(xùn)練。

2)將驗(yàn)證樣本的低維特征向量輸入概率神經(jīng)網(wǎng)絡(luò)進(jìn)行分類，并對存在疵點(diǎn)的子窗口進(jìn)行標(biāo)注。

2 實(shí)驗(yàn)結(jié)果與分析

在相同條件下，將5種方法進(jìn)行對比實(shí)驗(yàn)研究，驗(yàn)證本文方法的合理性和有效性。這5種方法如下所示。

1)本文提出的基于Gabor小波與Isomap方法相結(jié)合的疵點(diǎn)檢測方法(簡稱Isomap方法)。

2)基于Gabor小波和灰度共生矩陣疵點(diǎn)檢測方法(簡稱Gabor方法)。該方法利用Gabor濾波簇進(jìn)行預(yù)處理的部分與方法1相同。不同之處在于：對每個(gè)子窗口提取4個(gè)方向(0°,45°,90°,135°)的對比度、能量、熵和相關(guān)性，共計(jì)16個(gè)特征。然后，將15張濾波圖像中相同位置子窗口獲得的16個(gè)特征依次連接組成高維特征向量，向量維數(shù)為240。

3)基于Gabor小波、灰度共生矩陣和主成分分析(PCA)相結(jié)合的檢測方法(簡稱G-PCA方法)。該方法采用PCA對方法2中高維特征向量進(jìn)行降維。降維過程中選取出特征值占總特征值0.85以上的相互正交的特征值。

4)基于Gabor小波和PCA相結(jié)合檢測方法(簡PCA方法)。該方法利用Gabor濾波簇進(jìn)行預(yù)處理的部分與方法1相同。不同之處在于：用PCA方法取代Isomap方法進(jìn)行降維。降維過程中選取出特征值占總特征值0.85以上的相互正交的特征值。

5)基于Gabor小波、灰度共生矩陣和Isomap方法相結(jié)合的檢測方法(簡稱G-Isomap方法)，該方法與方法3近似，不同之處在于：利用Isomap方法取代PCA方法對方法2中高維特征向量進(jìn)行降維。

研究樣本來源于自建的Fabric Defects Database of WTU數(shù)據(jù)庫中E1子數(shù)據(jù)庫(平紋組織；緯密為256根/10 cm；經(jīng)密為256根/10 cm；線密度為15.63 tex；采集分辨率為0.044 mm/像素)和E2子數(shù)據(jù)庫(平紋組織；緯密為264根/10 cm；經(jīng)密為340根/10 cm；線密度為14.08 tex；采集分辨率為0.089 mm/像素)。2個(gè)子數(shù)據(jù)庫中包含雜經(jīng)、缺經(jīng)、破洞、油污、粗緯和竹節(jié)等疵點(diǎn)類型，疵點(diǎn)灰度圖尺寸為300像素×300像素。分別從子數(shù)據(jù)庫中選取40張圖片作為訓(xùn)練樣本集，每類疵點(diǎn)圖片數(shù)量為10張；20張圖片作為驗(yàn)證樣本集，每類疵點(diǎn)圖片數(shù)量為5張。

利用式(5)計(jì)算驗(yàn)證樣本集的疵點(diǎn)分類準(zhǔn)確率。

(5)

式中：R為準(zhǔn)確識(shí)別率；Pr為準(zhǔn)確識(shí)別子窗口個(gè)數(shù)；PA為識(shí)別子窗口總數(shù)。

使用5種方法對E1、E2數(shù)據(jù)庫中織物疵點(diǎn)進(jìn)行疵點(diǎn)對比檢測，疵點(diǎn)識(shí)別過程中分別選用尺寸為25像素×25像素、30像素×30像素和50像素×50像素的子窗口，研究子窗口尺寸對檢測精度的影響。

表1示出5種方法對E1數(shù)據(jù)庫的疵點(diǎn)檢測準(zhǔn)確率。表2示出5種方法對E2數(shù)據(jù)庫的疵點(diǎn)檢測準(zhǔn)確率。圖2示出在子窗口大小為30像素×30像素時(shí)，5種方法對E1數(shù)據(jù)庫的部分檢測結(jié)果。圖3示出在子窗口大小為30像素×30像素時(shí)，5種方法對E2數(shù)據(jù)庫的部分檢測結(jié)果。

觀察表1、2發(fā)現(xiàn)：在不同尺度下，PCA方法對E1和E2數(shù)據(jù)庫的檢測準(zhǔn)確率大約在70%～77%之間；而Isomap方法對E1和E2數(shù)據(jù)庫的檢測準(zhǔn)確率在94%～98%之間，明顯高于PCA方法。2種方法的高維特征向量提取方法相同，造成檢測準(zhǔn)確率相差較大的原因可能是：由灰度值構(gòu)成的高維特征向量在高維空間不是分布在一個(gè)超平面上，而是分布在一個(gè)曲面流形上。PCA方法是一種線性降維方法，只能對分布在超平面上的數(shù)據(jù)進(jìn)行有效降維；對分布在曲面流形上的數(shù)據(jù)點(diǎn)，降維過程中可能將它們映射到低維空間中同一個(gè)點(diǎn)上，降低了類的可分性；而Isomap方法屬于非線性降維方法，在降維過程中可保持流形的全局幾何結(jié)構(gòu)，并且利用Ncut準(zhǔn)則進(jìn)行參數(shù)優(yōu)化，增加了低維特征向量在低維空間的類可分性。因此，使Isomap方法的識(shí)別精度高于PCA方法。Gabor方法、G-PCA方法和G-Isomap方法對2個(gè)數(shù)據(jù)庫在不同尺度下的識(shí)別準(zhǔn)確率都在90%以上，其中G-Isomap方法的檢測準(zhǔn)確率最高，在94.2%～98%之間。G-PCA方法最低，在90%～95%之間，Gabor方法居中，在93%～96.7%之間。3種方法檢測準(zhǔn)確率相差不大，主要原因在于Gabor簇對疵點(diǎn)圖片進(jìn)行不同尺度和方向?yàn)V波后，增加了疵點(diǎn)和正常區(qū)域的對比度，消除光照不均的影響，提高基于灰度共生矩陣特征值(對比度、能量、熵和相關(guān)性)的靈敏度。由這些特征值組成的高維特征向量自身對疵點(diǎn)區(qū)域具有較高靈敏度，且特征維數(shù)固定，不受子窗口尺度的影響，不會(huì)強(qiáng)化PNN神經(jīng)網(wǎng)絡(luò)輸入層和樣本層的泛化能力，降低樣本層和競爭層間的擬合精度，導(dǎo)致分類器分類精度下降。Gabor方法檢測準(zhǔn)確率比G-PCA方法高的原因可能是：高維數(shù)據(jù)在數(shù)據(jù)空間呈弱非線性，PCA在降維過程中沒有完整保留高維特征向量的空間幾何結(jié)構(gòu)和有效信息，使得低維特征向量靈敏度下降。G-Isomap方法使用Isomap方法降維，可克服高維特征向量非線性影響，并有效地從高維數(shù)據(jù)中剔除冗余信息，降低了PNN神經(jīng)網(wǎng)絡(luò)輸入層和樣本層間的計(jì)算難度，提高了樣本層和競爭層間的擬合精度，增加了分類器的分類準(zhǔn)確性。G-Isomap方法和Isomap方法對2個(gè)數(shù)據(jù)庫的檢測準(zhǔn)確率相差不大，都在94%～98%之間，說明Isomap方法降維能力很強(qiáng)，對特征向量維數(shù)不敏感。其次，G-Isomap方法中存在提取基于灰度共生矩陣特征值的環(huán)節(jié)，在計(jì)算時(shí)間上肯定要比Isomap方法長，不利于疵點(diǎn)的快速檢測。

表1 數(shù)據(jù)庫E1辨識(shí)正檢率Tab.1 Recognition results positive detection rate of database E1

表2 數(shù)據(jù)庫E2辨識(shí)正檢率Tab.2 Recognition results positive detection rate of database E2

圖2、3分別示出數(shù)據(jù)庫E1(針對竹節(jié)疵點(diǎn))、E2(針對粗緯疵點(diǎn))的匹辨識(shí)效果。由圖可知：子窗口尺寸在30像素×30像素情況下，PCA方法正常紋理的誤判位置過多，導(dǎo)致該方法的檢測準(zhǔn)確率較低；而Isomap方法誤檢位置數(shù)量明顯強(qiáng)于PCA方法。進(jìn)一步說明PCA方法不適用于基于灰度值的高維特征向量降維。在相同子窗口尺寸下，在Gabor、G-PCA和G-Isomap方法對疵點(diǎn)檢測結(jié)果中發(fā)現(xiàn)：3種方法均存在少量的誤檢及漏檢位置，其中G-PCA方法中的誤檢數(shù)量略多于Gabor方法和G-Isomap方法。而G-Isomap方法誤檢數(shù)量比Gabor方法少，說明Isomap在對灰度共生矩陣特征向量降維更為有效，其低維特征向量具有更高的類可分性；與Isomap方法相比，PCA在對疵點(diǎn)特征降維過程中并不具備優(yōu)勢。G-Isomap方法和Isomap方法的誤檢和錯(cuò)檢數(shù)目基本相同，均可實(shí)現(xiàn)準(zhǔn)確檢測，與表1和表2的統(tǒng)計(jì)結(jié)果一致。

圖2 數(shù)據(jù)庫E1針對竹節(jié)疵點(diǎn)的匹辨識(shí)效果Fig.2 Slab yarn identification effect of database E1

圖3 數(shù)據(jù)庫E2針對粗緯疵點(diǎn)的匹辨識(shí)效果Fig.3 Coarse pick identification effect of database E2

3 結(jié) 論

基于Gabor濾波器簇提出的高維特征向量呈現(xiàn)非線性，Isomap方法可克服非線性影響，剔除高維特征中冗余信息，弱化分類器的泛化性能，強(qiáng)化分類器的擬合性能，提高分類器的分類精度和疵點(diǎn)的檢測精度。利用Ncut進(jìn)行Isomap方法的參數(shù)選擇，可增加低維特征向量的靈敏度和類可分性，有助于疵點(diǎn)檢測精度的提高。基于濾波圖像子窗口灰度值構(gòu)造的高維特征經(jīng)過Isomap方法降維后獲得的低維特征，在疵點(diǎn)檢測中具有較高的靈敏度。

FZXB

[1] 管聲啟, 石秀華. 基于頻域?yàn)V波的織物疵點(diǎn)檢測[J].計(jì)算機(jī)應(yīng)用, 2008, 28(10):2673-2675. GUAN Shengqi，SHI Xiuhua. Fabric defects detection based on frequency domain filtering[J]. Journal of Computer Application, 2008, 28(10):2673-2675.

[2] 楊曉波. 基于自適應(yīng)離散小波變換的混合特征畸變織物疵點(diǎn)識(shí)別[J]. 紡織學(xué)報(bào)，2013, 34(1):133-137. YANG Xiaobo. Research of mixture feature aberrance fabric defect recognition based on self-adaptive discrete wavelet transform[J]. Journal of Textile Research, 2013, 34(1):133-137.

[3] 李立輕. 基于計(jì)算機(jī)視覺的織物疵點(diǎn)自動(dòng)檢測研究[D]. 上海：東華大學(xué), 2003:2-4. LI Liqing. Study on computer vision-based automatic detection of fabric defects[D]. Shanghai: Donghua University, 2003:2-4.

[4] 槐向兵, 厲征鑫, 劉建立, 等. 基于輪廓波變換的織物疵點(diǎn)圖像消噪新方法[J]. 計(jì)算機(jī)工程與應(yīng)用, 2014, 50(19):143-146. HUAI Xiangbing, LI Zhengxin, LIU Jianli, et al．New fabric defect image denoising method based on contourlet transform[J]．Computer Engineering and Applications, 2014,50(19):143-146.

[5] 屈博, 盧朝陽, 李靜, 等. 一種改進(jìn)的多通道Gabor濾波器布匹瑕疵檢測方法[J].紡織學(xué)報(bào)，2009, 30(12):37-40. QU Bo, LU Zhaoyang, LI Jing, et al. An improved multichannel Gabor filter algorithm for fabric defect detection[J]. Journal of Textile Research, 2009, 30(12):37-40.

[6] JING Junfeng, YANG Panpan, LI Pengfei, et al．Supervised defect detection on textile fabrics via optimal Gabor filter[J]. Journal of Industrial Textiles, 2013, 44(1):40-57

[7] ZHENG Zhonglong，F(xiàn)AN Yang，TAN Wenan, et al. Gabor feature-based face recognition using supervised locality preserving projection[J]. Signal Processing, 2007, 87(10):2473-2483

[8] YAN Deqin, LIU Shenglan, LI Yanyan. An embedding dimension reduction algorithm based on sparse anal-ysis[J]. Acta Automatica Sinica, 2011, 37(11): 1306-1312.

[9] AJAY Kumar, GRANTHAM Pang. Defect detection in textured materials using Gabor filters[C]. Industry Applications IEEE Transactions on, 2002, 38(2):425-440.

[10] WANG L，ZHANG Y，F(xiàn)ENG J. On the Euclidean distance of images[J]. IEEE Trans on Pattern Analysis and Machine Intelligence,2005,27(8):1334-1339.

[11] LIU Yan, GU Xueping. Skeleton-network reconfiguration based on topological characteristics of scale-free networks and discrete particle swarm optimization [J]. IEEE Transactions on Power Systems, 2007, 22(3):1267-1274.

[12] ALBERTO García-Villoria, RAFAEL Pastor. Introducing dynamic diversity into a discrete particle

swarm optimization [J]. Computer and Operation Research, 2009, 36(3):951-966.

[13] 胡峰, 蘇訊, 劉偉, 等. 基于改進(jìn)局部線性嵌入算法的故障特征提取方法[J]. 振動(dòng)與沖擊，2015, 34(15):201-204. HU Feng, SU Xun, LIU Wei, et al. Fault feature extraction based on improved locally linear embe-dding [J]. Journal of Vibration and Shock, 2015, 34(15):201-204.

[14] 胡峰, 王傳桐, 吳雨川, 等. 基于改進(jìn)監(jiān)督LLE算法的故障特征提取方法[J]. 振動(dòng)與沖擊，2015(21):119-123. HU Feng, WANG Chuantong, WU Yuchuan, et al. Fault features extraction based on improved supervised locally linear embedding [J] Journal of Vibration and Shock, 2015(21):119-123.

[15] LIANG Dong, YANG Jie, ZHENG Zhonglong. A facial expression recognition system based on supervised locally linear embedding[J]. Pattern Recognition Letters, 2005, 26(15):2374-2389.

[16] 史峰, 王小川, 郁磊, 等. MatLab神經(jīng)網(wǎng)絡(luò)30個(gè)案例分析[M]. 北京:北京航空航天大學(xué)出版社,2010：197-182. SHI Feng, WANG Xiaochuan,YU Lei, et al. 30 Case Analysis of MatLab Neural Network[M]. Beijing: Beijing University of Aeronautics and Astronautics Press, 2010：197-182.

Detection of fabric defects based on Gabor filters and Isomap

WANG Chuantong, HU Feng， XU Qiyong， WU Yuchuan， YU Lianqing

(SchoolofMechanicalEngineeringandAutomation,WuhanTextileUniversity,Wuhan,Hubei430074,China)

In order to improve the correction rate of fabric defects, Gabor filters and Isomap were used to detect the fabric defect. Firstly, the images of fabric defect were filtered by 15 Gabor filters with 3 orientations and 5 scales, which contributed to overcome the effect of uneven illumination and low contrast. Then, the filtered images were divided into non-overlapping rectangular patches and high-dimensional features were extracted. Simultaneously, Isomap algorithm was applied to reduce the dimensionality of feature and eliminate the redundant information. Besides, a mapping model of new samples was proposed to detect the low dimensional embedding results. Finally, the performance of proposed algorithm was estimated off-line by two sets of fabric defect images. The theoretical argument is supported by experimental results.

fabric defect detection； Gabor filter； Isomap； image processing

2016-03-25

2016-12-13

國家自然科學(xué)基金資助項(xiàng)目(51205294；61271008；51275363)

王傳桐(1991—)，男，碩士生。主要研究方向?yàn)榧徔椘焚|(zhì)量檢測和模式識(shí)別。胡峰，通信作者，E-mail：hufeng@wtu.edu.cn。

10.13475/j.fzxb.20160305306

TP 311.131； TS 101.97

A