清空記錄
歷史記錄
取消
清空記錄
歷史記錄
1、引言
草原生態(tài)系統(tǒng)對于促進(jìn)經(jīng)濟(jì)社會可持續(xù)發(fā)展、筑牢國家生態(tài)安全屏障、鑄牢中華民族共同體意識和鞏固邊疆少數(shù)民族地區(qū)穩(wěn)定等方面 具有戰(zhàn)略性作用。十八大以來,草原生態(tài)系統(tǒng)保護(hù)修復(fù)工作成效顯 著,部分地區(qū)草原生態(tài)恢復(fù)明顯,計(jì)劃到2025年和2035年 草原綜合植被蓋度要分別達(dá)到57%和60%以 上,但目前部分區(qū) 域的草原生態(tài)系統(tǒng)承受力和抵御干擾的耐受力極度脆弱,有向逐漸 裸露化或荒漠化發(fā)展的趨勢。草原荒漠化主要表現(xiàn)為植被覆蓋度 (FVC)減少和裸地面積增加等,傳統(tǒng)的地面小尺度覆蓋度采集 受天氣、時(shí)間和地形地貌的影響較大,且成本高、勞動(dòng)強(qiáng)度大、費(fèi) 時(shí)費(fèi)力;而衛(wèi)星雖然可以大尺度采集覆蓋度信息,但是收集速度緩 慢且獲得圖像的細(xì)節(jié)水平下降。無人機(jī)(UAV)+ 高 光 譜 可 獲 得中尺度區(qū)域的高空 - 譜分辨率圖像,兼顧了效率和精度要求,在 相對較小區(qū)域的草地退化監(jiān)測中將發(fā)揮重要作用,正成為傳統(tǒng)地面 監(jiān)測和航空、衛(wèi)星遙感的優(yōu)越補(bǔ)充。FVC定義為統(tǒng)計(jì)范圍內(nèi)植被冠層的垂直投影面積占土地面積的比例,是研究生態(tài)系統(tǒng)平衡、 土 地退化和土壤侵蝕等的關(guān)鍵參數(shù)之一,也是評價(jià)草地退化和荒漠化的敏感指標(biāo)。
目前低空UAV 遙感FVC反演法可以分為回歸模型法、 混合像元分解法和深度學(xué)習(xí)法等。UAV高光譜遙感滿足了當(dāng)前遙感朝著高空- 譜分辨率、多時(shí)相的發(fā)展方向,關(guān)鍵在于能否選擇或建立一種適用性強(qiáng)、耗時(shí)少、準(zhǔn)確性高的覆蓋度估算方法。利用光譜分辨率高、光譜波段信息豐富的高光譜成像儀結(jié)合低空UAV組建 UAV 高光譜遙感系統(tǒng),對試驗(yàn)區(qū)荒漠草原的高光譜影像進(jìn)行數(shù)據(jù)采集和數(shù)據(jù)預(yù)處理,并利用回歸模型法、深度學(xué)習(xí)法ResNet18模型及其改進(jìn)的3D-ResNet18模型對影像中的植 被和裸土等地物進(jìn)行分類,并在分類基礎(chǔ)上開展基于UAV高光譜遙感信息的荒漠草原FVC計(jì)算方法研究。旨在實(shí)現(xiàn)利用高光譜遙感系統(tǒng)對荒漠草原從數(shù)據(jù)獲取、數(shù)據(jù)分析、理論方法及驗(yàn)證等完整流程的嘗試,為草原退化評價(jià)及治理提供依據(jù),對于草原畜牧養(yǎng)殖和生態(tài)環(huán)境保護(hù)均具有重要的現(xiàn)實(shí)意義,符合國家和地區(qū) “助力山水林田湖草沙生命共同體建設(shè) ”的發(fā)展需求。
2、數(shù)據(jù)采集與分析
2.1 研究區(qū)概況
研究區(qū)位于烏蘭察布市四子王旗格根塔拉草原 (41°47′33″N,111°53′58″E),詳 見 圖1。
圖1 研究區(qū)位置圖
海拔高度介于1100—1200m之間,降水量稀少,水資源匱乏,具有荒漠草原典型的地域代表性,是典型草原與荒漠的緩沖地帶。研究區(qū)植被草層稀 疏低矮、交錯(cuò)覆蓋、種類匱乏,屬短花針茅荒漠草原地帶。植被 平均高度為8cm,平均蓋度為17%—25%。
2.2 試驗(yàn)設(shè)計(jì)
2.2.1野外調(diào)查及樣方布置
野外地面調(diào)查包括記錄GPS信息、樣方編號,實(shí)測樣方FVC。為了樣方框堅(jiān)固耐用且避免被風(fēng)刮走,由白色PVC管制作,用 U 型鐵絲將其固定于地面上,隨機(jī)布置1m×1m的樣方框60個(gè)。FVC 的確定采用照相法。樣方冠層照片通過在樣方中心的垂直上方2.2m 處拍攝全范圍俯視照實(shí)現(xiàn)。選用 RTK 測定試驗(yàn)區(qū)和樣方的具體位置,收集了所有研究樣方框的 GPS信息,選定地面控制點(diǎn)對圖像進(jìn)行幾何校正。
2.2.2空中數(shù)據(jù)采集
采用UAV 懸停方式采集試驗(yàn)區(qū)內(nèi)地物高光譜遙感影像,結(jié)合2021年草原氣候特點(diǎn)和牧草生長期特性,在牧草的生長茂盛期———2021年7月27日— 8月8日時(shí)采集數(shù)據(jù)。為保證采集質(zhì)量,選擇晴朗、無云、光照條件較好,無卷云、濃積云等,風(fēng)力較小時(shí)進(jìn)行測量,采集時(shí)間為10∶00—14∶00,且要求每10—20min使用標(biāo)準(zhǔn)白板進(jìn)行1次校正,以便消除光照強(qiáng)度變化對UAV高光譜圖像產(chǎn)生的影響。UAV搭載高光譜儀垂直于地面的方向30m高度測量研究地塊植被的冠層、群落的高光譜數(shù)據(jù)。每個(gè)懸停點(diǎn)采集2幅高光譜遙感影像,單幅圖像采集時(shí)間耗時(shí)7s,包含懸停點(diǎn)間飛行時(shí)間,1個(gè)架次約可采集86幅遙感影像,即43個(gè)懸停點(diǎn)。
2.2.3 數(shù)據(jù)預(yù)處理
首先,通過人工檢查去除過曝、欠曝、彎曲、抖動(dòng)的遙感影像,選出成像質(zhì)量最好的遙感影像。其次使用光譜專業(yè)軟件進(jìn)行反射率校正,并識別真正的反射率值和感興趣的特征。最后,使用方法進(jìn)行數(shù)據(jù)降維,計(jì)算式如下:
(1)式中,X為張量,r?yàn)楦吖庾V圖像中的行數(shù),c為高光譜圖像中的列數(shù),b為高光譜圖像中的維數(shù)。式(1)計(jì)算的高光譜圖像F-norm2在 0—67波段(波長為397—548.2nm)的F-norm2的斜率較小,即表示該波段范圍內(nèi)包含的信息量較少;而在 214—256波段(波 長 為 909.7—1019.4nm)的F-norm2的斜率較大,即表示該波段范圍內(nèi)存在明顯的噪聲干擾,因此經(jīng)行值降維后保留的68—213波段(波長為550.5—906.8nm),單個(gè)圖像所占空間大小也隨之大幅減少 ,有效提高了數(shù)據(jù)后處理的效率。
2.3 驗(yàn)證值的確定
驗(yàn)證值的確定一方面是依據(jù)實(shí)地設(shè)立1m×1m 的樣方框通過照相 法確定FVC,但是該方法存在一定的主觀判斷性而且只代表了局部的FVC。研究表明,基于照相法和全圖人工目視解譯相結(jié)合的方法對FVC的提取精度較高。為了雙重保險(xiǎn),后續(xù)又結(jié)合人工目視解譯方法對整張UAV高光譜圖像中植物群落和非植物群落(枯草和裸土)進(jìn)行覆蓋度提取,其結(jié)果作為地表 FVC 的 驗(yàn) 證 值,對各植被指數(shù)提取的草地 FVC 的精度進(jìn)行驗(yàn)證。本文從1080 組數(shù)據(jù)中選?。玻敖M有代表性的數(shù)據(jù)作為示例進(jìn)行研究。20 組數(shù)據(jù)分別命名為T1、T2、T3、…、T19、 T20。由于圖像較多,將其中的T1、T8、T16組高光譜數(shù)據(jù)設(shè)為展示圖像,T1、T8、T16組數(shù)據(jù)覆蓋度的驗(yàn)證值詳見表1。
表1植物群落和非植物群落驗(yàn)證值
3、研究方法
3.1 回歸模型法
回歸模型法是通過對遙感數(shù)據(jù)的某一波段、波段組合或利用遙感數(shù)據(jù)計(jì)算出的植被指數(shù)與FVC進(jìn)行回歸分析,建立經(jīng)驗(yàn)估算模型。通過回歸模型法構(gòu)建荒漠草原FVC估算模型,選擇了2種傳統(tǒng)的適合于低密度覆蓋綠色植被區(qū)域的植被指數(shù),分別為歸一化植 被指 數(shù)(NDVI)和土壤調(diào)整植被指數(shù) (SAVI),在分析植被和土壤光譜曲線差異的基礎(chǔ) 上,通過簡單波段自相關(guān)選擇法及編程篩選最敏感的特征波段及組合,確定了植被指數(shù)相關(guān)系數(shù)較高的波段組合為525—600nm、620—780nm 和 780—900nm,相關(guān)系數(shù)均達(dá)0.8以上。對基于傳統(tǒng)植被指數(shù)的像元二分模型進(jìn)行了FVC估算,并選取 20組數(shù)據(jù)進(jìn)行精度驗(yàn)證。選擇SAVI和 NDVI進(jìn)行 連續(xù)統(tǒng)去除及光譜增強(qiáng),并確定特征波段和波段最佳組合,對高光譜圖像進(jìn)行連續(xù)統(tǒng)去除后,峰值點(diǎn)對應(yīng)的反射率為1,而其他點(diǎn)的反射率均小于1,植物群落在可見光范圍中的500nm 和680nm 左右的吸收谷特征被放大,綠光550nm 附近反射峰更加明顯, 紅邊斜率增加明顯,這些變化有利于對光譜吸收特征波段及參數(shù)的提取,因此提出了綠光連續(xù)值土壤調(diào)整植被指數(shù)(G_CR_SAVI)和綠光連續(xù)值歸一化 植被指數(shù)(G_CR_NDVI),原始及改進(jìn)的植被指數(shù)詳見表2。
表2原始及改進(jìn)的植被指數(shù)
3.2 ResNet模型
該模型將殘差學(xué)習(xí)的思想引入到網(wǎng)絡(luò)結(jié)構(gòu)中,通過跨層鏈接將殘差塊的輸入與輸出進(jìn)行疊加求和,此項(xiàng)操作提升了反向傳播的效率,有效解決了深層網(wǎng)絡(luò)中出現(xiàn)的梯度爆炸及消失問題,提高了神經(jīng)網(wǎng)絡(luò)訓(xùn)練的速度,實(shí)現(xiàn)了利用深層網(wǎng)絡(luò) 結(jié)構(gòu)提取更細(xì)的特征。利用 Pytorch模塊搭建 ResNet18CNN,該模型 共有20層,包括17個(gè) 卷 積 層、2個(gè)池化層和1個(gè)全連接層,模型結(jié)構(gòu)圖詳見圖2。
圖2 ResNet深度學(xué)習(xí)模型結(jié)構(gòu)圖
每殘差塊中進(jìn)行2次卷積層運(yùn)算,并與卷積前的參數(shù)進(jìn)行疊加,隨后利用激活函數(shù)ReLU 激活進(jìn)入下一個(gè)殘差塊。卷積層的卷積核數(shù)最少為64個(gè),最多為512個(gè)。
3.3改進(jìn)3D-ResNet18模型
3D卷積核為 一 個(gè) 立 方 體,由3D卷積核組成的CNN 稱為3D 卷積神經(jīng)網(wǎng)絡(luò),通過3D卷積核 可以同時(shí)提取高光譜圖像lines、samples和bands方 向上的特征。公式如下:
(2) 式中,f為激活函數(shù),l、w、h 分別代表卷積核立方體 的長、寬、高,m、n 分別代表上一層和本層的卷積核數(shù),ulwhkn為圖像上(l,w,h)位 置上的值與上一層第m個(gè)卷積核的計(jì)算值,uδερkmn為 圖 像(δ,ε,ρ)位置上的值與卷積核的第n層、第k個(gè)卷積核的計(jì)算值。利用F-norm2降低噪聲干擾和高光譜數(shù)據(jù)的維數(shù)。在典型的深度學(xué)習(xí)模型 ResNet18的基礎(chǔ)上,將其2D卷積核3D卷積核,即3D-ResNet18。改進(jìn)后的3D-ResNet18模 型 與 ResNet18模型具有相同的深度,共有20層,包括17個(gè)卷積層、2個(gè)池化層和1個(gè)全連接層,17個(gè)卷積層被分成1個(gè)卷積層和4個(gè)殘差塊。將卷積層中的2D卷積核改進(jìn)為3D卷積核,在圖3中用Conv3×3×3-N表 示,N為卷積核數(shù)量,分別為64、64、128、256、512。每殘差塊中進(jìn)行兩次卷積層運(yùn)算,并與卷積前的參數(shù)進(jìn)行疊加,隨后利用激 活 函 數(shù) ReLU 激活進(jìn)入下一個(gè) 殘差塊。為了提高處理速度,程序用高光譜影像裁剪為500lines×500samples×146bands,總 共250000個(gè) 像素。為了識別特征,隨機(jī)選擇60%的標(biāo)記樣本 作 為訓(xùn)練數(shù)據(jù),其余的40%作為測試數(shù)據(jù)。
4、結(jié)果分析
4.1 評價(jià)指標(biāo)
在進(jìn)行覆蓋度估算模型的精度分析時(shí),選 用以下3個(gè)指標(biāo)進(jìn)行驗(yàn)證 分 析:估算值和實(shí)測值間的決定系數(shù)(R2)、估算值和實(shí)測值間的均方根誤差(RMSE )和估算 精 度(EA )。當(dāng) R2 值越 接近1、RMSE的值越小且估算精度越 時(shí),說明估測方程所得到的估算值與實(shí)測值的擬合效果就越好,模型的精度就越高。 R2、RMSE、EA 的計(jì)算式如下:
式中,為第i個(gè)樣本的估算,yi 為第i個(gè)樣本的實(shí) 測值,:為實(shí)測值的平均值,x 為樣本數(shù)量。
式中,Rm為均方根誤差,Me 為實(shí)測值的均值,EA 為估算精度。
4.2回歸模型法
改進(jìn)植被指數(shù)中提取的 FVC 與驗(yàn)證值最為接近,其RMSE和EA 分別為0.018和 95.42%;其次為 G_CR_NDVI,RMSE 和EA 分別為 0.018和95.36%,SAVI提取的FVC與驗(yàn)證值較為接近 ,其RMSE和EA分別為0.032和91.64%,詳見表3。
表3 植被覆蓋度與各方法對應(yīng)的模型精度
從估算值與驗(yàn)證值之間的散點(diǎn)圖來看,G_CR_SAVI和 G_CR_NDVI估算值緊密分 布在1∶1線的兩側(cè),其中 G_CR_SAVI估算值與驗(yàn)證值之間的R2高達(dá)0.942,G_CR_NDVI估算值與驗(yàn)證值之間的R2高達(dá)0.939?;诟倪M(jìn)植被指數(shù)與傳統(tǒng)植被指數(shù)的荒漠草原UAV高光譜圖像FVC估算精度大小依次為 G_CR_SAVI>G_CR_NDVI> SAVI>NDVI。分析表明:通過光譜增強(qiáng)及最佳波段選擇的改進(jìn) G_CR_SAVI和 G_CR_NDVI植被指數(shù)優(yōu)于傳統(tǒng)的 SAVI、NDVI植被指數(shù),FVC 估 算 精 度較傳統(tǒng)的植被指數(shù) 高出約4%,G_CR_SAVI在4個(gè)植被指數(shù)中的提取效果最好,最適合于荒漠草原UAV 高光譜影像 FVC的提取。最佳波段組合的植被指數(shù)計(jì)算結(jié)果詳見圖4。
圖4 植被指數(shù)變換分布圖
4.3 深度學(xué)習(xí)法
為實(shí)現(xiàn)基于UAV 高光譜遙感的荒漠草原覆蓋 度人工智能化估算,制作了荒漠草原 UAV 高光譜覆 蓋度數(shù)據(jù)集,通過 ResNet18經(jīng)典深度學(xué)習(xí)網(wǎng)絡(luò)模型對覆蓋度數(shù)據(jù)集進(jìn)行分類,發(fā)現(xiàn)其對覆蓋度數(shù)據(jù)集取得了較理想的總體分類精度,并對覆蓋度數(shù)據(jù)集中的植被和裸土取得較高的單體分類精度。為進(jìn)一步實(shí)現(xiàn)對荒漠草原UAV 高光譜數(shù)據(jù)中的覆蓋度高精度估算,基 于 ResNet模 型 建立3D卷積模型3D- ResNet,并進(jìn)行模型結(jié)構(gòu)與超參數(shù)優(yōu)化,通過對卷積核數(shù)、卷積核尺寸、Batchsize等參數(shù)優(yōu)化與對比,發(fā)現(xiàn)對覆蓋度數(shù)據(jù)集分類性能最佳的模型為3D-Res- Net,總體估算精度達(dá)97.56%,覆蓋度估算的RMSE 為0.016,詳見表3。
3D卷積網(wǎng)絡(luò)模型對覆蓋度數(shù)據(jù)集具有更好的提 取性能。相比于2D卷積僅能同時(shí)提取覆蓋度數(shù)據(jù)集中的2D光譜 信 息,3D 卷積可以同時(shí)提取光譜- 空間聯(lián)合信息,充分利用了高光譜覆蓋度數(shù)據(jù)集的數(shù)據(jù)特性,對細(xì)小特征的提取展現(xiàn)出良好潛力,得益于3D卷積核可以同時(shí)提取高光譜數(shù)據(jù)中3個(gè)維度的光譜信息,相比2D卷積核同時(shí)提取 X、Y2個(gè)維度的光譜信息,3D卷積核對地物的空間特征提取能力更強(qiáng)。通過對模型結(jié)構(gòu)與超參數(shù)優(yōu)化逐步提高了模型 性能,同時(shí)發(fā)現(xiàn)具有更卷積結(jié)構(gòu)的3D-ResNet模 型展現(xiàn)出更好的性能。部分?jǐn)?shù)據(jù)可視化結(jié)果如圖5 所示。
圖5覆蓋度估算模型可視化結(jié)果
3D-ResNet模型實(shí)現(xiàn)了對覆蓋度數(shù)據(jù)集地物的高精度、高效率、智能化識別,為基于UAV高光譜 遙感的荒漠化草原覆蓋度的快速、高效、精準(zhǔn)統(tǒng)計(jì)奠定基礎(chǔ)。
5、總結(jié)
針對荒漠草原覆蓋度提取效率和精度較低的現(xiàn)狀,基于人工智能技術(shù)和3D 卷積深度學(xué)習(xí)方法,深挖了高光譜數(shù)據(jù)中的空間特征,建立了高效準(zhǔn)確的高光譜遙感圖像的覆蓋度估算模型,驗(yàn)證了3D-Res- Net模型在計(jì)算FVC為指標(biāo)的草地退化評價(jià)的可行性和精確性。就估算精度而言,深度學(xué)習(xí)法優(yōu)于回歸模型法,尤其是經(jīng)過3D卷積核改進(jìn)的深度學(xué)習(xí)模型最佳,總體估算精度達(dá)97.56%,可以實(shí)現(xiàn)智能化、 高精度、準(zhǔn)確地荒漠草原 FVC提取。
推薦:
便攜式地物光譜儀iSpecField-NIR/WNIR
專門用于野外遙感測量、土壤環(huán)境、礦物地質(zhì)勘探等領(lǐng)域的最新明星產(chǎn)品,由于其操作靈活、便攜方便、光譜測試速度快、光譜數(shù)據(jù)準(zhǔn)確是一款真正意義上便攜式地物光譜儀。
無人機(jī)機(jī)載高光譜成像系統(tǒng)iSpecHyper-VM100
一款基于小型多旋翼無人機(jī)機(jī)載高光譜成像系統(tǒng),該系統(tǒng)由高光譜成像相機(jī)、穩(wěn)定云臺、機(jī)載控制與數(shù)據(jù)采集模塊、機(jī)載供電模塊等部分組成。無人機(jī)機(jī)載高光譜成像系統(tǒng)通過獨(dú)特的內(nèi)置式或外部掃描和穩(wěn)定控制,有效地解決了在微型無人機(jī)搭載推掃式高光譜照相機(jī)時(shí),由于振動(dòng)引起的圖像質(zhì)量較差的問題,并具備較高的光譜分辨率和良好的成像性能。
便攜式高光譜成像系統(tǒng)iSpecHyper-VS1000
專門用于公安刑偵、物證鑒定、醫(yī)學(xué)醫(yī)療、精準(zhǔn)農(nóng)業(yè)、礦物地質(zhì)勘探等領(lǐng)域的最新產(chǎn)品,主要優(yōu)勢具有體積小、幀率高、高光譜分辨率高、高像質(zhì)等性價(jià)比特點(diǎn)采用了透射光柵內(nèi)推掃原理高光譜成像,系統(tǒng)集成高性能數(shù)據(jù)采集與分析處理系統(tǒng),高速USB3.0接口傳輸,全靶面高成像質(zhì)量光學(xué)設(shè)計(jì),物鏡接口為標(biāo)準(zhǔn)C-Mount,可根據(jù)用戶需求更換物鏡。