ArcGIS中的土地利用变化分析

╰+攻爆jí腚メ 2022-06-08 01:21 562阅读 0赞

### 文章目录 ###

*   *   *  一、不同时期的土地利用矢量数据，如何分析其图形及属性变化？
         *  二、如何统计不同土地利用类型的面积？
         *  三、如何统计不同区域各类土地利用类型的面积？
         *  四、如何将行政区划的信息追加到土地利用数据中？
         *  五、如何制作土地利用变化转移矩阵？

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

统计的操作视频，点击[这里][Link 1]观看。  
用地变化分析操作视频，点击[这里][Link 2]观看。  
转移矩阵操作视频，点击[这里][Link 3]观看。

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

![format_png][]  
土地利用图（左图为1993年，右图为2003年）

**思路如下：**

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

**操作步骤如下：**  
   
1、将两个图层的属性放到一个图层里。

使用GP工具Union（联合）。注意JoinAttribute，选择All，也就是结果中包含两个图层的所有字段。  
![format_png 1][]

结果属性表如下：  
![format_png 2][]

Unoin可以将数据的属性合并到一起，同时，图斑不一样的区域也可以拆分开。见下图：  
![70][]

2、新建一个短整型字段。  
![format_png 3][]

3、使用字段计算器，比较两个字段的内容。如果属性没变，就赋值为1，如果属性变化了就赋值为0。这里我们定义isSame()函数，在函数中比较字段值，相同就赋值1，不相同就赋值0。  
![format_png 4][]

Python代码如下：

Pre-Logic Script Code中写：

def isSame(x, y):
        fieldA = str(x)
        fieldB = str(y)
        if fieldA == fieldB:
            return 1
        else:
            return 0

变化中写：（比较的是Class\_Name和Class\_Name1两个字段）

isSame( !Class_Name! , !Class_Name_1! )

> **简单方法：**  
> Python语法中比较字段值是否相等，可以用 == 。注意，中间不能有空格。所以表达式可以简化为：  
> !Class\_Name! == !Class\_Name\_1!

![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2NlaWJha2U_size_16_color_FFFFFF_t_70_pic_center]  
结果如下：  
![format_png 5][]

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

结果如下：  
![format_png 7][]

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

使用GP工具[Summary Statistics（汇总统计数据）][Summary Statistics]，统计各类的总面积。  
![format_png 8][]

结果是一张统计表，如下：  
![format_png 9][]

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

可以使用GP工具[Tabulate Intersection（交集制表）][Tabulate Intersection]进行统计。  
![format_png 10][]

结果如下：  
![format_png 11][]

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

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

一个图斑可能会跨区域，如果图斑需要切分，一部分属于行政区划A，另一部分属于B，可以使用GP工具Intersect（相交），参数中的JoinAttribute，选择All。结果如下：  
![70 2][]

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

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

![70 3][]

结果如下：

![70 4][]

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

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

**思路如下：**

> 1、对两期土地利用数据分别进行融合Dissolve（融合）操作，使每个土地利用类型仅有一条记录；  
> 2、再对融合的结果进行Intersect（相交）操作，分析土地利用变化；  
> 3、转出为Excel表；  
> 4、在Excel中采用透视表进行统计。

**模型如下：**

针对上面的操作思路，制作了一个Model模型，转出结果为Excel。模型内容如下：  
模型下载网址，链接:https://pan.baidu.com/s/1pKNHDeb 密码:n9hy

![format_png 12][]

双击工具箱中此工具，界面如下：  
![format_png 13][]

关于土地利用转移矩阵模型，适用ArcMap 10.2及以上版本。如何打开模型，编辑模型，见下图操作。大家可以参考模型中的工具，自己来制作这个模型。模型构建参考帮助文档：[在模型构建器中执行工具][Link 4]。

![这里写图片描述][70 5]

结果为一个Excel表格，**制作透视图方法如下：**

1、点击数据透视图，然后选中所有区域的数据；  
![format_png 14][]

2、先勾选第一期的土地利用类型字段；  
![format_png 15][]

3、再将第二期的土地利用类型字段拖拽到透视图的第一行位置；  
![format_png 16][]

4、最后将面积字段拖拽到透视图的中间空白位置；  
![format_png 17][]

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

关于制图及出图，可以参考[ArcGIS制图及出图小技巧——以土地利用图为例][ArcGIS]。

更新了[土地整理中的小案例][Link 5]，欢迎查看。

欢迎关注公众号，之后会发操作视频哦。

![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2NlaWJha2U_size_16_color_FFFFFF_t_70_pic_center 1]

[Link 1]: https://mp.weixin.qq.com/s?__biz=MzU2MjY1NzIzNg==&mid=2247484250&idx=2&sn=a368baaf164264ccf3fc83c18842722f&chksm=fc6761e1cb10e8f78f9276ffe55a5bd88fadf6d6210b74db6c957d65087271d1a3637f01e8cc&token=551350061&lang=zh_CN#rd
[Link 2]: https://mp.weixin.qq.com/s?__biz=MzU2MjY1NzIzNg==&mid=2247484292&idx=2&sn=405fa71a1dd1bec4f89d712bccd0ef52&chksm=fc67613fcb10e82916984cd878c7a6eacb18e08cf1ae388bbb5c011881244ccbc13540854c01&token=1172524420&lang=zh_CN#rd
[Link 3]: https://mp.weixin.qq.com/s?__biz=MzU2MjY1NzIzNg==&mid=2247484292&idx=3&sn=e893c361c64218d82e588fede8606b68&chksm=fc67613fcb10e829d41c0717ecac0d9648088f726c481cc7f7c1e774a5ac74861e83ef4d01f6&token=1172524420&lang=zh_CN#rd
[format_png]: /images/20220608/a01c31357c344342aa54fc53f149f97a.png
[format_png 1]: /images/20220608/b62262511a604c2b9b835246a280499a.png
[format_png 2]: /images/20220608/66e372d868ce4a0fb26b0ec2e0c42acc.png
[70]: /images/20220608/4b9709c7fc114210ba630c07ab9e3621.png
[format_png 3]: /images/20220608/ab96459f28e74de9845c22b6e7786e13.png
[format_png 4]: /images/20220608/4adcee6641c7401f9e5e6e5e84341292.png
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2NlaWJha2U_size_16_color_FFFFFF_t_70_pic_center]: /images/20220608/4f9bdff851a24c7c997e2068479963f2.png
[format_png 5]: /images/20220608/8e27e79836cc45caac8d6e9be0a814d8.png
[format_png 6]: /images/20220608/1827008b37e64426b0095d336fcced20.png
[format_png 7]: /images/20220608/ee05fc0a108b4cf3bbfd7601746f0d43.png
[Summary Statistics]: http://desktop.arcgis.com/zh-cn/arcmap/latest/tools/analysis-toolbox/summary-statistics.htm
[format_png 8]: /images/20220608/c20fd8168e2b48e79e1a85403e9c72e8.png
[format_png 9]: /images/20220608/9360a77cf201485bb55cde4a85da381e.png
[Tabulate Intersection]: http://desktop.arcgis.com/zh-cn/arcmap/latest/tools/analysis-toolbox/tabulate-intersection.htm
[format_png 10]: /images/20220608/5bd380464e6c4e69b69c574cc71aea9f.png
[format_png 11]: /images/20220608/2daa7291a6b146508494d0c25a0ef5d8.png
[70 1]: /images/20220608/a789ea28459642ec9e45417aaf05eaf6.png
[70 2]: /images/20220608/f7b5a17661c94ba192e2d44f930a3da8.png
[70 3]: /images/20220608/50a993c6faef4285a55a61a59c306f15.png
[70 4]: /images/20220608/13519304da5b416caf49820ce102170a.png
[format_png 12]: /images/20220608/31655ba6fb3b45c9989fe9d24953a07c.png
[format_png 13]: /images/20220608/a733901a4b98498aa856faf230fb67b5.png
[Link 4]: http://desktop.arcgis.com/zh-cn/arcmap/10.3/analyze/modelbuilder/executing-tools-in-modelbuilder-tutorial.htm
[70 5]: /images/20220608/497ce3a44bf84959a60654eeeaa91cdf.png
[format_png 14]: /images/20220608/c6d6ef0f181e41c49dbe84101cb27862.png
[format_png 15]: /images/20220608/5f1f2b12e53c414e99fd2e2eeb39c4b6.png
[format_png 16]: /images/20220608/6d947e9c3cba4ffba220c128a1be7116.png
[format_png 17]: /images/20220608/06f31a16973e4cb9b47aa27ac100d681.png
[ArcGIS]: https://blog.csdn.net/ceibake/article/details/79843986
[Link 5]: https://blog.csdn.net/ceibake/article/details/102916255
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2NlaWJha2U_size_16_color_FFFFFF_t_70_pic_center 1]: /images/20220608/ca8828b43b1f46dfabfd599a3e2b80fc.png