我想用150个栅格对象的频率写一个df。
我知道我可以使用
读取单个栅格对象raster_x <- raster::raster()
我可以通过
进一步获得频率raster_freq_y <- raster::freq(raster_x)
之后，我可以使用
将多个栅格对象的频率输出绑定到df上。cbind.data.frame(raster_freq_x, raster_freq_y)
我不知道如何一次性处理150个栅格对象？
我应该使用循环吗？如果是这样，什么样的循环才有意义？
感谢您的帮助。
答案 0 :(得分：1)
如果RasterLayer对象具有相同的范围和分辨率，则可以将它们合并到RasterStack中。以下示例来自?freq
示例数据：
library(raster)
r <- raster(nrow=18, ncol=36)
r[] <- runif(ncell(r))
r[1:5] <- NA
s <- stack(r, r*2, r*3)
解决方案：
freq(s, merge=TRUE)
如果RasterLayer对象具有相同的范围和分辨率，则可以将它们放在一起并使用lapply
ss <- list(r, r*2, r*3)
lapply(ss, freq)