孫 鈺 脫小倩 蔣 琦 張海燕 陳志泊 宗世祥 駱有慶

(1. 北京林業(yè)大學(xué)信息學(xué)院 北京 100083； 2. 北京林業(yè)大學(xué)林學(xué)院 北京 100083；3. 北京航空航天大學(xué)網(wǎng)絡(luò)空間安全學(xué)院 北京 100191)

害蟲識(shí)別是蟲害預(yù)測(cè)預(yù)報(bào)的首要工作，對(duì)害蟲綜合治理具有重要意義(陳梅香等， 2015)。鉆蛀性害蟲以林木為寄主蛀食成孔洞或隧道，直接危害林木主干和主梢生長(zhǎng)，輕則阻礙林木養(yǎng)分、水分傳導(dǎo)，重則導(dǎo)致林木枯萎死亡(馮國(guó)民, 2011; 高曉兵, 2010; 王曉園, 2011)。目前，蟲害監(jiān)測(cè)一般利用人工或圖像識(shí)別技術(shù)監(jiān)測(cè)成蟲(Yaoetal.， 2012; 孫鈺等， 2018; 張怡, 2017)，然而鉆蛀性害蟲幼蟲在取食階段對(duì)寄主的危害更早也更嚴(yán)重，傳統(tǒng)監(jiān)測(cè)方法難以發(fā)現(xiàn)隱蔽生活的幼蟲，早期預(yù)警能力有限。隨著聲音識(shí)別技術(shù)不斷發(fā)展，利用幼蟲鉆蛀振動(dòng)信號(hào)進(jìn)行蟲害監(jiān)測(cè)，具有高效、簡(jiǎn)單、成本低、預(yù)警時(shí)間早等優(yōu)勢(shì)，應(yīng)用前景廣闊(卜宇飛等， 2017; 韋雪青等， 2010; 祁驍杰, 2016)。

國(guó)內(nèi)外基于聲音識(shí)別技術(shù)對(duì)害蟲進(jìn)行偵聽監(jiān)測(cè)，研究對(duì)象多聚焦于木材檢疫害蟲(許小芳等, 2011)、倉(cāng)儲(chǔ)害蟲(郭敏等， 2001; Njorogeetal.， 2016)、水果害蟲(Hansenetal.， 1988)、林業(yè)鉆蛀性害蟲(Mankinetal.， 2008)等隱蔽性活動(dòng)的害蟲。在偵聽監(jiān)測(cè)過程中，由于鉆蛀振動(dòng)信號(hào)在樹干-空氣界面會(huì)大幅衰減，因此不宜使用麥克風(fēng)直接采集傳播至空氣中的聲音(婁定風(fēng)等， 2013)，而需將壓電式振動(dòng)傳感器嵌入樹干直接偵聽(Mankinetal.， 2011)。目前，采集鉆蛀振動(dòng)信號(hào)的典型儀器是美國(guó)AEC公司搭配SP-1L接觸式探頭的便攜式聲音探測(cè)儀AED-2000/2010L。在國(guó)外，Herriick等(2012)采用AED探測(cè)儀對(duì)密閉和開放環(huán)境下的紅棕象甲(Rhyncophorusferrugineuss)幼蟲鉆蛀振動(dòng)進(jìn)行研究，發(fā)現(xiàn)在至少5～10 cm的距離內(nèi)能檢測(cè)到所有齡期的幼蟲；Dosunmu等(2014)采用AED探測(cè)儀利用Raven和DAVIS程序分析紅棕象甲鉆蛀振動(dòng)的時(shí)頻特征，實(shí)現(xiàn)了紅棕象甲幼蟲的識(shí)別；Mankin等(2016)采用AED探測(cè)儀在果園環(huán)境下對(duì)紅棕象甲和蛀犀金龜(Orycteselegans)幼蟲的鉆蛀振動(dòng)進(jìn)行研究，發(fā)現(xiàn)其脈沖特征相差較大，實(shí)現(xiàn)了2種幼蟲的識(shí)別；Njoroge 等(2017)在裝滿谷物的密閉玻璃罐中采用AED探測(cè)儀錄制米象(Sitophilusoryzae)成蟲的鉆蛀振動(dòng)，研究了密封貯藏對(duì)害蟲活動(dòng)和死亡率的影響；Mankin等(2018)采用AED探測(cè)儀在疑似被南美硬木鋸天牛(Mallodondasystomus)侵害的鱷梨(Perseaamericana)樹上進(jìn)行測(cè)試,從11株樹中正確識(shí)別所有4株被侵害樹木，僅誤報(bào)1株健康樹木。在國(guó)內(nèi)， 祁驍杰(2016)對(duì)楊(Populus)樹木段中不同數(shù)量的光肩星天牛(Anoplophoraglabripennis)幼蟲脈沖進(jìn)行研究，發(fā)現(xiàn)脈沖個(gè)數(shù)與幼蟲數(shù)量存在線性關(guān)系； 卜宇飛等(2016)采用AED探測(cè)儀對(duì)7種蛀干害蟲鉆蛀振動(dòng)進(jìn)行研究，總結(jié)了害蟲時(shí)域、頻域特征的基本規(guī)律，證明了構(gòu)建聲音偵聽系統(tǒng)的實(shí)用性； 卜宇飛等(2017)還采用AED探測(cè)儀對(duì)2種天牛的4類行為特征進(jìn)行探索，發(fā)現(xiàn)取食鉆蛀振動(dòng)脈沖時(shí)間短、振幅大且頻率高，可用于實(shí)際的聲音偵聽中。也有研究者使用美國(guó)農(nóng)業(yè)部Mankin團(tuán)隊(duì)錄制的昆蟲聲音庫(kù)，實(shí)現(xiàn)了基于梅爾倒譜系數(shù)和矢量量化的昆蟲聲音自動(dòng)鑒別識(shí)別(竺樂慶等， 2010)，且探索了基于梅爾倒譜系數(shù)和混合高斯模型的昆蟲聲音自動(dòng)識(shí)別方法(竺樂亦等， 2012)，均獲得較高識(shí)別精度。

鉆蛀振動(dòng)信號(hào)經(jīng)傳感器錄制后被保存為音頻格式，可使用聲音識(shí)別技術(shù)進(jìn)行分類。聲音識(shí)別的一個(gè)重要任務(wù)是關(guān)鍵詞檢測(cè)，即對(duì)聲音信號(hào)中目標(biāo)關(guān)鍵詞進(jìn)行識(shí)別。關(guān)鍵詞檢測(cè)與害蟲鉆蛀振動(dòng)識(shí)別類似，均以音頻中的短脈沖為識(shí)別對(duì)象。傳統(tǒng)關(guān)鍵詞檢測(cè)是基于隱馬爾可夫模型實(shí)現(xiàn)的(陳玉平等， 2008)，隨著以深度學(xué)習(xí)為代表的人工智能技術(shù)迅速發(fā)展，將其應(yīng)用于聲音識(shí)別領(lǐng)域也取得了重大突破。Chen等(2014)使用深度神經(jīng)網(wǎng)絡(luò)進(jìn)行關(guān)鍵詞檢測(cè)，與隱馬爾可夫模型相比識(shí)別性能提高了45%；Sainath等(2015)和Sun等(2016)分別將卷積神經(jīng)網(wǎng)絡(luò)和長(zhǎng)短期記憶網(wǎng)絡(luò)用于關(guān)鍵詞檢測(cè)，相比深度神經(jīng)網(wǎng)絡(luò)識(shí)別性能分別提高27%～44%和67.6%；利用谷歌公開的語(yǔ)音識(shí)別數(shù)據(jù)集(Google speech commands dataset)，Tang等(2018)研究深度殘差學(xué)習(xí)和空洞卷積在關(guān)鍵詞檢測(cè)的應(yīng)用，設(shè)計(jì)的ResNet15變體網(wǎng)絡(luò)得到了95.8%的精度； Zeng等(2019)使用DenseNet網(wǎng)絡(luò)提取局部特征和時(shí)間序列特征，得到了96.6%為精度。

本研究針對(duì)園林場(chǎng)景中側(cè)柏(Platycladusorientalis)和臭椿(Ailanthusaltissima)混交種植的情況，使用壓電式傳感器監(jiān)測(cè)主要蛀干害蟲雙條杉天牛(Semanotusbifasciatus)和臭椿溝眶象(Eucryptorrhynchusbrandti)幼蟲蛀干取食發(fā)出的振動(dòng)，使用聲音識(shí)別技術(shù)構(gòu)建鉆蛀振動(dòng)識(shí)別模型，設(shè)計(jì)并訓(xùn)練輕量級(jí)神經(jīng)網(wǎng)絡(luò)，自動(dòng)識(shí)別雙條杉天牛鉆蛀和臭椿溝眶象鉆蛀振動(dòng)，以期為提高鉆蛀性害蟲的早期預(yù)警能力提供技術(shù)支撐。

1 數(shù)據(jù)與處理

1.1 鉆蛀振動(dòng)錄音

害蟲鉆蛀振動(dòng)信號(hào)采集在實(shí)驗(yàn)室隔音箱內(nèi)進(jìn)行，采集工具為搭配SP-1L探頭的AED-2010L便攜式聲音探測(cè)儀(圖1)。SP-1L探頭為壓電式傳感器探頭(壓電晶體諧振頻率40 kHz，前置放大器帶寬1～50 kHz，放大倍數(shù)40 dB)，通過直徑6 mm金屬探針嵌入樹干，直接耦合鉆蛀振動(dòng)。AED-2010L主機(jī)連接錄音筆,設(shè)置采樣頻率為44.1 kHz，采樣精度為16 bit，錄制害蟲鉆蛀振動(dòng)信號(hào)。

圖1 鉆蛀振動(dòng)采集環(huán)境與工具Fig.1 The environment and tools for the boring vibration collection

從林場(chǎng)采集25段長(zhǎng)約30 cm的健康側(cè)柏木段，分成3組并于3月下旬分別接入不同數(shù)量的雙條杉天牛雌雄成蟲，以飼養(yǎng)初孵幼蟲；采集6段長(zhǎng)約1 m的健康臭椿木段，分成2組并于7月下旬分別接入不同數(shù)量的臭椿溝眶象幼蟲。于木段中間位置鉆孔，每隔5～7天在適宜溫度的時(shí)間段內(nèi)(9、10、11時(shí))利用AED-2010L便攜式聲音探測(cè)儀采集害蟲鉆蛀振動(dòng)信號(hào)(卜宇飛, 2016; 祁驍杰, 2016)。每投音頻時(shí)長(zhǎng)約90 s，保存為.wav格式。當(dāng)幼蟲進(jìn)入老熟階段停止進(jìn)食，害蟲鉆蛀振動(dòng)信號(hào)采集結(jié)束。害蟲鉆蛀實(shí)驗(yàn)室環(huán)境下，害蟲鉆蛀振動(dòng)信號(hào)主要分布在10 kHz內(nèi)(圖2)，具有持續(xù)時(shí)間短、能量集中等特點(diǎn)(卜宇飛等， 2017)。

1.2 數(shù)據(jù)集建立

整理各木段音頻，將其分為雙條杉天牛鉆蛀振動(dòng)、臭椿溝眶象鉆蛀振動(dòng)和無(wú)鉆蛀振動(dòng)3類。雙條杉天牛和臭椿溝眶象鉆蛀振動(dòng)的訓(xùn)練、測(cè)試集以各木段為單位進(jìn)行劃分，選取具有代表性的3段雙條杉天牛和2段臭椿溝眶象木段，將其所有音頻作為測(cè)試集，其余木段的音頻作為訓(xùn)練集，所有木段采集到的無(wú)鉆蛀振動(dòng)音頻以3∶2比例隨機(jī)劃分為訓(xùn)練集和測(cè)試集。最終數(shù)據(jù)集音頻的劃分?jǐn)?shù)量及提取用于訓(xùn)練的對(duì)數(shù)梅爾聲譜數(shù)量如表1所示，其中測(cè)試集中無(wú)鉆蛀振動(dòng)的對(duì)數(shù)梅爾聲譜數(shù)量指進(jìn)入第2分支檢測(cè)出的數(shù)量。

表1 數(shù)據(jù)集劃分?jǐn)?shù)量Tab.1 Number of data sets divided

1.3 鉆蛀振動(dòng)信號(hào)預(yù)處理

提取對(duì)數(shù)梅爾聲譜作為神經(jīng)網(wǎng)絡(luò)學(xué)習(xí)的音頻特征(Kahletal.， 2018)。對(duì)害蟲鉆蛀振動(dòng)信號(hào)的預(yù)處理包括端點(diǎn)檢測(cè)、時(shí)間規(guī)整和提取對(duì)數(shù)梅爾聲譜。害蟲鉆蛀振動(dòng)信號(hào)短時(shí)能量較高，具有間隔、短促、尖銳(卜宇飛等， 2016)的特性，首先對(duì)害蟲鉆蛀振動(dòng)信號(hào)進(jìn)行端點(diǎn)檢測(cè)獲取鉆蛀振動(dòng)信號(hào)片段，然后利用時(shí)間規(guī)整算法將信號(hào)片段規(guī)整至統(tǒng)一時(shí)長(zhǎng)，最后經(jīng)過時(shí)頻變換和對(duì)數(shù)操作提取對(duì)數(shù)梅爾聲譜。

圖2 害蟲鉆蛀振動(dòng)信號(hào)波形和聲譜Fig.2 Waveform and spectrogram of insect boring vibration signals

圖3 預(yù)處理示意Fig.3 The diagram of pretreatment

1.3.1 端點(diǎn)檢測(cè) 端點(diǎn)檢測(cè)是指從一段包含有效片段的連續(xù)音頻中找到能量較高的有效片段起止端點(diǎn)，本研究采用基于短時(shí)能量的單參數(shù)雙門限方法(黎煊等， 2018)對(duì)采集到的害蟲鉆蛀振動(dòng)信號(hào)進(jìn)行端點(diǎn)檢測(cè)。以200個(gè)采樣點(diǎn)(4.54 ms)為1幀，計(jì)算每幀的短時(shí)能量，與設(shè)定的低門限閾值(T1)和高門限閾(T2)進(jìn)行比較，確定有效音頻片段的端點(diǎn)(邢亞從, 2011)。在端點(diǎn)檢測(cè)時(shí)，一般選取音頻樣本前導(dǎo)無(wú)話段計(jì)算高、低門限閾值，本研究隨機(jī)選取5 s無(wú)鉆蛀振動(dòng)的音頻信號(hào)代替前導(dǎo)無(wú)話段進(jìn)行計(jì)算。如圖 3端點(diǎn)檢測(cè)部分所示，紅色豎線為害蟲鉆蛀振動(dòng)脈沖的起始位置，綠色豎線為該脈沖的結(jié)束位置，單參數(shù)雙門限端點(diǎn)檢測(cè)可有效檢測(cè)出害蟲鉆蛀脈沖。

1.3.2 時(shí)間規(guī)整 端點(diǎn)檢測(cè)后的信號(hào)長(zhǎng)度不一，無(wú)法滿足神經(jīng)網(wǎng)絡(luò)輸入層要求，需要利用時(shí)間規(guī)整算法將信號(hào)片段規(guī)整至統(tǒng)一時(shí)長(zhǎng)。經(jīng)試驗(yàn)探索，本研究將長(zhǎng)短不一的有效音頻片段規(guī)整至0.5 s。以端點(diǎn)檢測(cè)到的起始點(diǎn)向前5幀位置作為時(shí)間規(guī)整的起始點(diǎn)，向后0.5 s位置作為時(shí)間規(guī)整的終點(diǎn)，端點(diǎn)檢測(cè)出的有效片段長(zhǎng)度多于0.5 s則拆分，少于0.5 s則前后延長(zhǎng)。如圖3時(shí)間規(guī)整部分所示，在端點(diǎn)檢測(cè)基礎(chǔ)上,將有效音頻片段規(guī)整至統(tǒng)一時(shí)長(zhǎng)片段。

1.3.3 提取對(duì)數(shù)梅爾聲譜 對(duì)規(guī)整后時(shí)長(zhǎng)0.5 s的信號(hào)進(jìn)行預(yù)加重(竺樂慶等， 2012)(預(yù)加重因子為0.95)，以200點(diǎn)為1幀、0.8的幀移重疊比例進(jìn)行分幀，對(duì)分幀信號(hào)加漢寧窗以增加幀左右兩端的連續(xù)性，再以1 024點(diǎn)的長(zhǎng)度對(duì)每幀信號(hào)作短時(shí)傅里葉變換，使用64組梅爾濾波器組濾波并進(jìn)行取對(duì)數(shù)操作，最后得到對(duì)數(shù)梅爾聲譜作為神經(jīng)網(wǎng)絡(luò)最終的輸入特征。對(duì)數(shù)梅爾聲譜的特征維度由信號(hào)的幀數(shù)和梅爾濾波器組的個(gè)數(shù)決定，如圖3提取對(duì)數(shù)梅爾聲譜部分所示，本研究得到的對(duì)數(shù)梅爾聲譜特征維度為109×64。信號(hào)的幀數(shù)計(jì)算方法如下：

Lhop=?Nfft×(1-Poverlap)」；

(1)

式中：?」表示向下取整；「?表示向上取整；Lhop為每次幀移的步幅；Nfft為短時(shí)傅里葉變換長(zhǎng)度；Povertap為幀移重疊比例；Nframe為幀數(shù)；Lsig為信號(hào)長(zhǎng)度。

2 鉆蛀性害蟲聲音識(shí)別

2.1 識(shí)別模型

設(shè)計(jì)面向2種蛀干害蟲鉆蛀振動(dòng)識(shí)別的輕量級(jí)卷積神經(jīng)網(wǎng)絡(luò)InsectFrames(圖4a)，以維度(109，64)的對(duì)數(shù)梅爾聲譜為輸入，主體部分由4層卷積組成，在全連接層前接全局平均池化以防止過擬合。

圖4 InsectFrames網(wǎng)絡(luò)結(jié)構(gòu)Fig.4 Network structure of InsectFrames

在卷積部分設(shè)計(jì)2種卷積塊結(jié)構(gòu)和2種中間層維度。卷積塊結(jié)構(gòu)為： 1) ConvBlock1(圖 4b)， 卷積核尺寸為3×3，步長(zhǎng)為1，每層卷積后引入ReLU激活函數(shù)實(shí)現(xiàn)網(wǎng)絡(luò)的非線性變換。連接批量歸一化層避免梯度消失，同時(shí)加快收斂速度，批量歸一化后接步長(zhǎng)為2的最大池化層實(shí)現(xiàn)降維。2) ConvBlock2(圖4c)， 同樣設(shè)置卷積核尺寸為3×3，每層卷積后引入ReLU激活函數(shù)，連接批量歸一化層。與ConvBlock1不同的是，ConvBlock2未采用最大池化策略，而是在卷積過程中設(shè)置步長(zhǎng)為2，直接進(jìn)行降維操作。為了對(duì)比中間層維度對(duì)識(shí)別結(jié)果的影響，每種卷積塊結(jié)構(gòu)分別對(duì)應(yīng)2種輸出維度： 1) 4層卷積的輸出維度n1—n4分別設(shè)置為[8,16,32,64]； 2) 4層卷積的輸出維度n1—n4分別設(shè)置為[16,32,64,128]。

綜上，本研究基于InsectFrames提出4種具有不同卷積結(jié)構(gòu)的網(wǎng)絡(luò)變體結(jié)構(gòu)InsectFrames_1—4(表2)，如InsectFrames_1卷積結(jié)構(gòu)采用ConvBlock1提取特征并降維，4層卷積對(duì)應(yīng)的輸出通道分別設(shè)為[8,16,32,64]； InsectFrames_4卷積結(jié)構(gòu)采用ConvBlock2提取特征并降維，4層卷積對(duì)應(yīng)的輸出通道分別設(shè)為[16,32,64,128]。

表2 InsectFrames 4種變體結(jié)構(gòu)Tab.2 Four different variant structures of InsectFrames

2.2 識(shí)別流程

本研究設(shè)計(jì)的識(shí)別模型首先經(jīng)過訓(xùn)練，通過深度學(xué)習(xí)有監(jiān)督地優(yōu)化模型參數(shù)至算法收斂，然后在測(cè)試集上評(píng)估4種模型準(zhǔn)確率。

圖5 害蟲鉆蛀振動(dòng)識(shí)別流程Fig.5 Recognition process of insect boring vibration

在訓(xùn)練階段(圖5a)，預(yù)處理提取訓(xùn)練集音頻的對(duì)數(shù)梅爾聲譜輸入到卷積神經(jīng)網(wǎng)絡(luò)中，通過前向傳播預(yù)測(cè)類別，利用損失函數(shù)計(jì)算預(yù)測(cè)類別與真實(shí)類別之間的誤差，再通過網(wǎng)絡(luò)反向傳播進(jìn)行參數(shù)更新，經(jīng)多次迭代完成模型訓(xùn)練。模型訓(xùn)練過程中，采用Adam算法進(jìn)行優(yōu)化； 設(shè)置初始學(xué)習(xí)率為0.001，并以余弦衰減方式進(jìn)行衰減，直至0.000 001停止； 以32個(gè)音頻為1個(gè)批次，最多迭代50次； 使用交叉熵函數(shù)(cross entropy function)作為損失函數(shù)。

在測(cè)試階段(圖5b)，將未參與訓(xùn)練的測(cè)試集音頻輸入到訓(xùn)練完畢的模型，測(cè)試模型在新數(shù)據(jù)上的分類準(zhǔn)確率。在測(cè)試過程中，為了節(jié)約計(jì)算資源并提高識(shí)別精度，進(jìn)行端點(diǎn)檢測(cè)時(shí)設(shè)計(jì)了2個(gè)分支： 1) 當(dāng)端點(diǎn)檢測(cè)未檢測(cè)出信號(hào)片段時(shí)(即當(dāng)前音頻的能量偏低)，音頻直接被判定為無(wú)鉆蛀振動(dòng)類； 2) 當(dāng)端點(diǎn)檢測(cè)檢測(cè)出信號(hào)片段時(shí)，將信號(hào)片段轉(zhuǎn)化為對(duì)數(shù)梅爾聲譜，并輸入到已訓(xùn)練好的卷積神經(jīng)網(wǎng)絡(luò)進(jìn)行分類預(yù)測(cè)，此時(shí)待識(shí)別音頻得到多個(gè)片段的預(yù)測(cè)分類結(jié)果，利用多數(shù)投票的方法進(jìn)行預(yù)測(cè)結(jié)果匯聚，得票高的類別被認(rèn)為是該音頻的預(yù)測(cè)類別。

3 試驗(yàn)與結(jié)果

3.1 試驗(yàn)環(huán)境

本研究訓(xùn)練模型采用的硬件平臺(tái)為Intel CoreTMi7-6700K CPU(32 GB內(nèi)存)以及GeForce GTX 1080 Ti(12 GB顯存)； 軟件為Ubuntu 16.04 LTS 64位操作系統(tǒng)，Lasagne (Raffeletal.， 2016)深度學(xué)習(xí)開源框架。

3.2 試驗(yàn)結(jié)果

害蟲識(shí)別一般以單音頻為單位進(jìn)行，為了避免CPU與GPU之間的數(shù)據(jù)傳輸，本研究以音頻分類的平均精度(average accuracy)以及在CPU上測(cè)試的平均識(shí)別時(shí)間(average recognition time of CPU)作為害蟲鉆蛀振動(dòng)識(shí)別的評(píng)價(jià)指標(biāo)。待測(cè)音頻輸入網(wǎng)絡(luò)模型前需經(jīng)預(yù)處理轉(zhuǎn)化為對(duì)數(shù)梅爾聲譜，該過程是音頻識(shí)別的基礎(chǔ)工作，與識(shí)別方法的選擇無(wú)關(guān)。

由表 3可知，單個(gè)音頻的平均預(yù)處理時(shí)間均為1.8 s左右，在整個(gè)害蟲識(shí)別過程中時(shí)間占比較大。在圖像分類問題中表現(xiàn)優(yōu)異的ResNet18網(wǎng)絡(luò)用于害蟲鉆蛀振動(dòng)識(shí)別能夠達(dá)到88.89%的平均識(shí)別準(zhǔn)確率，但由于ResNet18具有相對(duì)復(fù)雜的網(wǎng)絡(luò)結(jié)構(gòu)，單個(gè)音頻的CPU平均識(shí)別時(shí)間為229.612 s。廣泛應(yīng)用于昆蟲聲音識(shí)別領(lǐng)域的混合高斯模型(Gaussian mixture model, GMM)時(shí)間效率極高，是4種網(wǎng)絡(luò)變體識(shí)別速度的15～60倍，但在此數(shù)據(jù)集上的識(shí)別精度僅為61.81%。本研究測(cè)試了4種變體結(jié)構(gòu)簡(jiǎn)單的卷積神經(jīng)網(wǎng)絡(luò)(InsectFrames_1—4)，其中在卷積過程中直接實(shí)現(xiàn)降維的InsectFrames_3、InsectFrames_4識(shí)別精度能夠達(dá)到90.28%和93.75%，與GMM、ResNet18相比識(shí)別精度有所提高； 改變特征降維方式，設(shè)置2×2最大池化步長(zhǎng)為2，對(duì)應(yīng)的InsectFrames_1、InsectFrames_2 的識(shí)別精度能夠達(dá)到92.36%和95.83%，與在卷積過程實(shí)現(xiàn)降維的策略相比均提高2.08%，具有更加明顯的識(shí)別優(yōu)勢(shì)。圖6 為識(shí)別性能最好的InsectFrames_2測(cè)試結(jié)果的混淆矩陣，其中有3個(gè)雙條杉天牛鉆蛀振動(dòng)音頻錯(cuò)分為臭椿溝眶象類，1個(gè)臭椿溝眶象鉆蛀振動(dòng)音頻錯(cuò)分為雙條杉天牛類，2個(gè)無(wú)鉆蛀振動(dòng)音頻被錯(cuò)分為雙條杉天牛類。

表3 識(shí)別結(jié)果Tab.3 Identification results

圖6 InsectFrames_2測(cè)試結(jié)果的混淆矩陣Fig.6 The confusion matrix of InsectFrames_20： 雙條杉天牛類 S. bifasciatus；1：臭椿溝眶象類E.brandti； 2：無(wú)鉆蛀振動(dòng)類 No boring vibration.

4 討論

當(dāng)前廣泛應(yīng)用于蟲聲識(shí)別領(lǐng)域的GMM以及在圖像分類領(lǐng)域性能較好的ResNet18用于害蟲鉆蛀振動(dòng)識(shí)別的精度僅為61.81%和88.89%，而且ResNet18運(yùn)算開銷過高，在CPU上平均識(shí)別時(shí)間是4種網(wǎng)絡(luò)變體的170～2 000多倍。而本研究實(shí)現(xiàn)的4種基于輕量級(jí)神經(jīng)網(wǎng)絡(luò)的鉆蛀振動(dòng)識(shí)別模型InsectFrames_1—4的測(cè)試集平均精度均能達(dá)到90%以上，采用最大池化降維且中間層特征維度較高的InsectFrames_2模型具有最佳識(shí)別性能，平均精度能夠達(dá)到95.83%，CPU平均識(shí)別時(shí)間為1.334 s，可利用鉆蛀振動(dòng)有效區(qū)分2種鉆蛀性幼蟲，為蟲情監(jiān)測(cè)預(yù)警提供技術(shù)支持。

傳統(tǒng)蟲聲識(shí)別領(lǐng)域廣泛應(yīng)用的GMM為無(wú)監(jiān)督學(xué)習(xí)，主要由數(shù)據(jù)驅(qū)動(dòng)實(shí)現(xiàn)聚類，在不同類型的昆蟲(儲(chǔ)糧害蟲、蛀干害蟲、土壤昆蟲、植物昆蟲等)識(shí)別中表現(xiàn)優(yōu)異；但雙條杉天牛、臭椿溝眶象均屬林業(yè)鉆蛀性害蟲，鉆蛀振動(dòng)信號(hào)特征相似度較高，利用GMM難以精準(zhǔn)區(qū)分。ResNet18采用殘差結(jié)構(gòu)優(yōu)化深層網(wǎng)絡(luò)，可解決網(wǎng)絡(luò)在深度增加過程中的退化問題，同時(shí)增加模型的表示能力，在相對(duì)復(fù)雜的圖像分類任務(wù)中性能優(yōu)異，對(duì)于鉆蛀振動(dòng)特征相似的2類害蟲也能取得較高識(shí)別精度。但鉆蛀振動(dòng)信號(hào)特征相對(duì)較少，ResNet18網(wǎng)絡(luò)結(jié)構(gòu)復(fù)雜且參數(shù)較多，所以時(shí)間效率偏低，無(wú)法真正應(yīng)用于蟲聲識(shí)別。本研究實(shí)現(xiàn)的4種輕量級(jí)神經(jīng)網(wǎng)絡(luò)模型很好解決了鉆蛀振動(dòng)信號(hào)特征相對(duì)較少且相似性高的問題。

由于最大池化策略能夠?qū)ι弦粚拥妮敵鲆刖植科揭撇蛔冃?，且可有效保留圖像紋理特征，采用最大池化策略的InsectFrames_1和InsectFrames_2網(wǎng)絡(luò)相較在卷積過程中直接實(shí)現(xiàn)降維策略的網(wǎng)絡(luò)具有更好的識(shí)別能力，但同時(shí)增加了運(yùn)算量，耗時(shí)增加約2倍，時(shí)間效率相對(duì)較低。4層卷積輸出維度為[16,32,64,128]的InsectFrames_2能夠獲取更多的中間特征用于學(xué)習(xí)，在試驗(yàn)中得到最好的識(shí)別結(jié)果。InsectFrames_2網(wǎng)絡(luò)結(jié)構(gòu)簡(jiǎn)單、泛化能力強(qiáng)，在CPU上單個(gè)音頻平均識(shí)別時(shí)間約1.3 s，能夠滿足正常害蟲監(jiān)測(cè)預(yù)警的時(shí)效要求，適用于蛀干害蟲的鉆蛀振動(dòng)識(shí)別。

5 結(jié)論

基于輕量級(jí)卷積神經(jīng)網(wǎng)絡(luò)的鉆蛀振動(dòng)識(shí)別方法，可高效識(shí)別出信號(hào)特征相似度較高的2種蛀干害蟲，將其應(yīng)用于林間蟲害監(jiān)測(cè)和防治領(lǐng)域，能夠提高監(jiān)測(cè)效率，降低監(jiān)測(cè)成本，為林業(yè)蛀干害蟲監(jiān)測(cè)提供新的研究思路和技術(shù)支撐。未來(lái)將探索在含噪森林環(huán)境中更多種類鉆蛀性幼蟲的自動(dòng)識(shí)別技術(shù)，提升對(duì)隱蔽鉆蛀性林業(yè)害蟲的早期預(yù)警能力。