0  引言

成像技術和光譜技術是傳統的光學技術的兩個重要方向，成像技術能夠獲得物體的影像，得到其空間信息；光譜技術能夠得到物體的光學信息，進而研究其物質屬性。20世紀70年代以前，成像技術和光譜技術是相互獨立的學科，隨著遙感技術的發展，成像光譜技術迅速發展起來，它是一種快速、無損的檢測技術，具有光譜分辨率高、多波段和圖譜合一的特點，能在大尺度范圍內識別地表并深入研究其地表物質的成分及結構。目前成像光譜技術已經成為遙感技術的發展趨勢之一，并在軍事偵察、海洋遙感、地質勘探、植被分析等領域得到越來越廣泛的應用。

隨著無人機技術的日益成熟，基于無人機平臺的新型遙感技術異軍突起，得到遙感工作人員的青睞，科研工作者更對其在行業上的應用前景予以眾望。目前無人機搭載成像光譜儀在農業上可用于診斷作物長勢、病蟲害情況以及土壤肥力等；在環境保護方面可用于檢測海洋、湖泊的化工原料污染、富營養化等；在林業領域可用于檢測林火、林業病害、林木存活率、林木種類區分等；在礦產資源領域上可用于石油、礦物礦產等領域的勘測；在考古領域可用于古村落的修復、文物遺址的勘查等；在通信部門可用于電纜絕緣子的勘察等。

我國是煙草種植大國，煙草在全國各省均有種植，其面積和產量目前均居世界*，稅負收入在我國稅收總量占有舉足輕重的地位，為我國經濟發展做出了相當貢獻。而我國煙葉總產量中烤煙占了80％以上，是生產卷煙的主要原料。長期以來，我國在觀察監測煙田煙草的生長長勢以及煙草品質等方面均采用耗時、耗力、成本高的傳統方法。無人機成像高光譜技術是高效的煙草種植管理手段，其能大范圍快速、準確地為決策者提供煙田煙草的生長狀況，為決策者施肥、灌溉、噴灑農藥提供有效的數據，有利于顯著提高農藥現代化水平，促進現代農業可持續發展。因此無人機高光譜技術在煙田監測方面具有較大的應用價值，為煙田管理的新趨勢。

1  材料與試驗部分

1.1  研究區域

野外試驗田在云南大理、楚雄、石林板橋鎮紅塔煙草(集團)有限責任公司的煙草試驗基地。

1.2  數據采集設備

本次試驗采用大疆六旋翼無人機M600 Pro(無人機凈重約4 kg，大載重約 10 kg)，在無人機遙感平臺上搭載四川雙利合譜科技有限公司自主研發的高光譜成像光譜儀GaiaSky-mini-2，該無人機遙感平臺采用的是無人機懸置空中，高光譜成像光譜儀采用內置推掃的方式獲取地面圖像(其主要參數見表1)。在云南大理和楚雄煙草試驗基地，無人機的飛行高度均為120m，分別采集了33景和67景圖像，每景圖像代表地面幅寬為48m*48m；在石林板橋鎮煙草試驗基地，無人機的飛行高度為250m，一共采集了86景圖像，每景圖像代表地面幅寬為100m*100m。無人機遙感平臺的實景圖，如圖1。

表1   GaiaSky-mini 機載成像高光譜儀系統參數

Table 1 GaiaSky-mini high airborne imaging spectrometer system parameters

圖1   無人機高光譜成像系統實景圖

1.3  無人機高光譜數據的預處理與分析

無人機高光譜圖像的預處理在四川雙利合譜科技有限公司自主研發的SpecView軟件中進行，包括鏡像變換、黑白幀校準（如公式1所示）。

         (1)

式中，Rref 是黑白校正過的圖像的反射率值，DNraw 是原始圖像的DN值，DNwhite為白板的白幀數據，DNdark 是相機的系統誤差DN值。

考慮到無人機飛到一定高度后，高光譜成像儀獲取的高光譜影像數據可能會受到大氣、水汽等因素的影響。為了消除這些因素的影響，我們在無人機起飛之前，在拍攝區域放置一塊經過國家計量院標定過的2m*2m灰布，在高光譜影像獲取的時候，只需要在其中的一景高光譜影像中覆蓋到灰布即可。消除大氣、水汽等因素影響的方法如公式2所示。

    (2)

式中，Rfixed 是消除大氣、水汽等因素后的圖像光譜反射率，Rref是經過黑白校正后的圖像反射率，Rstandard是經過國家計量院標定的灰布的光譜反射率，Rgrayref 是經過黑白校正后圖像中灰布的光譜反射率。

1.4  無人機高光譜影像拼接

無人機高光譜影像的拼接采用四川雙利合譜科技有限公司自主研發的無人機高光譜拼接軟件HiSpectralStitcher進行拼接，該拼接軟件有圖像篩選、拼接預覽、投影方式選擇、拼接方法選擇、重采樣方法選擇、是否勻色、拼接結果格式輸出選擇等功能。拼接軟件界面如圖2所示。

圖2  無人機高光譜拼接軟件HiSpectralStitcher

2  結果與分析

2.1  拼接結果預覽

圖3為利用無人機高光譜拼接軟件HiSpectralStitcher對云南大理、楚雄、石林板橋鎮紅塔煙草(集團)有限責任公司的煙草試驗基地的無人機高光譜影像的三波段拼接效果預覽圖（RGB分別代表750 nm/650 nm/550 nm鄰近波長）。從拼接結果來看，除楚雄煙草試驗基地，由于航線規劃過程沒規劃成功，導致部分煙草田塊沒有拼接成功外，大理和石橋板鎮煙草基地的無人機拼接結果較好。

圖3.1 大理煙草試驗基地無人機高光譜拼接預覽圖

圖3.2 楚雄煙草試驗基地無人機高光譜拼接預覽圖

圖3.3 石橋板鎮煙草試驗基地無人機高光譜拼接預覽圖

2.2  無人機高光譜影像煙草提取方法

無人機高光譜影像不僅包含有煙草，還有其他的作物、雜草、土壤、道路、建筑物等，從影像中提取感興趣的目標物，通常的方法有非監督分類、監督分類、決策樹、面向對象分類等，本研究采用監督分類的大似然法對拼接好的無人機高光譜影像進行分類，分別提取了煙草和其他植物做植被指數分析，葉綠素、總氮、煙堿等生化指標的反演。其中大理和楚雄煙草試驗基地只提取煙草做生化指標反演，石板橋鎮煙草基地將煙草和其他作物一起參與生化指標的反演，觀察煙草與其他作物在生化指標的差異。

2.3  煙草與其他地物的光譜分析

   圖4為無人機高光譜影像上煙草、雜草、作物以及土壤在400-1000 nm范圍的光譜曲線。從圖中可知，煙草的光譜反射率高于雜草的光譜反射率，而雜草的光譜反射率則高于作物的光譜反射率，土壤的光譜反射率在400-760nm范圍內保持穩定上升的趨勢，在760-900nm范圍內光譜差異不大，在900-1000nm范圍內有一峰谷，這是因為水汽吸收的原因。

圖4 煙草與其他地物的光譜反射率曲線

2.4  植被指數分析

植被指數是一類具有一定生化意義的不同波段光譜值的組合，通常有比值植被指數、線性組合植被指數、修正植被指數、差值植被指數等。不同波段組合的植被指數對于不同指標預測效果不同。在農業上，基于光譜技術檢測作物生理指數的波段范圍一般在 400～2 500 nm 之間，涉及到色素（葉綠素、類胡蘿卜素等）、氮、水分等吸收和葉片細胞的內部結構。在 400-740 nm  可見光波段，葉綠素在 480、650、670-680、740 nm  處有吸收峰，類胡蘿卜素在420、425、440、450、470、480 nm 均有吸收峰，葉黃素在  425、445、475 nm  有吸收峰。而在  740-1 300 nm  近紅外波段由于健康的葉肉細胞反射作用，其反射率急劇升高；作物水分的吸收峰主要集中在 970、1 450、1 944 nm  處。因此當作物受到脅迫作用時，相應的氮、色素、酶等發生變化，通過應用各種植被指數監測這些生理指標變化，可判斷作物脅迫情況、生長狀況以及產量情況。圖5為大理、楚雄、石板橋鎮煙草基地煙草的NDVI分布圖（NDVI采用的綠光波長與紅光波長的），NDVI值越大則說明長勢越好。

圖5.1 大理煙草基地煙草NDVI分布圖

圖5.2 楚雄煙草基地煙草NDVI分布圖

圖5.3 石板橋鎮煙草基地煙草、作物、雜草的NDVI分布圖

2.5  煙草葉綠素分析

葉綠素含量是植被生長狀態的良好指示器，它與植被的光合作用能力，發育狀況，脅迫

程度及氮素含量等密切相關，又是植被與外界發生物質能量交換的重要條件，因此估測作物葉綠素含量已成為評價其長勢的一種有效手段。王強等以棉花冠層葉綠素密度及冠層高光譜反射率為數據源，分析了葉綠素密度同原始光譜反射率(R)，一階導數光譜反射率(DR)，已有光譜指數及全波段組合指數的相關性，認為由一階導數比值為自變量估測棉花冠層葉綠素密度的模型效果好。丁永軍等和宋開山等分別利用高光譜建立了預測番茄、玉米等葉綠素含量的模型，效果較好。徐新剛等研究了高光譜曲線特征同葉綠素的相關性，發現可見光波段520-740nm之間的歸一化光譜反射率同葉綠素含量呈良好的負相關關系；綠峰反射光譜曲線特征邊的變化速率及綠峰兩邊的夾角等變量與水稻葉片葉綠素含量具有良好的相關性。楊巧明等采用高光譜分析和線性回歸技術分析了3個橡膠樹品種葉片的室內反射光譜數據和葉綠素含量，認為橡膠樹葉片的葉綠素含量敏感波段主要指藍邊范圍。

針對煙草葉綠素的監測研究，本研究利用云南農業大學研究的模型進行反演，圖6分別對大理、楚雄和石板橋鎮煙草基地的煙草進行葉綠素a和葉綠素b進行反演。

圖6.1 大理煙草基地煙草葉綠素a和葉綠素b的反演圖

圖6.2 楚雄煙草基地煙草葉綠素a和葉綠素b的反演圖

圖6.3 石橋板鎮煙草基地煙草、作物、雜草葉綠素a和葉綠素b的反演圖

2.6  煙草總氮分析

國內外對多種作物的營養元素、多種生化組分反射光譜特性及其與品質和產量的關系研究有較多的研究和報道。李佛琳等對可見近紅外波段(350nm-1650nm)單波段光譜和100 種光譜指數共兩類光譜參量進行了與生化組分之間線性函數、冪函數、指數函數共3種形式相關分析和基于決定系數的篩選，結果表明，對于煙堿、總氮、總糖方法分別是在 1135nm反射率倒數對數二階微分的線性擬合(R2=0.20)、在 666nm 反射率倒數對數一階微分的冪函數擬合(R2=0.44)、在532nm反射率一階微分線性擬合(R2=0.54)。牛錚等以小麥活體葉片為研究對象，利用高光譜數據，探索了利用成像光譜遙感預測小麥生化組分的基礎性研究，采用多元逐步回歸方法，分析了小麥葉片7種生化組分含與特征光譜參量的關系。張喜杰等在溫室內，利用自然光照反射光譜研究了黃瓜葉片的含氮量預測模型，發現利用單一敏感波長520 nm就可獲得理想模型；原始光譜及其一階微分光譜都可用于黃瓜葉片含氮量預測。Thomas 等對 7 種不同的作物葉片在不同氮素營養水平下的光譜特性，結果證明在氮素營養缺乏條件下，7 種作物葉片的可見光波段的反射率表現為增加，但不同植物的光譜反射率增加幅度不一樣。蘇永士等研究了不同施肥條件下煙草冠層光譜反射率差異顯著，經篩選綠波段植被指數(GNDVI)與葉綠素含量和葉面積指數關系顯著，建立的回歸方程的 R2 分別為 0.436 和 0.568，均達到顯著水平。

本研究利用現有的煙草總氮非線性模型，通過算法分析，反演大理、楚雄、石板橋鎮煙草基地煙草的總氮含量，如圖7所示。

圖7.1 大理煙草基地煙草總氮反演圖

圖7.2 楚雄煙草基地煙草總氮反演圖

圖7.3 石板橋鎮煙草基地煙草、作物、雜草總氮反演圖

2.7  煙草煙堿分析

煙草是我國一項非常重要的經濟作物，中國是世界上烤煙種植面積和產量大的國家，煙草除了可為卷煙制品提供原料外，同時在食品和醫藥領域還有很多潛在的應用價值。煙堿是絕大多數煙草屬作物葉片的重要的生物堿，其含量是評價煙葉品質的重要指標。雖然煙堿對人類健康和環境具有危害性，但是美國食品和藥物管理局近年來證實了它可用來治療一些疾病，也可被用來制作戒煙產品和殺蟲劑，同時，煙堿和煙草中的氮及色素含量密切相關。因此，快速地監測遮蔭條件下煙草葉片中煙堿含量的變化有助于采取有效的農藝管理措施。李向陽等以不同成熟度烤煙葉片為研究對象，開展了室內光譜紅邊參數研究，研究表明紅邊位置是預測烤煙葉片成熟度精度較高的指標，紅邊位置位于 693-695nnl。劉國順等利用冠層光譜估測煙草葉面積指數和地上生物量，篩選出相應的特征變量為 Rg／Dr。蔣錦鋒等以煙草為研究對象，應用近紅外光譜技術建立了煙草主要化學成分的快速、無損檢測方法。

本文利用研究者利用光譜設備研究煙草觀察煙堿的模型反演到大尺度無人機高光譜影像上，從而分析得到每一田塊煙草煙堿的含量分布情況，如圖8所示。

圖8.1 大理煙草基地煙草煙堿含量反演圖

圖8.2 楚雄煙草基地煙草煙堿含量反演圖

圖8.3 石板橋鎮煙草基地煙草、作物、雜草煙堿含量反演圖

3  結論

從煙草的NDVI和葉綠素a/b、總氮和煙堿反演分布圖來看，石橋板鎮煙草基地的煙草與其他作物和雜草的含量差異顯著，特別是煙堿作為煙草*的屬性，其與其他作物、雜草差異尤為顯著，因此利用無人機高光譜可大范圍快速識別煙草，可幫助決策者判斷煙草的種植面積。同時根據大理、楚雄和石板橋鎮煙草基地的煙草NDVI、和葉綠素a/b、總氮和煙堿反演分布圖來看，利用無人機高光譜技術可以判別煙草的長勢及監測生化參數的含量，為決策者定量施肥、灌溉、噴灑農藥等提供技術支持，同時根據煙草煙堿的分布圖，可以判斷煙草的品質。

因此利用高光譜技術可以對煙草主要生化成分進行很好的監測，而且具有廣闊的應用前景。然而煙草學者們雖然已經提出的一系列的煙草生化成分的監測和估算模型，但每個模型都有特定的研究方法和適用條件，并不適用于所有品種的煙草，所以很難找到通用的模型。同時高光譜遙感數據有一些缺陷，比如景觀異質性、大氣噪音、太陽位置等的干擾都會影響高光譜遙感技術在實際的應用能力。所以要使煙草光譜資料更加完備，今后的研究方向：一是煙草生化參數，要一步步地上升到較為確的模型水平 ，進一步對其進行細化和確定。為了降低實際生產中煙草的品種類型 、生態條件多樣和栽培管理等方面的差異，需要建立更加全面和更具規模的樣本參數，以便進行修正。二是遙感數據的分析方法和遙感信息的信噪比要進一步提高，在今后的研究中，要注重于完善和擴充煙草的光譜數據庫，加強煙草高光譜數據的采集和挖掘，結合3S技術的應用 ，進一步推動高光譜遙感技術在煙草生化成分診斷和監測中的應用。