基于ArcGIS统计格网中不同土地利用类型占比
- 数据准备
- ArcGIS操作步骤
-
- [1、创建渔网(Create Fishnet)](#1、创建渔网(Create Fishnet))
- 2、建立唯一标识
- 3、选择格网
- 4、提取不同类别土地利用类型
- 5、各类用地面积计算
- 参考
- 另:可能出现的问题总结
-
- [Q1:ArcGIS获取唯一值(Get Unique Value)失败](#Q1:ArcGIS获取唯一值(Get Unique Value)失败)
数据准备
1、研究区矢量边界
2、土地利用数据
数据下载-The 30 m annual land cover datasets and its dynamics in China from 1985 to 2022
- 空间位置:中国
- 地理坐标系: D_WGS_1984
- 投影坐标系: Albers_Conic_Equal_Area
- 数据格式:栅格(.tif)
- 空间分辨率:30*30米
土地利用类型: 9类
编号 | 类别 |
---|---|
1 | Cropland |
2 | Forest |
3 | Shrub |
4 | Grassland |
5 | Water |
6 | Snow/Ice |
7 | Barren |
8 | Impervious |
9 | Wetland |
ArcGIS操作步骤
1、创建渔网(Create Fishnet)
通过绘制渔网我们可以统计出网格所占要素的多少,并符号化显示。
运行ArcToolbox,选择【数据管理工具】,点击【样例】中的【创建渔网】工具,进行参数设置。
范围选择流域研究区;像元宽度与高度选择1 km×1 km;几何类型选择面
渔网创建完成如下:
2、建立唯一标识
创建完渔网后,为了确保每个网格有一个固定标识,将其生成的ID字段进行计算,利用字段计算输入下图所示公式:
3、选择格网
图中有多余的网格中并没有我们需要计算的土地利用类型要素,所以应该提取对应图层范围的渔网。
点击界面的【选择】---【按位置选择(Select By Location)】,打开此属性面板
【按位置选择(Select By Location)】,选择与研究区边界相交 的渔网。
目标图层选择渔网,而来源图层我们选择研究区,之后应用就可以看到需要的渔网范围已经被选中,之后导出渔网数据即可
右键点击渔网,在调出的面板中选中【导出】---【导出数据】,弹出的面板中保存输出路径,存为shp文件,ok即可
渔网范围被提取出来显示如下,
修改符号化。点击提取渔网,将其颜色改为无色即可。
4、提取不同类别土地利用类型
运行ArcToolbox,打开【空间分析工具】,选择【提取】工具的【按属性提取】工具,调出面板后进行参数设置。输入栅格选择研究区;打开语法建立面板,在里面进行语法设置。这里以提取耕地(Cropland)为例!("编号"= 1)
举例1:输入SQL代码即可:Value IN ( 1 , 3 , 5 , 10 , 12 , 17),含义:提取Value等于1、3、5、10、12、17的对象。
提取结果如下:
其它各类土壤类型提取原则相同。
5、各类用地面积计算
运行ArcToolbox,打开【空间分析工具】,选择【邻域分析】里面的【表格显示统计方法】调出界面,进行参数设置。输入栅格或邻域要素数据选择提取的范围渔网,邻域字段选择具有标识的FID,输入用途栅格为提取的耕地要素;设置输出路径,应用即可
耕地面积统计表导出成功,我们打开查看属性表的数据,显示如下,每一个FID字段都有相应的面积、数量等要素统计值。
需要说明的是统计得到FID与渔网数并不相同,当渔网内数据为空时,不做统计。
使用同样的方法计算其它土地利用类型要素的统计表。
参考
1、知乎-【ArcGIS教程】(66)利用渔网统计不同土地利用类型面积
另:可能出现的问题总结
Q1:ArcGIS获取唯一值(Get Unique Value)失败
按属性提取时,,SQL代码咋写Verify都是错误的,界面如下:
推测可能是土地利用数据有点问题吧。土地利用分类对应的应该是直方图,IRasterBand.Histogram,IRasterHistogram.Counts属性会获得一个256长的数组,里面保存了每个值的count。因此可见这里只能针对0-255的整形栅格,浮点栅格就比较难了,唯一值太多arcmap也崩溃。
**解决方案:**将土地利用类型数据更改为唯一值,再执行按属性提取时,可用