r语言 - 栅格砖中所选图层的子分配

Question

我想通过将栅格乘以另一个栅格来修改栅格砖中的图层子集。

例如，如果我们有一个名为"r.brick"的光栅砖，我们尝试将其图层 2：4 乘以具有相同行和列尺寸的栅格"r.mult"：

• r.brick[[2:4]]按预期返回层 2：4

• r.brick[[2:4]] * r.mult如预期的那样成功地增加了这些层

但是，如果我尝试将结果分配回子集图层，则会出现错误

r.brick[[2:4]] = r.brick[[2:4]] * r.mult
# Error in value[] <- val : 
#   incompatible types (from S4 to double) in subassignment type fix

错误消息表明分配正在尝试分配栅格值，而不是栅格本身。 但是如果我尝试使用 getValues 进行作业，则会出现不同的错误：

r.brick[[2:4]] = getValues(r.brick[[2:4]] * r.mult)
# Error in .local(x, values, ...) : length(values) is not equal to ncell(x)

正确的方法是什么？

一些可重复的数据：

library(raster)
r.list = vector("list", 20)
set.seed(123)
for (i in 1:20) {
  r.list[[i]] = raster(vals=runif(100), nrows=10, ncols=10, ext=extent(c(0,25,0,25)))
}
r.brick = brick(r.list)
r.mult = raster(vals=sample(2,100,T), nrows=10, ncols=10, ext=extent(c(0,25,0,25)))

Answer 1

然后我只为您提供了一个解决方法（使用循环）：

layers <- 2:4
for(i in layers) {
  r.brick[[i]]  <- r.brick[[i]] * r.mult
}

注意：显然，具有[]子集的任务仅适用于单个图层。

Answer 2

我认为这是一个缺失的功能。感谢您指出这一点。我会按照 maRtin 进行循环（也许在创建 RasterStack 之后，这可能会更有效）。如果数据集不是太大，你可以

# example data
library(raster)
b <- brick(nrows=2, ncols=2, nl=6)
values(b) <- rep(1:4, 6) 
r.mult <- raster(vals=10, nrows=2, ncols=2)

values(b)[,3:4] <- values(b[[3:4]] * r.mult)
# values(b)

Answer 3

以下代码包含具有 r.mult 值的 r.brick 栅格值的多个图层 2-4，并将结果分配给 r.brick 图层 2-4。

attr(attr(r.brick, 'data'), 'values')[,2:4] = attr(attr(r.brick, 'data'), 'values')[,2:4] * attr(attr(r.mult, 'data'), 'values')

或

attr(attr(r.brick, 'data'), 'values')[,2:4] = getValues(r.brick[[2:4]] * r.mult)