vlambda博客
学习文章列表

R语言-以表格显示分区统计

      在做区域统计时,我们经常需要统计某个区域内的均值或者中位数等信息,这时候可能会有不同的选择,如ArcGIS可以选择“Zonal Statistics as Tabl”,可以实现以矢量面或整型栅格来统计栅格信息。


      但当数据量大或者栅格为浮点型时,比如统计我们想统计降水在300-600mm之间气温的均值、中位数等特征时,或许你会想到首先重分类,之后再分区统计,这固然是一种方法,但当数据量很大时,似乎有些重复,那么用R如何首先呢?
       首先,读取栅格数据,这里我们以2015年全球降水和气温为例:
数据源:( 链接:https://pan.baidu.com/s/1LF9aSJtMGaSuEirfWMOupw   提取码:mj27


R语言-以表格显示分区统计


第一种方法:直接统计

R语言-以表格显示分区统计


第二种方法:结合条件语句判断统计

R语言-以表格显示分区统计


第三种方法:通过重分类来进行统计

R语言-以表格显示分区统计


第四种方法:通过将栅格转化为数据框格式来处理,适合数据量较大的栅格,速度较快。

R语言-以表格显示分区统计


第五种方法:结合自编函数与重分类后的降水栅格数据进行分区统计


同时,以上方法均适用于栅格数据的批处理,只需要将气温数据替换为栅格栈数据即可,以最后一种方法举例:

以上均是自己在学习中用到的一点点小代码,欢迎大家留言区补充指正,感谢!