在R中更新foreach循环中的栅格值？(光栅时间序列NA估算)

# 安装和加载必要的包
if (!require("foreach")) install.packages("foreach")
if (!require("doParallel")) install.packages("doParallel")
library(foreach)
library(doParallel)

# 创建一个示例栅格数据
set.seed(123)
raster_data <- matrix(runif(100), nrow = 10)
raster_data[sample(1:100, 20)] <- NA

# 转换为Raster对象
library(raster)
r <- raster(raster_data)

# 定义一个函数来估算NA值
estimate_na <- function(x) {
  if (is.na(x)) {
    return(mean(x, na.rm = TRUE))
  } else {
    return(x)
  }
}

# 注册并行后端
cl <- makeCluster(detectCores() - 1)
registerDoParallel(cl)

# 使用foreach循环更新栅格值
r_new <- foreach(i = 1:ncell(r), .combine = 'c') %dopar% {
  cell_value <- getValues(r)[i]
  estimate_na(cell_value)
}

# 关闭并行后端
stopCluster(cl)

# 将新的值赋回栅格对象
values(r) <- r_new

# 查看更新后的栅格数据
print(r)