【ArcGIS教程】土地利用转移矩阵及土地利用数据获取

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

第一步 准备数据

我们使用的数据为河南省2000年与2010年土地利用类型栅格数据（已经与土地利用分类表格连接。PS：如何连接分类等级表格详见此教程======）

（1）加载2000年河南省土地利用类型数据

对2000年数据符号化显示。打开图层属性，在【符号表示法】中按【唯一值表示】，点击选择Value Filed字段为【名称】，显示即可

2000年河南省土地利用类型一级分类显示如下

查看属性表，这里我们选用【名称】字段进行下序的操作

（2）加载2010年河南省土地利用类型数据

打开2000年土地利用类型图层属性，按下面显示顺序操作（这样可以导入如2000年图层符号化一致的颜色显示）

符号化显示如下

查看2010年河南省土地利用类型数据属性表

第二步 栅格转面

（1）对2000年土地利用类型数据栅格转矢量

运行ArcToolbox，打开【转换工具】，选择【从栅格转出】里面的【栅格转面工具】，调出面板进行参数设置。输入栅格选择裁剪的2000年土地利用类型河南省栅格数据；字段选择【名称】；不要勾选【可选】参数（效果较好）

2000年河南省土地利用栅格数据转化成矢量数据成功，显示如下

查看属性表，发现被分成了很多个图层（需要融合属性）

（2）对2000年土地利用类型数据栅格转矢量

使用同样方法对2010年河南省栅格土地利用类型数据进行转化矢量

查看其属性表，显示和2000年矢量数据一样，都被分成了很多个图层

第三步 融合

（1）2000年矢量土地利用类型属性融合

运行ArcToolbox，打开【数据管理工具】，选择【制图综合】里面的【融合工具】，调出面板进行参数设置。输入要素选择2000年土地利用类型河南省矢量数据；融合字段选择【名称】；勾选【创建多边形要素】参数（这样的融合效果较好）

2000年矢量土地利用类型融合完成，如下（属性表只有六大类，但是有些字段显示不太全，我们可以新创建一个字段NAME）

修改字段(可选)

符号化显示

(2)2010年矢量土地利用类型属性融合

方法如2000年数据一致

符号化显示

第四步 变化操作

运行ArcToolbox，打开【分析工具】，选择【叠加分析】里面的【相加工具】，调出面板进行参数设置。输入要素选择2000年的河南省矢量数据与2010年的河南省矢量数据；其它参数默认即可

【相交】操作完成，显示如下

查看属性表

这里我们新建一个字段【面积】，面积类型为双精度。从而进行变化面积的统计

打开几何计算面板进行计算面积，单位选择平方千米

变化面积计算成功，显示如下

但是有的时候计算面积时显示面积计算被禁用（需要为图层选择投影）
打开图层属性，在坐标系系统中选择 阿尔伯特或者 兰伯特投影即可（中心经度与第一纬度、第二纬度可根据地区经纬度自行选取合适数值）

第五步 表转Excel

运行ArcToolbox，打开【转化工具】，选择【Excel】里面的【表转Excel】工具，调出面板进行参数设置。输入表选择进行相交操作的数据；其它参数默认即可

之后我们在保存位置打开刚才的文件Excel

字段数据显示如下

可以删除一些不需要的字段要素

选中数据，点击Excel中的【插入】按钮，选中【数据透视表】，之后点击确定即可

在数据透视表字段面板中对所有要素打上√，将较前的年份数据拖至【列】下面

点击值调出值字段设置面板，选择求和操作，确定即可

表格显示如下。那些右对角线数据表示此地类型没有发生变化

第六步 查看变化

添加一个【变化】字段，类型为文本（用于显示土地利用类型变化）

在字段计算器中进行计算

表格显示如下，可以明显看出由什么地类变为什么地类

对【Change】字段符号化显示

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

1.地理科学生态网

 网站地址www.csdn.store

最近更新了全国30m分辨率数据1980-2020年连续年份土地利用数据，分类为一级分类（耕地林地草地水体建设用地未利用地等）；据同学下载实验，和landsat数据进行了对比，发现精度还是不错的。

并且有各种地理科学数据，竟然发现了很多宝藏数据。比如道路网数据、NPP数据、土地利用数据（二级分类）、乡镇边界、NDVI数据、地铁站点、气象数据（降雨量、气温、蒸散量、辐射、湿度、日照时数、风速）、径流量数据、30米土地利用二级分类数据、12.5米高精度DEM高程、2020年行政区划、10、30米NDVI数据、30m精度NPP净初级生产力、夜间灯光数据、统计年鉴、道路网、水系、景区医院等poi、三级流域矢量边界、地质灾害分布数据、30m土壤理化性质数据集、30mGDP柵格数据、30m人口柵格数据、30m精度植被类型数据、30m精度连续年份土地利用数据、GPP初级生产力数据、农田作物类型分布数据、自然保护区分布数据、高精度遥感影像数据、1:10万沙漠沙地分布矢量数据、蒸散发数据、降雨量气温蒸发数据、地表径流量、水汽压、土壤侵蚀数据等等！

2.地理空间数据云

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

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

3.地理遥感生态网

 网站地址www.gisrs.cn

该网站同样更新了很多有关地理的数据，30米土地利用二级分类数据、12.5米高精度DEM高程、2020年行政区划珍贵的乡镇边界、村边界、气象数据(降雨量、气温)、径流量数据等。

4.马里兰大学数据集

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

5.MCD12Q1.006 MODIS Land Cover Type

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

​

6.IGBP（国际地圈生物圈计划）

在大多的重要的、普遍使用的陆地数据集、土地覆盖数据不能满足一致的、长期的。如今的全球土地覆盖数据不管是在满足全球变化研究的须要方面,还是在满足各国或国际可行的资源管理规划方面都是不足的。在1992年法国图卢兹召开的IGBP工作会议上，讨论了IGBP核心科学计划对全球数据的需求，在IGBP-DIS的领导下对眼下的数据进行了评估，最后将LUCC作为其几个核心研究计划之首。与此同一时候。做出了用AVHRR数据制作全球1公里土地覆盖特征数据的建议，IGBP-DIS成立了土地覆盖工作组(LCWG)。启动了基于AVHRR的全球1公里土地覆盖数据项目，项目由U.S. Geological Survey (USGS),University of Nebraska-Lincoln (UNL), 内布拉斯加州林肯大学和欧洲联合研究中心承担。其目标是：採用USGS的方法，利用1992年四月到1992年三月的AVHRR数据开发出1km分辨率的全球土地覆盖数据集。分类系统採取IGBP制定的分类系统，把全球分为17类。

7.清华大学10m/30m数据下载

网站地址http://data.ess.tsinghua.edu.cn/