基于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。
-----提取完成
