基于landsat数据的水体提取流程

1输入：

一幅landsat光学影像、鄱阳湖湖面矢量图

2过程：

2.1导入数据

1、解压landsat光学遥感影像包。
2、打开envi
3、在envi里，打开为–光学传感器–Landsat–带有元数据的GeoTIFF–定位到解压好的landsat光学遥感影像包，导入XXXXMTL.txt的数据。

2.2数据预处理

2.2.1辐射定标

------点击右侧的Toolbox–辐射校正（Radiometric Correction）–辐射定标（Radiometric Calibration）

------点击应用FLAASH按键，设置文件输出的路径。

2.2.2辐射校正

------点击右侧的Toolbox–辐射校正（Radiometric Correction）–大气校正模块（Atmospheric Correction Module）–FLAASH大气校正（FLAASH Atmospheric Correction）

------第一项选择上一步辐射定标的数据

------选择第二个

------设置输出路径（整个路径不能出现中文字）

------将Sensor Type的Multispectral设置为Landsat-8 OLI

------将Ground Elevation (km)设置为0.121

------设置Atmospheric Model

------按照以下的表格，根据月份和纬度数据进行选择大气模型。具体如下：
如果是1月2月选择Jan列、3月4月选择March列、5月6月选择May列、7月8月选择July列、9月10月选择Sept列、11月12月选择Nov列。

再根据纬度选择行，定位到Atmospheric Model。

（比如示例的这一期影像的日期为1月，纬度为28°51′52.63″，四舍五入一下lat为30。那么根据下表可以定位Atmospheric Model选择MID-LATITUDE SUMMER模式。）

Lat	Jan(1月)	March(3月)	May(5月)	July(7月)	Sept(9月)	Nov (11月)
80	SUB-ARCTIC WINTER	SUB-ARCTIC WINTER	SUB-ARCTIC WINTER	MID-LATITUDE WINTER	MID-LATITUDE WINTER	SUB-ARCTIC WINTER
70	SUB-ARCTIC WINTER	SUB-ARCTIC WINTER	MID-LATITUDE WINTER	MID-LATITUDE WINTER	MID-LATITUDE WINTER	SUB-ARCTIC WINTER
60	MID-LATITUDE WINTER	MID-LATITUDE WINTER	MID-LATITUDE WINTER	SUB-ARCTIC SUMMER	SUB-ARCTIC SUMMER	MID-LATITUDE WINTER
50	MID-LATITUDE WINTER	MID-LATITUDE WINTER	SUB-ARCTIC SUMMER	SUB-ARCTIC SUMMER	SUB-ARCTIC SUMMER	SUB-ARCTIC SUMMER
40	SUB-ARCTIC SUMMER	SUB-ARCTIC SUMMER	SUB-ARCTIC SUMMER	MID-LATITUDE SUMMER	MID-LATITUDE SUMMER	SUB-ARCTIC SUMMER
30	MID-LATITUDE SUMMER	MID-LATITUDE SUMMER	MID-LATITUDE SUMMER	TROPICAL	TROPICAL	MID-LATITUDE SUMMER
20	TROPICAL	TROPICAL	TROPICAL	TROPICAL	TROPICAL	TROPICAL
10	TROPICAL	TROPICAL	TROPICAL	TROPICAL	TROPICAL	TROPICAL
0	TROPICAL	TROPICAL	TROPICAL	TROPICAL	TROPICAL	TROPICAL
-10	TROPICAL	TROPICAL	TROPICAL	TROPICAL	TROPICAL	TROPICAL
-20	TROPICAL	TROPICAL	TROPICAL	MID-LATITUDE SUMMER	MID-LATITUDE SUMMER	TROPICAL
-30	MID-LATITUDE SUMMER	MID-LATITUDE SUMMER	MID-LATITUDE SUMMER	MID-LATITUDE SUMMER	MID-LATITUDE SUMMER	MID-LATITUDE SUMMER
-40	SUB-ARCTIC SUMMER	SUB-ARCTIC SUMMER	SUB-ARCTIC SUMMER	SUB-ARCTIC SUMMER	SUB-ARCTIC SUMMER	SUB-ARCTIC SUMMER
-50	SUB-ARCTIC SUMMER	SUB-ARCTIC SUMMER	SUB-ARCTIC SUMMER	MID-LATITUDE WINTER	MID-LATITUDE WINTER	SUB-ARCTIC SUMMER
-60	MID-LATITUDE WINTER	MID-LATITUDE WINTER	MID-LATITUDE WINTER	MID-LATITUDE WINTER	MID-LATITUDE WINTER	MID-LATITUDE WINTER
-70	MID-LATITUDE WINTER	MID-LATITUDE WINTER	MID-LATITUDE WINTER	MID-LATITUDE WINTER	MID-LATITUDE WINTER	MID-LATITUDE WINTER
-80	MID-LATITUDE WINTER	MID-LATITUDE WINTER	MID-LATITUDE WINTER	SUB-ARCTIC WINTER	MID-LATITUDE	WINTER MID-LATITUDE WINTER

注：月份信息和纬度信息在选择了输入文件时会自动带入。

------点击Multispectral Settings 将Defaults设置为Over-Land Retrieval standard(660:2100nm)

------点击ok再点击apply，开始FLAASH大气校正。
------校正完成，数据预处理完成。

2.2.3水体提取

------使用envi toolbox–Band Algebra–Band Math

-----输入以下的公式： (b3-b6)*1.0/(b3+b6)*1.0

-----B3对应第三项的波段、B6对应第六项的波段。
然后得到一个ndwi的图

-----然后右击生成的ndwi的数据

-----这两个数值都设为0。

-----提取完成