使用循环计算.nc环境数据集月均值与SD时遇递归索引错误

问题背景 #

尝试从2018-2019年每日海表温度（SST）.nc数据集中，按月份计算变量的均值和标准差，但运行循环时触发错误。

错误信息 #

Error in h(simpleError(msg, call)) : 
  error in evaluating the argument 'x' in selecting a method for function 'mean': recursive indexing failed at level 2

原始脚本 #

# Load required libraries (install the required libraries using the Packages tab, if necessary)
library(raster)
library(ncdf4)

#Open the .nc file with the environmental data 

ENV = nc_open("SST.nc")
ENV

#create an index of the month for every (daily) capture from 2018 to 2019 (in this dataset)

m_index = c()
for (y in 2018:2019) {
  # if bisestile year (do not apply for this data but in case a larger year set is used)
  if (y%%4==0) { m_index = c(m_index, rep(1:12 , times = c(31,29,31,30,31,30,31,31,30,31,30,31))) }
  # if non-bisestile year
  else { m_index = c(m_index, rep(1:12 , times = c(31,28,31,30,31,30,31,31,30,31,30,31))) }
}

length(m_index) #  expected length (730)
table(m_index) # expected number of records assigned to each of the twelve months

# computing of monthly mean and standard deviation.
# We first create two empty raster stack...
SST_MM = stack() # this stack will contain the twelve average SST (one per month)
SST_MSD = stack() # this stack will contain the twelve SST st. dev. (one per month)

# We run the following loop (this can take a while)
for (m in 1:12) { # for every month
  
  print(m) # print current month to track the progress of the loop...
  
  sstMean = mean(ENV[[which(m_index==m)]], na.rm=T) # calculate the mean SST for all the records of the current month
  sstSd = calc(ENV[[which(m_index==m)]], sd, na.rm=T) # calculate the st. dev. of SST for all the records of the current month
  
  # add the monthly records to the stacks
  
  SST_MM = stack(SST_MM, sstMean)
  SST_MSD = stack(SST_MSD, sstSd)
  
}

问题根源 #

nc_open()返回的是ncdf4对象，而非raster包可直接处理的RasterStack/RasterBrick对象。你试图用raster包的索引方式（[[）和函数（mean()、calc()）操作ncdf4对象，导致索引逻辑冲突，触发递归索引错误。

修正后的解决方案 #

用raster包直接读取nc文件生成raster对象，同时简化月份索引的生成逻辑，确保所有操作都基于raster兼容的对象：

# 加载依赖包
library(raster)

# 用raster包的brick()函数读取nc文件，生成RasterBrick（适合多图层的nc数据）
ENV <- brick("SST.nc")

# 生成2018-2019年每日的日期序列，从中提取月份索引
dates <- seq.Date(as.Date("2018-01-01"), as.Date("2019-12-31"), by = "day")
m_index <- as.integer(format(dates, "%m"))

# 初始化存储结果的空栈
SST_MM <- stack()
SST_MSD <- stack()

# 循环计算每月的均值和标准差
for (m in 1:12) {
  print(m)
  # 筛选当前月份对应的所有SST图层
  month_layers <- ENV[[which(m_index == m)]]
  # 计算月度均值
  sstMean <- mean(month_layers, na.rm = TRUE)
  # 计算月度标准差
  sstSd <- calc(month_layers, fun = sd, na.rm = TRUE)
  # 将结果添加到栈中
  SST_MM <- stack(SST_MM, sstMean)
  SST_MSD <- stack(SST_MSD, sstSd)
}

# 可选：给结果图层命名，方便识别
names(SST_MM) <- paste0("SST_mean_", sprintf("%02d", 1:12))
names(SST_MSD) <- paste0("SST_sd_", sprintf("%02d", 1:12))

关键修正说明 #

1. 读取方式修正：用brick()替代nc_open()，直接生成raster包可处理的对象，支持图层索引和统计计算。
2. 月份索引简化：通过日期序列自动提取月份，避免手动计算天数可能出现的错误（如闰年判断疏漏）。
3. 对象类型匹配：确保所有统计计算的输入都是raster图层集合，让mean()和calc()函数正常工作。

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