查看原文
其他

ArcGIS中的土地利用变化分析详解

ceibake GIS前沿 2021-06-14


本篇主要是针对矢量数据的分析。


一、不同时期的土地利用矢量数据,如何分析其图形及属性变化?


土地利用图(左图为1993年,右图为2003年)


思路如下:


可以先对2个图层进行Union操作,然后在结果中新建一个字段,分析属性的变化。如果属性没变,就赋值为1,如果属性变化了就赋值为0。



操作步骤如下:

1、将两个图层的属性放到一个图层里。


使用GP工具Union(联合)。注意JoinAttribute,选择All,也就是结果中包含两个图层的所有字段。


结果属性表如下:


Unoin可以将数据的属性合并到一起,同时,图斑不一样的区域也可以拆分开。更多干货敬请关注:GIS前沿  见下图:


2、新建一个短整型字段。


3、使用字段计算器,比较两个字段的内容。如果属性没变,就赋值为1,如果属性变化了就赋值为0。


Python代码如下:


1 def isSame(x, y):
2    fieldA = str(x)
3   fieldB = str(y)
4   if fieldA == fieldB:
5      return 1
6 else:
7     return 0


变化中写:(比较的是Class_Name和Class_Name1两个字段)


isSame( !Class_Name! , !Class_Name_1! )


如果想知道变化的图斑是从什么类型转为什么类型,可以新建一个文本型字段,将土地利用类型的两个字段赋值过去。


结果如下:


二、如何统计不同土地利用类型的面积?


使用GP工具Summary Statistics(汇总统计数据),统计各类的总面积。


结果是一张统计表,如下:


三、如何统计不同区域各类土地利用类型的面积?


可以使用GP工具Tabulate Intersection(交集制表)进行统计。


结果如下:


四、如何将行政区划的信息追加到土地利用数据中?


有一份行政区划矢量数据,希望知道每个图斑隶属于哪个行政区划。数据可以是行政区域,也可以是坡度及其他数据。


一个图斑可能会跨区域,如果图斑需要切分,一部分属于行政区划A,另一部分属于B,可以使用GP工具Intersect(相交),参数中的JoinAttribute,选择All。结果如下:


如果不希望图斑被切分,在属性表中记录属于A、B、C区,可以使用GP工具Spatial Join(空间连接)。


在区域名字段上右键,合并规则选择连接,其他设置默认就可以了。



结果如下:



五、如何制作土地利用变化转移矩阵?


对于分类后的土地利用栅格数据,在ENVI里很容易计算土地利用变化矩阵。针对矢量数据,我们可以在ArcGIS中统计数据,然后在Excel中制作矩阵。


思路如下:


1、对两期土地利用数据分别进行融合Dissolve(融合)操作,使每个土地利用类型仅有一条记录;

2、再对融合的结果进行Intersect(相交)操作,分析土地利用变化;

3、转出为Excel表;

4、在Excel中采用透视表进行统计。


模型如下:


针对上面的操作思路,制作了一个Model模型,转出结果为Excel。模型内容如下:

模型下载网址,链接:https://pan.baidu.com/s/1pKNHDeb 密码:n9hy



双击工具箱中此工具,界面如下:


关于土地利用转移矩阵模型,适用ArcMap 10.2及以上版本。如何打开模型,编辑模型,见下图操作。大家可以参考模型中的工具,自己来制作这个模型。




结果为一个Excel表格,制作透视图方法如下:


1、点击数据透视图,然后选中所有区域的数据;


2、先勾选第一期的土地利用类型字段;


3、再将第二期的土地利用类型字段拖拽到透视图的第一行位置;


4、最后将面积字段拖拽到透视图的中间空白位置;


再修改下字头内容,添加下边框,土地利用转移矩阵就做好了。


- END -

一步步详解ArcGIS中制作专题地图(附练习数据下载)

我把谷歌地图、高德地图、天地图、Bing....都装进了ArcGIS
ArcGIS将SHP的等高线数据生成CASS标准的高程点
公开课—ArcGIS坐标系统与投影变换(附课件PPT下载)
推荐 | ArcGIS地图制图教程——超详细(附练习数据下载)

    您可能也对以下帖子感兴趣

    文章有问题?点此查看未经处理的缓存