無人機高光譜在樹種分類識別上的應用研究
0 引言
高光譜遙感具有光譜分辨率高、波段範圍窄、圖譜合一、連續成像等特點,能夠區分出地物光譜的細微差別,探測到其他寬波段遙感無法探測的信息。因此,高光譜遙感在生態、大(da)氣(qi)和(he)海(hai)洋(yang)等(deng)諸(zhu)多(duo)應(ying)用(yong)領(ling)域(yu)具(ju)有(you)很(hen)大(da)優(you)勢(shi)。近(jin)年(nian)來(lai),高(gao)光(guang)譜(pu)遙(yao)感(gan)在(zai)林(lin)業(ye)方(fang)麵(mian)的(de)一(yi)個(ge)重(zhong)要(yao)應(ying)用(yong)是(shi)對(dui)森(sen)林(lin)樹(shu)種(zhong)類(lei)型(xing)進(jin)行(xing)識(shi)別(bie)。森(sen)林(lin)樹(shu)種(zhong)類(lei)型(xing)識(shi)別(bie)的(de)主(zhu)要(yao)目(mu)的(de)是(shi)提(ti)取(qu)樹(shu)種(zhong)的(de)專(zhuan)題(ti)信(xin)息(xi),為(wei)劃(hua)分(fen)森(sen)林(lin)類(lei)型(xing)、繪製林相圖和清查森林資源提供基礎和依據。目前,國內外利用高光譜遙感進行樹種識別主要是從葉片、冠層和高光譜影像3個研究尺度開展。基於葉片的樹種識別主要是對葉片反射率及其變換形式運用統計方法、遺傳算法等進行分析,以樹種識別的可行性分析與識別潛力為主要研究內容;基於冠層的樹種識別主要運用光譜信息散度法、光譜角填圖法等基於光譜信息的遙感圖像分類方法,並利用地物光譜儀獲取的林分冠層反射率曲線,進行樹種分類;基於高光譜影像的樹種識別主要通過對影像進行去噪降維等預處理後,運用監督或非監督分類的方法進行樹種識別。
國外已有很多學者對高光譜樹種識別進行了研究。Gong等利用ANN分類法對光譜數據進行判別,區分出1種闊葉樹種和6種針葉樹種,分類精度大於90%;Martin等利用AVIRIS高光譜數據與樹種葉片化學成分之間的關係,鑒別出11種樹種類型,可有效進行樹種分類;Petropoulos等分別采用支持向量機和基於對象的分類方法,對Hyperion高光譜影像進行土地覆蓋類型分類,雖2 zhongfenleixiaoguojunjiaohao,danjiyuduixiangdefenleifangfajingdugenggao。guoneiyeyouyuelaiyueduodexuezhejinxingsenlinshuzhongshibietanjiu。tongqingxidengliyongguangpuboxingpipeisuanfaduiMAIS高光譜影像進行植被類型識別,獲得了潘陽湖典型濕地的植被分類圖,對高光譜的樹種識別提供了實用依據;王圓圓等采用隨機子空間法對OMIS高光譜遙感數據進行識別,並利用遺傳算法來提高分類精度;劉秀英等使用地物光譜儀測得的光譜數據,采用逐步判別分析方法、特征波段選擇等方法識別出4種zhong樹shu種zhong。綜zong上shang所suo述shu,國guo內nei外wai利li用yong高gao光guang譜pu數shu據ju進jin行xing樹shu種zhong識shi別bie的de研yan究jiu已yi取qu得de階jie段duan性xing的de進jin展zhan,從cong研yan究jiu方fang法fa看kan,主zhu要yao基ji於yu不bu同tong樹shu種zhong具ju有you不bu同tong光guang譜pu特te征zheng的de原yuan理li,通tong過guo特te征zheng波bo段duan的de合he理li選xuan擇ze,或huo者zhe通tong過guo波bo段duan變bian換huan對dui高gao光guang譜pu數shu據ju進jin行xing降jiang維wei來lai識shi別bie樹shu種zhong。
gaoguangpuyingxiangboduanshuduo,xinxiliangda,weidiwudejingxishibietigongyoushidetongshi,yedailaileshujuliangduo,boduanjianxiangguanxingda,chulijingduhexiaolvxiajiangdewenti。benwenyishanghaijiaotongdaxuezhiwuyuanweiyanjiuqu,liyongwurenjigaoguangpushuju,yunyongzuijiazhishuboduanxuanzefahesee5.0決策樹自動分類相結合進行樹種識別,有利於無人機高光譜數據分類識別精度的進一步提高。
1 研究區及數據源
研究區位於上海交通大學閔行校區植物園,如圖1所示。

圖1 研究區域位置
2017年7月28日,運用大疆的無人機M600 Pro搭載四川雙利合譜科技有限公司自足研發的GaiaSky-mini2在上海交通大學閔行校區植物園上空進行數據采集。數據獲取當天晴朗無(少)雲,飛行高度為200米,采用無人機懸停高光譜相機內置推掃的方式獲取高光譜數據。獲取的高光譜影像數據具體參數見表1,無人機高光譜相機如圖2所示。
表1 高光譜影像數據具體參數
|
參數 |
光譜範圍/nm |
波段個數 |
圖像分辨率 |
空間分辨率 |
單景幅寬/m |
光譜分辨率/nm |
鏡頭焦距/mm |
|
GaiaSky |
400-1000 |
176 |
1920*1400 |
10cm |
95*95 |
3.5 |
18.5 |

圖2 無人機高光譜起飛示意圖
無人機搭載的高光譜相機獲取的高光譜數據為數字量化值(簡稱DN值),無物理意義,需轉化為具有物理意義的反射率數據。具體轉化方式如公式1所示。

(1)
其中,Reftarget為目標物的反射率,DNtarget為目標物的DN值,DNdark為相機本身的暗點流DN值,DNwhite為參考板的DN值,Refwhite為參考板的反射率。
降噪的目的主要是為了突出圖像的特征信息,提高圖像的信噪比。本研究利用ENVI5.3自帶的curve smoothing對圖像進行降噪處理,處理結果如圖3。
.jpg)
圖3光譜降噪前後的反射率(左為降噪前,右為降噪後)
圖4分別列舉了四種不同樹種、雜草、樹蔭、水體和裸土的光譜反射率曲線,從圖4可知,水體、luotuheshuyindefanshelvguangpuquxianyuzhiwudefanshelvguangpuquxianchayijiaoda,keliyongboduanyuzhidengfangfajiangqiyuzhiwuqufenkai,danduiyuzacaohebutongshuzhongzhijiandefenleishibiezebunengtongguojiandandeboduanyuzhijinxingqufen。

圖4不同樹種與地麵其他地物的光譜反射率
目前應用比較廣泛的最佳波段選取方法有各波段信息量的比較、波段間相關性比較、最佳指數法(O IF)、各波段數據的信息熵和聯合熵、協方差矩陣特征值法、波段指數法等。
在各種方法中,由美國查維茨提出的最佳指數法( OIF)綜合考慮單波段圖像的信息量及各波段間的相關性,更接近於波段選擇的原則,且計算簡單,易於實現,得到廣泛的應用。OIF指數的計算公式如下:
(2)
其中:Si為第i個波段的標準差,Rij為i、j 兩波段的相關係數。對n波段圖像,先統計計算單波段圖像的標準差,計算各波段間的相關係數矩陣,再分別求出所有可能的波段組合對應的OIF指數,根據該指數大小來判斷各種波段組合的優劣。指數越大,則相應組合影像所包含的信息量就越大。對OIF指數從大到小進行排序,最大O IF指數對應的波段組合即為最佳波段組合。以圖5為例,利用OIF篩選高光譜數據的特別波段,表2分別列舉了前10個特征波段組合及對應的OIF指數。從表2可知利用最佳指數法篩選的特征波段前10個波長組合對應的波長位置相差不大,且OIF指數值也十分接近。
.jpg)
圖5 無人機高光譜采集的單景影像(RGB彩色合成)
表2 OIF指數篩選的前10個波段組合
|
組合排列 |
波段組合(nm) |
OIF指數 |
|
1 |
485.6、784.2、878 |
0.29319887 |
|
2 |
485.6、780.7、878 |
0.29317102 |
|
3 |
485.6、784.2、881.6 |
0.29310915 |
|
4 |
485.6、784.2、874.3 |
0.29310489 |
|
5 |
485.6、780.7、881.6 |
0.29308453 |
|
6 |
485.6、780.7、878 |
0.29307423 |
|
7 |
482.4、784.2、878 |
0.29307401 |
|
8 |
482.4、780.7、878 |
0.29304684 |
|
9 |
485.6、863.4、870.7 |
0.29303626 |
|
10 |
485.6、784.2、874.3 |
0.29300225 |
see5.0機器學習規則軟件是美國USGS在完成國家土地覆蓋製圖(NLCD)項目中開發的一個自動提取分類規則的數據挖掘工具。表3為利用全波段和特征波段進行樹種識別的總體分類精度、Kappa係數和運行時間,從表中可知,利用特征波段運行時間更短,且精度與利用全波段進行分析的精度接近。
表3 不同組合分類精度對比分析
|
分類方法 |
總體分類精度 |
Kappa係數 |
運行時間 |
|
全波段+see5.0 |
96.7% |
0.94 |
38min |
|
特征波段+see5.0 |
95.3% |
0.93 |
1min |
考慮到用戶科研數據的保密性,本研究僅用單景高光譜影像數據進行不同樹種的分類識別,分類識別結果如圖6所示。用戶可利用無人機高光譜相機獲取的多組數據進行拚接,然後再進行樹種的分類,步驟和算法均是相同的。

圖6 利用特征波段+see5.0的分類效果圖
本研究利用四川雙利合譜科技有限公司自主研發的無人機高光譜相機GaiaSky-mini2,通過鏡像處理、黑白幀校正和降噪等預處理,並采用最佳指數法篩選特征波長,選用see5.0分類方法進行分類識別不同的樹種和其他地物,得出如下結論:
1) 該無人機高光譜數據共176個波段,采用最佳指數法獲取的三個特征波長,可以很好地刻畫區分各個類別,實現樹種級的識別;
2) 采用特征波長和全波段進行地物的區分,其分類精度差別不大;
3) 經實地驗證,各類別範圍和分布區域較準,特征波段+see5.0的總體分類精度達到95.3%,Kappa係數為 0.93。
4) 由you於yu數shu據ju的de保bao密mi性xing,未wei對dui拚pin接jie好hao的de無wu人ren機ji高gao光guang譜pu影ying像xiang數shu據ju進jin行xing不bu同tong樹shu種zhong的de分fen類lei識shi別bie。拚pin接jie好hao的de高gao光guang譜pu影ying像xiang數shu據ju覆fu蓋gai的de地di麵mian範fan圍wei更geng廣guang,樹shu種zhong更geng多duo,影ying像xiang中zhong的de其qi他ta地di物wu也ye更geng多duo,不bu同tong地di物wu之zhi間jian可ke能neng存cun在zai“同物異譜,異物同譜”的現象,對於不同樹種之間分分類識別難度更大,可能需要進一步改進波段選擇的方法和分類方法,方能更好地區分不同的樹種。
地址:北京市海澱區中關村大街19號新中關B座北翼1701-1706室
電話:
傳真:
郵箱: