gis地形起伏计算(8影像计算鄱阳湖面积变化)
gis地形起伏计算(8影像计算鄱阳湖面积变化)非监督分类结果点击运行后得到2017年影像非监督分类结果,如下图所示。调用Iso聚类非监督分类工具在显示的Iso聚类非监督分类对话框内,输入栅格波段为2017年的影像;因为主要是获取鄱阳湖面积,不需要设置太多的类数,所以类数设置为4;设置输出分类的栅格数据,后面的参数默认,如下图所示。Iso聚类非监督分类
概述鄱阳湖是我国第一大淡水湖,其面积在不同的季节差异比较大,这里我们通过ArcGIS Pro使用Landsat 8影像来计算鄱阳湖在不同季节的面积,关于如何下载Landsat 8影像,可以参考文章《如何在USGS下载Landsat 8-9影像(2022年版)》【点击了解】。
对比地图分别将2017年11月1日和2022年4月21日的影像加载到ArcGIS Pro内,通过卷帘工具可以对比两个时间点的差异,如下图所示。
对比图
非监督分类在工具箱中点击“Spatial Analyst工具\多元分析\Iso聚类非监督分类”,调用Iso聚类非监督分类工具,如下图所示。
调用Iso聚类非监督分类工具
在显示的Iso聚类非监督分类对话框内,输入栅格波段为2017年的影像;因为主要是获取鄱阳湖面积,不需要设置太多的类数,所以类数设置为4;设置输出分类的栅格数据,后面的参数默认,如下图所示。
Iso聚类非监督分类
点击运行后得到2017年影像非监督分类结果,如下图所示。
非监督分类结果
同样的方法对2022年的影像进行Iso聚类非监督分类,分类结果如下图所示。
非监督分类结果
调整样式分别对2017年和2022年影像非监督分类结果进行调整,除了水系之外均调整为透明,如下图所示。
调整后图
清除分类上面分类后可能存在一些被分类为水体但不属于鄱阳湖的小型孤立像素,比如小型池塘或水体,所以需要尽可能的将其删除,以免影响后面的面积计算。
在工具箱中点击“Spatial Analyst工具\栅格综合\众数滤波”,调用众数滤波工具,如下图所示。
调用众数滤波工具
在显示的众数滤波对话框内,输入栅格为2017年Iso聚类非监督分类影像,设置输出栅格名称,要使用的相邻像元数设置为四,替换阈值为半数,这样更加平滑,如下图所示。
众数滤波设置
点击运行完成众数滤波调整,调整后的影像如下图所示。
调整后影像
同样的方法调整2022年影像,调整后效果如下图所示。
调整后效果
清理边界上一步清除了细小的像素,这一步将处理不同分类之间的边界,通常会删除个别像素并用它们周围的像素值替换它们,与众数滤波工具能够实现相似的效果,但将侧重于类边界,让边界变得更加平滑。
在工具箱中点击“Spatial Analyst工具\栅格综合\边界清理”,调用边界清理工具,如下图所示。
调用边界清理工具
在显示的边界清理对话框内,输入栅格为上一步得到的清除分类后的2017年影像,设置输出栅格,后面默认即可,如下图所示。
边界清理设置
点击运行完成2017年影像边界清理,清理后图如下图所示。
清理后影像
同样的方法对2022年影像进行边界清理,清理后影像如下图所示。
清理后影像
计算面积在图层属性内,可以看到每个像素代表30*30,考虑到坐标系为UTM,所以单位是米,即每个像素代表的面积是30*30=900㎡。
图层属性
在属性表内新建一个字段area,用于存放计算得到的面积,如下图所示。
新建字段
在字段计算器内,输入公式!Count! * 900 / 1000000,如下图所示。
计算字段设置
点击应用后得到计算结果,如下图所示。
计算结果
同样的方法为2022年影像计算面积,计算结果如下图所示。
计算结果
从上面的计算结果即可得出2022年鄱阳湖面积比2017年鄱阳湖面积增加了4771.735- 4172.919=598.816平方公里。
结语以上就是如何使用Landsat 8影像计算鄱阳湖面积变化的详细说明,主要包括了对比地图、非监督分类、调整样式、清除分类、清理边界和计算面积等功能。
最后,如果您在工作中遇到与GIS相关的任何难题,请通过关注“水经注GIS”并在后台留言告诉我们,我们将针对比较有典型性和代表性的问题提供解决方案,从而将我们“为您提供丰富的地图数据与专业的GIS服务”之宗旨落到实处!