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

Posted 2021-04-18 栅格空间分析

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

      但当数据量大或者栅格为浮点型时，比如统计我们想统计降水在300-600mm之间气温的均值、中位数等特征时，或许你会想到首先重分类，之后再分区统计，这固然是一种方法，但当数据量很大时，似乎有些重复，那么用R如何首先呢？

       首先，读取栅格数据，这里我们以2015年全球降水和气温为例：

数据源：（ 链接：https://pan.baidu.com/s/1LF9aSJtMGaSuEirfWMOupw   提取码：mj27）

第一种方法：直接统计

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

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

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

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

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

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