一、引言

草原生態系統是陸地生態系統的重要組成部分。我國草地面積約為4億hm2，每年平均可為畜牧業發展提供3.5億噸優質牧草。但由于自然因素與人為因素影響，全球草原生態系統出現退化問題，逐漸形成荒漠草原。據統計，我國超過90%的草原出現不同程度退化現象。其中，內蒙古自治區草原總面積約為8.8×107hm2，退化面積約占73.5%。草原退化對農牧民生活及農牧業經濟發展造成嚴重威脅。荒漠草原地表微斑塊主要包括各種群植被、裸土、鼠洞。草原退化主要表現為植被矮小、植被整體覆蓋度逐漸減少、鼠洞數量增加等。因此，地表微斑塊精確識別與分類是荒漠草原研究的重要指標。傳統荒漠草原研究主要依靠人工實地監測，精確度較高但費時費力。為滿足大空間尺度觀測需求，遙感技術應用于草原研究。遙感采用非接觸、遠距離探測技術，可實時獲取圖像數據，與人工調查相比具有成本低、實驗不受人主觀因素影響等優點。常用遙感方式為衛星遙感，分辨率在1m左右，應用于大面積植被動態監測。荒漠草原使用衛星遙感監測具有局限性。

近年來，低空無人機遙感系統逐漸興起，彌補了衛星遙感與人工調查間尺度空缺。以無人機搭載高光譜成像儀組建無人機高光譜遙感系統拍攝圖像，不僅分辨率高、數據源豐富，而且方便快捷、不損壞原有實驗環境，可重復進行實驗。

隨著深度學習的發展，卷積神經網絡不斷應用于高光譜圖像分類處理任務。CNN是能夠進行深度學習的人工神經網絡，能做到自動提取圖像特征。本文通過構建無人機高光譜遙感系統，采集荒漠草原地表微斑塊高光譜圖像，建立輕量化卷積神經網絡模型―高分辨率圖像卷積神經網，用于地表微斑塊識別與分類任務，為荒漠草原退化監測提供理論基礎。

二、材料與方法

2.1 研究區概況

試驗區域如圖1，地處內蒙古自治區烏蘭察布市格根塔拉草原(N41.78o、E111.88o)，草原類型為荒漠草原。植被種群稀疏低矮，平均高度10cm。草類主要為短花針茅、無芒隱子草等。鼠類主要有長爪沙鼠、達烏爾黃鼠等。

圖1 試驗區域衛星圖

2.2數據采集

高光譜圖像數據采集時間：2022年的7月至8月，中午11:00~14:00。采集條件：光照充足、無云霧遮擋，少風或微風。試驗區域總面積為2.5hm2。在試驗區域內，隨機設置50組樣方，樣方尺寸為1m×1m。樣方內部含有植被、裸土、鼠洞，外部使用實驗輔助物（地墊與小旗）進行標記。此外，實驗輔助標記物還有充當草原環境中因人為活動導致外來垃圾的作用。無人機飛行高度30m，空間分辨率1.73cm，采集地物樣方光譜圖像尺寸大小為696lines×775samples×256bands，且對每個樣方拍攝3次，以提高數據的可用性。

2.3 反射率矯正

無人機飛行拍攝過程中因受環境因素影響導致部分拍攝圖像發生扭曲變形，故首先將變形及拍攝質量較差圖像進行人工去除，其次，將初步篩選數據（包括黑幀、白幀、樣方原始圖像）導入SpecView軟件進行輻射校正。校正計算如公式：

式中：R為反射率校正后圖像，Iraw為原始圖像，Iwhite為白幀圖，Idark為黑幀圖。

圖像經反射率校正后，使用ENVI5.3軟件提取每類地物純凈像元各100個進行空間統計（因鼠洞像元較少，故提取20個進行研究），并使用OriginPro2021軟件繪制每類地物反射率曲線如圖2所示。

圖2 地物反射率校正圖

從圖2中可看出每類地物反射率差異明顯，其中，地墊反射率波動較大。小旗在460Band后急劇上升，560Band后增長趨于平穩。裸土與植被反射率差異處為543~687Band，植被存在明顯‘波峰’與‘波谷’。鼠洞反射率增長最為緩慢。每類地物同波段反射率變化不一為后續地物分類提供可能。

2.4 數據預處理

高光譜圖像具有豐富的光譜和空間信息，通常包含幾十到數百個連續光譜波段，但波段間存在較高的相關性和冗余度，這將使模型在訓練階段容易出現Hughes現象。因此，高光譜數據使用前應做數據降維處理。本文采用廣泛應用于高光譜數據降維的主成分分析算法作為特征提取的方法，有效剔除冗余信息。同時為有效防止圖像邊緣大量扭曲混合像元對實驗影響，對數據尺寸進行剪裁，最終獲取高光譜圖像數據。

根據研究目的，本次實驗將草原地物分為植被、裸土、地墊、小旗、鼠洞共五類。通過ENVI5.3軟件對圖像進行ROI感興趣區域標簽處理，每類樣本數量分別為31600、27419、686、179、32個，共計59916個樣本。本次實驗將50%的地物樣本劃分為訓練集，為保證實驗可靠性將訓練集中50%的數據劃分為驗證集。在圖像識別與分類處理任務中，深度學習效果嚴重依賴于樣本數量，當訓練樣本較少時，易使網絡陷入局部最優。本實驗為增加訓練樣本數量，通過圖像隨機翻轉方法進行數據增強，使訓練集樣本數量增加至14979個。此外，為一定程度上解決小樣本數據集過擬合問題，對樣本數據進行隨機打亂處理。樣本數據如表1所示。

表1 樣本數據表

三、草原物種分類模型

卷積神經網絡特征提取方法是影響模型優劣的關鍵因素。本文針對荒漠草原地表微斑塊分類問題，提出高分辨率圖像卷積神經網絡模型。該模型主要特點包括使用帶有不同大小卷積核的深度可分離卷積作為各分支結構，以降低計算量和參數量。通過逐步使用不同感受野的卷積核，將圖像底層特征提取為高層語義特征。同時，引入殘差結構將底層特征與高層語義特征進行融合，以防止卷積運算中的特征丟失和網絡模型退化問題的發生。在卷積層最后，使用轉置卷積將提取的高層語義特征圖映射為高分辨率特征圖像，從而降低圖像中混合像元對最終分類效果的影響，并提高模型分類精度。該模型屬于低參數量的輕量化網絡模型，共包含兩層分支結構與一層組合結構。分支結構內部包括輸入層、三個深度可分離卷積層、兩個殘差結構、最大池化層與轉置卷積層。組合結構中包括維度拼接層、全連接層與分類層。此外，每層卷積運算結束后，均進行批量歸一化(BatchNormalization)與激活函數(Relu)操作，避免梯度消失與爆炸，加快模型收斂速度。具體網絡模型結構如圖3所示，網絡模型運行參數如表2所示。

圖3 高分辨率圖像卷積神經網絡

表2網絡模型參數

結果與分析

4.1 實驗結果

本模型在訓練實驗過程中，除模型訓練初期驗證集損失出現震蕩外，訓練集與驗證集整體呈下降趨勢。在第20個epoch后趨于平穩，此時共花費144s，訓練集損失值穩定變化在0.023~0.025，驗證集損失值穩定變化在0.062~0.065。為保證訓練模型穩定，30個epoch全部訓練完畢共花費214s。訓練結果如圖4所示。PLSR模型效果最佳。而在兩種草地類型下的N和K含量以及溫性草原下P含量構建的RF模型均表現最優。

圖4 HR-CNN模型訓練損失曲線

4.2 實驗分析

為驗證本模型有效性，在采用相同數據預處理及實驗環境基礎上，采用ResNet34、GoogLeNet兩種常用高光譜圖像分類算法進行對比。同時，為驗證深度可分卷積與常規卷積所構模型區別，將深度可分卷積改為常規卷積，構建常規2D卷積神經網絡模型(2D-CNN)進行對比驗證，因二者具有相同模型結構，故在生成模型所占內存及運行參數量方面具有較大差異。

本文在不同模型分類混淆矩陣（如圖5所示）基礎上，通過對其進行運算，分別以整體分類精度(Overallaccuracy,OA)(公式)、F1得分(F1-score,F1)(公式)及Kappa系數(公式)作為評價指標。

其中：TP為預測正樣本，TN為預測負樣本，FP為真實正樣本，FN為真實負樣本，N為樣本總數量。

為便于模型間對比驗證，本文對Kappa系數進行放大100倍處理。各模型F1分類精度如表3所示，OA及Kappa系數分類精度對比結果如圖6所示。此外，在相同數據預處理及圖像分類預測結果條件下，模型構建優劣與訓練時間、模型運行參數量所占內存有關。訓練時間越短且內存占有量越小，構建模型越好。四類模型訓練時間及內存占有量如表4所示。

圖5 混淆矩陣對比

由表3可知，四類模型對高光譜圖像中具有較大數據量的植被、裸土均具有較高識別精度。但對于小樣本數據集（地墊、小旗、鼠洞），四類模型圖像識別精度具有較大差異，但HR-CNN模型各地物分類精度均表現優異。其中，ResNet34、GoogLeNet網絡模型對于鼠洞識別精度較低，相較HR-CNN模型精度分別降低13.40%、18.10%。

表3 F1分類精度對比表

由圖6可知，(c)2D-CNN與(d)HR-CNN模型OA與Kappa值相差不大，但(a)ResNet34相較HR-CNN模型OA與Kappa值分別降低1.91%、0.97，(b)GoogLeNet相較HR-CNN模型OA與Kappa值分別降低0.47%、0.93。

圖6 OA及Kappa系數分類精度

實驗中，四類模型在訓練階段模型損失逐漸降低并趨于平穩（損失平穩即為模型構建完畢）。各模型構建完成時間（表4），其中HR-CNN模型訓練時間最短為144s，僅占ResNet34、GoogLeNet、2D-CNN模型訓練時間的34.12%、34.29%、86.23%。HR-CNN模型內存占有量最小為19.5MB，相比ResNet34、GoogLeNet、2D-CNN分別降低92.11%、79.21%、43.64%。HR-CNN模型運行參數所占內存最小為6.06MB，較ResNet34、GoogLeNet、2D-CNN分別降低92.57%、79.03%、46.18%。

表4模型訓練時間及內存占有量對比表

為驗證模型分類效果，對高光譜圖像進行可視化分類研究。由圖7可看出，HR-CNN模型對草原地物分類表現良好。

圖7高光譜圖像原圖與模型分類預測圖

結論

本文旨在通過構建無人機高光譜遙感系統，并結合深度學習中的卷積神經網絡，提出一種高分辨率圖像卷積神經網絡（HR-CNN）用于荒漠草原地表微斑塊分類識別研究。實驗證明該網絡模型在此類任務中具有較高的識別精度和計算效率。該模型通過引入殘差結構，在高光譜圖像中使用不同大小的卷積核提取各地物影像特征。轉置卷積用于重組特征，以生成高分辨率的特征圖，并通過通道拼接的方式減少特征丟失。此外，引入了深度可分離卷積以減少網絡參數量，提升模型性能。與2D-CNN、ResNet34和GoogLeNet網絡模型相比，HR-CNN在訓練時間和模型內存占用方面有顯著改進。本模型的輕量化設計對于荒漠草原地表微斑塊的精細分類研究具有重要意義，為進一步深入研究荒漠草原地物分類、識別和反演提供了新思路。

推薦：

便攜式高光譜成像系統iSpecHyper-VS1000

專門用于公安刑偵、物證鑒定、醫學醫療、精準農業、礦物地質勘探等領域的最新產品，主要優勢具有體積小、幀率高、高光譜分辨率高、高像質等性價比特點采用了透射光柵內推掃原理高光譜成像，系統集成高性能數據采集與分析處理系統，高速USB3.0接口傳輸，全靶面高成像質量光學設計，物鏡接口為標準C-Mount，可根據用戶需求更換物鏡。

審核編輯 黃宇