董本志，聶麗酈，景維鵬，崔 航

東北林業(yè)大學(xué) 信息與計算機(jī)工程學(xué)院，哈爾濱 150040

1 引言

榆樹是東北地區(qū)常見樹種，也是重要的經(jīng)濟(jì)樹種和觀賞樹種，但在其生長過程中經(jīng)常受到害蟲侵?jǐn)_。準(zhǔn)確識別害蟲判別災(zāi)害情況，對有效治理蟲災(zāi)有重要意義，利用圖像處理方法可有效對其進(jìn)行識別。

目前利用圖像處理對昆蟲進(jìn)行識別的方法為：先使用圖像去噪[1]、圖像分割等[2]方法對圖像進(jìn)行預(yù)處理，然后采用灰度直方圖[3-4]、隨機(jī)森林算法[5]、方向梯度直方圖（Histogram of Oriented Gradient，HOG）[6]、詞袋算法（Bag of Words，BOW）[7-8]等算法提取圖像的特征，最后將提取到的特征送入支持向量機(jī)（Support Vector Machine，SVM）[9]、前饋神經(jīng)網(wǎng)絡(luò)（Back Propagation，BP）[10]、自組織映射網(wǎng)絡(luò)（Self Organizing Map，SOM）[11]等分類器中，對提取到的特征信息進(jìn)行分類表達(dá)。但上述識別方法均需人為參與到特征提取模板和特征提取算法的設(shè)計中，往往加入了先驗知識，具有很強(qiáng)的主觀性，影響分類器的判斷。同時，針對昆蟲這種目標(biāo)小、紋理特征多樣、結(jié)構(gòu)豐富、姿態(tài)多樣、種間相似度高的目標(biāo)[12]，很難人工設(shè)計出有效描述昆蟲典型特征的模板與算法。

深度學(xué)習(xí)模型由數(shù)據(jù)直接驅(qū)動特征的提取，可以利用像素之間的位置特征，自主學(xué)習(xí)出數(shù)據(jù)間許多容易被人忽視的潛在特征，有效解決了特征提取模板針對性不強(qiáng)的問題。深度學(xué)習(xí)模型中卷積神經(jīng)網(wǎng)絡(luò)在圖像處理方面有較廣泛的應(yīng)用，如漢字識別驗證碼[13]、道路中車輛識別[14]等。目前最成熟且應(yīng)用最廣泛的深度卷積神經(jīng)網(wǎng)絡(luò)模型Faster R-CNN[15]將目標(biāo)檢測的候選區(qū)域生成、特征提取、分類、位置精修四個步驟統(tǒng)一到一個深度網(wǎng)絡(luò)架構(gòu)之內(nèi)，形成了統(tǒng)一的體系結(jié)構(gòu)。識別復(fù)雜背景中的目標(biāo)時具有很好的魯棒性和很強(qiáng)的辨別待識別物體間的細(xì)微差別的能力。同時，利用數(shù)據(jù)自主驅(qū)動模型提取特征，有效克服了傳統(tǒng)識別分類方法中依賴于人工設(shè)計特征提取模板的局限性，可較好地識別出自然環(huán)境中復(fù)雜背景下的目標(biāo)，因此本文采用Faster R-CNN網(wǎng)絡(luò)對榆紫葉甲蟲進(jìn)行框定。但是標(biāo)準(zhǔn)Faster R-CNN網(wǎng)絡(luò)模型是針對標(biāo)準(zhǔn)數(shù)據(jù)集VOC2007的20分類任務(wù)設(shè)計的，生成的初始候選框為固定的三種尺寸、三種比例，用其識別榆紫葉甲蟲時，初始候選框的長寬比不符合榆紫葉甲蟲的長寬比形態(tài)學(xué)特征，容易造成候選框冗余過大。加之榆紫葉甲蟲的甲殼反光，在框定榆紫葉甲蟲目標(biāo)時會出現(xiàn)誤差，在榆紫葉甲蟲和榆樹葉片豁口或孔洞相鄰時或者兩只榆紫葉甲相鄰時框定誤差尤其嚴(yán)重。因此，本文使用Faster R-CNN模型識別榆紫葉甲蟲時，對初始候選框生成網(wǎng)絡(luò)進(jìn)行改進(jìn)，使其生成的初始候選框更加貼合榆紫葉甲蟲本身的特征，減少周圍復(fù)雜環(huán)境造成的影響，以提高識別精度。

2 基于K-means聚類算法的Faster R-CNN網(wǎng)絡(luò)模型

為得到符合榆紫葉甲蟲長寬比形態(tài)學(xué)特征的初始候選框，可以用聚類算法對樣本的長寬值進(jìn)行統(tǒng)計，得出最理想的候選框長寬比。目前常用的幾種聚類算法有K-means聚類、層次聚類、SOM聚類、FCM聚類。其中K-means聚類算法是基于距離的聚類算法，使聚類后每個子類內(nèi)的點到當(dāng)前類中心點的距離之和最小，符合本文的聚類需要。本文以最小化子類內(nèi)的點到當(dāng)前類中心點距離之和為目標(biāo)函數(shù)，通過計算評價聚類結(jié)果中各子類的類內(nèi)緊密程度與類間分離程度關(guān)系的BWP（Between-Within Proportion）指標(biāo)來確定最佳聚類中心點個數(shù)kopt，對n個榆紫葉甲蟲訓(xùn)練樣本矩形標(biāo)簽的長寬比值X{x1,x2,…,xn}進(jìn)行聚類。這里kopt的計算方法見公式（1）。

令聚類中心點個數(shù)k在[2,n)內(nèi)循環(huán)，計算各k值對應(yīng)的所有樣本的BWP值的平均值avgBWP(k)，kopt的值取令avgBWP(k)最大值時對應(yīng)的k值，計算方法見公式（2）。

其中BWP(j,i)是最小類間距離和類內(nèi)距離的衡量指標(biāo)，其計算方法見公式（3）。

公式（3）中b(j,i)為最小類間距離，是第 j類中第i個樣本到其他每個類中樣本的平均距離最小值，其計算方法見公式（4）；w(j,i)為類內(nèi)距離，是第 j類中第i個樣本到其第 j類中其他所有樣本的平均距離，其計算見公式（5）。

其中 j∈[1,ki]，i∈[1,nj]。nj為第 j類的樣本數(shù)；m和j表示類標(biāo)；表示第 j類的第i個樣本；表示第m類的第 p個樣本；nm表示第m類的樣本數(shù)；表示平方歐氏距離。

其中nj表示第 j類的樣本數(shù)，表示第 j類的第q個樣本。

通過計算得到kopt的值后，可同時得到kopt個聚類中心點，可以用集合A={A1,A2,…,Akopt}表示。A中的元素可以用來代替標(biāo)準(zhǔn)網(wǎng)絡(luò)生成的初始候選框的長寬比值0.5、1、2，生成更加符合榆紫葉甲蟲長寬比形態(tài)學(xué)特征的初始候選框，減少框定誤差。

圖1 改進(jìn)的Faster R-CNN網(wǎng)絡(luò)示意圖

調(diào)整了初始候選框的生成比例后，進(jìn)一步對初始候選框的生成尺寸進(jìn)行調(diào)整。標(biāo)準(zhǔn)網(wǎng)絡(luò)生成的初始候選框的三種尺寸固定為1282，2562，5122，但是由于榆紫葉甲蟲的尺寸相對于VOC2007標(biāo)準(zhǔn)數(shù)據(jù)集中的飛機(jī)、汽車、馬等目標(biāo)類別的尺寸來說占據(jù)整張圖像比例很小，屬于小目標(biāo)。而5122這種尺寸對于榆紫葉甲蟲來說冗余過大，導(dǎo)致位置精修時的初始候選框邊框平移量過大，容易造成框定不準(zhǔn)。因此，本文舍棄這種尺寸的候選框，以提高檢測的準(zhǔn)確率。

改進(jìn)的Faster R-CNN網(wǎng)絡(luò)示意圖如圖1所示。

改進(jìn)主要涉及到初始候選框生成文件generate_anchor.py和候選框?qū)游募roposal Layer.py。首先，在proposal Layer.py文件中的類proposal Layer（Caffe.Layer）里根據(jù)公式（1）～（5）定義一個用來對訓(xùn)練數(shù)據(jù)集標(biāo)簽長寬比值進(jìn)行聚類的函數(shù)K_means，返回值為kopt和{A1,A2,…,Akopt}。然后，將兩個返回值傳給generate_anchor.py文件中的generate_anchors函數(shù)，用{A1,A2,…,Akopt}代替標(biāo)準(zhǔn)網(wǎng)絡(luò)中候選框生成比例0.5、1、2。其中g(shù)enerate_anchors函數(shù)根據(jù)傳入的參數(shù)生成生成conv5特征圖上的特征點數(shù)×kopt×2個初始候選框。在對網(wǎng)絡(luò)進(jìn)行訓(xùn)練前，先對訓(xùn)練樣本矩形標(biāo)簽的長寬比值進(jìn)行聚類，并對生成初始候選框的面積加以調(diào)整。這樣在Faster R-CNN網(wǎng)絡(luò)后面的全連接層部分，在分類層得分大于0.6的初始候選框會進(jìn)入后面的位置回歸層進(jìn)行候選框的四個邊框的精修。由于改進(jìn)后的Faster R-CNN網(wǎng)絡(luò)模型生成的初始候選框更加符合榆紫葉甲蟲本身的形態(tài)學(xué)特征，所以初始候選框四周冗余較標(biāo)準(zhǔn)網(wǎng)絡(luò)生成的初始候選框少。對初始候選框四邊框精修時，四個邊框從初始位置平移到標(biāo)準(zhǔn)位置時的平移量較少，使得復(fù)雜背景對框定的干擾較少，從而能夠更加準(zhǔn)確地判定出邊框精修的終止位置，達(dá)到更為準(zhǔn)確的對目標(biāo)榆紫葉甲蟲進(jìn)行框定的目的。

3 實驗

3.1 實驗環(huán)境

本文實驗軟件平臺采用linux Ubuntu 16.04 LTS系統(tǒng)、label Image 1.3.2、pycharm 2.7.3、openCV 2.4.6，硬件環(huán)境是4 GB內(nèi)存，Intel?Core?i5-5450M CPU@2.50 GHz，主頻3.4 GHz的計算機(jī)。

3.2 實驗數(shù)據(jù)

實驗數(shù)據(jù)采集地點為東北林業(yè)大學(xué)帽兒山實驗林場，以其30～40年生的大葉榆、大果榆、裂葉榆、春榆等榆樹上榆紫葉甲蟲為研究對象。拍攝到3 000張像素大小為4 000×6 000的圖片數(shù)據(jù)。對數(shù)據(jù)集進(jìn)行篩選以避免錯誤、重復(fù)和模糊的圖像并人工標(biāo)注，同時按照7∶1∶2的比例隨機(jī)抽取訓(xùn)練數(shù)據(jù)集、驗證數(shù)據(jù)集、測試數(shù)據(jù)集。

3.3 實驗結(jié)果及分析

根據(jù)本文樣本數(shù)據(jù)集得出聚類中心數(shù)kopt=3，中心點分別為A1=0.64、A2=0.99、A3=1.27。用三個新的生成初始候選框的長寬比值代替標(biāo)準(zhǔn)網(wǎng)絡(luò)的生成初始候選框的長寬比值0.5、1、2，網(wǎng)絡(luò)改進(jìn)前后的初始候選框生成的尺寸和比例框圖如圖2所示。其中，圖2（a）為按照標(biāo)準(zhǔn)初始候選框生成比例和尺寸規(guī)則生成的候選框情況，圖2（b）為加入聚類策略和調(diào)整了初始候選框的尺寸的網(wǎng)絡(luò)的生成候選框情況。兩圖均為以O(shè)點為中心生成的候選框，圖2（a）中的三個藍(lán)色框為最大尺寸的三種比例候選框，三個紅色框為中間尺寸的三種比例候選框，三個黃色框為最小尺寸的三種比例的候選框。

圖2 網(wǎng)絡(luò)改進(jìn)前后的初始候選框生成比例和尺寸框圖

從圖2中可以看出，藍(lán)色候選框尺寸對于榆紫葉甲蟲來說冗余過大，在后面全連接層精修回歸過程中易受復(fù)雜背景中諸多特征的干擾，造成誤圈。紅色框和黃色框的長寬比例不符合榆紫葉甲蟲本身的長寬比形態(tài)學(xué)特征，在邊框精修時四個邊框平移的位移過大，易受復(fù)雜背景的干擾。在某一局部區(qū)域計算的損失小于規(guī)定閾值時就會停止回歸，造成邊框平移中斷或者平移過大，出現(xiàn)框定過大或者過小的情況。圖2（b）中去掉了藍(lán)色尺寸的三種比例的候選框。同時，針對紅色框和黃色框的情況，在相同尺寸下利用聚類生成新的初始候選框的長寬比，使其更加符合榆紫葉甲蟲本身的生物學(xué)特性。在后面的回歸過程受復(fù)雜的背景影響較小，可排除大部分復(fù)雜背景干擾圈定更為準(zhǔn)確。

榆紫葉甲蟲和葉片豁口或孔洞相鄰時，標(biāo)準(zhǔn)網(wǎng)絡(luò)模型與加入聚類策略和調(diào)整了生成候選框尺寸Faster RCNN網(wǎng)絡(luò)輸出結(jié)果對比圖如圖3所示。其中，圖3（a）～（e）為標(biāo)準(zhǔn)網(wǎng)絡(luò)的檢測效果圖，圖3（f）～（j）為改進(jìn)后的網(wǎng)絡(luò)輸出結(jié)果圖。圖3（a）與（f）、（b）與（g）、（c）與（h）、（d）與（i）、（e）與（j）分別為同一張圖像的兩個模型的輸出結(jié)果。圖中黃色圓圈部分為和榆紫葉甲蟲相鄰的葉片豁口或孔洞。

從圖3（a）～（e）可以看出，由于榆紫葉甲蟲的形狀、顏色等特征與葉片豁口相似，會出現(xiàn)框定榆紫葉甲蟲范圍時將臨近的葉片豁口也框定進(jìn)去，造成框定范圍冗余的現(xiàn)象。并且由于榆紫葉甲蟲的甲殼反光，標(biāo)準(zhǔn)模型會將一部分榆紫葉甲蟲識別為背景，出現(xiàn)框定不全的現(xiàn)象。而改進(jìn)后的網(wǎng)絡(luò)模型，減少了葉片豁口對榆紫葉甲蟲框定結(jié)果的影響，框定范圍比較準(zhǔn)確，如圖3（f）～3（j）所示。實驗結(jié)果證明，經(jīng)過對生成初始候選框網(wǎng)絡(luò)的改進(jìn)，在框定有葉片豁口干擾的榆紫葉甲蟲圖片時，取得了優(yōu)于標(biāo)準(zhǔn)候選框的效果。

相鄰榆紫葉甲蟲框定結(jié)果對比，如圖4所示。其中，圖4（a）～（e）為標(biāo)準(zhǔn)網(wǎng)絡(luò)的檢測效果圖，圖4（f）～（j）為改進(jìn)候選框生成網(wǎng)絡(luò)的輸出結(jié)果圖。圖4（a）與（f）、（b）與（g）、（c）與（h）、（d）與（i）、（e）與（j）分別為同一張圖像的兩個模型的輸出結(jié)果。

從圖4中可以看出，圖4（a）將兩只榆紫葉甲蟲框定到一起；圖4（b）中將一只榆紫葉甲蟲包含于另一只的范圍內(nèi)；圖4（c）中兩只榆紫葉甲蟲范圍均未框定完全；圖4（d）中框定的榆紫葉甲蟲范圍冗余過大；圖4（e）中的兩只榆紫葉甲蟲范圍均未圈定完全。這是由于兩只榆紫葉甲蟲特征相似且位置相臨，相似的特征對框定結(jié)果造成干擾，而改進(jìn)候選框生成網(wǎng)絡(luò)模型有效的解決了標(biāo)準(zhǔn)網(wǎng)絡(luò)模型出現(xiàn)的問題，如圖4（f）～（j）。由此可見，在兩只榆紫葉甲蟲相鄰的情況中，改進(jìn)后的網(wǎng)絡(luò)比標(biāo)準(zhǔn)網(wǎng)絡(luò)框定效果更好。

圖3 榆紫葉甲蟲和葉片豁口或孔洞相鄰時網(wǎng)絡(luò)模型改進(jìn)前后輸出效果對比圖

圖4 相鄰榆紫葉甲蟲框定結(jié)果對比圖

圖5 候選框生成網(wǎng)絡(luò)改進(jìn)前后的榆紫葉甲蟲和與其特征相似的昆蟲識別效果對比

由于榆樹花甲蟲和盾瘤胸葉甲蟲也是榆樹常見害蟲之一，且和榆紫葉甲蟲形態(tài)類似，在采集圖像時很容易采集到這兩種蟲，造成錯誤框定，如圖5所示。其中，圖5（a）～（e）為標(biāo)準(zhǔn)網(wǎng)絡(luò)的檢測效果圖，圖5（f）～（j）為改進(jìn)候選框生成網(wǎng)絡(luò)輸出結(jié)果圖。圖5（a）與（f）、（b）與（g）、（c）與（h）、（d）與（i）、（e）與（j）分別為同一張圖像的兩個模型的輸出結(jié)果，圖5（a）與（f）、（b）與（g）、（c）與（h）為盾瘤胸葉甲，圖5（d）與（i）、（e）與（j）為榆樹花甲蟲。由實驗結(jié)果可以看出，標(biāo)準(zhǔn)網(wǎng)絡(luò)模型存在誤將與榆紫葉甲蟲形態(tài)相似的昆蟲框定出來的現(xiàn)象，而改進(jìn)候選框生成網(wǎng)絡(luò)后的模型不會輸出這樣的結(jié)果圖，只將榆紫葉甲蟲框定出來，其余種類的昆蟲將不會框定出來。

分類實驗中，對于結(jié)果的處理，一般僅用一種指標(biāo)很難得到對算法的正確評估的。所以，一般用精準(zhǔn)率（Precision，P），召回率（Recall，R）來共同對算法進(jìn)行評估。衡量的最終指標(biāo)是識別的平均精度值A(chǔ)P。一般將以召回率和準(zhǔn)確率為橫縱軸的PR曲線與兩坐標(biāo)軸圍成的面積作為衡量指標(biāo)。

準(zhǔn)確率的計算公式如下：

召回率的計算公式如下：

TP（True Positives）是模型預(yù)測為正樣本實際為正樣本的特征數(shù)量，F(xiàn)P（False Positives）是模型預(yù)測為正樣本但實際是負(fù)樣本的特征數(shù)量。TN（True Nagetives）是模型預(yù)測為負(fù)樣本實際也為負(fù)樣本的特征數(shù)量，F(xiàn)N（False Nagetives）是模型預(yù)測為負(fù)樣本但實際為正樣本的特征數(shù)。其中，準(zhǔn)確率P指檢索出的榆紫葉甲蟲數(shù)量與此次檢索得到的檢測框的數(shù)量的比值，衡量的是查準(zhǔn)率。召回率R指檢索出的榆紫葉甲蟲數(shù)量與數(shù)據(jù)集中所有榆紫葉甲的數(shù)量的比值，衡量的是查全率。因為準(zhǔn)確率和召回率不可兼得，本文選擇AP作為衡量實驗效果好壞的指標(biāo)。將曲線與橫縱坐標(biāo)圍成的面積作為AP的值。

改進(jìn)前后模型的PR（Precision-Recall）曲線如圖6所示。

圖6 改進(jìn)候選框生成網(wǎng)絡(luò)模型和標(biāo)準(zhǔn)網(wǎng)絡(luò)模型的PR曲線

從圖6中可以看出，改進(jìn)后的網(wǎng)絡(luò)框架測試結(jié)果的AP曲線與橫縱坐標(biāo)軸的面積大于標(biāo)準(zhǔn)框架輸出的測試結(jié)果與坐標(biāo)軸的面積，即改進(jìn)后的網(wǎng)絡(luò)框架的輸出平均精度高于標(biāo)準(zhǔn)網(wǎng)絡(luò)框架。同時通過圖3～圖5實驗輸出效果對比圖可以看出，改進(jìn)后的模型輸出的識別效果圖圈定的目標(biāo)范圍更加精準(zhǔn)，且不會出現(xiàn)誤識別其他種類昆蟲的情況，識別精度從90.58%提升至94.73%。

4 結(jié)束語

本文提出了一種基于聚類算法的自適應(yīng)Faster R-CNN網(wǎng)絡(luò)模型用以榆紫葉甲蟲的識別：通過利用K-means算法結(jié)合BWP指標(biāo)對訓(xùn)練數(shù)據(jù)的矩形標(biāo)簽的長寬比值進(jìn)行聚類，利用聚類中心點代替標(biāo)準(zhǔn)網(wǎng)絡(luò)中生成初始候選框的長寬比例，使生成的初始候選框更加符合榆紫葉甲蟲長寬比形態(tài)學(xué)特征，減少邊框精修時的平移量，提高了識別精度。實驗結(jié)果表明，本文提出的基于K-means的Faster R-CNN算法有效解決了傳統(tǒng)識別算法中特征提取模板的局限性，并針對初始候選框不貼合待識別目標(biāo)造成的誤差加以改進(jìn)，提高了識別準(zhǔn)確率。改進(jìn)后的網(wǎng)絡(luò)在框定單只榆紫葉甲蟲、榆紫葉甲蟲與葉片豁口或孔洞相鄰、榆紫葉甲蟲相鄰時和其他與榆紫葉甲特征類似種類的昆蟲的框定效果均優(yōu)于標(biāo)準(zhǔn)網(wǎng)絡(luò)。測試相同數(shù)據(jù)集時的準(zhǔn)確率優(yōu)于標(biāo)準(zhǔn)網(wǎng)絡(luò)模型，證明了改進(jìn)后的網(wǎng)絡(luò)的優(yōu)越性。