您好我在这里给出了一个重复示例。我有三个栅格,我把它们堆叠起来。 我想用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) = >="
有任何帮助吗?谢谢
答案 0 :(得分:5)
> 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