如何通过条件将所有值分配给NA？

Question

library(raster)
r1 <- r2 <- r3 <- raster(ncol=10, nrow=10)
r1[] <- runif(ncell(r1))
r2[] <- runif(ncell(r2)) / 2
r3[] <- runif(ncell(r3)) * 1.5
s <- stack(r1, r2, r3)


 r11 <- r22 <- r33 <- raster(ncol=10, nrow=10)
 r11[] <- runif(ncell(r11))
 r22[] <- runif(ncell(r22)) / 2
 r33[] <- runif(ncell(r33)) * 1.5
 g <- stack(r11, r22, r33)

如果g中的任何像素是＆gt; = 5，则可以在s = na中进行相应的像素，

可以这样完成：

  s[g >= 5]= NA

，但此像素的所有时间序列s = na

Answer 1

让我们用“ Terra”（替换“栅格”）来做到这一点。

示例数据

library(terra)
set.seed(0)
s <- rast(ncol=10, nrow=10, nlyr=3, vals=rep(1:100, 3))
g <- rast(ncol=10, nrow=10, nlyr=3, vals=sample(8, 300, replace=TRUE))

将s中的所有值设置为na，如果g中的单元格大于或等于其任何层中的5个。

g5 <- any(g >= 5)  
x <- mask(s, g5, maskvalue=1)