基于R语言按日期从NetCDF提取对应位点的SST值需求

首先精简依赖包，terra 包可以高效处理NetCDF和空间数据，配合dplyr处理数据框会更顺畅：

library(terra)
library(dplyr)

第一步：整理NetCDF文件并匹配年月 #

先从文件名中提取对应的年月信息，方便后续和鲸类记录匹配：

# 获取所有nc文件路径
filenames <- list.files('SST/Ocean_ColorSST_2002_2022', pattern = '\\.nc$', full.names = TRUE)

# 从文件名提取年月（假设文件名包含类似"200207"的格式，根据实际文件名调整正则）
file_info <- tibble(file_path = filenames) %>%
  mutate(year_month = stringr::str_extract(file_path, "\\d{6}"), # 提取6位数字的年月
         year = substr(year_month, 1, 4) %>% as.integer(),
         month = substr(year_month, 5, 6) %>% as.integer())

注意：如果你的文件名格式不是包含6位数字年月，需要修改str_extract的正则表达式，确保能正确提取每个文件对应的年份和月份。

第二步：预处理鲸类出现记录 #

筛选出SST数据覆盖时间段（2002-07至2022-07）的记录，并创建年月匹配标识：

# 假设你的鲸类数据框名为whale_data
whale_data <- whale_data %>%
  # 筛选在SST数据时间范围内的记录
  filter(year >= 2002 & (year > 2002 | month >=7),
         year <=2022 & (year <2022 | month <=7)) %>%
  # 统一年月格式为YYYYMM，方便匹配
  mutate(year_month = sprintf("%04d%02d", year, month))

第三步：批量提取SST值 #

为了避免重复打开/关闭NetCDF文件，按年月分组处理更高效：

# 定义提取函数：输入年月和对应经纬度，返回SST值
extract_sst <- function(year_month, lon, lat) {
  # 找到对应年月的nc文件路径
  file_path <- file_info$file_path[file_info$year_month == year_month][1]
  if (is.na(file_path)) return(NA)
  
  # 读取nc文件为SpatRaster
  sst_raster <- rast(file_path, varname = "sst")
  
  # 提取对应经纬度的SST值
  extract(sst_raster, data.frame(x = lon, y = lat))[,1]
}

# 按年月分组提取
whale_data_with_sst <- whale_data %>%
  group_by(year_month) %>%
  mutate(SST = extract_sst(year_month, longitude, latitude)) %>%
  ungroup()

关键注意事项 #

• 确保鲸类数据的经度、纬度字段名和代码中的longitude、latitude一致，不一致的话要替换成实际字段名。
• 如果NetCDF文件的坐标系统和鲸类记录的坐标系统不一致，需要先做坐标转换（用terra::project），不过一般海洋数据都是WGS84（EPSG:4326），和鲸类记录的经纬度一致的话可以跳过。
• 如果部分经纬度位点不在SST数据的范围内，提取结果会是NA，后续可以根据需求处理这些缺失值。

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