基于無(wú)人機(jī)高光譜不同高度的地物快速識(shí)別研究

基于無(wú)人機(jī)高光譜不同高度的地物快速識(shí)別研究

摘要：無(wú)人機(jī)遙感以其操作簡(jiǎn)便、高時(shí)效、高分辨率、低成本、低風(fēng)險(xiǎn)、低損耗等傳統(tǒng)

遙感無(wú)法比擬的優(yōu)勢(shì)，給高空目標(biāo)識(shí)別、地物分類等研究提供了新的平臺(tái)。為了研究不同高度下地物的光譜特征變化以及不同時(shí)期、不同高度下地物的有效快速識(shí)別方法，本研究以2015年8月7日200 m高空高光譜影像和2015年9月25日100 m、200 m、300 m高空高光譜影像共4景影像為研究對(duì)象，研究不同高度下，地物的高光譜特征變化以及不同時(shí)期、不同高度下地物的分類精度。結(jié)果表明不同高度下，植被的光譜反射率差異顯著，隨著高度的升高，植被特有的特征如“綠峰”降低、“紅谷”升高、“紅邊位置”出現(xiàn)“紅移”，在近紅外范圍內(nèi)，光譜反射率降低；在綜合考慮人工參與程度、處理時(shí)間和分類精度等方面進(jìn)行比較可知，基于ISODATA法可實(shí)現(xiàn)不同時(shí)期、不同高度下地物的快速識(shí)別研究。研究結(jié)果為利用無(wú)人機(jī)高光譜遙感在其他領(lǐng)域的應(yīng)用奠定了基礎(chǔ)。

關(guān)鍵詞：無(wú)人機(jī)遙感；高光譜；不同高度；地物識(shí)別；ISODATA

無(wú)人機(jī)(unmanned aerial vehicle，UAV），是一種通過(guò)無(wú)線遙控或規(guī)劃航線飛行的無(wú)人駕駛飛機(jī)，它一般有動(dòng)力系統(tǒng)、飛控系統(tǒng)、無(wú)線通訊遙控系統(tǒng)、有效載荷（武器、偵查設(shè)備）等部分組成[1]。無(wú)人機(jī)與遙感技術(shù)的結(jié)合，即無(wú)人機(jī)遙感，克服了近地面小范圍的作物遙感監(jiān)測(cè)的極限，同時(shí)也克服了衛(wèi)星影像受時(shí)間分辨率、空間分辨率的影響。無(wú)人機(jī)遙感以其全天時(shí)、實(shí)時(shí)化、高分辨率、靈活機(jī)動(dòng)、高性價(jià)比等優(yōu)勢(shì)，在農(nóng)業(yè)、生態(tài)環(huán)境、新農(nóng)村建設(shè)規(guī)劃、自然災(zāi)害監(jiān)測(cè)、公共安全、水利、礦產(chǎn)資源勘探、測(cè)繪等國(guó)民經(jīng)濟(jì)及社會(huì)發(fā)展各個(gè)領(lǐng)域發(fā)揮著越來(lái)越重要的作用，成為繼衛(wèi)星遙感和有人通用航空遙感技術(shù)之后的新興發(fā)展方向[2]。

目前，國(guó)內(nèi)外學(xué)者通過(guò)無(wú)人機(jī)遙感平臺(tái)搭載各種傳感器獲取地面圖像光譜信息，取得了較大的研究進(jìn)展。如Sugiura 等[3]將成像傳感器搭載在無(wú)人機(jī)上，獲取小面積農(nóng)田信息，劃分出研究區(qū)域內(nèi)的作物葉面積指數(shù)分布情況。李冰等[4] 利用低空無(wú)人機(jī)遙感監(jiān)測(cè)冬小麥不同生長(zhǎng)時(shí)期覆蓋度的變化。Mitch Bryson 等[5]使用可見(jiàn)光低空拍攝無(wú)人機(jī)影像基于植被顏色及紋理對(duì)不同植被類型進(jìn)行分類。Francisco 等[6]使用無(wú)人機(jī)影像對(duì)植被黃龍病進(jìn)行監(jiān)測(cè)。Córcoles 等[7]利用旋翼無(wú)人機(jī)實(shí)現(xiàn)了洋蔥郁閉度的無(wú)損測(cè)量，并建模分析了郁閉度與葉面積指數(shù)之間的關(guān)系。劉峰等[8]設(shè)計(jì)構(gòu)建基于無(wú)人直升機(jī)平臺(tái)的低空遙感系統(tǒng)，對(duì)板栗生育期內(nèi)植被覆蓋度變化情況實(shí)施了監(jiān)測(cè)。高林等[9] 以多旋翼無(wú)人機(jī)為平臺(tái)同步搭載高清數(shù)碼相機(jī)和多光譜傳感器組成的無(wú)人機(jī)農(nóng)情監(jiān)測(cè)系統(tǒng)對(duì)研究大豆葉面積指數(shù)反演，取得了較高的精度。

然而，國(guó)內(nèi)外大部分學(xué)者利用無(wú)人機(jī)遙感平臺(tái)運(yùn)用于作物的農(nóng)情監(jiān)測(cè)、礦產(chǎn)資源探測(cè)等領(lǐng)域較多[10-12]，然而很少有人利用無(wú)人機(jī)平臺(tái)搭載成像高光譜傳感器研究不同高度下同一地物光譜反射率的變化以及不同高度下不同時(shí)期地物的快速識(shí)別研究。因此本研究利用無(wú)人機(jī)搭載成像高光譜儀，以期分析不同高度下地物的光譜反射率差異和不同高度不同時(shí)期地物的快速精準(zhǔn)識(shí)別，為利用無(wú)人機(jī)高光譜遙感在其他領(lǐng)域的應(yīng)用奠定基礎(chǔ)。

2.1  研究區(qū)域

野外實(shí)驗(yàn)分別于2015年8月7日和9月25日在北京市通州區(qū)馬駒橋鎮(zhèn)某農(nóng)田開(kāi)展。

2.2  數(shù)據(jù)采集設(shè)備

本次試驗(yàn)采用大疆八旋翼無(wú)人機(jī)S1000 (無(wú)人機(jī)凈重約4 kg，最大載重約 6 kg)，在無(wú)人機(jī)遙感平臺(tái)上搭載江蘇雙利合譜科技有限公司自主研發(fā)的成像光譜儀GaiaSky-mini，該無(wú)人機(jī)遙感平臺(tái)采用的是無(wú)人機(jī)懸置空中，成像光譜儀內(nèi)置推掃獲取地面圖像(其主要參數(shù)見(jiàn)表1)。8月7日僅飛行了200 m高度，地面范圍約80 m×80 m；9月25日飛行了100 m、200 m、300 m三個(gè)不同高度，地面范圍分別約為40 m×40 m、80 m×80 m、120 m×120 m，無(wú)人機(jī)遙感平臺(tái)的實(shí)景圖，如圖1。

表1   GaiaSky-mini 機(jī)載成像高光譜儀系統(tǒng)參數(shù)

Table 1 GaiaSky-mini high airborne imaging spectrometer system parameters

圖1   無(wú)人機(jī)遙感平臺(tái)實(shí)景圖

Fig. 1 UAV remote sensing platform picture

2.3  數(shù)據(jù)的預(yù)處理與分析

無(wú)人機(jī)高光譜圖像的預(yù)處理在江蘇雙利合譜科技有限公司自主研發(fā)的SpecView軟件中進(jìn)行，包括鏡像變換[13]、黑白幀校準(zhǔn)[14-15]（如公式1所示）。高光譜圖像噪聲去除采用最小噪聲分離變換（Minimum Noise Fraction Rotation，MNF Rotation）[16-18]，圖像的感興趣區(qū)域提取、圖像分類方法均在ENVI/IDL 5.3上完成。

   (1)

式中，Rref 是校正過(guò)的圖像，DNraw 是原始圖像，DNwhite為白板校正圖像，DNdark 是黑板校正圖像。

高光譜遙感圖像分類是高光譜遙感圖像分析和應(yīng)用的重要內(nèi)容，根據(jù)是否使用類別的先驗(yàn)知識(shí)，可分為監(jiān)督分類和非監(jiān)督分類，常用的監(jiān)督分類方法有最小距離分類法、馬氏距離分類法、最大似然分類法、光譜匹配法和神經(jīng)網(wǎng)絡(luò)法等，非監(jiān)督分類常用的方法有K-Means和ISODATA法[19]。本研究運(yùn)用監(jiān)督分類方法中的馬氏距離分類法、最大似然分類法、神經(jīng)網(wǎng)絡(luò)法以及非監(jiān)督分類的ISODATA法，比較這四種分類方法在不同高度下不同時(shí)期的無(wú)人機(jī)高光譜影像中地物的分類識(shí)別精度，本研究利用總體分類精度和Kappa系數(shù)對(duì)分類精度進(jìn)行評(píng)價(jià)[20]。

3.1  同一研究區(qū)域不同高度下光譜曲線的差異分析

    圖2為9月25日，天氣為晴，基于無(wú)人機(jī)獲取的三個(gè)不同高度的高光譜圖像的（RGB真彩色合成圖），光譜曲線分別為同一研究區(qū)域，100 m、200 m、300 m三個(gè)高度的光譜曲線差異。以植被為研究對(duì)象，分析同一研究區(qū)域不同高度下光譜曲線的差異。以一個(gè)具體的植被象元為例，從圖2中可知，不同高度下植被光譜反射率差異顯著，隨著高度的升高，植被特有的特征“綠峰”降低、“紅谷”升高、“紅邊位置”出現(xiàn)“紅移”。且在400-510 nm范圍內(nèi)光譜反射率隨無(wú)人機(jī)飛行高度的升高而升高；在近紅外范圍內(nèi)，光譜反射率隨無(wú)人機(jī)飛行高度的升高而降低，該這是因?yàn)殡S著高度升高，空間分辨率降低，圖像中一個(gè)像元代表地面的范圍加大，出現(xiàn)混合像元的概率增加。從圖2中可知，隨著無(wú)人機(jī)高度的升高，混合像元中土壤所占比例增加，從而造成“綠峰”降低、“紅谷”升高，紅邊位置發(fā)生了“紅移”等現(xiàn)象。

圖2 同一研究區(qū)域不同高度下光譜曲線

         Fig. 2 The spectral curve at different heights in the same study area

3.2  機(jī)載影像中不同地物的光譜曲線差異分析

  圖3為8月7日天氣為晴，利用無(wú)人機(jī)獲取的近地面200 m高度的高光譜圖像及不同作物與土壤的光譜反射率曲線。從圖中可知，土壤的光譜曲線與三種作物的光譜曲線差異顯著，在400-515 nm和580 - 690 nm范圍內(nèi)，土壤的光譜反射率高于三種作物的光譜反射率。在可見(jiàn)光區(qū)域，茄子苗、西瓜苗、葡萄樹(shù)三種作物的光譜反射率差異也較為顯著，其中西瓜苗的光譜反射率最高，其次是茄子苗，最后是葡萄樹(shù)；在近紅外區(qū)域，西瓜苗和茄子苗的光譜反射率十分相近，總體而言，茄子苗的光譜反射率最高，其次是西瓜苗，最后是葡萄樹(shù)。

圖3  8月7日不同作物與土壤的光譜反射率曲線

Fig. 3 The spectral reflectance curves of different crops and soil on August 7

    圖4為9月25日天氣為晴，利用無(wú)人機(jī)獲取的近地面300 m高度的高光譜圖像及不同地物的光譜反射率曲線。從中可知，植物、土壤、公路、大棚植物等地物光譜曲線差異顯著。在400-920 nm范圍內(nèi)公路的光譜反射率曲線保持較平穩(wěn)的趨勢(shì)，土壤的光譜反射率在此范圍內(nèi)呈上升趨勢(shì)；400-700 nm范圍內(nèi)，大棚植物的光譜反射率高于裸露植物，而在700-100 nm范圍內(nèi)，則是裸露植物的光譜反射率高于大棚植物。

圖4  9月25日不同地物的光譜反射率曲線

Fig. 4 The spectral reflectance curves of different ground objects on September 25

3.3  不同地物的快速識(shí)別分類研究

分別運(yùn)用馬氏距離法、最大似然法、神經(jīng)網(wǎng)絡(luò)、ISODATA（K-Means分類效果不佳，在此不列舉）四種分類方法對(duì)8月7日的200 m高光譜圖像和9月25日的300 m高光譜圖像進(jìn)行分類研究。完成分類之后，使用目視判讀的方法從圖像上選取不同地物各1000像元計(jì)算混淆矩陣進(jìn)行驗(yàn)證。結(jié)合驗(yàn)證數(shù)據(jù)，采用混淆矩陣的方法對(duì) 4 類方法進(jìn)行精度評(píng)價(jià)，結(jié)果如表 2和表3 所示。從表2和表3可知，神經(jīng)網(wǎng)絡(luò)法的分類精度最高，其次是馬氏距離分類法，再次是ISODATA法，最低的最大似然分類法。但從處理時(shí)間來(lái)看，神經(jīng)網(wǎng)絡(luò)的處理時(shí)間較長(zhǎng)，分別為35.14 min和38.26 min，其次是最大似然分類法，分別為15.61 min和16.75min，再次是馬氏距離法，分別是12.54 min和14.55 min，處理時(shí)間最短的是ISODATA法，處理時(shí)間分別為3.64 min和3.97 min。

表2  8月7日的200 m高光譜圖像分類精度評(píng)價(jià)

Table 2 200 m hyperspectral image classification accuracy evaluation on August 7

表3  9月25日的300 m高光譜圖像分類精度評(píng)價(jià)

Table 3 300 m hyperspectral image classification accuracy evaluation on September 25

4種分類方法的原理各不同，從人工參與程度、處理效率和分類精度等方面進(jìn)行比較分析可知，ISODATA法的分類精度雖然低于神經(jīng)網(wǎng)絡(luò)和馬氏距離法，但是不需要太多的人工參與，處理卻大大提高，且分類精度高于最大似然分類法。因此，它同時(shí)具備了監(jiān)督分類的高精度和非監(jiān)督分類的低時(shí)間的特點(diǎn)，在無(wú)人機(jī)低空遙感研究中，能夠以較低的工作成本快速獲取高精度的地物識(shí)別。圖5和圖6分別為8月7日200 m和9月25日300 m高光譜圖像基于ISODATA法的分類效果圖。從圖5可知，少量的雜草被錯(cuò)分為葡萄樹(shù)、西瓜苗、茄子苗，但是從大的區(qū)域來(lái)看，原味葡萄樹(shù)、西瓜苗、茄子苗的像素并未存在錯(cuò)分和漏分現(xiàn)象。在圖6中，少量的土壤，可能是沙土，被錯(cuò)分為大棚或公路斑馬線，其他地物不存在錯(cuò)分和漏分現(xiàn)象?？傮w而言，ISODATA法可實(shí)現(xiàn)無(wú)人機(jī)高光譜遙感對(duì)地面地物的快速有效的識(shí)別。

圖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.  結(jié)論

本研究基于無(wú)人機(jī)平臺(tái)搭載成像高光譜相機(jī)，分別獲取了兩個(gè)時(shí)間段不同高度的高光譜圖像，分析了不同高度同一地物光譜曲線的差異和基于不同高度下不同分類方法的分類精度評(píng)價(jià)。研究表明隨著無(wú)人機(jī)高度的升高，混合像元中土壤所占比例增加，從而造成“綠峰”降低、“紅谷”升高，紅邊位置發(fā)生了“紅移”等現(xiàn)象，研究結(jié)果與地面光譜儀ASD的研究結(jié)果一致[21-22]。ISODATA法的分類精度雖然不及神經(jīng)網(wǎng)絡(luò)和馬氏距離分類法，但是在綜合考慮人工參與程度、處理效率和分類精度時(shí)可知，ISODATA法對(duì)于快速實(shí)現(xiàn)地物的識(shí)別效果最佳，研究結(jié)果為利用無(wú)人機(jī)高光譜遙感進(jìn)行軍事偽裝快速識(shí)別、海洋油污面積估測(cè)、災(zāi)害損失評(píng)估等領(lǐng)域具有重大意義。

參考文獻(xiàn)

[1]  鄒湘伏，何清華，賀繼林.  無(wú)人機(jī)發(fā)展現(xiàn)狀及相關(guān)技術(shù)[J]. 飛航導(dǎo)彈，2006(10): 9 -14.

[2]  周清波. 國(guó)內(nèi)外農(nóng)情遙感現(xiàn)狀與發(fā)展趨勢(shì)[J]. 中國(guó)農(nóng)業(yè)資源與區(qū)劃，2004，25(5): 9 -14.

[3]  Sugiura R, Noguchi N, Ishii K. Remote sensing technology for vegetation monitoring using an unmanned helicopter[J]. Biosystems Engineering, 2005, 90(4): 369 -379.

[4]  李冰，劉镕源，劉素紅，等. 基于低空無(wú)人機(jī)遙感的冬小麥覆蓋度變化監(jiān)測(cè)[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2012，28(13): 160-165.

[5]  Bryson Mitch, Reid Alistair, Ramos Fabio, et al. Airborne vision-based mapping and  classification of large farmland environments[J]. Journal of Field Robotics, 2010, 27(5): 632-655.

[6]  Garcia Ruiz Francisco, Sankaran Sindhuja, Maja Joe Mari, et al. Comparison of two aerial  imaging platforms for identification of Huanglongbing-infected citrus trees[J]. Computers  and Electronics in Agriculture, 2013, 91: 106-115.

[7]  Córcoles J I, Ortega J F, Hernández D, et al. Estimation of leaf area index in onion (Allium cepa L.) using an unmanned aerial vehicle[J]. Biosystems Engineering, 2013, 115(1): 31- 42.

[8]  劉峰，劉素紅，向陽(yáng). 園地植被覆蓋度的無(wú)人機(jī)遙感監(jiān)測(cè)研究[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào)，2014，45(11): 250 - 257.   

[9]  高林，楊貴軍，王寶山等. 基于無(wú)人機(jī)遙感影像的大豆葉面積指數(shù)反演研究[J]. 中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào)，2015，23(7): 868-876.

[10]  田振坤，傅鶯鶯，劉素紅等. 基于無(wú)人機(jī)低空遙感的農(nóng)作物快速分類方法[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2013，29(7): 109-116.

[11]  王利民，劉佳，楊玲波等. 基于無(wú)人機(jī)影像的農(nóng)情遙感監(jiān)測(cè)應(yīng)用[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2013，29(18): 136 -145.

[12]  邴媛媛. 無(wú)人機(jī)遙感在某鐵礦礦區(qū)資源監(jiān)測(cè)中的應(yīng)用[M]. 遼寧，遼寧工程技術(shù)大學(xué)， 2008.

[13]  賽恒吉雅. 數(shù)字圖像鏡像變換方法的實(shí)現(xiàn)[J]. 信息科技，2010，18(2)：142-143.

[14]  張東彥. 基于高光譜成像技術(shù)的作物葉綠素信息診斷機(jī)理及方法研究[M]. 浙江，浙江大學(xué)，2012.

[15]  Liu Jiangui, Wu Changshan, Zhang Bing, et al. Spectral reflectance retrival of image data by imaging spectrometer PHI[J]. Journal of Remote Sensing, 1999, 3(4): 290-294.

[16]  Numata Izaya, Roberts Dar A, Chadwick Oliver A, et al. Evaluation of hyperspectral data for pasture estimate in the Brazilian Amazonusing field and imaging spectrometers[J]. Remote Sensing of Environment, 2008, 112:1569-1583.

[17]  Petropoulos George P, Arvanitis Kostas, Sigrimis Nick. Hyperion hyperspectral imagery analysis combined with machine learning classier for land use/cover mapping[J]. Expert Systems with Applications, 2012, 39:3800-3809.

[18]  Wang JianJun, Zhang Yun, Bussink Coen. Unsupervised multiple end member spectral mixture analysis based detection of opium poppy ?elds from an EO-1 hyperion image in Helmand, Afghanistan[J]. Science of the Total Environment, 2014, 476:1-6.

[19]  梅安新，彭望碌，秦其明，劉慧平. 遙感導(dǎo)論[M]. 北京，高等教育出版社，2005.

[20]  趙英時(shí). 遙感應(yīng)用分析原理與方法[M]．北京，科技出版社，2003.

[21]  儲(chǔ)旭. 不同覆蓋度條件下水稻葉層氮素營(yíng)養(yǎng)的高光譜監(jiān)測(cè)研究[M]，南京，南京農(nóng)業(yè)大學(xué)，2013.

[22]  Yao X, Ren H, Cao Z, et al. Detecting leaf nitrogen content in wheat with canopy hyperspectrum under different soil backgrounds. International Journal of Applied Earth Observation and Geoinformation, 2014, 32:114 -124.