董哲，王凌，朱西存，馮文斌，張美霞

(山東農(nóng)業(yè)大學(xué) 資源與環(huán)境學(xué)院，山東 泰安 271000)

0 引言

山區(qū)地表水體影響著山區(qū)的生態(tài)系統(tǒng)和人類的生產(chǎn)勞動，快速、大范圍地獲取山區(qū)水體的分布情況對于山區(qū)的災(zāi)害評估、用水安全保障、資源分配有著重要的意義[1-2]。山區(qū)水體分布零散，面積相對較小，且易與山地陰影相混淆，在水體的提取上存在一定的難度。遙感作為采集地球數(shù)據(jù)以及變化信息的手段，已經(jīng)廣泛應(yīng)用于國民經(jīng)濟和社會發(fā)展的各個領(lǐng)域[3]，也成為獲取水體真實范圍的一個可行手段[4]。中低分辨率遙感數(shù)據(jù)很難用于高精度山區(qū)地表水體的提取，而高分二號(GF-2)作為我國第一顆分辨率為亞米級的光學(xué)遙感衛(wèi)星[5]，為獲取更真實、更豐富的地表水體信息提供了可能。

針對遙感影像的水體提取，國內(nèi)外學(xué)者進行了大量的研究。目前遙感影像水體提取最常用的方法是水體指數(shù)法[6]，通過分析水體的光譜特征選取相關(guān)波段構(gòu)建水體指數(shù)，并給定相關(guān)閾值來進行水體提取[7]。Mcfeeters[8]根據(jù)植被和水體的反射不同構(gòu)建歸一化差異水體指數(shù)(normalized difference water index，NDWI)，該指數(shù)突出了水體信息，可以很好地提取水體，但會造成陰影與水體的混淆，形成的噪聲降低了水體提取的精度。馬吉晶等[9]通過濁水指數(shù)對蘭州市安寧區(qū)濕地公園進行了渾濁水體提取，證明濁水指數(shù)對區(qū)分渾濁水體的效果較好。鄭逢斌等[10]利用決策樹法，結(jié)合最大似然法、單波段閾值法和城區(qū)水體陰影指數(shù)(urban shadow water index，USWI)提取了安陽縣彰武水庫的水體，結(jié)果表明該方法的水體提取精度比NDWI高。

面向?qū)ο蠓ㄒ彩且环N常見的遙感信息提取方法，它能夠利用更加豐富的地物信息，同時還具有良好的抗噪能力。相較于水體指數(shù)法，面向?qū)ο蠓梢詼p少“椒鹽現(xiàn)象”的出現(xiàn)[11]，能更好地利用高分辨率影像豐富的信息，使得面向?qū)ο蠓椒ㄌ崛∷w信息更加符合實際，很好地彌補水體指數(shù)法的缺陷。付勇勇等[12]基于面向?qū)ο蠓?，采用分離閾值法建立提取規(guī)則，實現(xiàn)杭嘉湖水網(wǎng)平原水體的提取，總體精度達到97.3%。王俊海等[13]通過統(tǒng)計簡單的光譜信息、紋理特征，利用面向?qū)ο蠓ㄌ崛P州市某區(qū)域的水體，精度為89.55%，很好地區(qū)分了水體和建筑物陰影，表明該方法在城區(qū)水體提取方面的有效性。

上述水體提取實例均基于GF-2影像。目前的GF-2水體提取研究多數(shù)針對平坦區(qū)域，而針對地勢起伏較大的山區(qū)進行的水體提取卻較少。面向?qū)ο蠓ㄌ崛∷w需要進行復(fù)雜的計算才能得到最優(yōu)的分割參數(shù)，用以保證水體邊界提取的正確性[14]。受限于GF-2僅有藍、綠、紅、近紅外4個波段，無法構(gòu)建出需要中波紅外波段的改進歸一化差異水體指數(shù)(modified normalized difference water index，MNDWI)、水體指數(shù)2015(water index 2015，WI2015),以及自動水體提取指數(shù)(automated water extraction index，AWEI)等提取效果較好的水體指數(shù)。

基于上述考慮，本文在NDWI指數(shù)的基礎(chǔ)上，提出差異水體光譜模型(difference water spectral model，DWSM)，通過與面向?qū)ο蠓ńY(jié)合，旨在為大范圍山區(qū)水體的提取提供一種更加準確的方法，簡化復(fù)雜計算。

1 研究區(qū)概況與數(shù)據(jù)來源

1.1 研究區(qū)概況

研究區(qū)為山東省泰安市泰山的主體區(qū)(117°0′E～117°10′E，36°12′N～36°20′N)，最高海拔1 532.7 m，高差達1 200 m以上，地勢起伏較大，屬于暖溫帶季風氣候，具有明顯的氣候垂直差異，山腳年均氣溫13.2 ℃，年均降水量678.5 mm；山頂年均氣溫5.9 ℃，年均降水量1 031.3 mm。水體分布較離散，包括橫嶺后、曹家莊、藥鄉(xiāng)附近面積較大的水庫及一些零散分布的坑塘。地表景觀主要包括大范圍的林地、分布在山腳四周的耕地及裸地、零散分布的水體及建筑用地。

1.2 數(shù)據(jù)來源

選用2017年11月8日分辨率為4 m的GF-2多光譜影像作為實驗影像。同時，為驗證本文方法的普適性，以2018年6月23日的GF-2多光譜影像作為驗證影像，兩期影像的云覆蓋率均為0%。利用Erdas Image對實驗影像和驗證影像進行預(yù)處理，包括正射校正、影像匹配、圖像裁剪，其中正射校正總體誤差小于0.3個像元。

2 研究方法

圖1 水體提取技術(shù)路線

本文技術(shù)路線如圖1所示，采用差異水體光譜模型結(jié)合面向?qū)ο蠓ㄗ鳛樗w提取方法，將陰影水體指數(shù)(shade water index，SWI)決策樹法[15]、改進的陰影水體指數(shù)(modified shade water index，MSWI)決策樹法[16]和支持向量機法(support vector machine，SVM)作為此方法的對照方法。

2.1 DWSM結(jié)合面向?qū)ο蠓?/h3>

NDWI利用水體在近紅外波段反射弱、在綠光波段反射強的特點，基于上述兩波段構(gòu)建歸一化差異指數(shù)來突出水體信息，表達如式(1)所示。

(1)

式中：Green為綠光波段DN值；NIR為近紅外波段DN值。NDWI可以抑制植被和土壤信息，但會混淆水體和陰影[17]，這種情況在山區(qū)尤為嚴重，地形起伏較大的山區(qū)不僅存在大面積的山體陰影，也存在零散分布的建筑物陰影，使得山區(qū)水體提取難度變得更大。為更加準確地提取水體，本文基于NDWI指數(shù)提取水體結(jié)果掩膜后的實驗影像，選取24 978個水體像元和27 956個陰影像元進行對比分析，其中水體主要包括研究區(qū)的水庫和坑塘，陰影包括建筑物陰影和山體陰影，分別統(tǒng)計水體和陰影DN值在4個波段像元的頻率，如圖2所示。從圖2可以看出，相較于綠光、紅光和近紅外波段，水體-陰影在藍光波段的DN值重疊范圍最小，這表明在選取合適閾值的前提下，水體-陰影在藍光波段有較好的區(qū)分度。為了進一步擴大水體-陰影的區(qū)分度，構(gòu)建差異水體光譜模型DWSM，表達如式(2)所示。

(2)

式中：Blue、Green、NIR分別為藍光波段、綠光波段、近紅外波段的DN值。

利用與圖2相同的方式在DWSM影像中選取水體像元和陰影像元，構(gòu)建出的水體-陰影對比頻率如圖3所示。從圖3可以看出，當分割閾值為119時，水體和陰影具有良好的區(qū)分度，但仍有1.14%(319個像元)的陰影被錯分為水體，將繼續(xù)利用面向?qū)ο蠓ㄌ蕹@部分陰影噪聲。

圖3 水體-陰影DWSM值頻率圖

面向?qū)ο笏w提取方法包括影像分割和水體提取規(guī)則建立兩個步驟。影像分割最終目的是滿足異質(zhì)性標準，規(guī)則集的建立主要依靠光譜特征、形狀特征以及紋理特征3個方面[18]。

圖2 4個波段水體-陰影DN值頻率圖

1)面向?qū)ο蟮亩喑叨确指?。根?jù)研究區(qū)水體形狀不規(guī)則、大小不一的特點，選取多尺度分割法[19]作為影像分割方法，用以提取水體分布信息。分割參數(shù)會影響圖像分割的質(zhì)量，進而直接影響到面向?qū)ο蟮姆诸惥取１疚膶WSM水體提取結(jié)果作為第5波段與原始GF-2 4個波段疊加，選取DWSM水體提取結(jié)果范圍內(nèi)的影像用于分割，每個波段的權(quán)重均為1。通常情況下，影像的最優(yōu)分割尺度難以把握[20]，需要進行繁雜的計算。而本文以DWSM所提取的水體區(qū)域作為影像分割范圍，在忽略背景0值的情況下，水體是獨立存在的，不與其他地物相鄰，要區(qū)分的類別只有水體和零散陰影兩類，且要剔除的陰影較NDWI法少，利用簡單的“試錯法”可以很快地找到合適的分割尺度。經(jīng)過實驗，當尺度參數(shù)為50，形狀因子為0.1，緊實度因子為0.3時，分割效果較好，不論是水庫還是坑塘，分割情況都符合實際，沒有出現(xiàn)水體與其相近的陰影區(qū)域分割到同一對象的情況，分割后的對象具有高度均質(zhì)性。分割后的影像及細節(jié)如圖4所示。

圖4 尺度為50的分割圖及細節(jié)(R：NIR；G：Red；B：Green)

2)規(guī)則建立。綜合考慮研究區(qū)內(nèi)水體分布較為零散、面積相對較小且輪廓不規(guī)則的特點，同時兼顧水體在不同波段的特性，經(jīng)多次實驗，本文選用水體強反射的藍波段、水體強吸收的近紅外波段和DWSM水體提取波段分別構(gòu)建比值率(ratio)、均值(mean)、亮度值(brightness)3種光譜特征值作為水體提取的規(guī)則。

在分割后的影像中盡可能均勻地選取水體和陰影訓(xùn)練樣本，創(chuàng)立訓(xùn)練樣本集，從而得到水體提取規(guī)則值域。實驗影像水體提取規(guī)則范圍如表1所示。

表1 實驗影像規(guī)則值域

2.2 對照方法

1)決策樹法。SWI決策樹法、MSWI決策樹法是基于GF-1影像提出的水體提取算法，這兩種決策樹法都是先利用水體在近紅外波段的強吸收性的特點，通過設(shè)定閾值在近紅外波段實現(xiàn)水體粗提取，再利用SWI、MSWI分別實現(xiàn)水體的精確提取。SWI、MSWI的計算如式(3)和式(4)所示。

SWI=Blue+Green-NIR

(3)

MSWI=(Blue-NIR)/NIR

(4)

式中：Blue、Green、NIR分別為藍光、綠光、近紅外波段的DN值。

2)支持向量機法。支持向量機法是建立在統(tǒng)計學(xué)習(xí)理論基礎(chǔ)上的機器學(xué)習(xí)方法[21]。SVM分類法的精度主要取決于核函數(shù)類型的選擇和核函數(shù)參數(shù)的設(shè)置[22]。研究表明，在選取合適的核函數(shù)類型和參數(shù)的前提下，即使選取樣本數(shù)量較少，SVM分類也能取得較高的精度[23]。本文借鑒段秋亞等[24]的研究，選用徑向基核函數(shù)(radial basis function，RBF)和對應(yīng)默認的參數(shù)(gamma為0.25；penalty為100；pyramid為0)進行GF-2影像水體信息提取。

2.3 精度評價

選取400個水體點位和400個非水體點位，利用混淆矩陣計算總體分類精度(overall accuracy，OA)、Kappa系數(shù)作為各種水體提取方式的精度評價指標。

3 結(jié)果分析

3.1 NDWI與DWSM提取結(jié)果對比

NDWI與DWSM水體提取結(jié)果對比如圖5所示。NDWI法混淆了水體和陰影，形成大量的噪聲；DWSM法可以很大程度上避免這種情況，在高低起伏的山區(qū)提取水體的效果更好，在基本剔除陰影的同時，保留完整的水體信息，但也存在一些陰影像元錯分，使得影像存在椒鹽現(xiàn)象，原因在于這些陰影和水體在4個波段的光譜特征極為相似，單純地利用波段之間的數(shù)學(xué)運算難以找出陰影和水體之間的區(qū)別，因而結(jié)合面向?qū)ο蠓ㄊ潜匾摹?/p>

圖5 NDWI法與DWSM法水體提取結(jié)果

3.2 本文方法與對照方法的結(jié)果對比

實驗影像和驗證影像水體提取結(jié)果分別如圖6和圖7所示，紅圈是錯分的山體陰影，藍圈是錯分的建筑物陰影。其中，圖7結(jié)果只展示精度較高的SVM法和DWSM結(jié)合面向?qū)ο蠓ā２煌椒ǖ乃w信息提取精度評價見表2。

表2 水體提取精度評價結(jié)果

由表2可知，無論是在實驗影像還是驗證影像中，4種水體提取方法的總體精度和Kappa系數(shù)由高到低均為DWSM結(jié)合面向?qū)ο蠓?、SVM分類法、SWI決策樹法和MSWI決策樹法。DWSM結(jié)合面向?qū)ο蠓ǖ腒appa系數(shù)和總體分類精度均高于其他算法，且水體提取精度較為穩(wěn)定。實驗影像和驗證影像的總體分類精度均大于95%，Kappa系數(shù)分別為0.953 3和0.933 4。DWSM結(jié)合面向?qū)ο蠓ㄔ趦煞跋裰械穆┓终`差分別為1.59%和1.41%，錯分誤差均為0%，總體誤差最低。該方法不僅可以很好地區(qū)分水體和其他易混地物，而且也能較為精準地提取細微水體，并未出現(xiàn)嚴重的水體漏分情況，在保證水體信息完整性的同時，減少椒鹽現(xiàn)象的發(fā)生(圖6、圖7)。

圖6 實驗GF-2影像水體提取結(jié)果及細節(jié)對比

圖7 驗證GF-2影像水體提取結(jié)果及細節(jié)對比

SVM法提取的水體輪廓較為符合實際情況，很好地保持了水體信息，在4種水體提取方法中的漏分誤差最低，但是伴有建筑物陰影以及部分山體陰影的誤提，兩幅影像中的總體精度較DWSM結(jié)合面向?qū)ο蠓ㄏ陆盗?.27%和7.01%。

MSWI決策樹法的精度在4種方法中最低，Kappa系數(shù)分別為0.615 6和0.616 1。與NDWI相似。MSWI決策樹法精度低主要在于光譜特征相似的山體陰影、植被陰影被誤分為水體。SWI決策樹法較MSWI決策樹法精度略有提高，總體分類精度分別提高了9.30%和6.61%，但依然存在著大量陰影錯分，同時漏分誤差也較高。無論是SWI決策樹法還是MSWI決策樹法，應(yīng)用于本研究區(qū)的水體提取效果均不理想，不能作為水體提取的最終結(jié)果。

4 結(jié)束語

基于GF-2影像，本文利用DWSM與面向?qū)ο蠓ㄏ嘟Y(jié)合的方法，在山區(qū)實現(xiàn)了水體信息的高精度提取，并與SWI決策樹法、MSWI決策樹法和SVM分類法進行對比，結(jié)論如下。

1)在GF-2影像中，DWSM能夠很好地減少由地形起伏造成的陰影噪聲，相較于NDWI，應(yīng)用于山區(qū)水體提取效果較好。

2)DWSM結(jié)合面向?qū)ο蠓ǖ奶崛【雀哂趯φ辗椒?，在兩個時相GF-2影像上都可以取得較高的精度；其次是SVM分類法，SVM法可以較為精確地提取水體輪廓，但會受到建筑物陰影的影響，降低了水體提取精度；而SWI決策樹法精度比前兩種方法低，存在大量陰影噪聲。

3)在波段較少的影像中，單純利用水體指數(shù)法提取山區(qū)水體的精度不高，如果要獲得更好的水體提取效果，水體指數(shù)與面向?qū)ο蠓ㄏ嘟Y(jié)合是一種有效的方法。

本文所提出的水體提取方法在山區(qū)取得了較高的精度，后續(xù)將進一步驗證其在平原地區(qū)的應(yīng)用。獲取合適閾值是光譜模型準確提取水體的前提，如何快速準確地在影像中確定水體閾值將是下一步的方向。