Arcgis 预处理和获取 30 米二级分类土地利用数据

最编程 2024-04-14 20:40:47

...

本篇以武汉市为例，主要介绍将土地利用数据转换成武汉市内各区土地利用详情的过程以及分区统计每个区内各地类面积情况，后面还有制作过程中遇到的面积制表后数据过小的解决方法以及一些相关的知识点：

示例数据下载链接：数据下载链接

（1）首先在Arcgis导入湖北30m精细地表土地利用数据，然后导入武汉各区矢量行政边界对土地利用数据进行掩膜提取。（ArcToolbox-Spatial Analyst-提取分析-按掩膜提取）

编辑

掩膜提取窗口第一项为原始土地利用数据，输入栅格集为武汉市各区行政区划图，输出栅格保存位置可以按实际需求放置，设置好后点击确定。提取结果如下：

编辑

（2）接下来利用提取好的武汉市土地利用图来进行各区面积统计。（ArcToolbox-Spatial Analyst-区域分析-面积制表）

编辑

面积制表窗口第一行为武汉市各区行政区划，区域字段选县级市，输入栅格数据为上面裁剪好的武汉市土地利用图，类字段选择vaule，输出位置*选择，处理像元大小默认。

编辑

然后按照实验需求导出表格就可以了，以上统计武汉市各行政区内各土地利用类型的面积就完成了，下面将介绍处理过程中遇到的问题：

（1）面积制表结果的数值过小，重采样显示数值过大和过小

编辑

这时候需要先对图像进行投影，然后进行重采样成30×30：

投影：ArcToolbox-数据管理-投影栅格

编辑

重采样：ArcToolbox-栅格-栅格处理-重采样

编辑

分辨率大小即设置X,Y,这个根据实验需求来设置，这里重采样为30×30

采样结果可以打开图层属性查看（像元大小一列变为30,30）

编辑

另附上获取30m土地利用数据获取来源：

1、地理遥感生态网

网站地址www.gisrs.cn

地理遥感生态网公布了1980年-2020年土地利用分类数据，精度有10m、30m、1000m。精度据官网介绍达到了92%精度。

编辑

（1）2000年全球土地覆盖计划（GLC2000）

（2）欧空局全球陆地覆盖数据（ESA GlobCover）

网站地址www.****.store

UMd基于AVHRR数据的5个波段及NDVI数据经过又一次组合建议数据矩阵，用分类树的方法进行了全球土地覆盖分类工作。其目的是希望建立一个比过去数据更高精度的数据集。

编辑

时间范围：2001-2020，MCD12Q1 V6 产品每年提供来自六种不同分类方案的全球土地覆盖类型。它是使用 MODIS Terra 和 Aqua 反射率数据的监督分类得出的。然后，受监督的分类经过额外的后处理，结合了先验知识和辅助信息，以进一步细化特定类别。

编辑