R语言：如何将栅格多值替换为NA并导出至指定目录？

现有代码是否正确？ #

不正确，存在两个明显问题：

• x1 <- clamp(s1, useValues=FALSE)这行完全多余：既没有指定clamp必需的上下限参数，后续代码也未用到x1变量；
• reclassify规则逻辑错误：cbind(-Inf, 65535, NA), right=FALSE会把所有小于65535的值替换成NA，和你要保留有效数值、仅替换65533/65534/65535的需求完全相反。

如何修改代码替换三个缺失值为NA？ #

推荐两种直观的实现方式：

方式1：直接索引赋值替换 #

利用栅格对象的索引特性，直接定位指定值并替换为NA：

# 读取文件并创建栅格栈（修正正则匹配，确保只选.tiff后缀文件）
r1 <- list.files(path="D:/Test/", pattern="\\.tiff$", full.names = TRUE)
s1 <- stack(r1)
# 替换指定缺失值为NA
y1 <- s1
y1[y1 %in% c(65533, 65534, 65535)] <- NA

方式2：用reclassify函数精准匹配 #

创建分类矩阵，明确指定每个需要替换的值：

r1 <- list.files(path="D:/Test/", pattern="\\.tiff$", full.names = TRUE)
s1 <- stack(r1)
# 构建分类规则：每行对应[要替换的值, 替换后的值]
reclass_mat <- rbind(c(65533, NA), c(65534, NA), c(65535, NA))
y1 <- reclassify(s1, reclass_mat)

如何导出处理后的文件到指定目录并保留原文件名？ #

需要遍历栅格栈中的每个图层，提取原文件名并拼接输出路径：

# 先创建输出目录（自动创建上级目录）
output_dir <- "D:/Test/Out"
if (!dir.exists(output_dir)) {
  dir.create(output_dir, recursive = TRUE)
}

# 逐个导出处理后的栅格
for (i in 1:nlayers(y1)) {
  # 获取原文件名（去掉路径，仅保留文件名）
  original_filename <- basename(r1[i])
  # 拼接输出路径
  output_path <- file.path(output_dir, original_filename)
  # 导出栅格（允许覆盖已有文件）
  writeRaster(y1[[i]], filename = output_path, format = "GTiff", overwrite = TRUE)
}

完整整合代码 #

# 定义路径
input_dir <- "D:/Test/"
output_dir <- "D:/Test/Out"

# 创建输出目录（若不存在）
if (!dir.exists(output_dir)) {
  dir.create(output_dir, recursive = TRUE)
}

# 读取所有.tiff文件
r1 <- list.files(path=input_dir, pattern="\\.tiff$", full.names = TRUE)
s1 <- stack(r1)

# 替换缺失值为NA
y1 <- s1
y1[y1 %in% c(65533, 65534, 65535)] <- NA

# 批量导出处理后的文件
for (i in 1:nlayers(y1)) {
  output_path <- file.path(output_dir, basename(r1[i]))
  writeRaster(y1[[i]], filename = output_path, format = "GTiff", overwrite = TRUE)
}

内容的提问来源于stack exchange，提问作者Alexia k Boston