r - 如何在堆栈中用 NA 替换值？

Question

您好，我在这里给出一个复制示例。我有三个栅格，我将它们堆叠起来。我想用 NA 替换任何高于 50 的值，但出现错误：

     r <- raster(nrows=10, ncols=10)
     r <- setValues(r, 1:ncell(r))
     r1 <- raster(nrows=10, ncols=10)
     r1 <- setValues(r1, 1:ncell(r))
     r2 <- raster(nrows=10, ncols=10)
     r2 <- setValues(r2, 1:ncell(r))
     St=stack(r,r1,r2)

     NAvalue(St) <- >50
   Error: unexpected '>=' in "NAvalue(St) =  >="

有什么帮助吗？谢谢

Answer 1

> St[9,9]
     layer.1 layer.2 layer.3
[1,]      89      89      89
> St[St >50] <- NA
> St[9,9]
     layer.1 layer.2 layer.3
[1,]      NA      NA      NA