如何将栅格中的多个值替换为 NA？

Question

我有多个.tiff图像，例如

20000101
20000113
20000124
.
.
.
20201224
(format is yyyymmdd)

栅格文件可能包含 65533、65534 和 65535 个值作为缺失数据。我想用 NA 替换这些值，因为稍后我需要计算这些文件的平均值，丢失的数据可能会导致问题。

我正在使用以下代码来实现相同的目的。代码是否正确？

r1 <- list.files(path="D:/Test/", pattern=".tiff",full.names = TRUE)
s1 <- stack(r1)
x1 <- clamp(s1, useValues=FALSE)
y1 <- reclassify(s1, cbind(-Inf, 65535, NA), right=FALSE)

如何在同一代码中将 65533 和 65534 值也转换为 NA。另外，如何导出缺失值标记为 NA 的新栅格？

我在R里做

Answer 1

我假设所有实际值都（低得多）。在那种情况下，你可以做

library(terra)
ff <- list.files(path="D:/Test/", pattern=".tiff$", full.names=TRUE)
r <- rast(ff)
x <- clamp(r, -Inf, 65500, values=FALSE)

否则，你可能会这样做

x <- classify(r, cbind(65533:65535, NA))