王永麗++宋希磊

摘要：以丹江口水庫流域水體為研究對象，根據(jù)其Landsat ETM+影像數(shù)據(jù)，利用ENVI4.7對其進行了處理分析。根據(jù)丹江口水庫庫區(qū)的水體特征，對單波段閾值法、譜間關(guān)系法、改進歸一化指數(shù)法等3種方法提取水體的優(yōu)缺點進行了比較和分析，最后選取適合大面積水體信息提取的譜間關(guān)系法進行提取。實踐表明：此方法不僅避免了閾值確定的麻煩，而且解決了陰影對水體提取的限制影響，實現(xiàn)了丹江口水庫大面積水體的精確提取。

關(guān)鍵詞：水體提??；陰影；單波段閾值法；譜間關(guān)系法；丹江口水庫

中圖分類號：TP722.6

文獻標識碼：A 文章編號：16749944（2017）10020803

1 引言

進入21世紀以來，數(shù)字流域的建設已逐漸成為我國水利現(xiàn)代化的重要組成部分[1]。遙感影像數(shù)據(jù)是自動提取水體的基礎，可大面積周期性實時地獲取水體信息數(shù)據(jù)，通過影像數(shù)據(jù)就可對研究區(qū)進行判讀和分析，大大提高了工作效率；在近20年的衛(wèi)星遙感發(fā)展過程中，利用衛(wèi)星遙感數(shù)據(jù)對資源和環(huán)境的調(diào)查研究得到了廣泛的應用，特別是在對水資源的研究。衛(wèi)星遙感影像記載了地物對電磁波的反射信息及地物本身的熱輻射信息，而由于不同的地物其自身的結(jié)構(gòu)屬性的不同，導致其對電磁波的反射及本身的熱輻射存在著差異[2]。與其他地物不同，水體在整個光譜范圍內(nèi)都呈現(xiàn)出較弱的反射率[3]；恰是因為這種特征，才能夠提取水體等其他地物的信息。

對丹江口水庫流域信息的提取研究最重要的是區(qū)別出水體和非水體的過程，即提取水體與其他地物的分界線，并對其水體面積進行統(tǒng)計；單波段閾值法和多波段法是現(xiàn)在通用的兩種水體提取方法，根據(jù)對研究的需求利用不同的提取方法。

單波段法主要選取對水陸界反應較好的Landsat-7的ETM+的第四或第五波段，即近紅外和中紅外波段，利用閾值法對水體進行提取[3]。多波段法根據(jù)不同波段的綜合特點進行影像的疊加和波段運算來實現(xiàn)水體信息的提取，由此可分為譜間關(guān)系法和水體指數(shù)法。以丹江口水庫的水體信息為例，選取2014年3月份LandsatETM+影像為數(shù)據(jù)源，用上述兩種方法對丹江口水庫的水體信息舉行提取，并對其提取成果展開分析，以選取適合該研究區(qū)的水體提取方法。

2 研究區(qū)概況和方法

2.1 研究區(qū)概況

丹江口水庫位于北緯32°36′～33°48′，東經(jīng)110°59′～111°49′之間；丹江口橫跨湖北和河南兩省，其主要流域位于河南省淅川縣和湖北省丹江口市，是國家南水北調(diào)工程的水源地，丹江口水庫大壩建成以后，蓄水以后庫容可達339.1億m3，丹江口水庫水域面積將達600 km2，其水域面積廣闊且較為連續(xù)，丹江口水庫是以供水、防洪、發(fā)電為主，兼具灌溉、養(yǎng)殖等功能的大型人工水庫。龐大的水域面積使其形成了自身的獨立的生態(tài)系統(tǒng)，使其周邊氣候適宜，良好的空氣質(zhì)量，達到飲用標準的水質(zhì)，周邊的地區(qū)以丹江口水庫獨特的性質(zhì)，建立了生態(tài)旅游、種植、養(yǎng)殖等多種經(jīng)濟社會功能；其是國家重點大型水利工程南水北調(diào)的水源地；該工程貫通后，將把丹江口水庫的水輸送到華北地區(qū)以供飲用，滿足用水需求。

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

為了盡量減少天氣對遙感影像數(shù)據(jù)的影響，根據(jù)研究區(qū)天氣預報情況及當天的云量分布情況，選用了2014-3-17條帶為WRS2 125/37的影像數(shù)據(jù)，該研究區(qū)的影像不存在地形陰影和云體陰影的影響。

2.3 研究方法

2.3.1 水體的影像特征

水體自身的屬性決定了水體的光譜特征，但不同的水態(tài)也會具有不同的光譜特征，使其與真實水體存在著一定的差異；相對于其他地物而言，水體在整個光譜范圍內(nèi)對入射的太陽能量具有強吸收性，在整個光譜范圍內(nèi)遙感都呈現(xiàn)出較弱的反射率，特別是在近紅外、中紅外及短紅外部分，水體幾乎吸收全部的入射能量，因此，水體在這些波段的反射能量很少[3]。而土壤、建筑物和植被等地物正好相反，具有較高的反射率，這使得水體與其他地物有明顯的區(qū)別。圖1所示為不同地物的光譜曲線。

2.3.2 單波段閾值法

單波段閾值法是比較簡單的一種水體自動提取方法，可以快速的提取水體，其跟據(jù)水體在整個光譜范圍內(nèi)具有較強的吸收能量，水體表現(xiàn)出較弱的反射率，特別是水體在近紅外和中紅外波段反射率幾乎為零，而非水體（如：植被、土壤、建筑物等地物）正好相反，具有較高反射率，這使得水體與其他地物有明顯的差異。利用水體信息和其他地物信息在遙感影像上的灰度值不同來確定閾值。閾值的確定很關(guān)鍵，必須經(jīng)過反復試驗得到，其大小都會影響水體提取的精度。

2.3.3 多波段法

（1）歸一化差異水體指數(shù)（NDWI）法。此研究方法是比較特殊的比值運算方法，在基本的比值運算中，如果分母中的波段的值比較?。ㄌ貏e是小于1的情況下），那么，比值的結(jié)果將夸大差異，使用歸一化指數(shù)可以避免這個問題，歸一化指數(shù)中像素值的變化是有約束的[4]。水體歸一化比值運算的原理就是在多光譜波段內(nèi)，把水體最強的反射率波段放在分子上，最弱的放在分母，通過兩者的相比運算從而進一步擴大兩者的差異，使水體的亮度達到最大值，而其他非水體地物則基本上都被抑制，從而突出研究水體信息的目的，即抑制建筑物、植物和土壤的信息提取水體[5]?？紤]水體的波譜特征是在綠光波段具有較強的透射性，在近紅外波段吸收性較強，可以更好地反應水體信息，因此Mcfeeters在1996年提出了歸一化差異水體指數(shù)NDWI[6]。

（2）譜間關(guān)系法。此研究方法利用水體具有獨特的譜間關(guān)系特征：（b2+b3）>（b4+b5），這一構(gòu)建模型來進行水體提取能將水體與陰影區(qū)分開來，特別適合山區(qū)水體的提取[7]。由圖1可知，只有水體的2、3波段相加大于其4、5波段之和，按照這種特性對水體進行提取。

3 研究結(jié)果

3.1 數(shù)據(jù)處理與分析

3.1.1 單波段閾值法

利用單波段閾值法對丹江口水庫進行水體信息的提取難點在于最佳閾值的確定，閾值大小的確定關(guān)系到了水體提取的精度；由于丹江口水庫庫區(qū)水陸界限清晰，范圍比較廣，因此為了避免選取閾值的麻煩，并不采用此方法進行水體提取。

3.1.2 多波段法

（1） 歸一化差異水體指數(shù)（NDWI）法。此方法最大程度地改善了誤提植被、建筑物和居民地等信息提取的問題，使水體信息得到了更好體現(xiàn)，提高了水體提取的準確度。但此指數(shù)研究方法在提取水體的時候忽略了山體、裸地以及山體陰影等信息對水體提取時候的影響，對提取結(jié)果的精確性造成了限制。

（2） 譜間關(guān)系法。依據(jù)丹江口水庫各波段的影像特征，選取了第2、3、4、5波段進行了波段運算，輸入運算公式（b2+b3）>（b4+b5）。利用此方法可以較為準確地提取丹江口流域的所有水體信息，包括影像中比較小的水體和流域，與上述方法相比，陰影對水體信息提取的限制也得到了很大的改善，但對于建筑物和居民地沒有更好的識別，存在著誤提的問題。

3.2 譜間關(guān)系法提取結(jié)果與分析

上述3種水體提取方法，都有各自的不足之處；單波段閾值法不能有效地精確提取細小的水體和流域，并且其水體的精確率在于閾值的確定，而閾值的確定到目前為止都是一個難點；NDWI法雖然改善了植被、建筑物和居民地等信息誤提的問題，但存在著誤把陰影和裸地提為水體的問題，直接限制了水體提取的精度。而研究范圍是丹江口水庫庫區(qū)，水域面積較大且連續(xù)性強，因此，選取譜間關(guān)系法來對研究區(qū)水體進行提?。淮颂崛》椒ú粌H可以對丹江口水庫庫區(qū)大面積水體進行有效提取，也對陰影的誤提有了很大的改進[7～10]；不用對閾值的確定進行反復試驗，避免了繁瑣的工作，實現(xiàn)水體快速有效的自動提取，提高了工作效率。采用譜間關(guān)系法對丹江口影像進行水體信息的提取，得到如圖2所示的利用譜間關(guān)系法提取水體的效果圖[11～13]。

4 結(jié)論與討論

4.1 主要結(jié)論

（1） 不同于上述兩種方法，譜間關(guān)系法不用進行閾值的確定，避免了最佳閾值選擇所帶來的誤差和麻煩，相比之下能更好地提取水體信息。

（2） 為了減少由于地形坡度、陰影等自然環(huán)境條件的限制，宜采用波段之間的比值運算進行水體信息的提取，從而提高水體識別的精確性。

（3） 通過幾種方法的比較和分析，最終根據(jù)丹江口水庫庫區(qū)本身的特征選取了譜間關(guān)系法提取水體信息，從得到的提取效果圖1中可以發(fā)現(xiàn)，陰影和建筑物基本上消除，水體輪廓清晰。

（4）丹江口水庫水體面積隨著季節(jié)而發(fā)生變化，受到降雨量、降水強度及降水時間的巨大影響。通常丹江口水庫水體面積在600 km2左右。通過對譜間關(guān)系法的提取效果圖進行判讀和分析，借助GIS相關(guān)軟件對丹江口水庫水體信息的提取面積進行計算為599.89 km2，同丹江口水庫實際水體面積相比，具有較高的提取率。

（5）最后采用譜間關(guān)系法對丹江口水庫進行了水體信息的提取，雖然有細小水體和水流的漏提，但相對于丹江口庫區(qū)大面積水體而言沒有影響，保證了對此研究區(qū)水體提取的精度。

4.2 問題探討

（1）通過對幾種水體信息提取方法的比較和分析發(fā)現(xiàn)，對于目前對水體提取的技術(shù)上，建筑物、陰影的誤提對水體信息提取有很大的干擾，因此今后應針對這兩個問題進行更深入的研究和分析，以找到更好、更精確的水體提取方法。

（2）要根據(jù)實地情況和需求，將遙感影像特征和研究區(qū)地域自然環(huán)境特征相結(jié)合，進而建立起適合不同地域環(huán)境的提取模型，以實現(xiàn)簡單快速精確的水體提取。

（3）波段閾值的確定一直是一個難點，它限制了水體提取的準確性，因此，應在閾值范圍選取及確定這方面加以深入研究，從而提高閾值確定的準確性。

（4）在遙感影像上，有些陰影與水體的光譜特征比較相似，這在遙感上稱為“異物同譜”現(xiàn)象；這種現(xiàn)象降低了提取水體的精度，因此，如何消除異物同譜現(xiàn)象是今后研究的重點。

此外，由于細小水體的面積小，光譜特征不明顯，容易漏提，因此，細小水體的精確提取也將成為今后的研究重點之一。

參考文獻：

[1]劉家宏，王光謙，王開.數(shù)字流域研究綜述[J].水利學報.2006，37 （2）：240～246.

[2]韋玉春，湯國安，楊 昕，等.遙感數(shù)字圖像處理教程[M].北京：科學出版社，2007.

[3]梅安新，彭望琭，秦齊明.遙感導論[M] .北京：高等教育出版社，2001.

[4]徐涵秋，利用改進的歸一化差異水體指數(shù)（MNDWI）提取水體信息的研究[J].遙感學報，2005，9（5）：589～595.

[5]楊存建，徐美.遙感信息機理的水體提取方法的探討[J].地理研究，1998，17（4）：86～89.

[6]McFEETERS S K . The Use of the Normalized Difference Water Index （NDWI） in the Delineation of Open Water Features [J]. International Journal of Remote Sensing， 1996， 17（7）： 1425.1432.

[7]玉素甫江·如素力，李蘭海. 基于Landsat ETM+的內(nèi)陸湖泊水體信息提取方法研究[J].西北農(nóng)林科技大學學報（自然科學版），2013（12）

[8]劉桂林，張落成，劉 劍，等. 基于Landsat TM影像的水體信息提取[J].中國科學院大學學報，2013，30（5）：644～650.

[9]程武學，李 悅，李惠君，等. 基于ETM+數(shù)據(jù)的水體識別方法對比研究[J].湖北農(nóng)業(yè)科學，2013，52（21）：5326～5328.

[10]蔣 怡，胡 濱，何政偉，等. 基于TM數(shù)據(jù)的水體信息提取試驗[J].地理空間信息，2013，11（6）：101～102.

[11]李 晶，李恩來，苗慧玲，等. 基于LandsatTM數(shù)據(jù)的濟寧市水體識別方法研究[J].環(huán)境科學與技術(shù)，2013（19）：175～179.

[12]杜先榮，楊建生，王飛紅.一種改進的基于TM影像的水體提取方法[J] .濟南大學學報（自然科學版），2014，28（6）：434～438.