鞏凌君

摘?要： 地理國情監(jiān)測是綜合利用全球衛(wèi)星導(dǎo)航定位技術(shù)（GNSS）、航空航天遙感技術(shù)（RS）、地理信息系統(tǒng)技術(shù)（GIS）等現(xiàn)代測繪技術(shù)，綜合各時期已有測繪成果檔案，對地形、水系、交通、地表覆蓋等要素進(jìn)行動態(tài)和定量化、空間化的監(jiān)測，并統(tǒng)計(jì)分析其變化量、變化頻率、分布特征、地域差異、變化趨勢等，形成反映各類資源、環(huán)境、生態(tài)、經(jīng)濟(jì)要素的空間分布及其發(fā)展變化規(guī)律的監(jiān)測數(shù)據(jù)、地圖圖形和研究報(bào)告。地理國情監(jiān)測通過對地理國情進(jìn)行動態(tài)地測繪、統(tǒng)計(jì)，從地理的角度來綜合分析和研究國情，為政府、企業(yè)和社會各方面提供真實(shí)可靠和準(zhǔn)確權(quán)威的地理國情信息。作為重要的地理國情，水體批量提取的工作，十分必要。針對幾種基于遙感影像的常見提取方法進(jìn)行對比分析，尋找較為完善的方案。

關(guān)鍵詞： ENVI;Landsat; ?水體提取

【中圖分類號】TP79?【文獻(xiàn)標(biāo)識碼】A?【DOI】10.12215/j.issn.1674-3733.2020.41.052

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

1、Landsat 8 是美國陸地衛(wèi)星計(jì)劃（Landsat）的第八顆衛(wèi)星，于2013年2月11號發(fā)射成功。Landsat 8上攜帶陸地成像儀（Operational Land Imager，OLI）和熱紅外傳感器（Thermal Infrared Sensor，TIRS）。OLI陸地成像儀共有9個波段，空間分辨率為30米。其中的短波紅外波段，通常應(yīng)用于分辨道路，裸露土壤，水，還能在不同植被之間有較好的對比度，并有較好的大氣、云霧分辨能力。

2、TM影像是指美國陸地衛(wèi)星Landsat4和Landsat5 攜帶的專題制圖儀（thematic mapper）所獲取的多波段掃描影像，包含7個波段。主要特點(diǎn)是具有較高空間分辨率、波譜分辨率、極為豐富的信息量和較高定位精度。美國的USGS網(wǎng)站提供最新的Landsat數(shù)據(jù)下載，產(chǎn)品幾何精度比較高，不需要再進(jìn)行幾何校正，但仍需要輻射定標(biāo)和大氣校正。

1?數(shù)據(jù)預(yù)處理

1.1?在ENVI中通過data manager 工具查看影像的相關(guān)屬性，確定研究區(qū)范圍，確定選擇的多光譜波段和區(qū)域平均高程。

1.2?輻射校正Radiometric Calibration

通過ENVI中的Toolboxs-Radiometric Correction-Radiometric Calibration。

1.3?大氣校正Atmospheric Correction

Toolboxs-Radiometric Correction-Atmospheric Correction Module-FLAASH Atmospheric Correction。

2?數(shù)據(jù)處理

利用單波段閾值法、譜間關(guān)系分析法、水體指數(shù)法對水體進(jìn)行提取。

2.1?單波段閾值法：

依據(jù)水體在不同波段光譜上的不同特征，以及其他地物與水體的區(qū)別，利用 landsat8OLI陸地成像儀的band 6 SWIR1短波紅外波段進(jìn)行水體的提取。具體步驟：

a）利用ENVI加載 band 6波段的全色圖像后，在Cursor Value中查看水體的dn值，發(fā)現(xiàn)基本都在700以內(nèi)，所以設(shè)置閾值為700。

b）利用 band math 添加公式：b1<700，并將b1賦以band 6波段。提取后黑色部分是水體提取的結(jié)果，白色高亮部分則是非水區(qū)域，其dn值統(tǒng)一為700。

2.2?譜間關(guān)系分析法：根據(jù)水體的（TM2+TM3）>（TM4+TM5）的特征，利用landsat8OLI陸地成像儀的band 3（green），band 4（red），band 5（near infrared），band 6（SWIR1）之間的波段計(jì)算，進(jìn)行水體提取。具體步驟：

a）ENVI-Basic Tools-band math-b1+b2>b3+b4;分別對應(yīng)TM3、TM4、TM5、TM6波段。

為了更好地提取水體，增加閾值，將公式修改為b1+b2-b3-b4>100。

2.3?水體指數(shù)法

水體指數(shù)法的基本原理是通過波段之間的比值運(yùn)算增強(qiáng)某些地物之間的反差。根據(jù)水體在TM2具有較強(qiáng)的透射性，TM5屬于水的吸收帶內(nèi)，歸一化水體指數(shù)（NDWI）的計(jì)算式為，NDWI =（TM2- TM5）/（TM2+TM5），當(dāng)滿足NDWI>C 時所提取的地物為水體。具體步驟：

a）ENVI-Basic Tools-band math-（b2-b4）/（b2+b4）-NDWI圖像。

b）確定水體的NDWI范圍： ENVI-Overlay-Region of Interest-ROI Tool-選擇樣本-確定范圍;max=0.49 min=0.22。

c）波段運(yùn)算輸出提取圖像：ENVI-Basic Tools-Band Math-b1<0.49，b1>0.22-NDWI的圖像。

3?統(tǒng)計(jì)分析

3.1?根據(jù)影像，繪制出各種地物的反射波譜曲線。

ENVI——Spectral——Spectral Libraries——Spectral Library Builder（工具構(gòu)建光譜庫）。

3.2?分析比較

（1）單波段閾值法方法簡單直接，但是無法區(qū)分陰影和水體，在山區(qū)密集地區(qū)，提取效果不佳。

（2）譜間關(guān)系分析法是國內(nèi)較為常用的的提取辦法，穩(wěn)定性強(qiáng)。相對于單波段閾值法，譜間關(guān)系分析法改進(jìn)了陰影區(qū)域的誤提取，但是還是會有少許陰影區(qū)域被誤提取，在城區(qū)附近，對于居民地也會存在一定的誤提取。

（3）水體指數(shù)法（NDWI）提取時，像建筑物陰影這些信息與水體在綠光波段和近紅外波段的特征較為相似，提取時會被誤提取，干擾水體提取效果。在山區(qū)提取水體時，也不容易區(qū)分山體陰影和水體。但是在實(shí)際操作時，可以通過用中紅外波段替換近紅外波段，計(jì)算出的明顯反差指數(shù)，降低陰影和居民地的影響，提高水提取的正確率，即改進(jìn)的歸一化差異水體指數(shù)（MNDWI）：MNDWI=（Green-MIR）/（Green+MIR）。

4?總結(jié)

單一的方法各有利弊，在實(shí)際操作中，針對各類問題，可以將各種方法有機(jī)結(jié)合，多次提取，提高水體的提取效果。例如在城市地區(qū)，可以利用改進(jìn)的歸一化差異水體指數(shù)（MNDWI）降低建筑居民地對水體的干擾，也可以先進(jìn)行提取，但是利用光譜特征曲線查看居民地的亮度值，通過閾值設(shè)置，降低居民地影響。對于陰影問題，可以利用譜間關(guān)系分析法降低陰影區(qū)域的影響，再通過閾值提高提取精度，或者先分類，在提取初步結(jié)果后，再聚類，增強(qiáng)效果。在復(fù)雜區(qū)域，可以分塊進(jìn)行提取，在非山區(qū)的植被區(qū)域用NDWI指數(shù)提取，在建筑居民區(qū)用MNDWI指數(shù)提取，唉喊去和農(nóng)林區(qū)域用譜間關(guān)系分析法，可以通過設(shè)置有效閾值來增強(qiáng)提取效果，真的注意的是，部分細(xì)小的水體不易區(qū)分，設(shè)置閾值時，較為困難，閾值過高會丟失細(xì)小水體，閾值過低，會誤提非水體信息。

參考文獻(xiàn)

[1]?吳文淵，沈曉華，雛樂君等.基于Landsat ETM+影像的水體信息綜合提取方法[J].科技通報(bào)，2008.3.

[2]?杜培軍.遙感原理與應(yīng)用[M].徐州：中國礦業(yè)大學(xué)出版社，2006.