R语言无数据栅格缓冲区加权均值为0的问题排查与解决

当使用raster::extract()并设置weights=TRUE和fun=mean, na.rm=TRUE时，如果多边形内所有栅格值都是NA，函数会返回0而非预期的NA。这是因为该函数在处理全NA场景时，内部逻辑会将空值的加权计算结果默认设为0，而非识别为缺失值，这和QGIS的处理逻辑存在差异。

解决方法1：手动计算加权均值（灵活可控） #

不直接在extract中指定fun=mean，而是提取栅格值与对应权重后手动计算，自定义全NA场景的返回值：

library(raster)
library(sf)

# 加载栅格数据
flooding_raster_data <- raster(raster_file)
# 确保缓冲区为sp对象（若原buffer是sf格式）
buffer_sp <- as(buffer, "Spatial")

# 提取值和权重（不指定fun，返回列表，每个元素是包含value和weight的data.frame）
extracted <- extract(flooding_raster_data, buffer_sp, weights=TRUE, cellnumbers=FALSE)

# 自定义加权均值计算函数
calc_weighted_mean <- function(x) {
  valid_idx <- !is.na(x$value)
  if (sum(valid_idx) == 0) {
    return(NA) # 全NA时返回NA
  } else {
    weighted.mean(x$value[valid_idx], x$weight[valid_idx], na.rm=TRUE)
  }
}

# 计算每个缓冲区的加权均值
raster_values_in_buffer_weight <- sapply(extracted, calc_weighted_mean)

# 后续写入shapefile的逻辑不变
buffer_with_mean <- st_as_sf(buffer)
buffer_with_mean$floodshare <- raster_values_in_buffer_weight
st_write(buffer_with_mean, output_file_path, append = FALSE)
rm(flooding_raster_data, buffer_with_mean, raster_values_in_buffer_weight, raster_file)

解决方法2：快速替换0为NA（适用于有效洪水值不为0的场景） #

如果你的洪水栅格有效数值不会是0，可以在提取结果后，将全NA缓冲区对应的0替换为NA：

# 原提取代码
raster_values_in_buffer_weight <- extract(flooding_raster_data, buffer, fun = mean, na.rm = TRUE, weights=TRUE)

# 提取原始栅格值，判断每个缓冲区是否全为NA
raw_values <- extract(flooding_raster_data, buffer)
is_all_na <- sapply(raw_values, function(x) all(is.na(x)))

# 替换全NA对应的结果为NA
raster_values_in_buffer_weight[is_all_na] <- NA

# 后续处理逻辑不变

解决方法3：改用terra包（推荐长期方案） #

raster包已进入维护状态，推荐使用替代包terra，其extract()函数在全NA场景下会正确返回NA：

library(terra)
library(sf)

# 加载栅格数据（terra用rast替代raster）
flooding_raster_data <- rast(raster_file)
# 将sf格式缓冲区转为terra的vect对象
buffer_vect <- vect(buffer)

# 提取加权均值
raster_values_in_buffer_weight <- extract(flooding_raster_data, buffer_vect, fun = mean, na.rm=TRUE, weights=TRUE)

# 后续处理
buffer_with_mean <- st_as_sf(buffer)
# terra返回的是data.frame，取第二列作为均值结果
buffer_with_mean$floodshare <- raster_values_in_buffer_weight[,2]
st_write(buffer_with_mean, output_file_path, append = FALSE)

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