黃亦其 李婕 趙建曄

摘要 以入侵植物薇甘菊高光譜圖像為研究對象，基于4種預(yù)處理方法對薇甘菊高光譜圖像進行降低噪聲處理，分別研究了基于主成分分析的特征提取方法和基于BP神經(jīng)網(wǎng)絡(luò)的分類模型，篩選出薇甘菊高光譜識別的最優(yōu)預(yù)處理方法，以實現(xiàn)薇甘菊的快速準(zhǔn)確識別。結(jié)果顯示，預(yù)處理方法為一階、二階微分的識別率分別為81.2%和76.92%;標(biāo)準(zhǔn)正態(tài)變量變換（SNV）和一階微分+SG平滑的識別率分別為89.74%和87.18%。多次試驗得到基于SNV預(yù)處理方法的識別率最穩(wěn)定，即得到最優(yōu)預(yù)處理方法為SNV。

關(guān)鍵詞 高光譜技術(shù);薇甘菊目標(biāo)識別;特征集選取;BP神經(jīng)網(wǎng)絡(luò)模型

中圖分類號 S126 ?文獻標(biāo)識碼 A

文章編號 0517-6611（2020）05-0246-04

doi：10.3969/j.issn.0517-6611.2020.05.069

開放科學(xué)（資源服務(wù)）標(biāo)識碼（OSID）：

Abstract With invasive plants Mikania mikania Kunt hyperspectral image as the research object，based on four kinds of pretreatment method of M.micrantha hyperspectral image noise reduction processing，we studied the feature extraction method based on principal component analysis and the classification model based on BP neural network.The optimal pretreatment method for M.micrantha hyperspectral identification was screened，in order to realize fast and exact recognition of M.micrantha.The experimental results showed that the recognition rates of first and second orders were 81.2% and 76.92%，respectively.The recognition rates of standard normal variable transformation （SNV） and first order differential +SG smoothing were 89.74% and 87.18%，respectively.Multiple experiments showed that the recognition rate of SNVbased pretreatment method was stable，in other words，the optimal pretreatment method was SNV.

Key words Hyperspectral technology;Mikania micrantha Kunt target recognition;Feature set selection;BP neural network

薇甘菊（Mikania micrantha Kunt）與附主植物在可見光下視覺辨識度低，分布無規(guī)律，受復(fù)雜多變的野外環(huán)境限制，一般檢測方法很難實現(xiàn)對薇甘菊實時精準(zhǔn)監(jiān)測。人工踏查[1]是比較普遍的針對薇甘菊的監(jiān)測方法之一，但其受地形影響，無法對薇甘菊的分布及危害進行準(zhǔn)確的定量評測。由于具有光譜分辨率高、波段多、信息量豐富等特點，高光譜監(jiān)測技術(shù)應(yīng)用于薇甘菊入侵的監(jiān)測上具有很大的優(yōu)勢。

目前，監(jiān)測方法主要有圖像識別法、高光譜監(jiān)測法、雷達數(shù)據(jù)輔助識別法和中低分辨率時序序列數(shù)據(jù)分析法等[2]。 Chance等[3]利用小型光譜成像儀（CASI）的圖像，在加拿大不列顛哥倫比亞省薩里的非森林植被環(huán)境中提取喜馬拉雅黑莓（Rubus armeniacus）、英國常春藤（Hedera helix）和其他物種的光譜，對所有的光譜通過光譜通道選擇算法處理，結(jié)果表明對英國常青藤和喜馬拉雅黑莓的檢測準(zhǔn)確率分別為80.0%和76.4%。而Narumalani等[4]對光譜圖像進行去噪處理，基于光譜角匹配法檢測了某區(qū)域內(nèi)4種入侵物種及其分布，總體準(zhǔn)確率為74%。

高光譜數(shù)據(jù)降維方法主要分為2種[5]：一種是基于變換的特征提取法，即通過數(shù)據(jù)變換產(chǎn)生一個新的特征空間，將全部波段的特征屬性完成從高維空間到低維空間轉(zhuǎn)換，從而保留盡可能高的分類精度。特征提取中典型的算法有包括主成分分析（principal components analysis，PCA）[6]、奇異值分解（singular value decomposition，SVD）[7]等。另一類是基于非變換的特征選擇，也稱局部特征波段選取。常見的波段選擇方法包括基于信息量的波段選擇以及基于類間可分性的波段選擇[5]。鑒于此，筆者使用主成分分析法完成特征提取，以較低的維數(shù)達到較好的識別結(jié)果，提高識別效率。

準(zhǔn)確地識別外來入侵植物是對其開展監(jiān)測工作的前提與基礎(chǔ)，任智偉等[5]通過將基于用PCA法和信息量改進PCA法對高光譜圖像數(shù)據(jù)降維的結(jié)果輸入最小距離分類器比較其分類精度來對比分析2種降維方法的效果。喻俊等[8]先利用S-Golay濾波剔除異常光譜數(shù)據(jù)，再采用光譜微分法與植被指數(shù)法構(gòu)建光譜特征最后應(yīng)用人工神經(jīng)網(wǎng)絡(luò)法結(jié)合因子分析法對典型植被進行分類識別。而吳培強等[9]以黃河口濕地為研究區(qū)分別建立支持向量機、人工神經(jīng)網(wǎng)絡(luò)和光譜角制圖3種分類模型開展利用全波段的和所選擇特征波段的分類對比試驗。筆者通過成像光譜儀獲取以薇甘菊為識別目標(biāo)的高光譜圖像，人工篩選獲取薇甘菊類和非薇甘菊類高光譜數(shù)據(jù)來用作識別薇甘菊目標(biāo)的研究，并基于主成分分析方法（PAC）對其做適當(dāng)?shù)慕稻S處理。

1 薇甘菊數(shù)據(jù)的采集

高光譜圖像的采集采用人工手持德國Cubert公司生產(chǎn)的S185高速成像光譜儀。共采集到薇甘菊高光譜圖像27張，分辨率為1 000×1 000。部分薇甘菊高光譜圖像樣本如圖1所示。觀察統(tǒng)計顯示，樣本中包括3類其他植物葉子以及其中1類植物的花，還有一些黑色背景。

Cube Ware是與S185配套的光譜分析軟件。通過易于操作的用戶界面來記錄和分析高光譜圖像。主要用到信號即圖像處理模塊，即從Cube Ware中獲取的高光譜數(shù)據(jù)樣本的預(yù)處理，數(shù)據(jù)降維、BP神經(jīng)網(wǎng)絡(luò)模型的構(gòu)建等都在Matlab 2018a中實現(xiàn)。

1.1 高光譜數(shù)據(jù)樣本獲取

通過Cube Ware軟件打開樣本集的薇甘菊高光譜圖像，隨機從18張圖像中以像素點為單位采集高光譜數(shù)據(jù)，每個像素點對應(yīng)1條高光譜曲線，每1條曲線所對應(yīng)的光譜數(shù)據(jù)都能以ASCII碼的格式保存（圖2）。共采集到高光譜數(shù)據(jù)樣本720個，該研究主要對薇甘菊進行識別，因此圖片中出現(xiàn)的其他類植物分別記為其他類植物1，其他類植物2，其他類植物3葉子和其他類植物3葉子花，具體的樣本采集數(shù)目見表1。

在主成分分析的基礎(chǔ)上還可以通過主成分得分圖初步地對樣本進行分析判別。由圖5可知，薇甘菊與非薇甘菊類（圖中的其他類）樣本均以坐標(biāo)原點為中心聚集分布，證明第一、二主成分對高光譜數(shù)據(jù)具有一定的聚類效果。比較發(fā)現(xiàn)，標(biāo)準(zhǔn)正態(tài)變量變換聚類效果較好。同時，各主成分得分圖中非薇甘菊類樣本比薇甘菊類聚集更緊密，聚類效果更明顯，但2類樣本間交叉重疊明顯，還需要基于支持向量機分類器對其進行精確分類。

4 BP神經(jīng)網(wǎng)絡(luò)建模

BP（back propagation）神經(jīng)網(wǎng)絡(luò)是1986年由Rumelhart和McClelland為首的科學(xué)家提出的概念，是一種按照誤差逆向傳播算法訓(xùn)練的多層前饋神經(jīng)網(wǎng)絡(luò)，是目前應(yīng)用最廣泛的神經(jīng)網(wǎng)絡(luò)[16]。BP神經(jīng)網(wǎng)絡(luò)是一種多層前饋，并通過誤差反傳訓(xùn)練的一種網(wǎng)絡(luò)，其主要方法是梯度下降法，使網(wǎng)絡(luò)的實際輸出值和期望輸出值的誤差均方差達到最小[16]。該研究使用了BP算法信號的前向傳播和誤差的反向傳播2個部分，計算誤差輸出時按從輸入到輸出的方向進行，而調(diào)整權(quán)值和閾值則從輸出到輸入的方向進行[17]。BP神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)如圖6所示。

將4種預(yù)處理后的主成分作為BP神經(jīng)網(wǎng)絡(luò)的輸入建立神經(jīng)網(wǎng)絡(luò)模型，并相互對比，確定最優(yōu)分類準(zhǔn)確率的預(yù)處理方法。BP神經(jīng)網(wǎng)絡(luò)建模后，進行參數(shù)調(diào)節(jié)，隱藏層均設(shè)置為10，輸出層神經(jīng)元個數(shù)均設(shè)置為2，而后調(diào)試隱藏層和輸出層學(xué)習(xí)率、訓(xùn)練次數(shù)等重要參數(shù)。對神經(jīng)網(wǎng)絡(luò)的參數(shù)進行反復(fù)調(diào)整直到獲得最佳分類準(zhǔn)確率?；?種預(yù)處理方法分別建立4個神經(jīng)網(wǎng)絡(luò)模型。建模時按2∶1分為訓(xùn)練集和測試集。

從表2可以看出，基于SNV和一階微分+SG預(yù)處理方法進行分類得到的分類準(zhǔn)確率較優(yōu)，但是就維度而言，SNV預(yù)處理方法的維度降到9，比一階微分+S-G卷積平滑預(yù)處理方法降維程度要好。通過比較得到二階微分對薇甘菊的識別效果最差。

從圖7和8可以看出，通過SNV方法處理過得薇甘菊數(shù)據(jù)，得到的識別程度相比于一階微分+S-G卷積平滑的更加穩(wěn)定。結(jié)合上述分析，得出基于SNV對薇甘菊的識別效果更好。并得到最優(yōu)的預(yù)處理方法為SNV，且主成分個數(shù)為9時，對薇甘菊的識別效果最佳。

5 結(jié)論

該研究主要以外來入侵植物薇甘菊為研究對象，對野外環(huán)境下高光譜數(shù)據(jù)進行分析研究。研究了高光譜數(shù)據(jù)的異常數(shù)據(jù)樣本剔除及預(yù)處理方法，基于高光譜數(shù)據(jù)降維特征提取方法，基于BP神經(jīng)網(wǎng)絡(luò)的模型建立方法。該研究得出了以下結(jié)論：

（1）采集野外薇甘菊高光譜圖像過程中，應(yīng)用馬氏距離剔除法剔除異常樣本，分別剔除了薇甘菊類9個異常樣本，以及非薇甘菊類9個異常樣本，得到用于后續(xù)研究的光譜數(shù)據(jù)樣本共702組，薇甘菊類351組，非薇甘菊類351組。

（2）通過4種預(yù)處理方法對高光譜數(shù)據(jù)預(yù)處理，從而達到了降噪、消除基線漂移等預(yù)期的效果，但其對識別模型建立的影響還需要進一步驗證分析。

（3）對上述所得4個特征集分別進行了基于主成分分析法的高光譜數(shù)據(jù)降維，得到最大主成分為前25個主成分用于識別薇甘菊的特征集。比較分析顯示，標(biāo)準(zhǔn)正態(tài)變量預(yù)處理方法使得高光譜數(shù)據(jù)聚類效果更好。

（4）對降維方法降維后特征集建立基于BP神經(jīng)網(wǎng)絡(luò)分類模型。在針對主成分分析降維后的4個特征集分別建立BP神經(jīng)網(wǎng)絡(luò)分類模型，確定了4個特征集中最終的主成分數(shù)目以達到特征集進一步降維的目的，并且得出所選主成分數(shù)目對應(yīng)的測試集識別準(zhǔn)確率。其中SNV-BP對應(yīng)的測試集不但識別準(zhǔn)確率高且其所建立的BP分類模型最穩(wěn)定，最終確定基于前9個主成分數(shù)目基礎(chǔ)上測試集識別準(zhǔn)確率最佳為89.74%，從而驗證了BP神經(jīng)網(wǎng)絡(luò)能夠?qū)σ巴廪备示者M行識別。

參考文獻

[1]柳帥，林輝，孫華，等.基于Pleiades-1衛(wèi)星數(shù)據(jù)薇甘菊信息提取[J].中南林業(yè)科技大學(xué)學(xué)報，2014，34（11）：116-119.

[2]孫玉芳，姜麗華，李剛，等.外來植物入侵遙感監(jiān)測預(yù)警研究進展[J].中國農(nóng)業(yè)資源與區(qū)劃，2016，37（8）：223-229.

[3]CHANCE C M，COOPS N C，CROSBY K，et al.Spectral wavelength selection and detection of two invasive plant Species in an urban area[J].Canadian journal of remote sensing，2016，42（1）：27-40.

[4] NARUMALANI S，MISHRA D R，WILSON R，et al.Detecting and mapping four invasive species along the floodplain of North Platte River，Nebraska[J].Weed technology，2009，23（1）：99-107.

[5]任智偉，吳玲達.基于信息量改進主成分分析的高光譜圖像特征提取方法[J].兵器裝備工程學(xué)報，2018，39（7）：151-154.

[6]張兵.高光譜圖像處理與信息提取前沿[J].遙感學(xué)報，2016，20（5）：1062-1090.

[7]VILLA A，CHANUSSOT J，JUTTEN C，et al.On the use of ICA for hyperspectral image analysis[C]//2009 IEE international geoscience and remote sensing symposium.Cape Town，South Africa：IEEE，2009.

[8]喻俊，李曉敏，張權(quán)，等.基于實測高光譜數(shù)據(jù)的太湖湖濱帶典型植被分類[J].江蘇農(nóng)業(yè)科學(xué)，2017，45（5）：240-244.

[9]吳培強，張杰，馬毅，等.基于地物光譜可分性的CHRIS高光譜影像波段選擇及其分類應(yīng)用[J].海洋科學(xué)，2015，39（2）：20-24.

[10]閔順耕，李寧，張明祥.近紅外光譜分析中異常值的判別與定量模型優(yōu)化[J].光譜學(xué)與光譜分析，2004，24（10）：1205-1209.

[11]FILZMOSER P，HRON K.Outlier detection for compositional data using robust methods[J].Mathematical geosciences，2008，40（3）：233-248.

[12]況潤元，曾帥，趙哲，等.基于實測高光譜數(shù)據(jù)的鄱陽湖濕地植被光譜差異波段提取[J].湖泊科學(xué)，2017，29（6）：1485-1490.

[13]陳彥兵，況潤元，曾帥.基于高光譜數(shù)據(jù)的鄱陽湖濕地典型植被識別分析[J].人民長江，2018，49（20）：19-23.

[14]顧紹紅，王永生，王光霞.主成分分析模型在數(shù)據(jù)處理中的應(yīng)用[J].測繪科學(xué)技術(shù)學(xué)報，2007，24（5）：387-390.

[15]唐紅，鄭文斌，李憲霞.主成分分析在光全散射特征波長選擇中的應(yīng)用[J].光學(xué)精密工程，2010，18（8）：1691-1698.

[16]邱晨，羅璟，趙朝文，等.基于BP神經(jīng)網(wǎng)絡(luò)的空氣質(zhì)量模型分類預(yù)測研究[J].軟件，2019，40（2）：129-132.

[17]來純曉，武振國，金松林，等.基于 BP 神經(jīng)網(wǎng)絡(luò)的小麥抗寒性模型構(gòu)建[J].河南科技學(xué)院學(xué)報（自然科學(xué)版），2019，47（3）：72-78.