01

概述

相信大家之前有看过很多通过电子地图提取水系的教程，该方法虽然可以提取水系，但是有一定的限制，如果在一些比较偏远的地方，电子地图上没有水系数据就无法提取，存在一定的局限性，之前已经为大家介绍了如何从卫星地图上提取道路的方法（详见《如何使用QGIS自动矢量化卫星图》【点击了解】），这里为大家介绍一下通过卫星地图提取水系数据的方法，希望能对大家有所帮助。

02

波段合成

从landsat 8的波段组合方式可以知道，NIR、SWIR1、Red（即B5 、B6、 B4）三个波段组合可以得到陆地和水域的效果，所以分别将这三个波段加载到ArcMap内。

在菜单栏上点击窗口，选择影像分析，如下图所示。

选择影像分析

在显示的影像分析对话框内，选中某一个波段图层，并点击添加函数，如下图所示。

点击添加函数

在显示的函数模板编辑器对话框内，找到图层，点击右键，选择“插入\波段合成函数”，如下图所示。

选择波段合成函数

在显示的栅格函数属性对话框内，将所有需要合成的波段添加进来，如下图所示。

添加所有波段

点击确定完成波段的合成，合成的影像如下图所示。

合成的影像

03

非监督分类

在ArcToolbox中点击“Spatial Analyst工具\多元分析\ Iso聚类非监督分类”，调用Iso聚类非监督分类工具，如下图所示。

调用Iso聚类非监督分类工具

在显示的Iso聚类非监督分类对话框内，输入栅格波段为上一步合成的影像，类目录设置为4，其余的默认，如下图所示。

Iso聚类非监督分类设置

点击确定完成Iso聚类非监督分类，分类结果如下图所示。

分类结果

04

影像重分类

在ArcToolbox中点击“Spatial Analyst工具\重分类\ 重分类”，调用重分类工具，如下图所示。

调用重分类工具

在显示的重分类对话框内，将旧值3和4设置为新值2，如下图所示。

重分类设置

点击确定完成重分类，重分类结果如下图所示。

重分类结果

05

矢量化栅格

分别新建一个线图层和一个面图层，如下图所示。

新建的图层

将线图层进入编辑状态，打开ArcScan工具条，在工具条上点击矢量化，选择矢量化设置，如下图所示。

选择矢量化设置

在显示的矢量化设置对话框内，点击样式，在显示的矢量化设置样式选择轮廓，如下图所示。

选择轮廓

确定之后，在ArcScan工具条上点击矢量化，选择生成要素，如下图所示。

选择生成要素

点击之后生成水系，生成的水系如下图所示。

生成的水系

保存编辑结果，调整样式并添加卫星图，最后效果如下图所示。

水系效果

06

结语

以上就是如何从卫星图中提取水系数据的详细说明，主要包括了波段合成、非监督分类、影像重分类和矢量化栅格等功能。