我想用150个光栅对象的频率写一个df。
我知道我可以阅读单个栅格对象
raster_x <- raster::raster()
我可以进一步获得频率
raster_freq_y <- raster::freq(raster_x)
之后，我可以将多个栅格对象的频率输出绑定到df
cbind.data.frame(raster_freq_x, raster_freq_y)
我不知道如何一次性为150个光栅对象做这个？
我应该使用循环吗？如果是这样，什么样的循环才有意义？
任何帮助表示赞赏。
如果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)