李興亮

（黔西南州水利電力勘測設計院，貴州 興義 562400）

1 引言

河流是指由一定區(qū)域內地表水和地下水補給，經?；蜷g歇地沿著狹長凹地流動的水流，河流湖泊作為陸地水圈的重要組成部分，在全球生態(tài)環(huán)境中扮演著重要角色。河流是地球上水文循環(huán)的重要路徑，參與自然界的水分循環(huán)，是泥沙、鹽類和化學元素等進入湖泊、海洋的通道。河流位置、水面寬度、水位及水儲量等信息的變化是區(qū)域甚至全球氣候變化的重要指示因子，是氣候變化、生態(tài)環(huán)境研究的主要對象之一[1-3]。

河流網作為流域地形特征的基本骨架在地理信息系統（GIS）中占據重要地位，也發(fā)揮著重要作用。河流在不同時期、不同環(huán)境中的狀態(tài)各不相同。遙感技術具有時效性強、數據獲取簡單、大面積同步觀測等優(yōu)勢，是目前河流湖泊面積、水位等研究的主要數據源，國內外學者基于此開展了許多有意義的工作[4]。

2 技術發(fā)展現狀

傳統的河流提取為人工數字化[5-6]，需要大量的野外調研與人力物力投入，時效性低。遙感憑借大面積同步觀測的優(yōu)勢，是目前河流湖泊等研究的主要數據源。水體在遙感影像中具有顯著的光譜特征，在可見光波段有較強的反射，在近紅外波段有較強的吸收，這種特征主要依賴所使用的光譜波段，與傳感器無關。

從遙感影像中提取湖泊邊界的單一波段閾值法采用了水體在近紅外波段的吸收特征[7]，如利用Landsat8 影像近紅外波段進行影像分割[8]，得出水體區(qū)域；波段比值法和水體指數法考慮了水體在可見光波段和近紅外波段的反射、吸收特性，以更好地分開水體與其他地物[9]，如徐涵秋等利用改進的歸一化水體指數進行水域提取，較容易提取出水體及剔除建筑陰影部分；影像分類方法也可較好提取水體信息，除基于水體光譜特征的方法，基于圖像紋理特征的方法也用在高空間分辨率全色遙感影像水體邊界的提取中。

目前高空間分辨率遙感影像已廣泛應用于河流的提取與變化分析中，如Landsat、SPOT，以及中低分辨率影像MODIS 等數據，本文基于近年興起的哨兵2多光譜影像（Sentinel-2）進行河流的遙感提取，以達到快速準確獲取河流信息的目的。

3 數據及解譯方法介紹

3.1 哨兵數據介紹

哨兵2 是高分辨率多光譜成像衛(wèi)星，攜帶一枚多光譜成像儀（MSI），分為2A 和2B 兩顆衛(wèi)星。哨兵2影像具有重訪周期短（5 ～10 天）、覆蓋范圍廣、波段多、分辨率高等優(yōu)點[10]，用于陸地監(jiān)測，可提供植被、土壤和水覆蓋、內陸水路及海岸區(qū)域等圖像。哨兵2 衛(wèi)星攜帶的高分辨率多光譜成像裝置主要用于土地環(huán)境監(jiān)測，可提供土壤覆蓋狀況、陸地植被生長、內河和沿海區(qū)域環(huán)境等信息，不僅為改善農林業(yè)種植、預測糧食產量提供宏觀信息，還可用于監(jiān)測洪水、火山噴發(fā)、山體滑坡等自然災害。哨兵2A 和2B 同時進入運行狀態(tài)后，每5 天可完成一次對地球赤道地區(qū)的完整成像，對于緯度較高的歐洲地區(qū)，這一周期僅需3 天。哨兵2 多光譜衛(wèi)星具有13 個波段，具體波段信息如表1 所示。

哨兵影像數據處理主要包括輻射定標、大氣校正、格式轉換等。哨兵數據的處理軟件主要為Sen2cor、SNAP、ENVI、ArcMap。

表1 哨兵2波段信息

3.2 數據預處理

下載哨兵2 多光譜影像2 級反射率產品L2A 數據，對原始數據進行解壓，然后在SNAP 軟件中的Sen2cor模塊進行影像波段重采樣，重采樣至10 米，輸出后導入ENVI 軟件，采用layer stacking 功能，將重采樣至10 米的單波段影像進行重組，最終得到哨兵2 多光譜影像，效果如圖1 所示。

圖1 哨兵2影像預處理效果

3.3 解譯方法

采用哨兵影像，基于分類與反演手段對天然河流及人工河流進行提取，馬別河流域數字線劃圖如圖2所示。采用遙感分類與反演方法進行河流信息提取，能避免測繪手段單一、容易產生錯誤識別、漏識別等缺點，具有融合兩種手段對水域信息進行計算等優(yōu)勢。

（1）遙感分類

遙感分類是利用計算機對遙感圖像中的各類地物進行光譜信息和空間信息分析，計算機自動選擇特征，然后將像元按照特定的規(guī)則或算法劃分類別，最終實現遙感圖像的分類。常用的分類方法為監(jiān)督分類與非監(jiān)督分類，相對于非監(jiān)督分類的聚類分析，監(jiān)督分類又稱訓練場地法，是以建立統計識別函數為基礎、依據典型樣本訓練方法進行分類的方法，根據提供的研究區(qū)樣本，選擇特征參數，求出特征參數作為分類規(guī)則，建立判別函數對影像進行分類。

本研究將影像影本劃分為水域、植被覆蓋區(qū)域、其他區(qū)域三類，并勾畫樣本。勾畫樣本時應在地物清晰位置，將樣本勾畫在地物中央，樣本內沒有混合像元。研究采用航天宏圖PIE 軟件面向對象多尺度分割進行水域提取，首先對影像進行多尺度分割，將相似像元分割成一類 ；在分割基礎上進行樣本勾選，勾選出水域、植被覆蓋區(qū)域、其他區(qū)域三類，然后將樣本帶到分類器中。研究采用監(jiān)督分類中精度最高的支持向量機分類法，對分類結果進行小斑塊剔除等后處理，最終得到高精度的河流信息。

（2）遙感反演

針對水域信息的提取，本文選取研究較多、應用較廣泛的歸一化差分水指數NDWI（Normalized Difference Water Index）進行水域提取，該指數是用遙感影像的特定波段進行歸一化差值處理，以凸顯影像中的水體信息。該指數是基于綠波段與近紅外波段的歸一化比值指數，指數在河流提取上也具有較好的應用。歸一化差分水體指數（NDWI）計算公式如下：

圖2 馬別河流域數字線劃圖

4 結果分析

4.1 分類結果分析

采用基于真彩色合成的哨兵2 多光譜影像勾畫的樣本，對影像進行分類，經過多次調整樣本數量以及樣本位置，最終的分類結果如圖3（a）所示，圖中紅色部分為水域、青色部分為植被覆蓋區(qū)域、藍色為其他區(qū)域。由圖3（a）可知，基于監(jiān)督分類的方法進行水域提取，效果較好，影像中的河流信息可完全分類識別出來，且未出現斷流現象，同時影像中較小且較多的湖泊水塘，也有很好的識別與分類。

對分類結果進行后處理，剔除分類過程中由混合像元引起的小斑塊，降低分類誤差，然后對后處理后的分類結果進行矢量化，矢量化后的結果如圖3（b）所示。對矢量化后的河流進行渲染，由圖可知河流信息整體被分類出來，并且河流信息符合自然分布規(guī)律。

圖3 遙感分類結果

4.2 反演結果分析

采用歸一化水體指數進行水體的反演與提取，效果如圖4 所示，基于反演手段的水域信息提取，相對于分類手段，以反演手段直接計算每個像元的信息，通過波段的線性與非線性組合，將水體信息在反演結果中放大并反映出來，相較于分類手段，反演手段較好地避免了混合像元等錯誤信息的引導。圖4（a）中藍色部分為反演顯示的水域信息，圖4（b）為通過反演手段與閾值提取的水體邊界信息。由圖可知，水體信息的反演整理效果很好，河流呈現連續(xù)分布，并且較小的湖泊水塘在反演的NWDI 圖上也有很好顯示。

圖4 水體反演效果

4.3 精度驗證

研究采用基于預處理后的影像采集地物點位，用這些點位驗證遙感分類與遙感反演的結果，結果如下：

（1）遙感分類驗證

針對遙感分類的驗證，基于預處理后的影像采集水域、植被覆蓋區(qū)域、其他區(qū)域三種類型的地物點位，每個地物類別采集10 個點位，用這些點位對分類結果進行驗證。驗證結果如表2 所示。

表2 遙感分類精度驗證

由表2 可知，遙感分類方法提取水域信息的精度驗證結果顯示，10 個點位均分類正確，即水體信息提取精度達到100%；植被覆蓋地區(qū)10 個點位驗證結果中1 個點位分類錯誤，驗證精度為90%；其他地區(qū)10 個點位驗證結果中1 個點位分類錯誤，驗證精度為90%；植被覆蓋區(qū)域與其他區(qū)域存在交互混合、過渡較明顯的形態(tài)分布，容易造成錯分現象，即植被覆蓋地區(qū)與其他地區(qū)識別顛倒，在遙感上稱為混合像元引起的分類錯誤。

（2）遙感反演驗證

采用散點圖的形式描述反演精度，由圖5 可知，遙感反演方法提取水域信息的精度驗證結果均在y=x線上，即每個驗證樣本均落在反演結果中，即反演方法提取水域的精度達到了100%。

圖5 遙感反演精度驗證

5 結論

本研究基于哨兵2 多光譜影像提取河流信息，提出了分類與反演相結合的方法對河流信息進行高精度提取，基于研究得出以下結論：

（1）哨兵2 多光譜影像具有重訪周期短、覆蓋范圍廣、影像波段信息豐富等優(yōu)點，為提取河流信息提供重要的數據支撐。

（2）基于遙感分類手段的河流提取精度很高，河流呈連續(xù)狀態(tài)分布，研究區(qū)水域的分類精度達100%，且植被區(qū)域與其他區(qū)域分類精度也較好，均為90%。

（3）基于遙感反演手段的河流提取同樣具有很高精度，反演手段從光譜信息中準確識別出水域信息并反映出來，河流信息符合自然規(guī)律且未出現斷流現象。

基于哨兵2 影像的河流信息提取雖然有較好的應用性，但仍有較多不足之處，如水域污染較嚴重的區(qū)域，光譜信息相差較大，遙感手段的分類與反演容易出現錯分漏分、錯誤識別的現象；同時遙感影像的像元分辨率的大小直接限制了遙感分類與反演能識別出的最小河流信息。因此，基于大力發(fā)展國產衛(wèi)星傳感器、提高衛(wèi)星影像分辨率的背景，今后在河流湖泊等水域信息的提取理論方面，仍有很長的科研道路要走。