NetCDF与CSV数据匹配报错求助：循环参数长度为零错误

错误根源 #

这个错误的核心是当grep未找到对应日期的SST文件时，r_sub变成空向量，调用raster()读取空对象触发内部判断错误。同时代码还存在几个隐性问题：

• 日期格式不匹配：CSV中的月/日是单数字（如6、1），但导出的ASCII文件名大概率是两位格式（如06、01），导致paste0(Y,M,D)拼接的字符串无法匹配文件名。
• 无效日期遍历：循环遍历了所有1-12月、1-31日，但很多日期（如2月30日）不存在对应的SST数据，徒增无效循环。
• 缺失ID列：CSV示例中没有ID列，但代码中用absence_dataset$ID进行匹配，会导致match函数出错。

修复步骤 #

• 统一日期格式：将CSV中的年/月/日转换为与SST文件名一致的格式（如YYYYMMDD的8位字符串），确保匹配逻辑正确。
• 添加空值判断：在读取栅格前检查目标文件是否存在，避免传入空对象给raster()。
• 移除无效循环：只遍历CSV中实际存在的日期，而非所有可能的年/月/日组合，提升效率。
• 补全ID列：给absence_dataset添加唯一ID，用于后续点位匹配。

修正后的完整代码 #

#### 修正版：SST数据与点位匹配脚本 ####

# 加载包
library(raster)
library(ncdf4)
library(rgdal)
library(sp)
library(readr)
library(dplyr)
library(stringr)

# 读取NetCDF并导出为ASCII（已完成可注释）
# temp <- raster::stack("cmems_mod_ibi_phy_my_0.083deg-3D_P1M-m_1694536161344.nc copy")
# writeRaster(temp, 
#             "/Users/lukeainsworth/Library/CloudStorage/OneDrive-UniversityofPlymouth/OneDrive - University of Plymouth/Ocean Science and Marine Conservation/Dissertation/Dissertation_R/Temperature/",
#             bylayer=T, format="ascii", overwrite=T)

# 读取点位数据
absence_dataset <- read.csv("absence_dataset.csv")

# 添加唯一ID列（用于后续匹配）
absence_dataset$ID <- seq(nrow(absence_dataset))

# 将点位日期转换为YYYYMMDD格式字符串，统一匹配规则
absence_dataset$date_str <- with(absence_dataset, 
                                 sprintf("%04d%02d%02d", Year, Month, Day))

# 读取所有SST的ASCII文件（使用full.names=T避免路径问题）
sst_files <- list.files(
  path = "/Users/lukeainsworth/Library/CloudStorage/OneDrive-UniversityofPlymouth/OneDrive - University of Plymouth/Ocean Science and Marine Conservation/Dissertation/Dissertation_R/Temperature",
  full.names = TRUE, 
  pattern = "\\.asc$"
)

# 提取文件名中的8位日期（根据实际文件名格式调整正则）
sst_dates <- str_extract(basename(sst_files), "\\d{8}")
names(sst_files) <- sst_dates

# 只遍历点位数据中存在的唯一日期，避免无效循环
unique_dates <- unique(absence_dataset$date_str)

# 初始化SST列
absence_dataset$sst <- NA

for(date_str in unique_dates) {
  # 跳过无对应SST文件的日期
  if(!date_str %in% names(sst_files)) next
  
  # 筛选当前日期的点位
  pts_sub <- subset(absence_dataset, date_str == date_str)
  
  # 读取对应日期的SST栅格
  r_sub <- raster(sst_files[date_str])
  
  # 仅当坐标系不匹配时进行投影转换
  if(!compareCRS(crs(r_sub), CRS("+proj=longlat +datum=WGS84 +no_defs"))) {
    r_sub <- projectRaster(r_sub, crs = CRS("+proj=longlat +datum=WGS84 +no_defs"))
  }
  
  # 提取点位对应的SST值
  pts_sub$sst <- extract(r_sub, pts_sub[, c("Longitude", "Latitude")])
  
  # 更新原数据集的SST值
  absence_dataset$sst[absence_dataset$ID %in% pts_sub$ID] <- pts_sub$sst
}

# 可选：删除临时的date_str列
# absence_dataset <- select(absence_dataset, -date_str)

额外说明 #

• 文件名匹配：如果你的ASCII文件名不是8位日期格式，需要调整sst_dates <- str_extract(...)中的正则表达式，确保能正确提取日期字符串。
• 投影优化：添加了compareCRS判断，避免重复执行投影转换，提升运行效率。
• 效率提升：通过只遍历存在的日期，减少了大量无效循环，同时用命名向量直接匹配文件，替代多次grep操作。

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