ArcGIS不同坡度植被覆盖率分析步骤
写在前面的话:这是一个栅格分析的入门级应用……好吧,算虾神我的实战操作科普系列的作品。怎么用ArcGIS做不同坡度的植被覆盖度分析,实际上在ArcGIS...
写在前面的话:这是一个栅格分析的入门级应用……好吧,算虾神我的实战操作科普系列的作品。
怎么用ArcGIS做不同坡度的植被覆盖度分析,实际上在ArcGIS里面做此类分析,已经是hello world级别的应用了,任何一个基础空间分析(栅格分析)的课程,都会有此类案例的操作……所以我本来认为一个百度就已经足够了,但是百度了半天,发现……居然没有任何一个完整的操作案例,所有人都认为这个东东是在太简单了,所以直接就给出几个关键字:坡度分析——重分类——叠加……OK打完收工(真是特简洁明了……)
好吧,要详细步骤,只能我自己来写一个了。
先来看看数据:
首先是30米DEM数据一套,用于做坡度分析。(下图山体阴影渲染效果,具体方法)见以前的文章:
然后就是用遥感分类手段提取出来的植被覆盖矢量数据一份(不要问怎么提取出来的,学计算机的虾神不懂遥感……那些神马NDVI,神马波段,神马监督非监督分类的,虾神统统不懂)。
OK,只要有这两份数据,就可以做不同坡度的植被覆盖的分析了。
首先,要有坡度,那么先利用ArcGIS空间分析工具集里面的坡度分析,来把地形进行坡度分级,坡度的算法原理非常简单,直接通过中学的一个三角函数就能直接算出来:
工具的位置如下:空间分析工具集——表面分析——坡度
工具非常简单,几乎没有任何需要设置的参数:
点击OK之后,就完成了坡度的分析,结果如下:
因为直接计算的坡度,数据太离散,无法进行分区统计,所以接下去我们进行重分类,把坡度分成若干要统计的类别,数据从0度到60度,所以我们按10度一个分类,分成6类,使用的工具,叫做重分类工具,如下:
然后直接设置分类:
设置为六类,每十度为一类
重分类结果如下:
重分类完成之后,就可以利用矢量数据的分析方法来进行叠加分析了,首先把栅格数据转换成矢量(当然,直接用栅格做叠加也行,就是把矢量的植被数据转换成栅格,反正二者之间,的数据维度要进行统一,为了简单起见,我直接用栅格转矢量来做)
在toolbox的转换工具集里面有个从栅格转出,选择转为面要素:
转换结果如下:
接下去,做叠加分析:注意,这里不能用相交,因为地形的数据要大于植被覆盖的数据,如果相交仅是取交集的话,就无法计算覆盖度了(当然,如果要做的不同坡度的植被覆盖面积计算,可以只做相交)。
标识工具位置:ArcToolbox——分析工具集——叠加分析——标识
标识工具是将一份数据的内容,进行空间提取,到另一份数据上面去,如下图所示:
执行标识工具,结果如下:
执行完成之后,数据就变成了这个样子:
如果该位置没有灌木,那么组合地这个字段,就会变成空值,而不是组合地为灌木的,就表示有被覆盖,下面就可以进行统计了:
利用分析工具里面的汇总统计工具来实现统计功能:
将面积字段进行累加,然后用gridcode(坡度分类)和组合地字段为统计汇总字段:
汇总结果如下:
有这个数据,就可以很轻易的做出各种分析图表出来了,比如采用Python进行可视化,结果如下:
打完收工。
开放原子开发者工作坊旨在鼓励更多人参与开源活动,与志同道合的开发者们相互交流开发经验、分享开发心得、获取前沿技术趋势。工作坊有多种形式的开发者活动,如meetup、训练营等,主打技术交流,干货满满,真诚地邀请各位开发者共同参与!
更多推荐
所有评论(0)