ArcGIS实战：利用LandSat8 提取水体并计算面积

ArcGIS实战：利用LandSat8 提取水体并计算面积

1.数据下载

LandSat 8 遥感影像属于多波段遥感影像，利用不同波段可以做一些分析，比如NDWI（水体提取），NDVI（绿地提取）等。
数据下载地址：地理空间数据云

（1）进入地理空间数据云下载数据

在地理空间数据云中选择【高级检索】，选择数据集为【LandSat系列数据】中的【LandSAT 8 OLI_TRIS卫星数字产品】，根据自身需求，可选择按行政区查询或经纬度或自定义范围查询，查询什么时间可自选。如果不会用ENVI做大气校正，辐射定标等操作的，最好选择云量少的数据。选择好后，点击【检索】，即可查看数据集，点击数据集，可查看影像覆盖范围

点击上图所示按钮，即可下载数据

2.提取水体

提取水体的公式为：NDWI =(p(Green)-p(NIR))/(p(Green)+p(NIR))
各波段介绍如图所示

实战：
运用【栅格计算器】，输入以下公式：

简单表示为：(FloatB3-FloatB6)/(FloatB3+floatB6)
输出后图像如图

此时，该图层的值变为-1到1，我们需要提取到水体层面，且对水体进行归一化运算。就需要设置水体的提取阈值，一般是0.2，但如需提取到更加准确的水体，需用识别按钮点击水体边缘，查看具体的值，例如，下图水体通过查看水体边缘值后，得知该边缘值为0.17…，所以设置阈值为0.17
使用【重分类】工具，输入上一步提取的水面栅格，如图进行分类

将小于0.17的部分设为NoData

点击【确定】，即可得到水面栅格，如图所示

此时，水体已经全部提取可对其进行矢量化
运用【栅格转面】工具，即可转出

3.计算面积

通过上一步提取的水体，即可计算水面面积，通过查看矢量水面的数据源，其坐标系为投影坐标系，可直接计算。
在矢量水面的属性中，新建一个双精度的字段，在字段上右键，计算几何，选择单位，即可得出水面面积。