time-series - 関心のあるポイントに最も近い複数のポイントから E-OBS から気候時系列を抽出する

Question

次のコードを更新して、E-OBS データベースの毎日の気候時系列を複数の (ユーザー定義の) 最も近いポイントに抽出したいと思います。

ncdf4 パッケージを使用して E-OBS データを読み取ります。

library(ncdf4)
obsdata <- nc_open("tg_ens_mean_0.1deg_reg_v19.0e.nc")

位置と WGS84 座標を含むデータフレームを作成します。

df <- data.frame(locality=c(1,2,3), lat=c(48,48,48), lon=c(17,18,19))

次に、これらのコードを実行して、最も近いポイントからデータを抽出し、場所の ID ごとに個別の csv ファイルに保存します。

for (i in 1:nrow(df)) {
  lat <- df[i,2] # longitude of location
  lon <- df[i,3] # latitude  of location
  
  
  # get dates
  obsdatadates <- as.Date(obsdata$dim$time$vals, origin = '1950-01-01')
  
  
  # get values at location lonlat
  obsoutput <- ncvar_get(obsdata, varid = "tg",
                         start= c(which.min(abs(obsdata$dim$longitude$vals - lon)), # look for closest long
                                  which.min(abs(obsdata$dim$latitude$vals - lat)),  # look for closest lat
                                  1),
                         count = c(1,1,-1)) #count '-1' means 'all values along that dimension'that dimension'
  # create dataframe
  datafinal <- data.frame(date= obsdatadates, value = obsoutput, latitude=lat, longitude=lon, variable="Tday")
  fn1<-paste0(df[i,1],"_","Tday",".csv")
  write.table(datafinal,fn1,sep=",", row.names=FALSE)
}

どんな助けでも大歓迎です！ありがとう！