基于圖像處理技術(shù)的儲糧害蟲快速識別

摘要 為了快速檢測、識別出儲糧蟲害圖像中的蟲體信息，為蟲害的綜合防治提供依據(jù)，提出了一種基于圖像處理技術(shù)的害蟲檢測與識別快速方法，該方法采用圖像預(yù)處理技術(shù)、形態(tài)學(xué)開重建以及基于改進Hausdorff距離的模板匹配算法。結(jié)果表明：該方法不僅可以正確檢測、識別出害蟲信息，并且在檢測與識別效率上基本達到了實時性檢測的要求。

關(guān)鍵詞 圖像處理；害蟲識別；形態(tài)學(xué)操作；模板匹配

中圖分類號S37；TP391文獻標識碼A文章編號0517-6611（2014）30-10784-04

基金項目國家自然科學(xué)基金項目（61261024）；海南省自然科學(xué)基金項目（614221）；海南省教育廳基金項目（Hjkj201314、HNKY201418）；海南大學(xué)青年基金項目（qnjj1186）。

作者簡介胡祝華（1979-），男，湖南桃江人，講師，在讀博士，從事信號處理與模式識別、無線通信研究。*通訊作者，講師，碩士，從事計算機視覺與模式識別研究。

利用計算機圖像處理技術(shù)實現(xiàn)對儲糧害蟲的分割和正確識別是近年來儲糧蟲害綜合防治中具有活力的交叉研究方向之一[1-3]。把害蟲數(shù)量控制在對儲糧損害的水平以下，既不會因為害蟲造成損失，也不會因盲目防治造成浪費，加重對糧食和環(huán)境的污染，因此及時準確地獲取蟲害信息至關(guān)重要。

進入21世紀以來國內(nèi)外學(xué)者基于圖像處理技術(shù)針對害蟲檢測和識別的研究有很多。2009年，胡濤等利用灰色理論進行儲糧害蟲的邊緣檢測[4]。2010年，張衛(wèi)芳等采用基于混合高斯模型對儲糧害蟲進行分割處理[5]。張建華等基于SVM（支持向量機）模型、遺傳算法對儲糧害蟲進行智能識別[6]。2011年，SolisSanchez等采用尺度不變特征的LOSS V2算法對害蟲進行檢測和識別[7]。2012年，Xia等基于計算機視覺的多重分形分析方法，用來檢測溫室中的微小害蟲（如粉虱等）[8]。陸光等采用基于SVM的多特征融合算法對天牛圖像進行識別[9]。

上述研究方法在特定場合下都取得了較好的分割和識別效果，但由于采用了比較復(fù)雜的智能算法，在檢測和識別效率上都具有一定的局限性，很難達到實時檢測的要求。為了提高儲糧害蟲的識別效率，在保證識別正確率的前提下，筆者提出了一種基于輕量級圖像處理技術(shù)的害蟲檢測與識別快速方法。通過該方法可以快速識別出儲糧中的害蟲，為下一步的綜合防治提供及時有效的科學(xué)依據(jù)。

1識別方法概述

1.1整體思路以儲糧害蟲為研究對象，運用圖像處理技術(shù)實現(xiàn)對儲糧害蟲的目標檢測與識別。該方法的整體結(jié)構(gòu)如圖1所示。圖1整體結(jié)構(gòu)獲取原始圖像后，首先需對原始圖像進行預(yù)處理和均衡操作，包括圖像格式轉(zhuǎn)換、彩色圖像灰度化、圖像灰度拉伸、二值化等操作；其次對預(yù)處理后的圖像進行形態(tài)學(xué)濾波和重建工作；接著對處理后的圖像進行目標害蟲的提取；最后將提取出的目標害蟲與庫中的信息進行配準，配準過程包括對目標蟲體的旋轉(zhuǎn)、縮放、模板匹配等操作，識別結(jié)束后給出害蟲的識別結(jié)果。

1.2試驗測試環(huán)境①硬件環(huán)境：高清數(shù)碼攝像機、PC[Inter（R），Core（TM）i53230M CPU @2.60GHZ，2.56GB內(nèi)存]。②軟件環(huán)境：Matlab 8.1、Visual C++ 6.0、WIN7操作系統(tǒng)。③試驗語言：C語言。

2圖像預(yù)處理

在拍攝以及傳輸圖像的過程中受噪聲、外部環(huán)境等因素的干擾，造成圖像質(zhì)量有所下降。為此，需要增強圖像中的有用信息，改善圖像效果，使圖像中感興趣區(qū)域的特征差別更加明顯。該研究主要通過對圖像的格式轉(zhuǎn)換、灰度化、灰度拉伸以及二值化操作對圖像進行預(yù)處理分析。分析結(jié)果為后續(xù)的處理流程提供好的輸入。

2.1圖像格式轉(zhuǎn)換與灰度化在儲糧蟲害圖像獲取的過程中，考慮到監(jiān)控設(shè)施存儲空間大小的限制，以及圖像傳輸?shù)呢摵?，最初的畫面存儲為JPEG格式。JPEG圖像占用空間小，但這種格式的圖像存儲方式較為復(fù)雜，涉及到Huffman編碼，對于圖像后期各方面的處理并不適合。為此，將JPEG圖像轉(zhuǎn)換為存儲方式簡單的BMP圖像，這樣對圖像的灰度化、二值化等處理比較方便，格式轉(zhuǎn)換后的效果如圖2（b）所示。

彩色圖像信息量很大，如果對其直接進行圖像處理將會很困難，處理的速度也會很慢，極大地降低檢測與識別的效率。因此首先需要將彩色圖像轉(zhuǎn)化成灰度圖像。該研究采用標準加權(quán)平均值法進行圖像的灰度化，如式（1）所示。

式中，wR、wG、wB分別為彩色圖像R、G、B分量的權(quán)值，經(jīng)過反復(fù)試驗測試獲取經(jīng)驗值為：wR=0.299，wG=0.587，wB=0114。

害蟲圖像灰度化后的效果如圖2（c）所示。

2.2圖像的均衡處理均衡處理主要包括對儲糧害蟲圖像進行灰度拉伸和圖像的二值化處理。

2.2.1灰度拉伸?；叶然蟮幕叶戎荡蟛糠旨性谝粋€小范圍的區(qū)域內(nèi)，給害蟲目標的分割閾值選擇帶來了一定的困難，因此需要將灰度拉伸到整個范圍，才能顯示出清晰的圖像，使得在該范圍內(nèi)的像素，亮的越亮，暗的越暗，達到增強對比度的目的。灰度拉伸示意圖如圖2（a）所示。圖中，L-1表示最大灰度級，r1和r2表示進行灰度拉伸的最小輸入灰度和最大輸入灰度，s1和s2是變換拉伸后的灰度輸出，a、b、c為3段直線的斜率，因為是對圖像中的重要區(qū)域進行對比度拉伸，所以斜率b必須大于1，而a、c的值往往相等。在該試驗中a、b、c的取值分別為：a=c=0.5，b=1.75?；叶壤旌蟮男Ч鐖D2（d）所示，從效果圖可知，灰度拉伸后害蟲圖像的對比度得到了明顯的增強，有利于二值化處理中閾值的選擇。

安徽農(nóng)業(yè)科學(xué)2014年圖2蟲害圖像預(yù)處理后的灰度拉伸效果2.2.2二值化處理?；叶葓D像數(shù)據(jù)量仍然較大，會影響到檢測與識別的效率。二值化是將具有256個灰度等級的圖像通過閾值處理轉(zhuǎn)換成只有2個灰度等級（0、255）的圖像，二值化處理后使圖像具有更小的數(shù)據(jù)和更強的目標對比，有利于對目標圖像的進一步處理分析。 二值化算法如式（2）所示。

3目標蟲體的形態(tài)學(xué)運算與提取

3.1形態(tài)學(xué)運算為了在目標蟲體提取時足夠精確，對二值化后的圖像進行了無損濾波操作，消除噪聲的同時填補蟲體目標上的空洞，采用的算法如下。

3.1.1 首先對二值圖像進行形態(tài)學(xué)開重建濾波運算。如式（3）所示。

3.1.2在第1步的基礎(chǔ)上進行形態(tài)學(xué)閉運算。開重建濾波后，再利用形態(tài)學(xué)閉運算可以彌合目標蟲體圖像中的較窄的間斷，消除小的空洞，填補圖像中細小斷裂。如式（5）所示。

3.2目標蟲體的提取為了提取出感興趣區(qū)域中的目標蟲體，并統(tǒng)計出目標蟲體的面積和周長，處理了如下兩個關(guān)鍵問題。

3.2.1判斷當前分析的目標點是否處于正在查找的某一目標區(qū)域。采用的算法如下：①找到一個合理點（如第1次遍歷時可以尋找第1個灰度為0的點，這個點肯定是某一區(qū)域的邊界點），將該點記為點A。②在A點的右、右下、下、左下4個鄰近點中至少有一個是邊界點（一個邊界點的上下左右4個相鄰點都是黑點），將其記為B。③從B開始查找，按照右、右下、下、左下、左、左上、上、右上的順序?qū)ふ蚁噜徑c中的邊界點，每尋找到一個邊界點將其記為點C，并記錄該點的位置信息。④判斷點C，如果C就是A點，則表明該區(qū)域已查找完成，否則重復(fù)執(zhí)行前面的操作，直到重新找到A為止。

3.2.2圖像中目標蟲體的提取。為了不漏掉任何一處區(qū)域，需要對圖像中的像素點逐個遍歷，將所有區(qū)域都提取出來。該研究采取的算法是：在判斷出某一蟲體目標區(qū)域后，將該區(qū)域的點寫入到新的圖像中，同時將原圖中這些像素點刪除（即將其灰度值置為255）。這樣在提取該區(qū)域的蟲體生成新圖像的同時，原圖中該部分區(qū)域也被刪除，使得查找余下蟲體區(qū)域更為方便。

4目標蟲體的識別

4.1識別前的預(yù)處理蟲體提取出來后，由于生成的新圖是按照原圖的大小進行操作的，所以為了方便后面的操作，減少操作的數(shù)據(jù)量，需要將新生成的蟲體圖像的空白無效區(qū)域裁剪掉。其算法操作步驟是：找出目標的最高點、最左點、最右點以及最下點，根據(jù)這些點的坐標求出蟲體的寬和高，然后生成新圖。

采用對蟲體姿態(tài)和形體特征的配準進行害蟲的識別。在配準的過程中，必須要建立一個配準的基準，該研究以目標蟲體的最長軸在水平方向上為基準，此時需要對目標蟲體進行旋轉(zhuǎn)，找到最長軸，旋轉(zhuǎn)到水平方向。另外，在識別中目標蟲體的分辨率大小需要與模板庫中樣本的大小保持一致，因此需要對目標蟲體參照樣本大小進行縮放。

4.2改進Hausdorff距離的模板匹配為了抵抗噪聲的干擾，采用了基于改進Hausdorff距離的模板匹配算法[12]。改進Hausdorff距離是用來描述兩組點集之間相似程度的一種度量方法，公式為：

5試驗結(jié)果

以拍到的米象害蟲圖像為例，圖3（a）為蟲害圖像進行預(yù)處理后的結(jié)果，圖3（e）為形態(tài)學(xué)操作后的結(jié)果。對目標蟲體進行提取，并去掉空白無效區(qū)域，得到如圖4（a）所示的效果。旋轉(zhuǎn)到基準位置后，因為無法預(yù)知是否和模板庫中樣本的水平朝向一致，因此對旋轉(zhuǎn)完成的圖像還需進行180°的旋轉(zhuǎn)，形成左、右兩個朝向的目標蟲體。另外，為了防止在配準中出現(xiàn)較大誤差，對兩個方向的蟲體圖像都進行角度的微調(diào)，這樣在兩個方向都有了與水平方向略有角度差異（正負2度）的3幅圖像。試驗效果如圖4（b）所示。

6結(jié)語

儲糧害蟲的檢測與識別是蟲害智能監(jiān)測與防治中的重要一環(huán)，該研究基于計算機圖像處理技術(shù)構(gòu)建了一套目標提取的快速方法，并將其應(yīng)用到蟲害圖像中蟲體的分割和提取上，然后再應(yīng)用快速而簡單的識別算法對蟲體進行識別。試驗結(jié)果表明，該研究提出的方法針對不是特別復(fù)雜場景中的害蟲識別具有較好的效果，并且識別的效率基本達到了實時性的要求。

參考文獻

[1] 胡麗華，郭敏，張景虎，等.儲糧害蟲檢測新技術(shù)及應(yīng)用現(xiàn)狀[J].農(nóng)業(yè)工程學(xué)報，2007（11）：286-290.

[2] 毛罕平，張紅濤.儲糧害蟲圖像識別的研究進展及展望[J].農(nóng)業(yè)機械學(xué)報，2008（4）：175-179，186.

[3] 王江寧，紀力強.昆蟲圖像分割方法及其應(yīng)用[J].昆蟲學(xué)報，2011（2）：211-217.

[4] 胡濤，徐磊.灰色理論在儲糧害蟲圖像處理中的應(yīng)用[J].湖北農(nóng)業(yè)科學(xué)，2009（3）：729-731.

[5] 張衛(wèi)芳，郭敏.基于圖割理論的儲糧害蟲圖像分割[J].科學(xué)技術(shù)與工程，2010（7）：1661-1664，1679.

[6] 張建華，朱春華.基于遺傳算法和支持向量機的儲糧害蟲圖像識別[J].安徽農(nóng)業(yè)科學(xué)，2010（17）：8833-8834.

[7] SOLISSANCHEZ L O，CASTANEDAMIRANDA R，GARCIAESCALANTE J J，et al.Scale invariant feature approach for insect monitoring[J].Computers and Electronics in Agriculture，2011，75（1）：92-99.

[8] XIA C，LEE J M，LI Y，et al.In situ detection of smallsize insect pests sampled on traps using multifractal analysis[J].Optical Engineering，2012，51（2）：27001.

[9] 陸光，滿慶麗，徐然.基于SVM的多特征融合的天牛圖像識別[J].森林工程，2012（4）：21-25.

[10] OTSU N.A threshold selection method from graylevel histogram [J].IEEE Transaction on Systems，Man and Cybernetics，1979，9（1）：62-66.

[11] GONZALEZ R C，WOODS R E.數(shù)字圖像處理[M].阮秋琦，阮宇智，譯.3版.北京：電子工業(yè)出版社，2011：402-411.