将栅格中的值设置为间隔内的NA

Question

假设我在R中有一个栅格，例如：

r <- raster(ncols=10, nrows=10)
r[] <- sample(50, 100, replace=T)

如何将[10,30]之间的值设置为NA？

我尝试过values = r[r<= c(10,30)] = NA但这不会删除10到30之间的值。

提前致谢。

Answer 1

您必须找到这样的值：

r[r>10 & r<30]<-NA 
r[]
  [1] NA 38 NA  1 NA NA  3 32 40 36 NA NA 30  5  2 38 47 NA 42 42  1 NA NA 32 43 NA  7 NA  8 35 NA NA NA 48 10
 [36] 32 49 33 NA 48 NA 37  2 45  9  7 37 42  2 42 NA  3 49 48 NA NA 48  5 NA 46 43 NA NA NA 42 39 41 NA 48 NA
 [71] NA NA 30 35 32 46  4 32 NA 48 40  2 44 45 NA NA NA  9 47 NA NA 41  2  4 42 30 NA 36 48 32

Answer 2

使用

r[r <= 10:30] <- NA

似乎很危险，实际上它发出警告：

警告消息：在getValues（e1）<= e2中：较长的对象长度不是较短的对象长度的倍数

对于大型栅格，它也不是内存友好的。

正确/最好的方法是使用raster::reclassify ：

r <- raster::reclassify(r, c(10, 30, NA))

，这也快得多。

也可以看看：

https://stackoverflow.com/a/49159943/6871135