吳文強(qiáng) 李 軍 陳 龍 金 勇 劉 斌 鐘科軍 毛偉俊 王藝斌 許文武

(1. 湖南中煙工業(yè)有限責(zé)任公司技術(shù)中心，湖南 長(zhǎng)沙 410007；2. 南京焦耳科技有限責(zé)任公司，江蘇 南京 210000)

煙絲結(jié)構(gòu)是影響卷煙物理指標(biāo)的主要因素之一，不同結(jié)構(gòu)的煙絲對(duì)卷煙的質(zhì)量和穩(wěn)定性都有一定的影響，對(duì)卷煙的煙支重量、吸阻[1]、端部落絲量等質(zhì)量指標(biāo)[2]也有顯著影響。目前煙草企業(yè)對(duì)煙絲結(jié)構(gòu)的檢測(cè)仍普遍采用煙絲振動(dòng)分選篩，這種方法采樣周期長(zhǎng)，數(shù)據(jù)實(shí)時(shí)性差，而且受人為影響較多。近年來有學(xué)者通過圖像識(shí)別技術(shù)來檢測(cè)煙絲結(jié)構(gòu)，楚晗等[3]通過對(duì)比同檔次短支卷煙和常規(guī)卷煙的煙絲結(jié)構(gòu)分布并對(duì)煙絲結(jié)構(gòu)與物理指標(biāo)進(jìn)行回歸分析，研究了煙絲結(jié)構(gòu)分布及其對(duì)物理質(zhì)量的影響差異；余娜等[4]采用分形理論，推導(dǎo)了煙絲篩下積累質(zhì)量百分比與分形維數(shù)的相關(guān)性方程；夏營(yíng)威等[5]通過計(jì)算機(jī)視覺系統(tǒng)，拍攝煙絲的可見光圖像，并對(duì)圖像進(jìn)行畸變校正，提出了煙絲切口位置識(shí)別算法，實(shí)現(xiàn)了無規(guī)則擺放煙絲的寬度測(cè)量；郭三刺等[6]采用中值濾波對(duì)煙絲圖像進(jìn)行預(yù)處理，然后再用邊緣提取方法計(jì)算煙絲的各個(gè)特征值，適用于離線檢測(cè)煙絲結(jié)構(gòu)；劉洋等[7]采用CCD掃描成像式對(duì)煙絲寬度進(jìn)行快速檢測(cè)；邸成良等[8]提出一種基于圖像處理的自適應(yīng)煙絲寬度測(cè)量方法，該方法具有很強(qiáng)的自適應(yīng)性和抗干擾性，可以借鑒到煙絲結(jié)構(gòu)的在線測(cè)量方面。但是這些研究都是針對(duì)煙絲結(jié)構(gòu)[9]離線檢測(cè)的，未見在線檢測(cè)的相關(guān)文獻(xiàn)報(bào)道。

研究擬采用風(fēng)選攤薄方式與圖像智能識(shí)別技術(shù)，設(shè)計(jì)開發(fā)在線煙絲結(jié)構(gòu)檢測(cè)系統(tǒng)，以期能夠準(zhǔn)確快速地檢測(cè)煙絲長(zhǎng)度、寬度、面積，并統(tǒng)計(jì)出煙絲中整絲、碎絲比例。

1 系統(tǒng)組成

如圖1所示，在線煙絲結(jié)構(gòu)檢測(cè)系統(tǒng)主要包括煙絲輸送皮帶機(jī)、風(fēng)選室、圖像采集系統(tǒng)、圖像處理系統(tǒng)和煙絲結(jié)構(gòu)在線檢測(cè)軟件。輸送皮帶機(jī)B1烘干后的煙絲輸送至風(fēng)選室，風(fēng)選室中有兩道氣動(dòng)吹散裝置，分別是氣動(dòng)吹散裝置A1和氣動(dòng)吹散裝置A2，實(shí)現(xiàn)煙絲的多次吹散。第一道氣動(dòng)吹散裝置A1把輸送皮帶機(jī)B1來料中的單獨(dú)的煙絲吹落到輸送皮帶機(jī)B2上，煙絲團(tuán)落到輸送皮帶機(jī)B3上，輸送皮帶機(jī)B3上的煙絲和煙絲團(tuán)經(jīng)過氣動(dòng)吹散裝置A2，將獨(dú)立的吹散的輕的煙絲落到輸送皮帶機(jī)B2上，煙絲團(tuán)落到輸送皮帶機(jī)B3上。經(jīng)過兩次吹散后，吹散的煙絲落到輸送皮帶機(jī)B2上，進(jìn)入下道工序。經(jīng)過兩次吹散后，依然落下來的是梗簽、焦油塊等雜物，進(jìn)入廢料收集盒。圖像采集系統(tǒng)采集煙絲在吹散狀態(tài)下的煙絲圖像(也就是輸送皮帶機(jī)B2上方的煙絲圖像)。如圖2所示，圖像處理系統(tǒng)對(duì)采集的煙絲圖像進(jìn)行預(yù)處理，通過煙絲結(jié)構(gòu)檢測(cè)軟件準(zhǔn)確計(jì)算出每根煙絲的長(zhǎng)度、面積，進(jìn)行質(zhì)量擬合，建立面積與質(zhì)量的關(guān)系模型；同時(shí)根據(jù)這些數(shù)據(jù)統(tǒng)計(jì)出煙絲中整絲、碎絲比例；計(jì)算出整絲率、碎絲率，得出煙絲結(jié)構(gòu)參數(shù)。

2 技術(shù)實(shí)現(xiàn)

在線煙絲結(jié)構(gòu)檢測(cè)系統(tǒng)中煙絲攤薄、煙絲圖像檢測(cè)所涉及的技術(shù)是研究的核心。

2.1 煙絲攤薄機(jī)構(gòu)設(shè)計(jì)

風(fēng)選室中氣動(dòng)吹散裝置A1將輸送皮帶機(jī)B1勻速定量輸送煙絲團(tuán)經(jīng)行吹散攤薄，形成相互獨(dú)立的狀態(tài)，重量輕的煙絲被吹落在輸送皮帶機(jī)B2上，部分煙絲、重量重的煙絲團(tuán)、梗簽及雜物、生產(chǎn)線上產(chǎn)生的焦油塊下落到輸送皮帶機(jī)B3上，輸送皮帶機(jī)B3運(yùn)行方向與輸送皮帶機(jī)B2相反，再經(jīng)由氣動(dòng)吹散裝置A2再次進(jìn)行吹散攤薄，部分煙絲團(tuán)被吹散至相互獨(dú)立的完全攤開的狀態(tài)，在風(fēng)力作用下最終落入輸送皮帶機(jī)B2，無法吹散的煙絲團(tuán)、有一定重量的梗簽及雜物、生產(chǎn)線上產(chǎn)生的焦油塊落入下方廢料收集盒。而煙絲的造碎情況可以通過后續(xù)試驗(yàn)重復(fù)測(cè)量同一批次煙絲的碎絲率進(jìn)行比較得到，如果設(shè)備在線計(jì)算的碎絲率與國標(biāo)方法統(tǒng)計(jì)結(jié)果差別不大，并且重復(fù)測(cè)量碎絲率的變異系數(shù)較小，則可以判定該設(shè)備的造碎情況較小。而煙絲的攤薄效果可以通過統(tǒng)計(jì)煙絲的成團(tuán)率得到，成團(tuán)率越小則煙絲的攤薄效果也就越好。

1.輸送皮帶機(jī)B1 2.風(fēng)選室 3.圖像采集系統(tǒng) 4.圖像處理軟件 5.光源 6.氣動(dòng)吹散裝A1 7.氣動(dòng)吹散裝置A2 8.輸送皮帶機(jī)B2 9.輸送皮帶機(jī)B3 10.廢料收集盒圖1 在線煙絲結(jié)構(gòu)檢測(cè)系統(tǒng)結(jié)構(gòu)圖Figure 1 Structure diagram of online cut tobacco structure detecting system

圖2 圖像處理系統(tǒng)Figure 2 Image processing system

2.2 檢測(cè)單元的設(shè)計(jì)及原理

煙絲圖像檢測(cè)單元主要包括光源、高速CCD線陣相機(jī)、采集卡和工控機(jī)。為了滿足大視場(chǎng)高分辨率的要求，采用DALSA P4-CC-08K050-00-R線陣CCD相機(jī)，該相機(jī)有8 192個(gè)感光單元，像元尺寸達(dá)到7.04 μm×7.04 μm，當(dāng)鏡頭放大倍率為0.7時(shí)，檢測(cè)裝置的最小測(cè)量值為0.01 mm，可滿足煙絲長(zhǎng)度的測(cè)量要求。圖像檢測(cè)單元集成了煙絲檢測(cè)識(shí)別算法、煙絲結(jié)構(gòu)檢測(cè)算法。該單元的工作流程：采集的煙絲在吹散狀態(tài)下的煙絲圖像(輸送皮帶機(jī)B2上方的煙絲圖像)經(jīng)檢測(cè)識(shí)別算法處理后，識(shí)別出的煙絲目標(biāo)區(qū)域，對(duì)確定的煙絲目標(biāo)區(qū)域經(jīng)煙絲結(jié)構(gòu)檢測(cè)算法計(jì)算出該目標(biāo)的長(zhǎng)度，并依據(jù)事先標(biāo)定出的煙絲質(zhì)量模型，計(jì)算出待測(cè)樣品中的整絲和碎絲的質(zhì)量，同時(shí)確定整絲率和碎絲率，得到煙絲的結(jié)構(gòu)參數(shù)。

2.3 圖像處理算法

由于煙絲的測(cè)量精度要求過高，采用大分辨率線陣CCD拍攝煙絲圖像時(shí)，得到煙絲圖像較大(80 mm×80 mm 的圖像，每幅圖像大小約為80 MB)，因此需要對(duì)煙絲圖像進(jìn)行處理后才能快速準(zhǔn)確地計(jì)算煙絲的長(zhǎng)度。處理過程為：首先將圖像進(jìn)行預(yù)處理，中值濾波去除噪聲，然后閾值化處理將彩色圖像轉(zhuǎn)換為二值圖像，進(jìn)行連通區(qū)域標(biāo)記并提取煙絲連通區(qū)域，得到單根煙絲的圖像，最后將單根煙絲進(jìn)行細(xì)化、去毛刺(分叉點(diǎn))測(cè)量煙絲的長(zhǎng)度。另外，根據(jù)標(biāo)記的連通區(qū)域的面積判定成團(tuán)煙絲，若煙絲判定為成團(tuán)，還需要對(duì)成團(tuán)煙絲進(jìn)行分割。

(1)圖像預(yù)處理。采集到煙絲圖像如圖3所示，由于煙絲的結(jié)構(gòu)不規(guī)則，需要對(duì)圖像進(jìn)行預(yù)處理去除噪聲。去噪方法選用中值濾波法[10]，采用3×3的濾波窗口，以盡可能地保存圖像的細(xì)節(jié)信息。由于煙絲圖像中的背景部分存在略微明顯的白色噪點(diǎn)，采用中值濾波可以取得較好的抑制效果，見圖4。

(2)利用閾值分割對(duì)煙絲進(jìn)行二值化，提取煙絲連通區(qū)域[11]。為標(biāo)記圖像中的不同煙絲，對(duì)分割后的二值圖像進(jìn)行連通提取和標(biāo)記操作。常見的連通域提取方法有：線標(biāo)記法[12]、區(qū)域增長(zhǎng)法[13]、像素標(biāo)記法[14]、基于特殊體系的并行標(biāo)記算法[15]等。采用廣度優(yōu)先的連通域檢測(cè)方法[16]獲取煙絲連通域，并使用不同的灰度值標(biāo)記連通操作得到的不同區(qū)域。標(biāo)記效果如圖5所示。

圖3 煙絲圖像圖Figure 3 Image of cut tobacco

圖4 圖像預(yù)處理Figure 4 Image pre-processing

圖5 連通域標(biāo)記Figure 5 Label the connected areas

(3)對(duì)于成團(tuán)的煙絲，通過提取每個(gè)煙絲連通區(qū)域內(nèi)的面積，因?yàn)橐话銌蝹€(gè)煙絲的面積都比較小，只有當(dāng)煙絲成團(tuán)時(shí)，輪廓區(qū)域面積較大，通過設(shè)定面積閾值，來判斷連通區(qū)域是否屬于成團(tuán)煙絲，見圖6。若被判定為成團(tuán)煙絲，就需要對(duì)煙絲進(jìn)行分割操作。通過基于預(yù)測(cè)密集的局部形狀描述子[17]，將這些描述子組裝成實(shí)例，通過卷積神經(jīng)網(wǎng)絡(luò)進(jìn)行訓(xùn)練預(yù)測(cè)，并對(duì)其進(jìn)行分割操作。

通過訓(xùn)練一個(gè)卷積神經(jīng)網(wǎng)絡(luò)來預(yù)測(cè)較為密集的局部形狀塊，從中在單通道中將所有的實(shí)例集合到一個(gè)圖像中。卷積神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)圖像中每個(gè)像素的密集塊，然后用它來找到塊大小內(nèi)每對(duì)像素的一致性。從中挑選出最符合的塊進(jìn)行連接。塊與塊的連接邊緣是通過連接塊一致性的賦值分?jǐn)?shù)確定的。卷積神經(jīng)網(wǎng)絡(luò)會(huì)在函數(shù)上產(chǎn)生一個(gè)估值p。

p*:Dom(I)×P→{0,1}，

(1)

(2)

式中：

Dom(I)——采集圖像I；

P——圖像中的密集偏移集；

Instance——實(shí)例；

x——圖像中的像素；

dx——像素x的變化量；

fg(I)——前景圖像。

對(duì)于每個(gè)像素x∈Rd在d維圖像I中的前景fg(I)中，當(dāng)x+dx和x屬于同一實(shí)例，則在密集偏移集P?Rd處的每個(gè)像素的x+dx是固定不變的。通過對(duì)每個(gè)密集塊與塊之間邊緣的確定，來完成對(duì)煙絲中成團(tuán)煙絲輪廓的分割。

利用訓(xùn)練深度卷積神經(jīng)網(wǎng)絡(luò)來預(yù)測(cè)函數(shù)p。它通過預(yù)測(cè)輸入圖像的每個(gè)像素的px(x+P)來實(shí)現(xiàn)。因此，集合P的基數(shù)決定了網(wǎng)絡(luò)輸出通道的數(shù)量。使用的是U-Net作為主干架構(gòu)，通過在所有輸出上訓(xùn)練網(wǎng)絡(luò)的標(biāo)準(zhǔn)交叉熵?fù)p失來計(jì)算平均值。為了方便預(yù)測(cè)形狀表示結(jié)果，在U-Net上層路徑中保持特征圖的數(shù)量固定，因此避免了僅從標(biāo)準(zhǔn)U-Net倒數(shù)第二層中的有限元特征圖預(yù)測(cè)到的高維像素輸出。

U-Net的輸入是一張經(jīng)過壓縮后的256×256圖像，網(wǎng)絡(luò)是由卷積和Max Pooling構(gòu)成的一系列降采樣操作的壓縮路徑。壓縮路徑由4個(gè)block組成，每個(gè)block使用了3個(gè)有效卷積和1個(gè)Max Pooling降采樣，每次降采樣之后特征圖的個(gè)數(shù)會(huì)乘2，因此特征圖的尺寸會(huì)發(fā)生變化，最終得到了尺寸為41×41的特征圖。如圖7所示，白色輪廓表示單個(gè)煙絲輪廓，同時(shí)含有斜線和灰點(diǎn)的煙絲輪廓代表了成團(tuán)的煙絲，而不同的輪廓煙絲代表了分割出來的不同煙絲。

圖6 部分成團(tuán)煙絲Figure 6 Part of the cut tobacco

圖7 分割成團(tuán)煙絲Figure 7 Cut into pieces

(4)在得到的煙絲連通域中，對(duì)圖像進(jìn)行細(xì)化，提取煙絲骨架[18]。細(xì)化是將圖像的線條從多像素寬度減少到單位像素寬度的過程，又叫骨架化，細(xì)化的結(jié)果是原連通域的中心線，其穩(wěn)定狀態(tài)是將圖像細(xì)化到最低限度相連且無斷點(diǎn)的線[19]。

得到細(xì)化后的煙絲圖像后，有時(shí)還需要對(duì)煙絲圖像去除毛刺，主要思路就是將所有分支與毛刺都看成不同長(zhǎng)度的分支，如果分支長(zhǎng)度過短小于一定長(zhǎng)度就認(rèn)為是毛刺被去掉。每一個(gè)分支的結(jié)尾點(diǎn)都是一個(gè)端點(diǎn)。首先，遍歷所有連通區(qū)所有點(diǎn)，找到所有端點(diǎn)，端點(diǎn)的定義為其8鄰域點(diǎn)只有一個(gè)像素點(diǎn)，如圖8所示，斜線像素點(diǎn)8鄰域內(nèi)只有黑色點(diǎn)一個(gè)像素點(diǎn)，那么判斷斜線點(diǎn)為端點(diǎn)。遍歷每一個(gè)端點(diǎn)，直到遍歷到該分支的起始點(diǎn)也就是分岔點(diǎn)，分岔點(diǎn)定義為：經(jīng)過細(xì)化后，8鄰域的像素點(diǎn)個(gè)數(shù)≥3的點(diǎn)，統(tǒng)計(jì)每個(gè)分支的長(zhǎng)度，長(zhǎng)度小于閾值T的分支被認(rèn)為是毛刺分支。之后將毛刺分支擦除，得到完好的煙絲細(xì)化圖像。圖9為分割后的部分煙絲圖像，圖10為去毛刺后的細(xì)化圖像。

圖8 判斷煙絲輪廓端點(diǎn)Figure 8 Judging the end point of cut tobacco contour

圖9 分割后的煙絲圖像Figure 9 Cut tobacco image after segmentation

圖10 細(xì)化后的煙絲圖像Figure 10 Cut tobacco image after skeletonization

統(tǒng)計(jì)骨架的長(zhǎng)度即為煙絲的長(zhǎng)度。根據(jù)《卷煙工藝規(guī)范》[20]中規(guī)定，若煙絲長(zhǎng)度>2.5 mm即為整絲，若煙絲長(zhǎng)度<1.0 mm即為碎絲。然后通過計(jì)算煙絲的面積，擬合出整絲和碎絲的重量，從而計(jì)算出煙絲的整絲率和碎絲率。

3 檢測(cè)效果

3.1 在線持續(xù)檢測(cè)時(shí)間的設(shè)定

為確定適合煙絲結(jié)構(gòu)在線持續(xù)檢測(cè)時(shí)間，把持續(xù)檢測(cè)時(shí)間設(shè)定為10個(gè)時(shí)間段(1～10 min)，每個(gè)時(shí)間段分別測(cè)量10次，計(jì)算煙絲結(jié)構(gòu)碎絲率的平均值和變異系數(shù)，結(jié)果見表1。由表1可知，持續(xù)檢測(cè)時(shí)間<5 min時(shí)，各組檢測(cè)結(jié)果變異系數(shù)>5%，持續(xù)檢測(cè)時(shí)間≥5 min時(shí)，變異系數(shù)<5%。因此，將在線持續(xù)檢測(cè)時(shí)間設(shè)定在5 min以上，能夠較為真實(shí)地反映制絲線煙絲結(jié)構(gòu)實(shí)際質(zhì)量狀況。

3.2 煙絲結(jié)構(gòu)裝置的攤薄效果分析

為驗(yàn)證在線煙絲結(jié)構(gòu)檢測(cè)裝置的攤薄效果，針對(duì)湖南某卷煙品牌A的成品煙絲，每間隔5 min統(tǒng)計(jì)一次成團(tuán)率，共統(tǒng)計(jì)5次(一次統(tǒng)計(jì)結(jié)束后立即開始下一次的統(tǒng)計(jì))，成團(tuán)率分別為4.28%，5.07%，4.51%，4.22%，4.86%，均值為4.59%。說明此設(shè)備對(duì)于煙絲的攤薄效果較好。

3.3 在線煙絲結(jié)構(gòu)裝置的長(zhǎng)度準(zhǔn)確性檢測(cè)分析

為了驗(yàn)證煙絲結(jié)構(gòu)檢測(cè)裝置的長(zhǎng)度的準(zhǔn)確性以及重復(fù)性，隨機(jī)選取了10根經(jīng)過游標(biāo)卡尺準(zhǔn)確測(cè)量的煙絲，放入設(shè)備中連續(xù)采集5次。由圖11可知，測(cè)量值和真實(shí)值的吻合度較高，相關(guān)系數(shù)達(dá)到了0.999 8。由表2可知，長(zhǎng)度測(cè)量結(jié)果的相對(duì)誤差<2.56%，變異系數(shù)<3.67%。說明系統(tǒng)準(zhǔn)確性和重復(fù)性均能滿足測(cè)量要求。

表1 不同持續(xù)檢測(cè)時(shí)間的煙絲結(jié)構(gòu)碎絲率Table 1 Broken silk rate of tobacco structure at different continuous detection time

圖11 樣品重復(fù)檢測(cè)結(jié)果Figure 11 Repeated test results of samples

3.4 與國標(biāo)振篩法對(duì)比

為驗(yàn)證在線煙絲結(jié)構(gòu)檢測(cè)裝置檢測(cè)效果，采用振篩法[21](取樣量為1 000 178 g/次，水分11%～14%)和在線煙絲結(jié)構(gòu)檢測(cè)裝置對(duì)湖南某品牌A同一批次的煙絲進(jìn)行離線和在線對(duì)比檢測(cè)，根據(jù)3.1節(jié)確定的在線連續(xù)取樣時(shí)間為5 min以上，兩種方法均連續(xù)取樣10次，結(jié)果如表3 所示，在相同的檢測(cè)任務(wù)下，振篩法耗時(shí)為在線煙絲結(jié)構(gòu)檢測(cè)裝置的4倍。相對(duì)于振篩法，在線煙絲結(jié)構(gòu)檢測(cè)裝置的標(biāo)準(zhǔn)偏差和變異系數(shù)明顯減小，且均小于5%，說明檢測(cè)裝置測(cè)量精度較高，穩(wěn)定性較好。

表2 在線煙絲結(jié)構(gòu)檢測(cè)裝置長(zhǎng)度測(cè)量結(jié)果Table 2 Length measurement results of online cut tobacco structure detection device

表3 不同方法下煙絲整絲率和碎絲率測(cè)量結(jié)果對(duì)比Table 3 Comparison of measurement results of cut tobacco whole silk rate and broken silk rate under different methods

4 結(jié)論

設(shè)計(jì)的煙絲結(jié)構(gòu)在線檢測(cè)裝置，利用風(fēng)選攤薄方式與圖像智能識(shí)別技術(shù)，實(shí)現(xiàn)了煙絲結(jié)構(gòu)在線檢測(cè)，解決了現(xiàn)有檢測(cè)設(shè)備只能進(jìn)行離線檢測(cè)的問題，有效提高了煙絲結(jié)構(gòu)的檢測(cè)精度和效率。該裝置僅對(duì)湖南中煙的部分品牌煙絲進(jìn)行了分析，后續(xù)需進(jìn)一步擴(kuò)展更多的煙絲品牌進(jìn)行分析，以期實(shí)現(xiàn)煙絲結(jié)構(gòu)在線檢測(cè)裝置的工業(yè)應(yīng)用。