杞麓湖與異龍湖水體提取及湖泊面積變化動(dòng)態(tài)監(jiān)測(cè)

肖茜+楊昆+曹瓅+洪亮

摘 要：該文以杞麓湖和異龍湖的近30年來的TM、ETM+和OLI遙感影像為數(shù)據(jù)源，采用NDWI、MNDWI、NWI、EWI和AWEI5種水體提取指數(shù)來提取杞麓湖和異龍湖的湖泊水體，并對(duì)比各種方法的提取精度。結(jié)果顯示，杞麓湖的面積在1985年和2000年呈增長趨勢(shì)，并在2000年達(dá)到最大值37.019km2，之后逐年減少，2015年面積僅為22.408km2。異龍湖在1985—1990年異龍湖湖面面積減少了7.16km2，到1995年面積大幅度回升到34.041km2，之后幾年的面積逐年減小，至2010年，湖面面積較1985年縮小了1/3，僅有24.483km2，至2015異龍湖的面積縮小到了12.84km2。分析其面積變化的影響因素發(fā)現(xiàn)，氣候因素與人類活動(dòng)對(duì)其影響較大。

關(guān)鍵詞：水體指數(shù)；湖泊面積；動(dòng)態(tài)監(jiān)測(cè)；降水量

中圖分類號(hào) P343.3 文獻(xiàn)標(biāo)識(shí)碼 A 文章編號(hào) 1007-7731（2017）24-0123-03

水是大自然賦予人類最寶貴的財(cái)富，它是生命的源泉，是人類賴以生存和發(fā)展不可缺少的最重要的物質(zhì)資源之一[1]。云南省高原湖泊是影響區(qū)域生態(tài)環(huán)境的關(guān)鍵因子之一，其在社會(huì)經(jīng)濟(jì)發(fā)展過程中占據(jù)了重要的戰(zhàn)略地位。云南省高原湖泊的流域生態(tài)環(huán)境十分脆弱[2]，受氣象因素及人類活動(dòng)影響較大。近年來，云南省持續(xù)嚴(yán)重旱災(zāi)，湖泊水質(zhì)下降，杞麓湖與異龍湖表面水體面積萎縮嚴(yán)重。因此，研究杞麓湖與異龍湖水體提取及其近幾十年的變化具有十分重要的科學(xué)意義。

隨著遙感（Remote sensing，RS）和地理信息技術(shù)（Geographic Information System，GIS）的發(fā)展，遙感數(shù)據(jù)逐漸成為各類地理研究中重要的數(shù)據(jù)源。近幾年來，一些學(xué)者應(yīng)用RS和GIS技術(shù)對(duì)表面水體面積進(jìn)行了提取和動(dòng)態(tài)檢測(cè)，最初，McFeeters在1996年提出了歸一化差值水體指數(shù)（NDWI）[3]；徐涵秋在NDWI的基礎(chǔ)上提出了改進(jìn)型歸一化差值水體指數(shù)MNDWI[4]；陳華芳、王金亮等在提取香格里拉水體時(shí)應(yīng)用差值法、多波段譜間關(guān)系法分別與閾值法結(jié)合能有效地消除山體陰影對(duì)水體提取的影響[5]；石振杰、溫興平等利用纓帽變換及密度分割技術(shù)快速準(zhǔn)確地提取了云南省撫仙湖水體信息，分析了撫仙湖2000—2015年的面積變化及其變化的影響因素[6]；此外還有一些學(xué)者利用RS與GIS技術(shù)對(duì)湖泊的水體表面面積的變化做了一些研究[7-9]。

目前，采用水體指數(shù)法提取水體已經(jīng)得到了廣泛的應(yīng)用，但對(duì)于杞麓湖與異龍湖的研究較少。本文以杞麓湖和異龍湖流域作為研究區(qū)，利用遙感和GIS技術(shù)，基于TM和OLI遙感影像，采用5種水體提取指數(shù)法提取杞麓湖和異龍湖的表面水體面積，并對(duì)1985—2015年近30年湖泊表面水體面積的變化進(jìn)行分析。

1 研究區(qū)概況及遙感數(shù)據(jù)

1.1 研究區(qū)概況 云南省地處祖國西南邊陲，屬于典型的山地高原地形。杞麓湖位于于云南省玉溪市通?？h，海拔約1796.62m，是玉溪市通?？h的重要水資源；異龍湖位于滇南地區(qū)紅河州石屏縣，是石屏人民的“母親河”。

1.2 數(shù)據(jù)選擇 本文采用的數(shù)據(jù)是從USGS網(wǎng)站上下載的杞麓湖與異龍湖的1985年、1990年、1995年、2000年、2005年、2010年和2015年的Landsat數(shù)據(jù)，所有數(shù)據(jù)云量不超過5%且沒有條帶。

2 水體提取方法

2.1 不同提取方法簡介

2.1.1 歸一化水體指數(shù)NDWI McFeeters提出了Normalized Difference Water Index（NDWI），即歸一化水體指數(shù)。

[NDWI=Green-NIRGreen+NIR] （1）

2.1.2 改進(jìn)歸一化水體指數(shù)MNDWI 改進(jìn)型歸一化水體指數(shù)是徐涵秋將NDWI中的近紅外波段換成了中紅外波段：

[MNDWI=Green-MIRGreen+MIR] （2）

2.1.3 增強(qiáng)型水體指數(shù)EWI 閆霈[10]等提出了增強(qiáng)型水體指數(shù)EWI，公式如下：

[EWI=Green-NIR-MIRGreen+NIR+MIR] （3）

2.1.4 新型水體指數(shù)NWI 丁鳳[11]提出新型水體指數(shù)NWI，其公式如下：

[NWI=Blue-（NIR+MIR+SWIR）Blue+（NIR+MIR+SWIR）*C] （4）

式中，C為一常數(shù)，目的是拉伸NWI的數(shù)值區(qū)間，在本實(shí)驗(yàn)中C取255。

2.1.5 自動(dòng)化水體提取指數(shù)AWEI Feyisa等[12]提出了Automated Water Extraction Index（AWEI），其公式有2個(gè)，一個(gè)是在理想狀態(tài)下（沒有陰影）提取水體，見公式（5）；另一個(gè)是在有陰影的情況下提取水體，見公式（6）。

[AWEInsh=4*（Green-MIR）-（0.25*NIR+2.75*SWIR）] （5）

[AWEIsh=Blue+2.5*Green-1.5*（NIR+MIR）-0.25*SWIR] （6）

注：上述式中，Green代表遙感影像的綠光波段，Blue代表藍(lán)光，波段NIR代表近中外波段，MIR代表中紅外波段，SWIR代表熱紅外波段。

2.2 水體提取精度對(duì)比分析 利用遙感高分辨率影像，在每個(gè)湖泊研究區(qū)選取80～100個(gè)測(cè)試樣本點(diǎn)，進(jìn)行混淆矩陣計(jì)算其精度，5種提取方法精度對(duì)比見表1。因此，杞麓湖選用NDWI法，異龍湖選用NWI法。

3 結(jié)果與分析

3.1 杞麓湖與異龍湖面積變化分析 杞麓湖和異龍湖1985—2015年的面積變化見圖1和圖2，杞麓湖的面積在1985年和2000年呈增長趨勢(shì)，增大了1.013km2，并在2000年達(dá)到最大值37.019km2，之后逐年減少，2000—2015年面積縮減14.611km2。異龍湖于1981年干涸，后因降雨湖水慢慢積蓄起來，到1985年達(dá)到蓄水量最大值。1985年至1990年異龍湖湖面面積減少了7.16km2，到1995年面積大幅度回升到34.041km2，之后幾年的面積逐年減小，至2005年，湖面面積較1985年縮小了1/3，僅有24.483km2，2010—2015年短短5年的時(shí)間，異龍湖面積縮小了50%，僅有12.84km2（圖2）。endprint

云南大旱導(dǎo)致降水量減少，工農(nóng)業(yè)用水緊張，湖泊的湖水量不能滿足蒸發(fā)和人類活動(dòng)的正常需要，多個(gè)水位下降至近30年來最低，其中，杞麓湖面積縮減了1/3，異龍湖面積縮減了2/3。但據(jù)氣象資料，2014—2015年降水增加，而杞麓湖和異龍湖面積仍然極速縮減，查閱相關(guān)資料得知，2個(gè)湖泊面積變化與人類活動(dòng)也息息相關(guān)。由于云南省2009年底開始的大旱，導(dǎo)致杞麓湖幾乎干涸，周邊村民將湖面干涸部分圍湖造田，使杞麓湖面積大大減少。而異龍湖面積持續(xù)下降的原因是云南連續(xù)干旱導(dǎo)致異龍湖的主要水源補(bǔ)給—赤瑞湖的干涸，加上周圍村民圍湖養(yǎng)魚、提水灌田，目前異龍湖已裂成兩半，西半部已成草原。因此，盡管2013—2015年降水量增加，但異龍湖的表面水體面積還是不可遏止地大幅度減少。

4 結(jié)論

本文選用5種不同的水體指數(shù)提取法提取了杞麓湖和異龍湖的水體，經(jīng)過精度對(duì)比，杞麓湖選用NDWI法，異龍湖選用NWI法。

近30年來，杞麓湖和異龍湖面積縮減均很嚴(yán)重。杞麓湖的面積在1985年和2000年呈增長趨勢(shì)，增大了1.013km2，并在2000年達(dá)到最大值37.019km2，之后逐年減少，2000—2015年面積縮減14.611km2。異龍湖在1985—1990年異龍湖湖面面積減少了7.16km2，到1995年面積大幅度回升到34.041km2，之后幾年的面積逐年減小，至2010年，湖面面積24.483km2，至2015年異龍湖的面積縮小到了12.84km2。至2015年杞麓湖面積僅為1985年面積的2/3，而異龍湖更是僅為1985年面積的1/3。分析其面積縮減的主要原因是云南省在2009—2012年長達(dá)4年的持續(xù)干旱及人類對(duì)湖泊環(huán)境的破壞。

參考文獻(xiàn)

[1]許素貞.淺談水資源的重要性[J].中國科技博覽，2012，34：370.

[2]董云仙，洪雪花.云南高原深水湖泊程海磷形態(tài)分布和變化[J].湖泊科學(xué)，2012，24（3）：341-346.

[3]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.

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

[5]陳華芳，王金亮，陳忠，等.山地高原地區(qū)TM影像水體信息提取方法比較[J].遙感技術(shù)與應(yīng)用，2014，6（19）：479-484.

[6]石振杰，溫興平，馬威.纓帽變換在撫仙湖水面變化分析中的應(yīng)用[J].安徽農(nóng)業(yè)科學(xué)，2015，43（25）：352-355.

[7]白潔，陳曦，李均力，等.1975—2007年中亞干旱區(qū)內(nèi)陸湖泊面積遙感分析[J].湖泊科學(xué)，2011，23（1）：80-88.

[8]Zhang B，Wu Y H and Zhu L P et a1.Estimation and trend detection of water storage at Nam Co Lake[J].Central Tibetan Plateau Journal of Hydrology，2011，405：161-170.

[9]邊多，邊巴次仁，拉巴，等.1975—2008年西藏色林錯(cuò)湖面變化對(duì)氣候變化的響應(yīng)[J].地理學(xué)報(bào)，2010，65（3）：313-319.

[10]閆霈，張友靜，張?jiān)?利用增強(qiáng)型水體指數(shù)EWI和GIS去噪音技術(shù)提取半干旱地區(qū)水系信息的研究[J].遙感應(yīng)用，2007，94（6）：62-67.

[11]丁鳳.一種基于遙感數(shù)據(jù)快速提取水體信息的新方法[J].遙感技術(shù)與應(yīng)用，2009，24（2）：167-171.

[12]Feyisa G L，Meilby H，F(xiàn)ensholt R，Proud S R.Automated water extraction index：A new technique for surface water mapping using Landsat imagery[J].Remote Sens.Environ，2014，140：23–35. （責(zé)編：張宏民）endprint