如何使用raster::aggregate函数按标准差聚合栅格数据集?
栅格聚合标准差实现方案
- 针对你的需求,
raster包的aggregate()函数中没有标准差对应的缩写参数,直接使用R基础函数sd作为fun的参数值即可,和mean的用法一致。 - 修改后的代码如下:
library(raster) merge1 <- raster("temp_merge1.tif") red1 <- aggregate(merge1, fact=16, fun=sd, expand=TRUE, na.rm=TRUE)
- 说明:
sd是R内置的标准差计算函数,aggregate()会自动调用它完成栅格聚合计算;na.rm=TRUE参数会确保计算时忽略栅格中的NA值,和你之前用mean时的逻辑保持一致。
内容的提问来源于stack exchange,提问作者Jessica
相关产品推荐
相关产品推荐

