遙感*的優勢,給高空目標識別、地物分類等研究提供了新的平臺。為了研究不同高度下地物的光譜特征變化以及不同時期、不同高度下地物的有效快速識別方法,本研究以2015年8月7日200 m高空高光譜影像和2015年9月25日100 m、200 m、300 m高空高光譜影像共4景影像為研究對象,研究不同高度下,地物的高光譜特征變化以及不同時期、不同高度下地物的分類精度。結果表明不同高度下,植被的光譜反射率差異顯著,隨著高度的升高,植被*的特征如“綠峰”降低、“紅谷”升高、“紅邊位置”出現“紅移”,在近紅外范圍內,光譜反射率降低;在綜合考慮人工參與程度、處理時間和分類精度等方面進行比較可知,基于ISODATA法可實現不同時期、不同高度下地物的快速識別研究。研究結果為利用無人機高光譜遙感在其他領域的應用奠定了基礎。
關鍵詞:無人機遙感;高光譜;不同高度;地物識別;ISODATA
無人機(unmanned aerial vehicle,UAV),是一種通過無線遙控或規劃航線飛行的無人駕駛飛機,它一般有動力系統、飛控系統、無線通訊遙控系統、有效載荷(武器、偵查設備)等部分組成[1]。無人機與遙感技術的結合,即無人機遙感,克服了近地面小范圍的作物遙感監測的極限,同時也克服了衛星影像受時間分辨率、空間分辨率的影響。無人機遙感以其全天時、實時化、高分辨率、靈活機動、高性價比等優勢,在農業、生態環境、新農村建設規劃、自然災害監測、公共安全、水利、礦產資源勘探、測繪等國民經濟及社會發展各個領域發揮著越來越重要的作用,成為繼衛星遙感和有人通用航空遙感技術之后的新興發展方向[2]。
目前,國內外學者通過無人機遙感平臺搭載各種傳感器獲取地面圖像光譜信息,取得了較大的研究進展。如Sugiura 等[3]將成像傳感器搭載在無人機上,獲取小面積農田信息,劃分出研究區域內的作物葉面積指數分布情況。李冰等[4] 利用低空無人機遙感監測冬小麥不同生長時期覆蓋度的變化。Mitch Bryson 等[5]使用可見光低空拍攝無人機影像基于植被顏色及紋理對不同植被類型進行分類。Francisco 等[6]使用無人機影像對植被黃龍病進行監測。Córcoles 等[7]利用旋翼無人機實現了洋蔥郁閉度的無損測量,并建模分析了郁閉度與葉面積指數之間的關系。劉峰等[8]設計構建基于無人直升機平臺的低空遙感系統,對板栗生育期內植被覆蓋度變化情況實施了監測。高林等[9] 以多旋翼無人機為平臺同步搭載高清數碼相機和多光譜傳感器組成的無人機農情監測系統對研究大豆葉面積指數反演,取得了較高的精度。
然而,國內外大部分學者利用無人機遙感平臺運用于作物的農情監測、礦產資源探測等領域較多[10-12],然而很少有人利用無人機平臺搭載成像高光譜傳感器研究不同高度下同一地物光譜反射率的變化以及不同高度下不同時期地物的快速識別研究。因此本研究利用無人機搭載成像高光譜儀,以期分析不同高度下地物的光譜反射率差異和不同高度不同時期地物的快速精準識別,為利用無人機高光譜遙感在其他領域的應用奠定基礎。
2.1 研究區域
野外實驗分別于2015年8月7日和9月25日在北京市通州區馬駒橋鎮某農田開展。
2.2 數據采集設備
本次試驗采用大疆八旋翼無人機S1000 (無人機凈重約4 kg,zui大載重約 6 kg),在無人機遙感平臺上搭載四川雙利合譜科技有限公司自主研發的成像光譜儀GaiaSky-mini,該無人機遙感平臺采用的是無人機懸置空中,成像光譜儀內置推掃獲取地面圖像(其主要參數見表1)。8月7日僅飛行了200 m高度,地面范圍約80 m×80 m;9月25日飛行了100 m、200 m、300 m三個不同高度,地面范圍分別約為40 m×40 m、80 m×80 m、120 m×120 m,無人機遙感平臺的實景圖,如圖1。
表1 GaiaSky-mini 機載成像高光譜儀系統參數
Table 1 GaiaSky-mini high airborne imaging spectrometer system parameters
序號 | 項目 | 參數 |
1 | 光譜掃描范圍/nm | 400~1000 |
2 | 光譜分辨率/nm | 4±0.5 |
3 | 成像鏡頭/mm | 23 |
4 | 光譜通道數 | 520 |
5 | 全幅像素 | 1392×1040 |
6 | 傳感器 | CCD Sony ICX285 |
7 | 重量 | 約1.0kg |
圖1 無人機遙感平臺實景圖
Fig. 1 UAV remote sensing platform picture
2.3 數據的預處理與分析
無人機高光譜圖像的預處理在四川雙利合譜科技有限公司自主研發的SpecView軟件中進行,包括鏡像變換[13]、黑白幀校準[14-15](如公式1所示)。高光譜圖像噪聲去除采用zui小噪聲分離變換(Minimum Noise Fraction Rotation,MNF Rotation)[16-18],圖像的感興趣區域提取、圖像分類方法均在ENVI/IDL 5.3上完成。
(1)
式中,Rref 是校正過的圖像,DNraw 是原始圖像,DNwhite為白板校正圖像,DNdark 是黑板校正圖像。
高光譜遙感圖像分類是高光譜遙感圖像分析和應用的重要內容,根據是否使用類別的先驗知識,可分為監督分類和非監督分類,常用的監督分類方法有zui小距離分類法、馬氏距離分類法、zui大似然分類法、光譜匹配法和神經網絡法等,非監督分類常用的方法有K-Means和ISODATA法[19]。本研究運用監督分類方法中的馬氏距離分類法、zui大似然分類法、神經網絡法以及非監督分類的ISODATA法,比較這四種分類方法在不同高度下不同時期的無人機高光譜影像中地物的分類識別精度,本研究利用總體分類精度和Kappa系數對分類精度進行評價[20]。
3.1 同一研究區域不同高度下光譜曲線的差異分析
圖2為9月25日,天氣為晴,基于無人機獲取的三個不同高度的高光譜圖像的(RGB真彩色合成圖),光譜曲線分別為同一研究區域,100 m、200 m、300 m三個高度的光譜曲線差異。以植被為研究對象,分析同一研究區域不同高度下光譜曲線的差異。以一個具體的植被象元為例,從圖2中可知,不同高度下植被光譜反射率差異顯著,隨著高度的升高,植被*的特征“綠峰”降低、“紅谷”升高、“紅邊位置”出現“紅移”。且在400-510 nm范圍內光譜反射率隨無人機飛行高度的升高而升高;在近紅外范圍內,光譜反射率隨無人機飛行高度的升高而降低,該這是因為隨著高度升高,空間分辨率降低,圖像中一個像元代表地面的范圍加大,出現混合像元的概率增加。從圖2中可知,隨著無人機高度的升高,混合像元中土壤所占比例增加,從而造成“綠峰”降低、“紅谷”升高,紅邊位置發生了“紅移”等現象。
圖2 同一研究區域不同高度下光譜曲線
Fig. 2 The spectral curve at different heights in the same study area
3.2 機載影像中不同地物的光譜曲線差異分析
圖3為8月7日天氣為晴,利用無人機獲取的近地面200 m高度的高光譜圖像及不同作物與土壤的光譜反射率曲線。從圖中可知,土壤的光譜曲線與三種作物的光譜曲線差異顯著,在400-515 nm和580 - 690 nm范圍內,土壤的光譜反射率高于三種作物的光譜反射率。在可見光區域,茄子苗、西瓜苗、葡萄樹三種作物的光譜反射率差異也較為顯著,其中西瓜苗的光譜反射率z高,其次是茄子苗,然后是葡萄樹;在近紅外區域,西瓜苗和茄子苗的光譜反射率十分相近,總體而言,茄子苗的光譜反射率zui高,其次是西瓜苗,然后是葡萄樹。
圖3 8月7日不同作物與土壤的光譜反射率曲線
Fig. 3 The spectral reflectance curves of different crops and soil on August 7
圖4為9月25日天氣為晴,利用無人機獲取的近地面300 m高度的高光譜圖像及不同地物的光譜反射率曲線。從中可知,植物、土壤、公路、大棚植物等地物光譜曲線差異顯著。在400-920 nm范圍內公路的光譜反射率曲線保持較平穩的趨勢,土壤的光譜反射率在此范圍內呈上升趨勢;400-700 nm范圍內,大棚植物的光譜反射率高于裸露植物,而在700-100 nm范圍內,則是裸露植物的光譜反射率高于大棚植物。
圖4 9月25日不同地物的光譜反射率曲線
Fig. 4 The spectral reflectance curves of different ground objects on September 25
3.3 不同地物的快速識別分類研究
分別運用馬氏距離法、zui大似然法、神經網絡、ISODATA(K-Means分類效果不佳,在此不列舉)四種分類方法對8月7日的200 m高光譜圖像和9月25日的300 m高光譜圖像進行分類研究。完成分類之后,使用目視判讀的方法從圖像上選取不同地物各1000像元計算混淆矩陣進行驗證。結合驗證數據,采用混淆矩陣的方法對 4 類方法進行精度評價,結果如表 2和表3 所示。從表2和表3可知,神經網絡法的分類精度z高,其次是馬氏距離分類法,再次是ISODATA法,zui低的zui大似然分類法。但從處理時間來看,神經網絡的處理時間較長,分別為35.14 min和38.26 min,其次是zui大似然分類法,分別為15.61 min和16.75min,再次是馬氏距離法,分別是12.54 min和14.55 min,處理時間z短的是ISODATA法,處理時間分別為3.64 min和3.97 min。
表2 8月7日的200 m高光譜圖像分類精度評價
Table 2 200 m hyperspectral image classification accuracy evaluation on August 7
分類方法 | 總體分類精度 | Kappa系數 | 處理時間/min |
馬氏距離 | 92.67% | 0.915 | 12.54 |
大似然法 | 86.65% | 0.844 | 15.61 |
神經網絡法 | 98.56% | 0.982 | 35.14 |
ISODAIA | 89.41% | 0.887 | 3.64 |
表3 9月25日的300 m高光譜圖像分類精度評價
Table 3 300 m hyperspectral image classification accuracy evaluation on September 25
分類方法 | 總體分類精度 | Kappa系數 | 處理時間/min |
馬氏距離 | 94.07% | 0.922 | 14.55 |
zui大似然法 | 89.45% | 0.863 | 16.75 |
神經網絡法 | 98.76% | 0.984 | 38.26 |
ISODAIA | 90.53% | 0.891 | 3.97 |
4種分類方法的原理各不同,從人工參與程度、處理效率和分類精度等方面進行比較分析可知,ISODATA法的分類精度雖然低于神經網絡和馬氏距離法,但是不需要太多的人工參與,處理卻大大提高,且分類精度高于zui大似然分類法。因此,它同時具備了監督分類的高精度和非監督分類的低時間的特點,在無人機低空遙感研究中,能夠以較低的工作成本快速獲取高精度的地物識別。圖5和圖6分別為8月7日200 m和9月25日300 m高光譜圖像基于ISODATA法的分類效果圖。從圖5可知,少量的雜草被錯分為葡萄樹、西瓜苗、茄子苗,但是從大的區域來看,原味葡萄樹、西瓜苗、茄子苗的像素并未存在錯分和漏分現象。在圖6中,少量的土壤,可能是沙土,被錯分為大棚或公路斑馬線,其他地物不存在錯分和漏分現象。總體而言,ISODATA法可實現無人機高光譜遙感對地面地物的快速有效的識別。
圖5 8月7日200 m高光譜圖像基于ISODATA法的分類效果圖
Fig. 5 The classification rendering of 200m hyperspectral image based on ISODATA method on August 7
圖6 9月25日300 m高光譜圖像基于ISODATA法的分類效果圖
Fig. 6 The classification rendering of 300m hyperspectral image based on ISODATA method on September 25
4. 結論
本研究基于無人機平臺搭載成像高光譜相機,分別獲取了兩個時間段不同高度的高光譜圖像,分析了不同高度同一地物光譜曲線的差異和基于不同高度下不同分類方法的分類精度評價。研究表明隨著無人機高度的升高,混合像元中土壤所占比例增加,從而造成“綠峰”降低、“紅谷”升高,紅邊位置發生了“紅移”等現象,研究結果與地面光譜儀ASD的研究結果一致[21-22]。ISODATA法的分類精度雖然不及神經網絡和馬氏距離分類法,但是在綜合考慮人工參與程度、處理效率和分類精度時可知,ISODATA法對于快速實現地物的識別效果更好,研究結果為利用無人機高光譜遙感進行軍事偽裝快速識別、海洋油污面積估測、災害損失評估等領域具有重大意義。