如何在R的RasterStack或SpatRaster中替换NA为指定数值？

需求：将terra包的SpatRaster或raster包的RasterStack中所有图层的NA值替换为-9999。

可复现测试数据 #

# 构建terra的SpatRaster
m <- terra::rast(ncol=100, nrow=100)
n <- terra::rast(ncol=100, nrow=100)
m[] <- runif(ncell(m))
n[] <- runif(ncell(n))
m[m < 0.2] <- NA
n[n < 0.2] <- NA
sr <- c(m,n) # 多层SpatRaster

# 转换为raster的RasterStack
rs <- raster::stack(sr)

原代码问题分析 #

你用for循环直接赋值时得到NULL，是因为for循环在R中默认返回NULL，没有将每次处理后的图层收集并保存到结果对象中。

解决方案 #

1. terra包（处理SpatRaster）

方法1：循环收集结果

初始化空的SpatRaster，逐图层处理后添加：

sr.na_fill <- terra::rast()
for (i in 1:terra::nlyr(sr)) {
  filled_layer <- terra::classify(sr[[i]], cbind(NA, -9999))
  sr.na_fill <- c(sr.na_fill, filled_layer)
}

方法2：批量处理（更高效，无需循环）

terra的classify函数支持直接对多层SpatRaster操作，无需逐个图层循环：

sr.na_fill <- terra::classify(sr, cbind(NA, -9999))

2. raster包（处理RasterStack）

方法1：循环收集结果

用列表暂存处理后的图层，最后转为RasterStack：

rs_list <- list()
for (i in 1:raster::nlayers(rs)) {
  rs_list[[i]] <- raster::reclassify(rs[[i]], cbind(NA, -9999))
}
rs.na_fill <- raster::stack(rs_list)

方法2：用lapply批量处理（更简洁）

rs.na_fill <- raster::stack(lapply(rs, function(x) raster::reclassify(x, cbind(NA, -9999))))

内容的提问来源于stack exchange，提问作者jlklein