倪衡　劉曉穎

摘 要：日益增多的不合理的圍填?；顒訉０稁Лh(huán)境造成了不可逆轉(zhuǎn)的破壞。隨著Landsat8的成功發(fā)射，我們需要一種針對海岸帶附近的用?；顒舆M(jìn)行及時(shí)有效的監(jiān)測的方法。根據(jù)圍填海活動由海向陸硬化過程不同階段代表地物的光譜特性，在如今大規(guī)模使用Landsat8遙感數(shù)據(jù)的條件下，根據(jù)其波段特性，提出新的水體指數(shù)NWI8對光譜進(jìn)行處理并比較其區(qū)分水體、灘涂和人工填海地物的能力。結(jié)果發(fā)現(xiàn)，基于Landsat8特性提出的新型水體指數(shù)可以充分滿足區(qū)分由海向路硬化過程的不同階段的需求，高效的完成快速監(jiān)測用?；顒拥娜蝿?wù)。

關(guān)鍵詞：水體指數(shù) 地物光譜 遙感 Landsat8

中圖分類號：TP79 文獻(xiàn)標(biāo)識碼：A 文章編號：1672-3791（2015）07（a）-0007-03

Abstract：Increasingly irrational reclamation activities are causing irreversible damage to the coastal environment.With the launch of Landsat8，we need a method to monitor the activities of using the sea quickly and effectively.According to the spectral characteristics of the objects which represent the various stages of the changing from sea to land in reclamation activities，we use a new method based on its bands which is called NWI8 to research the ability to distinguish water，bench and artificial reclamation features.Finally we find this new method can complete the task well.At the end，this new method based on the bands of Landsat8 can do the job which monitor the sea using.

Key Words：Water Index；Spectrum；Remote Sensing；Landsat8

在土地資源越發(fā)緊張的情況下，易開發(fā)的濱海地區(qū)成為了時(shí)下熱門的新型開活動場所。濱海地區(qū)的海岸帶資源是我國重要的自然資源，不合理的違法開發(fā)利用會對濱海地區(qū)的生態(tài)環(huán)境造成不可逆轉(zhuǎn)的破壞。因此，為了保護(hù)海岸帶生態(tài)系統(tǒng)的平衡，促進(jìn)海岸帶自然資源的合理高效利用，實(shí)現(xiàn)海洋資源的可持續(xù)發(fā)展，海洋管理部門就有必要準(zhǔn)確、及時(shí)、高效地監(jiān)視監(jiān)測海岸帶開發(fā)利用活動。而海岸帶開發(fā)利用活動以圍填海為主，圍填海的過程也是由海向陸的硬化過程。對此，為了提取水體信息而創(chuàng)建的水體指數(shù)可以有效的反映這一硬化過程。目前常見的水體指數(shù)有：NDWI、MNDWI、EWI和NWI。在Landsat8遙感數(shù)據(jù)開始廣泛應(yīng)用的今天，該文基于其數(shù)據(jù)新的特性創(chuàng)建了New Water Index 8（NWI8），并針對由海向陸硬化過程不同階段提出不同參考指標(biāo)。

1 數(shù)據(jù)源

1.1 Landsat 8波段特性

2013年2月新一代的陸地資源衛(wèi)星Landsat8成功發(fā)射并在隨后開始向地面?zhèn)鬏敂?shù)據(jù)。Landsat8搭載的OLI傳感器在改進(jìn)了Landsat7衛(wèi)星ETM+傳感器的7個(gè)波段的基礎(chǔ)上并增加了2個(gè)新的波段——觀測海岸帶氣溶膠的波段1和觀測卷云的波段9[1]。

1.2 樣本數(shù)據(jù)

2012年7月，在連云港地區(qū)選取了多種地物進(jìn)行地物光譜測定，該文選取水體、灘涂和填海地物等三種作為圍填海活動由海向陸硬化過程的三個(gè)階段。

該文所選取的地物光譜數(shù)據(jù)波段范圍從330.48nm-2389.74nm共1560個(gè)波值，并且在可見光波段的分布密度大于1nm。

水體地物光譜包括普通海水、堿液和鹵水三種類型，共41個(gè)樣本，每種類型選取兩個(gè)樣本，其對應(yīng)的地物光譜曲線見圖1。

對于清水，在藍(lán)綠波段反射率4%～5%，0.6以后的紅光部分反射率降到2%～3%[6]。由于水體在近紅外及隨后的中紅外波段內(nèi)（740～2500nm）所具有的強(qiáng)吸收的特點(diǎn)，導(dǎo)致了清澈水在這一波長范圍內(nèi)幾乎無反射率，因此，這一波長范圍通常被用來研究水陸分界、圈定水體范圍[7]。

灘涂地物光譜包括灘涂、淤泥和泥混合三種類型，共27個(gè)樣本，每種類型選取兩個(gè)樣本，其對應(yīng)的地物光譜曲線見圖2。

如圖所示，灘涂各樣本的光譜曲線相似，反射率的高低主要是由于不同樣本的含水量的差別。上圖顯示灘涂樣本反射率曲線的“峰-谷”變化較弱，曲線的形態(tài)沒有水體樣本光譜曲線變化明顯。[6]總體看來，在可見光和近紅外波段，灘涂的反射率一般都是隨著波長的增加而增加；在中紅外波段范圍內(nèi)，灘涂的反射率會因?yàn)楦髯缘慕M成成分與含水量的高低而發(fā)生不同程度的波動。

填海是圍填?；顒拥淖罱K階段，是由海向陸硬化過程的最重要的一環(huán)，其標(biāo)志性的地物就是各種人工地物及各種填海材料。該文選取的填海地物包括混凝土、瀝青、砂石、土壤和巖石五種類型，共207個(gè)樣本，在5種類型中選取一個(gè)樣本，其對應(yīng)的地物光譜曲線見圖3。

五種填海地物的反射率整體上隨著波長的增加而變大。其中，瀝青、混凝土和砂石等人工地物的光譜曲線較為平緩，而土壤和巖石等填海材料的光譜曲線在中紅外波段處有較大波動。這是因?yàn)槠渌惖匚锏暮亢艿?，而土壤和巖石在填海過程中會附著一定程度的水分，所以在中紅外波段上有較大波動。

2 水體指數(shù)

2.1 常見水體指數(shù)

目前，區(qū)分水體與其他地物信息的方法主要有單波段閾值法、譜間關(guān)系法、指數(shù)模型法、監(jiān)督分類法和決策樹法，其中水體指數(shù)法結(jié)合閾值分割提取水體較為快捷可靠，且已經(jīng)成為最普遍的提取水體的方法。因?yàn)樗w指數(shù)是利用各地物在不同波段的反射率不同的特性所組成的歸一化差值方法，而各遙感器的波段特性不一，所以目前存在著多種水體指數(shù)，常見的幾種水體指數(shù)如下：

（1）McFeeters（1996）借鑒NDVI的構(gòu)思，提出了歸一化差異水體指數(shù)NDWI（Normalized Difference Water Index）[2]：

（2）徐涵秋（2005）在針對TM數(shù)據(jù)中城市、水體及其背景地物的反射特點(diǎn)分析的基礎(chǔ)上，提出了改進(jìn)歸一化差異水體指數(shù)MNDWI（Modified NDWI）[3]：

（3）閆霈利用ETM+影像的Green、NIR和MIR構(gòu)建了增強(qiáng)型水體指數(shù)EWI（Enhanced Water Index）[4]：

（4）丁鳳（2009）結(jié)合水體在NIR和MIR波段同時(shí)具有強(qiáng)吸收的特點(diǎn)，提出了新型水體指數(shù)NWI（New Water Index）[5]：

2.2 新型水體指數(shù)

通過仔細(xì)分析水體與其它地物的光譜特征發(fā)現(xiàn)，水體尚有以下2個(gè)在以往研究中未被充分挖掘和利用的光譜特征：（1）水體在藍(lán)光波段（Band 1）的反射要高于綠光波段（Band 2），以EWI水體指數(shù)構(gòu)建為例，若采用藍(lán)光波段代替綠光波段，可以進(jìn)一步擴(kuò)大水體與其它地物的反差，且使得大部分水體的亮度值為正；（2）在眾多地物類型中，唯有水體具有如下特點(diǎn)，即水體在OLI影像的第6、7波段同時(shí)具有強(qiáng)烈的吸收，而植被、干土壤和建筑物等非水體地物在這幾個(gè)波段范圍內(nèi)的反射則均高于水體。引入水體的上述兩個(gè)典型特征，本文提出了一種可用于提取水體信息的新型水體指數(shù)NWI8，其公式如下：

3 數(shù)據(jù)處理分析

為了測試新型水體指數(shù)對由海向陸硬化過程不同階段的區(qū)分能力，將41個(gè)水體樣本、27個(gè)灘涂樣本和207個(gè)人工填海地物樣本的地物光譜數(shù)據(jù)按照不同水體指數(shù)方法求值。這些求得的水體指數(shù)值在一定范圍內(nèi)有集中分布的趨勢，在（-1，1）的區(qū)間上，以0.02為組距，分布頻率為縱坐標(biāo)，各地物光譜經(jīng)水體指數(shù)求出的值呈現(xiàn)出圖4的分布形式。這些水體指數(shù)值的總體特征參數(shù)取四位有效數(shù)字見表2。將地物值域相交區(qū)域的樣本個(gè)數(shù)與各自樣本總數(shù)相比可得到表2中的樣本混合比例。

從上圖和上表中，可以看出在新型水體指數(shù)方法下，按照Landsat 8波段特性計(jì)算出來的水體指數(shù)值，水體顯著區(qū)別于灘涂和人工填海地物，說明可以通過在設(shè)定一定的閾值下，有效地將水體信息提取出來。

對于灘涂和人工填海地物，新型水體指數(shù)都無法將它們完全區(qū)分開，NWI8對于它們的區(qū)分能力較好，分別達(dá)到了62.97%和54.11%。

4 結(jié)語

按照Landsat8波段特性，提出的新型水體指數(shù)NWI8在對實(shí)測地物光譜處理后，發(fā)現(xiàn)其可以較好地提取海岸帶地物信息，對于由海向陸硬化過程不同階段代表地物類型水體、灘涂和人工填海地物的區(qū)分也能到達(dá)較好水平，所以在Landsat8遙感數(shù)據(jù)條件下對海岸帶的圍填?；顒舆M(jìn)行監(jiān)測時(shí)，可以選用NWI8。

參考文獻(xiàn)

[1] 徐涵秋，唐菲.新一代Landsat系列衛(wèi)星：Landsat 8遙感影像新增特征及其生態(tài)環(huán)境意義[J].生態(tài)學(xué)報(bào)，2013，33（11）：3249-3257.

[2] McFeeters S K.The Use of Normalized Difference Water Index（NDWI）in the Delineation of Open Water Features [J].International Journal of Remote Sensing，1996，17（7）：1425-1432.

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

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

[5] 丁鳳.基于新型水體指數(shù)（NWI）進(jìn)行水體信息提取的實(shí)驗(yàn)研究[J].測繪科學(xué)，2009，39（4）：155-157.

[6] 趙英時(shí).遙感應(yīng)用分析原理與方法[M].北京：科學(xué)出版社，2003.

[7] Jensen J R. Introductory Digital Image Processing： A Remote Sensing Perspective 3rd [M]. New Jersey： Prentice Hall，2004.

[8] 徐涵秋.從增強(qiáng)型水體指數(shù)分析遙感水體指數(shù)的創(chuàng)建[J]. 地球信息科學(xué)，2008，10（6） ：776-780.

[9] 王志輝，易善楨.不同指數(shù)模型法在水體遙感提取中的比較研究[J].科學(xué)技術(shù)與工程，2007，7（4）：534-537.