周國慶，吳　錫，袁汝旺，蔣秀明（天津工業(yè)大學(xué)天津市現(xiàn)代機(jī)電裝備技術(shù)重點(diǎn)實(shí)驗(yàn)室，天津　300387）

?

基于線陣CCD的紗線核直徑提取算法

周國慶，吳錫，袁汝旺，蔣秀明
（天津工業(yè)大學(xué)天津市現(xiàn)代機(jī)電裝備技術(shù)重點(diǎn)實(shí)驗(yàn)室，天津300387）

摘要：針對紗線直徑測量過程中核直徑提取算法的不足，提出了基于線陣CCD測量系統(tǒng)的紗線核直徑提取的斜率閾值算法，分析了斜率閾值大小和擬合點(diǎn)數(shù)對算法準(zhǔn)確性和快速性的影響，并且建立了適用于紗線在線測量斜率閾值算法中的最佳擬合模型.仿真及實(shí)驗(yàn)結(jié)果表明：擬合點(diǎn)數(shù)n = 3、擬合組不同點(diǎn)數(shù)s = 3、斜率閾值KT= 45時得到的紗線平均直徑與理論直徑較為接近，得到的紗線變異系數(shù)與烏斯特公報基本吻合.該提取算法適用于紗線在線檢測，為以后紗線質(zhì)量在線檢測提供了新的思路.

關(guān)鍵詞：紗線核直徑；線陣CCD；灰度值；斜率閾值

紗線質(zhì)量的高低直接影響到紡織品的外觀、透氣性和強(qiáng)度，同時對織物對染料的吸附力產(chǎn)生影響.所以紗線質(zhì)量指標(biāo)準(zhǔn)確、快速的測定對于紗線生產(chǎn)企業(yè)的工藝調(diào)控、質(zhì)量控制有著重要的現(xiàn)實(shí)意義[1].其中紗線直徑是影響紗線外觀質(zhì)量的一個重要因素，為準(zhǔn)確提取紗線直徑，盛國俊等[2]提出了K階矩濾波和最小熵濾波法來提取紗線直徑信號，避免了紗線毛羽的影響，但特征信號在提取的過程中產(chǎn)生了一定的損失，影響測量的準(zhǔn)確性；Ohta和Ozkaya等[3-4]利用數(shù)字圖像處理技術(shù)分析不同實(shí)驗(yàn)條件對測量效果的影響；Carvalho[5]開發(fā)了紗疵評價系統(tǒng)，獲取紗線直徑、線密度等參數(shù)，但該系統(tǒng)的設(shè)備較為昂貴，算法較復(fù)雜，周期較長.同時袁汝旺等[6]提出了紗線線輪廓信號的遞歸平均濾波法，避免了紗線毛羽的干擾；程立超等[7]提出了輪廓坐標(biāo)跟蹤法，通過設(shè)定浮動閾值線范圍，可以準(zhǔn)確測量紗線直徑.

隨著線性CCD技術(shù)的發(fā)展，線性陣列測量技術(shù)已逐步用于紗線外觀直徑的檢測[8].針對以前檢測系統(tǒng)周期長、算法難實(shí)現(xiàn)等問題，本文依據(jù)線性陣列的測量原理，提出了斜率閾值算法，分析了斜率閾值大小和擬合模型對算法準(zhǔn)確性和快速性的影響，通過仿真和實(shí)驗(yàn)，驗(yàn)證了該提取算法在紗線在線檢測中的可行性，為以后紗線質(zhì)量在線檢測提供了理論支撐.

1　紗線直徑測量原理

紗線結(jié)構(gòu)成分包括纖維和纖維間的空氣間隙，紗線在受到外部摩擦?xí)r，其表層纖維就有可能伸出從而成為毛羽，通常所測的紗線直徑是在一定膨松條件下的直徑值[9-10].圖1為顯微鏡下紗線外觀結(jié)構(gòu)，沿著紗線長度方向做垂直線ab，本文定義紗線ab部分為紗線核直徑，ab外側(cè)突出部分為紗線毛羽.

圖1　紗線外觀結(jié)構(gòu)Fig.1 Appearance of yarn structure

圖2　線陣CCD測量原理Fig.2 Principle of linear CCD measuring

依據(jù)紗線核直徑和紗線毛羽透光度的不同，線陣CCD傳感器檢測紗線原理圖如圖2所示.

通過光學(xué)系統(tǒng)將紗線直徑尺寸以一定的、準(zhǔn)確的倍率成像于線性陣列光敏面上[11-12]，對線性陣列輸出信號進(jìn)行A/D轉(zhuǎn)換，送往計算機(jī)，實(shí)現(xiàn)相應(yīng)的測量，式（1）給出了紗線直徑測量的基本模型[13].

式中：d為紗線直徑（mm）；h為單個線性陣列像素寬度（mm）；ΔP為紗線遮擋的線性陣列像素個數(shù)；β為光學(xué)系統(tǒng)的放大倍數(shù).

2　紗線核直徑提取算法

2．1斜率閾值法原理

圖3描述了紗線直徑線輪廓特征像素與灰度值的關(guān)系，由紗線直徑測量的模型，核直徑提取的關(guān)鍵在于準(zhǔn)確快速的獲取紗線遮擋線陣像素的個數(shù)ΔP.其中灰度值表示接收到光線的明暗程度；i為斜率閾值法得出最低點(diǎn)左側(cè)的第i個點(diǎn)；j為斜率閾值法得出最低點(diǎn)右側(cè)的第j個點(diǎn).

基于紗線信號的邊緣特性，紗線信號陰影邊界呈現(xiàn)斜坡狀態(tài)[14]，本文提出斜率閾值法.斜率閾值法的核心是基于紗線直徑線輪廓特征，首先設(shè)定斜率閾值的范圍，其次從紗線線輪廓最低點(diǎn)沿線輪廓左右兩側(cè)進(jìn)行尋點(diǎn)擬合，最終當(dāng)擬合直線的斜率達(dá)到設(shè)定的范圍時，確定紗線遮擋線陣像素的個數(shù)△P.

圖3　紗線直徑線輪廓特征曲線Fig.3 Characteristic curve of line profile on yarn diameter

設(shè)紗線線輪廓最低點(diǎn)對應(yīng)的坐標(biāo)為（Pmin，Dmin）；每組尋點(diǎn)擬合的點(diǎn)數(shù)為n；各擬合組之間不同點(diǎn)數(shù)為s （s＞0）；Ki、Kj分別為最低點(diǎn)左右兩側(cè)擬合直線的斜率，i，j?。?，s，2s，3s‥）.由最小二乘法得出Ki、Kj的表達(dá)式為：

2．2斜率閾值法算法實(shí)現(xiàn)

利用斜率閾值法提取紗線直徑的具體步驟如下：

（1）設(shè)定紗線斜率閾值KT；

（2）獲取紗線直徑線輪廓信號；

（3）找出線輪廓信號中灰度值最小點(diǎn)對應(yīng)的坐標(biāo)（Pmin，Dmin）；

（4）從（Pmin，Dmin）點(diǎn)沿紗線輪廓向左右兩側(cè)進(jìn)行尋點(diǎn)擬合；

（5）當(dāng)擬合直線斜率大于等于斜率閾值時，尋找停止；

（6）確定紗線遮擋線陣像元的個數(shù)△P，進(jìn)而得出紗線直徑值.

下面以從最低點(diǎn)沿線輪廓左側(cè)為例進(jìn)行說明斜率閾值法的流程：

設(shè)Pmin左側(cè)的像素坐標(biāo)依次為Pmin-1，…，Pmin-i，Pmin-i-1對應(yīng)的灰度值為Dmin-1，…，Dmin-i，Dmin-i-1.從Pmin點(diǎn)向左沿輪廓曲線尋點(diǎn)，第一組選取Pmin，Pmin-1，Pmin-2，…點(diǎn)，得出第一組擬合斜率K1，然后選取Pmin-n+1，Pmin-n，Pmin-n-1，…點(diǎn)得出其擬合斜率為K2，依次類推第1組斜率Kl，當(dāng)Kl≥KT時尋點(diǎn)停止.此時按照邊界提取原理，其第l組擬合點(diǎn)數(shù)為（Pmin-i+1，Pmin-i，…）取左側(cè)邊界PL=Pmin-i.同理，沿Pmin向右尋點(diǎn)時設(shè)Pmin右側(cè)的像素坐標(biāo)依次為Pmin+1，…，Pmin+j，Pmin+j+1，找到第r組斜率Kr≥KT時尋找停止.取PR= Pmin+j，則△P = PR- PL，然后由式（1）得紗線直徑：

3　實(shí)驗(yàn)研究

3.1實(shí)驗(yàn)條件

紗線直徑測量系統(tǒng)主要有歐姆龍ZX-GT型智能激光傳感器、樣本喂入速度張力控制系統(tǒng)、數(shù)據(jù)分析系統(tǒng)等組成.

選用的實(shí)驗(yàn)研究樣本為31.1 tex的環(huán)錠紡紗線，光學(xué)放大倍數(shù)β= 1，采樣時間間隔為10 ms，張力T = 50 mN，紗線的喂入速度v = 30 m/min.

由線密度Tt和紗線理論直徑D0的關(guān)系式[15]：

式中：δ為紗線體積質(zhì)量（g/cm3）.計算出31.1 tex的環(huán)錠紡紗線理論直徑D0= 0.234 mm.

3.2實(shí)驗(yàn)討論

在從線輪廓最低點(diǎn)沿紗線線輪廓左右兩側(cè)進(jìn)行尋點(diǎn)擬合的過程中，由式（2）得出，每組擬合點(diǎn)數(shù)n和各擬合組之間不同點(diǎn)數(shù)s的大小都會對紗線直徑的準(zhǔn)確性和算法的快速性產(chǎn)生影響.

其中圖4描述的是s = 0和s = 3時直線擬合的示意圖，圖5描述了n = 5時進(jìn)行直線擬合的示意圖.

當(dāng)n取值較大時，其達(dá)到設(shè)定閾值斜率的速度較快，但在選取遮擋邊界時產(chǎn)生的誤差也較大.當(dāng)s取值較小時其提取的精度很高，但同時其達(dá)到斜率閾值的速度較慢.

表1和表2分別為n和s取不同值時所對應(yīng)的紗線直徑.

圖4　s=0和s=3擬合方式示意圖Fig.4 Fitting method of s=0 and s=3

圖5　五點(diǎn)擬合方式示意圖Fig.5 Fitting method of using five points

表1　不同擬合點(diǎn)數(shù)n對應(yīng)的直徑值Tab.1 Corresponding value of n of different fitting points

表2　擬合組不同點(diǎn)s對應(yīng)的直徑值Tab.2 Value of s in different points of fitting group

為了滿足紗線在線檢測的要求，通過與理論直徑之間的對比，最終選擇n = 3、s = 3以保證算法精度和運(yùn)算速度.

為了研究斜率閾值的合理范圍，選用n = s = 3，得出不同斜率閾值對應(yīng)的紗線直徑如表3所示.

由表3中可以看出，隨著閾值斜率KT的增加，紗線直徑在逐漸增大，當(dāng)設(shè)定閾值KT＞60時，其直徑會發(fā)生較大變化，通過多組數(shù)據(jù)對比，結(jié)合31.1 tex環(huán)錠紡紗線的理論直徑，最終取KT= 45.

表3　不同斜率閾值KT對應(yīng)的直徑值Tab.3 Diameters at different slope threshold

3.3實(shí)驗(yàn)驗(yàn)證

為了驗(yàn)證選用的擬合模型和斜率閾值的設(shè)定范圍是否符合紗線在線檢測的要求，基于系統(tǒng)不失真的條件是紗線樣本不規(guī)則的最小波長為線陣接收單元寬度的10倍[16].按照10個樣本點(diǎn)為1個樣本子片段，10個子片段組成1個大樣本，圖6表示的是整個采樣過程中紗線直徑像素的變化情況，表4為子片段直徑測量情況.

圖6　紗線直徑測量信號Fig.6 Measured signal of yarn diameter

表4　樣本子片段直徑測量結(jié)果Tab.4 Measuring result on diameter of sample sub-segment

通過計算得出，斜率閾值法測得紗線平均直徑為0.241 mm，變異系數(shù)為15.38%.通過查烏斯特公報2012數(shù)據(jù)，所提出的斜率閾值法在滿足紗線在線檢測速度的同時，符合紗線檢測精度的要求.

4　結(jié)論

（1）通過實(shí)驗(yàn)和仿真分析，選取擬合點(diǎn)數(shù)n = 3，擬合組不同點(diǎn)數(shù)s = 3，斜率閾值KT= 45的擬合方式得到子片段直徑的變化比較規(guī)則，滿足紗線在線檢測快速性、準(zhǔn)確性的要求.

（2）通過仿真和實(shí)驗(yàn)結(jié)果表明斜率閾值算法測得的紗線平均核直徑為0.241 mm和紗線變異系數(shù)為15.38%，與烏斯特公報數(shù)據(jù)基本吻合.

（3）基于線陣CCD測量系統(tǒng)的斜率閾值算法算法較為簡單，運(yùn)算周期較短，提高了檢測的效率，為以后紗線質(zhì)量在線檢測提供了新的思路.

參考文獻(xiàn)：

[1] TàPIAS M，RALLó M，ESCOFET J. Automatic measurements of partial cover factors and yarn diameters in fabrics using image processing[J]. Textile Research Journal，2011，81（2）：173-186.

[2]盛國俊，董永貴.紗線質(zhì)量光電信號的消噪處理和異常檢測[J].清華大學(xué)學(xué)報：自然科學(xué)版，2010，50（2）：229-231. SHENG G J，DONG Y G. The optical signal quality yarn denoising and anomaly detection[J]. Journal of Tsinghua University：Natural Science Edition，2010，50（2）：229-231（in Chinese）.

[3] OHTA K，HATA Y，OHTA C，et al. Measuring the appearance diameter of yarns using visual information [J]. Sen -i Gakkaishi，2000，56（11）：518-523.

[4] OZKAYA Y A，ACAR M，JACKSON M R . Yarn twist measurement using digital imaging [J]. Journal of the Textile Institute，2008，101（2）：91-100.

[5] CARVALHO V H，CARDOSO P J，BELSLEY M S，et al. Yarn hairiness characterization using two orthogonal directions [J]. Instrumentation and Measurement，IEEE Transactions on，2009，58（3）：594-601.

[6]袁汝旺，蔣秀明，周國慶，等.基于線性陣列的紗線直徑與毛羽測量方法[J].紡織學(xué)報，2013，34（8）：132-137. YUAN R W，JIANG X M，ZHOU G Q，et al. The yarn linear array and the diameter of hairinessmeasurement method based on line array CCD[J]. Journal of Textile Research，2013，34 （8）：132-137（in Chinese）.

[7]程立超，蔣秀明，袁汝旺，等.基于線陣電荷耦合元件的紗線直徑測量方法[J].紡織學(xué)報，2015，36（6）：124-128. CHENG L C，JIANG X M，YUAN R W，et al. Measuring method of yarn diameter based on line arraycharge coupled device [J]. Journal of Textile Research，2015，36（6）：132-137 （in Chinese）.

[8]涂文驍，左保齊.線狀物直徑檢測方法[J].現(xiàn)代絲綢科學(xué)與技術(shù)，2013，28（2）：72-76. TU W X，ZUO B Q. Method for measuring diameter of thread [J]. Modern Silk Science and Technology，2013，28（2）：72-76 （in Chinese）.

[9] KOGANEYA K，YUKISHITA Y，F(xiàn)UJISAKI H，et al. Instrument for measuring the appearance width of running double ply staple yarn using lightness threshold for yarn image[J]. Sen-I Gakkaishi，2012，68（4）：98-105.

[10] HALEEM N，WANG X. Recent research and developments on yarn hairiness[J]. Textile Research Journal，2015，85（2）：211-224.

[11]隋修武，劉蕾，杜玉紅，等.基于線陣CCD的可吸收縫合線線徑在線檢測系統(tǒng)[J].傳感技術(shù)學(xué)報，2014，27（8）：1154-1158. SUI X W，LIU L，DU Y H，et al. Online measurement system of absorbable suture diameter based on linear CCD[J]. Chinese Journal of Sensors and Actuators，2014，27（8）：1154-1158 （in Chinese）.

[12]于之靖，劉月林，諸葛晶昌.高精度雙線陣CCD非接觸直徑測量系統(tǒng)[J].傳感器與微系統(tǒng)，2014，33（1）：120-122. YU Z J，LIU Y L，ZHUGE J J. Non-contactdiameter measurement system of high precision doublelinear array CCD[J]. Sensor and Micro System，2014（1）：120-122（in Chinese）.

[13]蘇波，王紀(jì)龍. CCD高精度測徑系統(tǒng)的研究[J].太原理工大學(xué)學(xué)報，2002，33（5）：506-509. SU B，WANG J L . High precision CCD diameter measurement system[J]. Journal of Taiyuan University of Technology，2002，33（5）：506-509（in Chinese）.

[14]蔣志年，易其順.基于CCD的非接觸式直徑測量實(shí)驗(yàn)的數(shù)據(jù)處理[J].廣西民族師范學(xué)院學(xué)報，2013，30（3）：22-24. JIANG Z N，YI Q S. Measurement of noncontact CCD diameter based on data processing[J]. Journal of Guangxi Normal University for Nationalities，2013，30（3）：22-24（in Chinese）.

[15]趙領(lǐng)航，蔡普寧，賈哲昆，等.金屬纖維及其混紡紗線直徑系數(shù)的計算及測試[J].產(chǎn)業(yè)用紡織品，2015，33（8）：41-44. ZHAO L H，CAI P N，JIA Z K，et al. Diameter coefficient of metal fiber and its blended yarncalculation and testing [J]. Industrial Textiles，2015，33（8）：41-44（in Chinese）.

[16]袁汝旺，蔣秀明，楊公源.纖維束條干外觀不勻的激光測量方法[J].中國激光，2012，39（7）：176-180. YUAN R W，JIANG X M，YANG G Y. Laser measurement method for appearance unevenness offiber budles[J]. China Laser，2012，39（7）：176-180（in Chinese）.

Algorithm of yarn core diameter extraction based on linear array CCD

ZHOU Guo-qing，WU Xi，YUAN Ru-wang，JIANG Xiu-ming
（Key Laboratoryof Advanced Mechatronics Equipment Technology，Tianjin Polytechnic University，Tianjin300387，China）

Abstract：Aiming at the shortage of yarn diameter extraction algorithm in the measurement process，a slope threshold algorithm based on linear array CCD measuring system is put forward to extract core diameter，and the influence of slope threshold degree and fitting points on accuracy and rapidity of the algorithm is analyzed，The best fitting model of the slope threshold algorithm for yarn online measuring is established. The simulation and experimental results show that the average diameter is closer to the theoretical diameter when the fitting points n is 3，fitting group with different number s is 3 and the slope threshold KTis 45，and the coefficient of variation is consistent with the USTER communique. The algorithm is suitable for on-line detection，which provides a new method for the online detection of yarn quality.

Key words：yarn core diameter；linear array CCD；grey level；slope threshold

通信作者：周國慶（1958—），男，教授，主要研究方向?yàn)樾滦图徔棛C(jī)械設(shè)計研究. E-mail：zhouguoqing@tjpu.edu.cn

基金項(xiàng)目：國家重點(diǎn)基礎(chǔ)研究發(fā)展973計劃（2010CB334711）；國家自然科學(xué)基金資助項(xiàng)目（51205288）；天津市科技計劃項(xiàng)目（13JCYBJC15900）

收稿日期：2015-12-07

DOI：10.3969/j.issn.1671-024x.2016.02.007

中圖分類號：TS 103.6；TS101.1

文獻(xiàn)標(biāo)志碼：A

文章編號：1671－024X（2016）02－0035－05