Question

我有一个非常大的栅格文件，其中包含值范围从1到25的像元。我想创建25个单独的栅格，它们仅包含原始值相同的像元。

我已经研究了r的“ raster”包，并认为RasterFromCells可能是我的答案，但是它需要像元编号，因为我的原始栅格有171,681,006个像元，每个像元都提供了每个像元的像元编号将非常耗时。

r中是否有一个函数可以用来自动创建一个包含所有值为1的像元的新栅格？

谢谢。

Answer 1

没有可复制的示例，这可能不是您想要的，但这有帮助吗？

library(raster)
logo <- raster(system.file("external/rlogo.grd", package="raster")) 
logo.1 <- (logo == 1)
plot(logo.1)