為了改善空氣質(zhì)量、防風(fēng)固沙、減少水土流失等環(huán)境問題，中國(guó)政府進(jìn)行了大規(guī)模的植樹造林活動(dòng)。臭椿、云南松和欒樹作為中國(guó)常見的綠化樹種，自然而然地成為植樹造林的首選。由于我國(guó)綠化造林面積巨大，林業(yè)病蟲害幾乎是不可避免的問題，而解決這一問題要花費(fèi)巨大的人力物力。高光譜激光雷達(dá)探測(cè)系統(tǒng)可以快速同步獲得地物目標(biāo)的測(cè)距信息和光譜信息，可應(yīng)用于目標(biāo)的探測(cè)與識(shí)別等。這為研究一種精準(zhǔn)高效的林業(yè)病蟲害監(jiān)測(cè)方法提供了保障，因此開展基于高光譜激光雷達(dá)的林業(yè)病蟲害監(jiān)測(cè)方法的研究具有重要的意義。

1實(shí)驗(yàn)系統(tǒng)與分類方法?

1.1?新型可調(diào)諧高光譜激光雷達(dá)系統(tǒng)?

?

新型可調(diào)諧高光譜激光雷達(dá)系統(tǒng)主要由發(fā)射系統(tǒng)、接收探測(cè)系統(tǒng)、掃?描控?制系?統(tǒng)三?個(gè)主?要?部分組成。發(fā)射系統(tǒng)主要由激光器和聲光可調(diào)諧濾波器（AOTF）組成；接收探測(cè)系統(tǒng)主要包括接收光學(xué)系統(tǒng)、光電探測(cè)器和數(shù)據(jù)采集系統(tǒng)；掃描控制系統(tǒng)負(fù)責(zé)對(duì)發(fā)射光束進(jìn)行空間位置和激光出射方向的控制，實(shí)現(xiàn)預(yù)定的空間掃描。圖?１?為?高?光譜激光雷達(dá)原理，表１為高光譜激光雷達(dá)系統(tǒng)的主要參數(shù)。

圖１?高光譜激光雷達(dá)原理示意圖

表１高光譜激光雷達(dá)的主要參數(shù)

高光譜激光雷達(dá)系統(tǒng)首先將超連續(xù)激光光束注入AOTF中，經(jīng)?AOTF調(diào)諧濾波后按照設(shè)置要求在不同?時(shí)刻發(fā)射不同波長(zhǎng)的單色光，激光經(jīng)準(zhǔn)直擴(kuò)散鏡和反?射鏡出射，出射光束有小部分投?射到光束采樣鏡后進(jìn)入主波探測(cè)器中并?傳遞?到?數(shù)據(jù)采集系統(tǒng)，大部分光束通過掃?描轉(zhuǎn)鏡后投射到目標(biāo)上，目標(biāo)的回波信號(hào)由光學(xué)接收系統(tǒng)接收，再經(jīng)光電探測(cè)器轉(zhuǎn)化為電信號(hào)，由數(shù)據(jù)?采集系統(tǒng)保存至上位機(jī)。對(duì)于采集到的數(shù)據(jù)首先進(jìn)行預(yù)處理，即將采集到的樣本回波強(qiáng)度信息轉(zhuǎn)化為樣本光譜反射率信息。為了減少外部環(huán)境因素的影響，保證同一次實(shí)驗(yàn)的數(shù)據(jù)之間可?以橫向比較，使用回波反射率對(duì)每個(gè)測(cè)量目標(biāo)進(jìn)行歸一化?；夭ǚ瓷渎视赏还庾V通道的目標(biāo)回波強(qiáng)度與參考板回波強(qiáng)度的歸一化處理得到。待測(cè)目標(biāo)的反射率可表示為

式中：Ｉtarget和Ｉreference分別為待測(cè)目標(biāo)和參考靶標(biāo)的激光回波強(qiáng)度；Ｒreference為參考靶標(biāo)的反射率，參考靶標(biāo)的反射率可通過地物光譜儀測(cè)量得到。起始波長(zhǎng)為650ｍ，終止波長(zhǎng)為1100ｎｍ，波長(zhǎng)間隔為5ｎｍ，共91個(gè)光譜通道，參考靶標(biāo)可以選用不同反射率的白色參考板。

1.2?分類方法

支持向量機(jī)是一種二類分類模型，其基本模型是定義在特征空間上間隔最大化的線性分類器，在分類問題中給定輸入數(shù)據(jù)和學(xué)習(xí)目標(biāo)ｘ＝｛ｘ1，ｘ2，ｘ3，…，ｘn｝和ｙ＝｛ｙ1，ｙ2，ｙ3，…，ｙｎ｝，其中輸入數(shù)據(jù)的每個(gè)樣本都包含多個(gè)特征?并由?此構(gòu)?成?特征空間，而學(xué)習(xí)目標(biāo)為二元變量并分別表示負(fù)類和正類。若輸入數(shù)據(jù)所在的特征空間存在作為決策邊界的超平面，則將學(xué)習(xí)目標(biāo)按正類和負(fù)類分開，并使任意樣本的點(diǎn)?Ｘ?到平面距離不小于最小距離，則稱該分類問題具有線性可分性，超平面的法向量ω?和截距ｂ可表示為

但是支持向量機(jī)還是局限于線性可分的數(shù)據(jù)，當(dāng)引入核函數(shù)之后，則使其成為實(shí)質(zhì)上的非線性分類器。最?早?Cortes等提?出?了?線?性?支?持?向?量機(jī)，之后Stephan通過引入核技巧提出了非線性支持向量機(jī)。支持向量機(jī)是一種常用的數(shù)據(jù)分類和回歸的機(jī)器學(xué)習(xí)方法，該方法已成功地應(yīng)用于遙感監(jiān)測(cè)領(lǐng)域。支持向量機(jī)的優(yōu)點(diǎn)是魯棒性強(qiáng)，對(duì)維數(shù)不敏感，對(duì)訓(xùn)練樣本數(shù)量要求較低，而且支持向量機(jī)具有良好?的?泛?化?性?能，具?有?較?強(qiáng)?的?統(tǒng)?計(jì)?學(xué)?習(xí)?理?論基礎(chǔ)。

2分類實(shí)驗(yàn)與數(shù)據(jù)分析?

2.1?實(shí)驗(yàn)樣本與樣本光譜反射率?

本實(shí)驗(yàn)采用的樣本包括標(biāo)準(zhǔn)反射率白板、健康的臭椿樣本、感染溝眶象的臭椿樣本、健康的云南松樣本、感染云南松切梢小蠹的云南松樣本、健康的欒樹樣本、感染潔長(zhǎng)棒長(zhǎng)蠹的欒樹樣本等。植被、土壤和水體等物質(zhì)是地球表面的主要組成部分，具有顯著不同的光譜特征。色素吸收決定著可見光波段的光譜反射率，細(xì)胞結(jié)構(gòu)決定近紅外波段的光譜反射率。一般情況下，植被在350～2500ｎｍ波長(zhǎng)范圍內(nèi)具有典型的反射光譜特征，而且植被分析常使用紅邊區(qū)域和近紅外區(qū)域的工作波段。因此，本實(shí)驗(yàn)的工作波長(zhǎng)選定為650～1100ｎｍ，光譜通道數(shù)量為91個(gè)，光譜分辨率為5ｎｍ，在吸收率100％的黑色參考板背景下對(duì)上述目標(biāo)在固定距離處分別采集回波強(qiáng)度信號(hào)。實(shí)驗(yàn)采集的部分樣本如圖２所示。

圖２?高光譜激光雷達(dá)數(shù)據(jù)采集的樣本實(shí)例

由于林木樣本的生長(zhǎng)環(huán)境和保存狀況一直處于變化之中，而且木質(zhì)樣本易受大氣介質(zhì)等因素的影響，所以林木樣本表面的成分隨著時(shí)間推移而不斷變化，因此激光回波強(qiáng)度受林木樣本表面特性的影響很大。首先對(duì)健康林木樣本表面與感染林木樣本表面的光譜反射率差異進(jìn)行初步分析。圖３為健康樣本表面與染病樣本表面的反射率趨勢(shì)。從圖３可以看到，除了健康臭椿與感染溝眶象臭椿的反射率在部分波長(zhǎng)范圍內(nèi)反射率差異較大，其他幾種林木樣本健康與染病的反射率都十分接近。

圖３?健康樣本表面與染病樣本表面的反射率。（ａ）臭椿；（ｂ）云南松；（ｃ）欒樹

2.2?支持向量機(jī)分類器模型的構(gòu)建及參數(shù)選擇

由于林木樣本的波長(zhǎng)與反射率是非線性的，為此本文采用了非線性支持向量機(jī)分類器。支持向量機(jī)可以配置不同的核函數(shù)來處理各種非線性決策邊界，目前常用的核函數(shù)有多項(xiàng)式核函數(shù)、高斯徑向基核函數(shù)、線性核函數(shù)。因?yàn)楦咚箯较蚧撕瘮?shù)具有平滑、穩(wěn)定、同步的特性，所以本研究采用的是高斯徑向基核函數(shù)，即

式中：ｘ和ｚ?為空間內(nèi)不同的樣本點(diǎn)；σ?為核半徑。關(guān)于實(shí)驗(yàn)數(shù)據(jù)處理，本文使用軟件對(duì)回波信號(hào)進(jìn)行處理、劃分樣本集、建模以及評(píng)價(jià)。通過數(shù)據(jù)預(yù)處理得到了樣本的光譜信息，之后將健康的臭椿樣本與感染溝眶象的臭椿樣本、健康的云南松樣本與感染云南松切?梢小?蠹的?云南松樣本、健康的欒樹樣本與感染潔長(zhǎng)棒長(zhǎng)蠹的欒樹樣本的表面反射率數(shù)據(jù)一分為二，一半作為支持向量機(jī)的訓(xùn)練集并用于構(gòu)建支持向量機(jī)模型，另一半作為測(cè)試集并用來驗(yàn)證支持向量機(jī)?分類?器模?型?的分?類效果。得到三組不同的支持向量機(jī)模型之后，調(diào)節(jié)模型參數(shù)使模型性能更加優(yōu)越。使用基于高斯徑向基核函數(shù)的?SVM?模型有兩個(gè)非常重要的參數(shù)?Ｃ?與γ。

Ｃ?為懲罰系數(shù)，即對(duì)誤差的寬容度，若Ｃ值越大，則說明越不能容忍出現(xiàn)誤差，容易過擬合；若?Ｃ值越小，則容易欠擬合，綜上Ｃ?值過大或過小，都會(huì)影響模型的泛化能力。γ?是?選?擇?高?斯?徑?向?基?函?數(shù)作為核函數(shù)后，該函數(shù)自帶的一個(gè)參數(shù)，其?隱?含地決定了數(shù)據(jù)映射到新的特征空間后的分布，γ?值越大，則支持向量越少，這會(huì)造成?SVM?只作用于支持向量樣本附近，則訓(xùn)練集上的準(zhǔn)確率很高，而測(cè)試集上的準(zhǔn)確率較低，存在過擬合的風(fēng)險(xiǎn)；γ值越小，則支持向量越多，則會(huì)造成平滑效應(yīng)太大，無法在訓(xùn)練集上得到特別高的準(zhǔn)確率，也會(huì)影響?測(cè)試集的準(zhǔn)確率。此外高斯徑向基核函數(shù)公式里面的σ?和γ?的關(guān)系，可以通過改變?chǔ)?值來影響γ，先根據(jù)經(jīng)驗(yàn)對(duì)其懲罰因子?Ｃ?和核函數(shù)的固有參數(shù)進(jìn)行初始化賦值，然后通過實(shí)驗(yàn)測(cè)試，再根據(jù)測(cè)試數(shù)據(jù)反復(fù)調(diào)整參數(shù)值，直到得到滿意的結(jié)果為止。因?為實(shí)驗(yàn)選用正類與負(fù)類的數(shù)據(jù)?數(shù)量一樣多，所以對(duì)懲罰系數(shù)在?０～64?之?間進(jìn)行選擇。根據(jù)臺(tái)灣大學(xué)林智仁教授的?研?究表明σ?值通?？梢栽冢?5～２5之間，所以對(duì)σ?在２-5～２5之間進(jìn)行選擇。圖４為不?同σ?值條?件下健康林木樣本與染病林木樣本的分類精度趨勢(shì)，圖?５?為不同Ｃ?值條件下健康林木樣本與染病林木樣本的分類精度趨勢(shì)。

圖４不同σ?值條件下健康林木樣本與染病林木樣本的分類精度趨勢(shì)。（ａ）臭椿；（ｂ）云南松；（ｃ）欒樹

σ?的取值集中在?０～１?之間，γ?的取值較大、支持向量減少，則展現(xiàn)出較高的?分?類?精度；當(dāng)σ?增大時(shí)，γ?值減小、支持向量增多，則樣本的分類精度呈現(xiàn)出明顯的下降趨勢(shì)。從圖４（ａ）可以看到，健康臭椿與感染溝眶象臭椿在σ＝２-2.7附近的分類精度可達(dá)0.9890。從圖４（ｂ）可以看到，健康云南松和感染切梢小?蠹云南?松?在σ＝２－２．１附?近?的?分?類?精?度?可?達(dá)0.9286。從圖４（ｃ）可以看到，健康欒樹與感染潔長(zhǎng)棒長(zhǎng)蠹欒樹在σ＝２-2.5附近的分類精度可達(dá)0.6703。對(duì)于懲罰系數(shù)?Ｃ?來說，隨著?Ｃ?值的增加，分類精度迅速增加，當(dāng)前取值超過一定范圍之后分類器模型的性能將不會(huì)有較大提升。通過分析發(fā)現(xiàn)，剛開始Ｃ?值增加，模型的復(fù)雜度增加，支持向量減少；而當(dāng)Ｃ?值足夠大之后，模型中邊界支持向量的數(shù)量為０，Ｃ?的變化就不會(huì)對(duì)模型的性能產(chǎn)生較大的影響。從圖５（ａ）可以看到，健康臭椿與感染溝眶象臭椿在參數(shù)Ｃ＝47附近的分類精度可達(dá)0.9945。從圖５（ｂ）可以看到，健康云南松和感染切梢小蠹云南松在?Ｃ＝18?附近的分類精度可達(dá)0.9286。從圖?５（ｃ）可以看到健康欒樹與感染潔長(zhǎng)棒長(zhǎng)蠹欒樹在?Ｃ＝48附近的分類精度可達(dá)0.6593。

圖５?不同Ｃ?值條件下健康林木樣本與染病林木樣本的分類精度趨勢(shì)。（ａ）臭椿；（ｂ）云南松；（ｃ）欒樹

2.2?測(cè)試集數(shù)據(jù)實(shí)驗(yàn)結(jié)果驗(yàn)證

完成支持向量機(jī)模型優(yōu)化之后，就得到了適合識(shí)別不同種類林木樣本健康或者染病的支持向量機(jī)模型。將測(cè)試集數(shù)據(jù)輸?入?支?持向量機(jī)分類器模型中，用來驗(yàn)證其分類效果。表２為樣本測(cè)試數(shù)據(jù)的分類精度。

表２測(cè)試集樣本數(shù)據(jù)的分類精度

通過測(cè)試集數(shù)據(jù)的驗(yàn)證，經(jīng)過參數(shù)優(yōu)化之后的支持向量機(jī)模型可以很好地分類識(shí)別各種林木樣本目標(biāo)，而且都達(dá)到了預(yù)期效果，并未出現(xiàn)過擬合的現(xiàn)象。從表２可以看到，臭椿和云南松的分類精度較高，可?以?達(dá)?到?90％?以?上；欒?樹?的?分?類?精?度?也?接?近70％。分?類?精?度?較?高?主?要?是?因?為?其?在?650～1100ｎｍ整個(gè)波長(zhǎng)范圍內(nèi)的反射率差異顯著，除了部分波段有混疊現(xiàn)象?外，其他?部?分差異?十?分明顯。而欒樹在650～1100ｎｍ?整個(gè)波長(zhǎng)范圍內(nèi)的反射率非常接近，這也是其分類精度較低的主要原因。綜合來看支持向量機(jī)對(duì)于健康與?染病?的多?種林木樣本，分類效果明顯，原因在于新型可調(diào)諧的高光譜激光雷達(dá)的91通道光譜信息為分類提供了豐富的光譜信息，而且具有主動(dòng)測(cè)量、占地面積小、抗干擾能力強(qiáng)等優(yōu)點(diǎn)，這就使測(cè)量精度得到保證。

2結(jié)??論

本研究利用自研的新型可調(diào)諧高光譜激光雷達(dá)系統(tǒng)，對(duì)固定距離下不同種類的健康林木樣本和染病林木樣本進(jìn)行了分類研究，通過分析健康林木樣本表面光譜與橫截面光譜的關(guān)系證實(shí)表面光譜用于分類的可行性；運(yùn)用支持向量機(jī)對(duì)健康和染病的林木樣本的表面光譜反射率展開了關(guān)于核函數(shù)選取、模型參數(shù)選擇、分類精度的相關(guān)研究；將實(shí)驗(yàn)獲得的數(shù)據(jù)輸入到優(yōu)化之后的支持向量機(jī)模型中并進(jìn)行驗(yàn)證，發(fā)現(xiàn)其對(duì)實(shí)驗(yàn)樣本有良好的分類效果。其中健康臭椿和感染溝眶象臭椿在參數(shù)σ＝２-2.7、Ｃ＝47?附近的分類精度可達(dá)96.98％；健康云南松和感染切梢小蠹云南松在參數(shù)σ＝２-2.1、Ｃ＝18??附近的分類精度可達(dá)?91.21％；健康欒樹與感染潔長(zhǎng)棒長(zhǎng)蠹欒樹在參數(shù)σ＝２-2.5、Ｃ＝48?附?近的分類精度?66.21％。

雖然欒樹的分類精度沒有臭椿和云南松高，但在日后的研究中可以通過對(duì)模型的調(diào)整優(yōu)化來進(jìn)一步提升分類精度。本文在有限的光譜通道條件下對(duì)健康和染病的林木樣本表面光譜反射率進(jìn)行分類識(shí)別，探索了一種適用于基于可調(diào)諧高光譜激光雷達(dá)光譜的病蟲害林木樣本分類方法，目前只獲得了初步結(jié)果，但具有一定的參考價(jià)值。本文的研究并未關(guān)注林木樣本的染病程度，只是定性地分析了林木樣本的健康染病與否，日后可以展開林木樣本染病程度的相關(guān)研究，實(shí)現(xiàn)更加精確的監(jiān)測(cè)。

萊森光學(xué)（深圳）有限公司是一家提供光機(jī)電一體化集成解決方案的高科技公司，我們專注于光譜傳感和光電應(yīng)用系統(tǒng)的研發(fā)、生產(chǎn)和銷售。