關於從多個nc文件的多個cordiante中提取數據並合並到單個柵格

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"))
}