arcgis中python坡度计算_ArcGIS不同坡度植被覆盖率分析步骤

原标题：ArcGIS不同坡度植被覆盖率分析步骤

写在前面的话：这是一个栅格分析的 入门级应用……好吧，算虾神我的实战操作科普系列的作品。

怎么用 ArcGIS做不同坡度的植被覆盖度分析，实际上在ArcGIS里面做此类分析，已经是hello world级别的应用了，任何一个基础空间分析（ 栅格分析）的课程，都会有此类案例的操作……所以我本来认为一个百度就已经足够了，但是百度了半天，发现……居然没有任何一个完整的操作案例，所有人都认为这个东东是在太简单了，所以直接就给出几个关键字： 坡度分析——重分类——叠加……OK打完收工（真是特简洁明了……）

好吧，要详细步骤，只能我自己来写一个了。

先来看看数据：

首先是 30米DEM数据一套，用于做坡度分析。（下图山体阴影渲染效果，具体方法）见以前的文章：

利用Python实现山体阴影效果渲染

然后就是用遥感分类手段提取出来的植被覆盖矢量数据一份（不要问怎么提取出来的，学计算机的虾神不懂遥感……那些神马NDVI，神马波段，神马监督非监督分类的， 虾神统统不懂

）。

OK，只要有这两份数据，就可以做不同坡度的植被覆盖的分析了。

首先，要有坡度，那么先利用ArcGIS空间分析工具集里面的坡度分析，来把地形进行坡度分级，坡度的算法原理非常简单，直接通过中学的一个三角函数就能直接算出来：

工具的位置如下： 空间分析工具集——表面分析——坡度

工具非常简单，几乎没有任何需要设置的参数：

点击OK之后，就完成了坡度的分析，结果如下：

因为直接计算的坡度，数据太离散，无法进行分区统计，所以接下去我们进行 重分类，把坡度分成若干要统计的类别，数据从0度到60度，所以我们按10度一个分类，分成6类，使用的工具，叫做重分类工具，如下：

然后直接设置分类：

设置为六类，每十度为一类

重分类结果如下：

重分类完成之后，就可以利用矢量数据的分析方法来进行 叠加分析了，首先把栅格数据转换成矢量（当然，直接用栅格做叠加也行，就是把矢量的植被数据转换成栅格，反正二者之间，的数据维度要进行统一，为了简单起见，我直接用栅格转矢量来做）

在toolbox的转换工具集里面有个从栅格转出，选择转为 面要素：

转换结果如下：

接下去，做叠加分析：注意，这里不能用相交，因为地形的数据要大于植被覆盖的数据，如果相交仅是取交集的话，就无法计算覆盖度了（当然，如果要做的不同坡度的植被覆盖面积计算，可以只做相交）。

标识工具位置： ArcToolbox——分析工具集——叠加分析——标识

标识工具是将一份数据的内容，进行空间提取，到另一份数据上面去，如下图所示：

执行标识工具，结果如下：

执行完成之后，数据就变成了这个样子：

如果该位置没有灌木，那么组合地这个字段，就会变成空值，而不是组合地为灌木的，就表示有被覆盖，下面就可以进行统计了：

利用分析工具里面的 汇总统计工具来实现统计功能：

将面积字段进行累加，然后用gridcode（坡度分类）和组合地字段为统计汇总字段：

汇总结果如下：

有这个数据，就可以很轻易的做出各种分析图表出来了，比如采用Python进行可视化，结果如下：

责任编辑：