鄭飛，李媛，劉德強(qiáng)，張大波

1.湖北中煙工業(yè)有限責(zé)任公司武漢卷煙廠，武漢市漢陽(yáng)區(qū)十升路22號(hào) 430051

2.中國(guó)煙草總公司鄭州煙草研究院煙草行業(yè)煙草工藝重點(diǎn)實(shí)驗(yàn)室，鄭州高新技術(shù)產(chǎn)業(yè)開發(fā)區(qū)楓楊街2號(hào) 450001

基于圖像處理的煙片加料均勻性評(píng)價(jià)方法

鄭飛1，李媛1，劉德強(qiáng)1，張大波*2

1.湖北中煙工業(yè)有限責(zé)任公司武漢卷煙廠，武漢市漢陽(yáng)區(qū)十升路22號(hào) 430051

2.中國(guó)煙草總公司鄭州煙草研究院煙草行業(yè)煙草工藝重點(diǎn)實(shí)驗(yàn)室，鄭州高新技術(shù)產(chǎn)業(yè)開發(fā)區(qū)楓楊街2號(hào) 450001

為量化評(píng)價(jià)加料工序單張煙片料液分布和加料過(guò)程均勻性，應(yīng)用數(shù)字圖像處理技術(shù)對(duì)料液附著點(diǎn)分布圖像進(jìn)行處理，建立了基于圖像處理的煙片加料均勻性評(píng)價(jià)方法，采用料液分布變異系數(shù)表征單張煙片料液分布均勻性，采用煙片間料液附著量變異系數(shù)表征加料過(guò)程均勻性。驗(yàn)證結(jié)果表明：基于圖像處理的加料均勻性評(píng)價(jià)方法檢測(cè)結(jié)果，與采用YC/T 353—2010行業(yè)標(biāo)準(zhǔn)檢測(cè)結(jié)果具有相同趨勢(shì)，檢測(cè)結(jié)果可有效表征整體加料過(guò)程的改善或劣化程度。

煙片加料；圖像處理；均勻性；評(píng)價(jià)方法

煙片加料作為制絲加工的重要工藝環(huán)節(jié)［1-3］，主要功能是按生產(chǎn)和配方的要求，將料液均勻、準(zhǔn)確地施加到煙片上。加料均勻與否對(duì)成品卷煙內(nèi)在品質(zhì)起著重要作用。為提高加料均勻性，煙草行業(yè)進(jìn)行了一些研究和改進(jìn)［4-5］，如改善噴嘴霧化狀態(tài)、調(diào)整引噴蒸汽壓力等。近年來(lái)開始采用的利用標(biāo)記物來(lái)表征加料均勻性的方法［6-8］，能夠?qū)恿线^(guò)程的均勻性進(jìn)行評(píng)價(jià)，但由于樣品量較大，制樣、檢測(cè)過(guò)程復(fù)雜，檢測(cè)結(jié)果滯后。數(shù)字圖像處理技術(shù)在煙草行業(yè)廣泛應(yīng)用于智能剔除煙草異物、片煙規(guī)格檢測(cè)等方面［9-11］，作者利用數(shù)字圖像處理技術(shù)［12］，通過(guò)對(duì)料液附著點(diǎn)分布圖像數(shù)據(jù)的提取和處理，建立了基于圖像處理的加料均勻性評(píng)價(jià)方法，旨在快捷高效地評(píng)價(jià)加料過(guò)程均勻性。

1 方法的建立

1.1 材料、試劑和儀器

40 mm×40 mm標(biāo)準(zhǔn)紙片作為煙片標(biāo)準(zhǔn)替代物［亞龍紙制品（昆山）有限公司，白色A4，定量80 g/cm2，厚度0.1 mm］。標(biāo)準(zhǔn)紙片單位面積質(zhì)量、尺寸與黃鶴樓軟藍(lán)生產(chǎn)線中加料過(guò)程煙片平均面積及質(zhì)量接近。

檸檬黃10%稀釋液（檸檬黃∶標(biāo)準(zhǔn)糖料=1∶9，質(zhì)量比）作為糖料標(biāo)準(zhǔn)替代物。其目的是提高辨識(shí)度，僅作為試驗(yàn)使用，不進(jìn)入正常產(chǎn)品。

奧林巴斯TG-3型相機(jī)。湖北中煙武漢卷煙廠制絲黃鶴樓軟藍(lán)生產(chǎn)線。

1.2 樣品制備及圖像處理

加料設(shè)備穩(wěn)定運(yùn)行后，在進(jìn)料振槽等時(shí)間間隔撒入120片標(biāo)準(zhǔn)紙片，在加料出口人工挑選出標(biāo)準(zhǔn)紙片，隨機(jī)挑選100片作為圖譜成像分析對(duì)象，并記錄加料精度。加料工序主要參數(shù)設(shè)置如下：物料流量2 500 kg/h、料液引噴蒸汽壓力0.18～0.22 MPa、滾筒轉(zhuǎn)速8 r/min、加料比例3.5%。

用相機(jī)定點(diǎn)拍攝樣品照片并轉(zhuǎn)換為JPG圖像格式，分別標(biāo)記為Aj，j∈（1，100）（如圖1a為A10樣品采集圖像）。采用MATLAB7.0對(duì)圖像進(jìn)行預(yù)處理，方法如下：①提取原始采集圖像有效區(qū)域：利用K均值聚類算法進(jìn)行背景檢測(cè)，保留樣品區(qū)域，將無(wú)效樣品區(qū)域標(biāo)記為紅色（A10樣品經(jīng)背景檢測(cè)處理后有效區(qū)域圖像如圖1b所示）。②統(tǒng)一樣品區(qū)域大小：將所有采集圖像以樣品區(qū)域?yàn)橹行闹夭蓸訛?00像素×400像素大小的方形區(qū)域（A10重采樣后的樣品圖像如圖1c所示）。③彩色圖像轉(zhuǎn)換為灰度圖像：將樣品圖像的RGB色彩模型轉(zhuǎn)化為HSV模型［13］，提取V分量作為樣品圖像各像素點(diǎn)的灰度值，從而得到樣品灰度圖像（A10樣品灰度圖像如圖1d所示），其像素點(diǎn)灰度值大小可表征該點(diǎn)料液施加量；單張紙片各像素點(diǎn)灰度值之和，可表征該樣品的料液施加量。④灰度圖像轉(zhuǎn)換為二值圖像：利用最大類間差法自動(dòng)選擇閾值，將樣品灰度圖像轉(zhuǎn)化為樣品二值圖像［14-16］（A10樣品二值圖像如圖1e所示），白色（像素值為1）為染色區(qū)域，黑色（像素值為0）為未染色區(qū)域。⑤二值圖像轉(zhuǎn)換為網(wǎng)格圖像：將樣品的二值圖像按區(qū)域大小均勻劃分為80像素×80像素的網(wǎng)格圖像（A10樣品網(wǎng)格圖像如圖1f所示）。

圖1 A10樣品圖像

1.3 加料均勻性評(píng)價(jià)模型的建立

單張紙片料液覆蓋率：指加料后，單張紙片上附著有料液的區(qū)域占整張紙片區(qū)域的比例，如式（1）所示：

式中：Kj—單張紙片料液覆蓋率，%；Nj—二值圖像統(tǒng)計(jì)出染色像素點(diǎn)的個(gè)數(shù)；j—第j個(gè)樣品；d—圖像重采樣的大小。

單張紙片料液分布均勻性：指加料后，單張紙片附著的料液在紙片每個(gè)網(wǎng)格區(qū)域分布的均勻性，如式（2）所示：

式中：Hj—單張紙片料液分布均勻性；i—網(wǎng)格圖像的第i個(gè)網(wǎng)格單元；Mj（i）—第j個(gè)樣品第i個(gè)網(wǎng)格單元附著有料液的像素點(diǎn)個(gè)數(shù)；RSD［Mj（i）］—第j個(gè)樣品Mj（i），i∈（1，5×5）的變異系數(shù)。

加料過(guò)程均勻性：指經(jīng)過(guò)加料過(guò)程后，所有紙片間附著料液量的分布均勻性。單張紙片料液附著量用各像素點(diǎn)灰度值的和表征，如式（3）所示；加料過(guò)程均勻性用加料均勻性系數(shù)來(lái)表征，如式（4）所示：

式中：Qj—單張紙片上料液附著量；Cj（x，y）—第j張樣品在（x，y）像素點(diǎn)的灰度值；PU—加料過(guò)程的均勻性；RSD（Qj）—100個(gè)樣品Qj，j∈（1，10）的變異系數(shù)。

2 方法的應(yīng)用

2.1 單張紙片上料液的覆蓋率和分布均勻性

分別檢測(cè)計(jì)算100片樣品的料液覆蓋率和分布均勻性。以A10樣品為例，根據(jù)單片料液覆蓋率和分布均勻性評(píng)價(jià)方法，利用公式（1）計(jì)算出A10樣品的單張紙片料液覆蓋率K10=29.8%，利用公式（2）計(jì)算出A10樣品單張紙片料液分布均勻性H10= 77.5%。

計(jì)算100片樣品的料液覆蓋率和分布均勻性，應(yīng)用SPSS17.0統(tǒng)計(jì)分析軟件［17］，分別對(duì)100片樣品的計(jì)算結(jié)果進(jìn)行Shapiro-Wilk的W統(tǒng)計(jì)量及矩陣法的正態(tài)分布檢驗(yàn)，檢驗(yàn)結(jié)果見表1。由表1可知，單張樣品料液覆蓋率和分布均勻性的W檢驗(yàn)確定的概率均符合P＞0.05，同時(shí)偏度和峰度的U＜1.96，說(shuō)明樣品的料液覆蓋率和分布均勻性的檢測(cè)結(jié)果服從正態(tài)分布［18］，加料試驗(yàn)過(guò)程穩(wěn)定。

表1 100片樣品正態(tài)分布檢驗(yàn)

2.2 加料過(guò)程均勻性評(píng)價(jià)

應(yīng)用公式（3），計(jì)算每片樣品上像素灰度值之和Qj；應(yīng)用公式（4），根據(jù)甲酸100片樣品的料液附著量的變異系數(shù)，計(jì)算得出加料過(guò)程均勻性PU值為68.3%。

3 方法的驗(yàn)證

行業(yè)標(biāo)準(zhǔn)方法（YC/T 353—2010［7］）是利用加料后煙葉中丙二醇含量的變異系數(shù)CU值表征加料均勻性。采用基于圖像處理的加料均勻性評(píng)價(jià)方法和文獻(xiàn)［7］中的方法，對(duì)加料過(guò)程進(jìn)行對(duì)比驗(yàn)證試驗(yàn)，驗(yàn)證試驗(yàn)步驟如下：①調(diào)配料液，檸檬黃∶1,2-丙二醇∶標(biāo)準(zhǔn)糖料=1∶1∶8（質(zhì)量比）。②加料試驗(yàn)，試驗(yàn)條件1設(shè)置同1.2小節(jié)；試驗(yàn)條件2更換加料噴嘴以增加料液的霧化面積；試驗(yàn)條件3更換加料噴嘴以減小料液的霧化面積。在3種試驗(yàn)條件下各進(jìn)行了5組試驗(yàn)。15組試驗(yàn)的結(jié)果見圖2。由圖2可以看出，應(yīng)用兩種方法檢測(cè)的加料均勻性評(píng)價(jià)結(jié)果具有相同的趨勢(shì)，說(shuō)明基于圖像處理的加料過(guò)程均勻性評(píng)價(jià)方法能有效地表征加料過(guò)程的改善或劣化程度，可用于對(duì)加料過(guò)程均勻性的快速評(píng)價(jià)。

圖2 YC/T 353方法與基于圖像處理的煙片加料均勻性評(píng)價(jià)方法對(duì)比結(jié)果

4 結(jié)論

（1）建立了基于圖像處理的加料均勻性評(píng)價(jià)方法，采用本方法能對(duì)單張煙片料液覆蓋率、單張煙片料液分布均勻性和加料過(guò)程均勻性進(jìn)行定量評(píng)價(jià)，檢測(cè)過(guò)程快捷高效。

（2）基于圖像處理的加料均勻性評(píng)價(jià)方法與YC/T 353—2010方法的檢測(cè)結(jié)果具有相同趨勢(shì)，檢測(cè)結(jié)果可有效表征整體加料過(guò)程的改善或劣化程度。

［1］國(guó)家煙草專賣局.卷煙工藝規(guī)范［M］.北京：中央文獻(xiàn)出版社，2003.

［2］陳良元.卷煙生產(chǎn)工藝技術(shù)［M］.鄭州：河南科學(xué)技術(shù)出版社，2002.

［3］黃嘉礽.卷煙工藝［M］.北京：北京出版社，2002.

［4］李向東.制絲工藝加料均勻性分析［J］.煙草科技，2003（4）：38-40.

［5］王軍，王聚奎，王雷.加料均勻性評(píng)價(jià)方法的應(yīng)用研究［J］.生物技術(shù)世界，2012（6）：14-16.

［6］廖惠云，張映，赫喜良.卷煙制絲加料工序均勻性的評(píng)價(jià)與控制［J］.煙草科技，2008（8）：5-7，21.

［7］YC/T 353—2010煙草制品和煙草工業(yè)設(shè)備卷煙加料均勻性的測(cè)定［S］.

［8］劉澤春，黃華發(fā)，洪偉齡，等.過(guò)程能力分析在卷煙加料均勻性評(píng)價(jià)中的應(yīng)用［J］.煙草科技，2010（10）：5-7，21.

［9］簡(jiǎn)慶龍，許興鹍.圖像處理在煙葉雜物剔除系統(tǒng)上的應(yīng)用研究［J］.微計(jì)算機(jī)信息，2005（5）：149-150.

［10］周文，韓力群.計(jì)算機(jī)圖像處理技術(shù)在烤煙煙葉形狀特征提取中的應(yīng)用［J］.煙草科技，2000（1）：12-14.

［11］張建平.計(jì)算機(jī)視覺(jué)在煙草行業(yè)中的研究及應(yīng)用展望［J］.煙草科技，1998（2）：22-23.

［12］Gonzalez R C，Woods R E.數(shù)字圖像處理［M］.2版.北京：電子工業(yè)出版社，2003.

［13］章衛(wèi)祥，周秉鋒.一種基于梯度域的彩色圖像轉(zhuǎn)灰度圖像的方法［J］.影像技術(shù)，2007（3）：20-22.

［14］田野，章錦文，趙廣州.基于二維Otsu閾值修正的圖像分割新方法［J］.計(jì)算機(jī)與數(shù)字工程，2012，40（5）：104-107.

［15］齊麗娜，張博，王戰(zhàn)凱.最大類間方差法在圖像處理中的應(yīng)用［J］.無(wú)線電工程，2006，37（7）：25-26.

［16］郭春秋，歐陽(yáng)松.一種快速、準(zhǔn)確的圖像閾值分割新方法［J］.湖南工業(yè)職業(yè)技術(shù)學(xué)院學(xué)報(bào)：自然科學(xué)版，2004（1）：9-11.

［17］劉慶武，胡志艷.如何用SPSS、SAS統(tǒng)計(jì)軟件進(jìn)行正態(tài)性檢驗(yàn)［J］.湖南學(xué)院學(xué)報(bào)：自然科學(xué)版，2005（3）：56-58.

［18］吳梅村.數(shù)理統(tǒng)計(jì)學(xué)基本原理和方法［M］.成都：西南財(cái)經(jīng)大學(xué)出版社，2006.

［19］楊玉波，程志遠(yuǎn).制絲線加香加料系統(tǒng)的改進(jìn)［C］//中國(guó)煙草學(xué)會(huì)工業(yè)專業(yè)委員會(huì)煙草工藝學(xué)術(shù)研討會(huì)論文集.延吉：中國(guó)煙草學(xué)會(huì)，2008：300-302.

責(zé)任編輯 周雅寧

Method for Evaluating Casing Uniformity of Tobacco Strips Based on Image Processing

ZHENG Fei1,LI Yuan1,LIU Deqiang1,and ZHANG Dabo*2
1.Wuhan Cigarette Factory,China Tobacco Hubei Industrial Co.,Ltd.,Wuhan 430051,China
2.Key Laboratory of Tobacco Processing Technology of CNTC,Zhengzhou Tobacco Research Institute of CNTC, Zhengzhou 450001,China

In order to quantitatively evaluate the distribution uniformity of casing on a single tobacco strip and the uniformity of casing process,the image of casing attached points on tobacco strip was processed by digital image processing technique,and a method for evaluating casing uniformity on tobacco strip based on image processing was established,wherein the variation coefficient of casing distribution was used to characterize the distribution uniformity of casing on a single strip and the variation coefficient of attached casing amount between strips was used to characterize the uniformity of casing process.The results of verification showed that the results evaluated by the developed method presented the same trend as did the results tested by industrial standard YC/T 353—2010,Casing Uniformity Measurement,and effectively characterized the variation of the casing process.

Tobacco strip casing；Image processing；Uniformity；Evaluation method

TS46

A

1002-0861（2015）11-0065-04

10.16135/j.issn1002-0861.20151112

2015-02-09

2015-07-10

鄭飛（1977—），本科，工程師，主要從事卷煙制絲工藝研究。E-mail：zhengfei@wh.hbtobacco.cn；*

張大波，E-mail：zhangdabo1979@126.com

鄭飛，李媛，劉德強(qiáng)，等.基于圖像處理的煙片加料均勻性評(píng)價(jià)方法［J］.煙草科技，2015，48（11）：65-68.

ZHENG Fei，LI Yuan，LIU Deqiang，et al.Method for evaluating casing uniformity of tobacco strips based on image processing［J］.Tobacco Science&Technology，2015，48（11）：65-68.