我想结合两个(或更多)光栅砖。这是我想做的一个例子
library(raster)
r <- raster(ncol=40, nrow=20)
r[] <- rnorm(n=ncell(r))
b1 <- brick(x=c(r, r*2, r))
b2 <- brick(x=c(r, r*3, r))
b <- brick(b1, b2)
此示例返回以下错误。
Error in (function (cl, name, valueClass) : assignment of an object of class “RasterBrick” is not valid for @‘nlayers’ in an object of class “.MultipleRasterData”; is(value, "integer") is not TRUE
如果我按照以下方法创建堆栈而不是砖块,则>]
b1 <- stack(x=c(r, r*2, r)) b2 <- stack(x=c(r, r*3, r)) b <- brick(b1, b2)
我收到以下错误消息
Error in if (values) { : argument is not interpretable as logical
我想结合两个(或更多)光栅砖。这是我要执行的示例library(raster)r
您需要stack
将具有相同范围和分辨率的多个栅格,堆栈或图块组合到单个多层对象中。然后,您可以转换为brick