关于从多个 nc 文件的多个坐标中提取数据并合并到单个栅格

Question

我一年有 365 个 .nc 文件，每个文件都包含一天的土壤水分信息。我想从三个不同的坐标中提取土壤水分信息，并将它们写入给定年份的三个 csv 文件。目前，我能够将所有 365 .nc 文件转换为 365 csv 文件，如附加代码中给出的。下一步怎么走？

rm(list = ls())
library(raster)
library(ncdf4)
ptf <- "D://SMOS_ECV_SM//SMOS_ECV_SM//ECV_SM_Data_1978_2010//1978"
setwd(ptf) # change your working directory
lf <- list.files(pattern="[.]nc$") # list of files ending in .nc
for(i in lf){
  nc.brick <- brick(i)
  nc.df <- as.data.frame(nc.brick[[1]], xy=T)
  write.csv(nc.df, sub("[.]nc$",".csv",i)) # write to the same file name substituting .nc to .csv
}

Answer 1

您应该能够执行以下操作

library(raster)
ptf <- "D://SMOS_ECV_SM//SMOS_ECV_SM//ECV_SM_Data_1978_2010//1978"
ff <- list.files(path=ptf, pattern="[.]nc$", full.names=TRUE)

s <- stack(ff)
points <- rbind(c(0,1), c(100,120), c(80,5))
v <- extract(s, points)

for (i in 1:ncol(v)) {
  write.csv(v[,i,drop=FALSE], paste0("file", i, ".csv"))
}