提取经纬度点的气候和土壤数据

Question

我想获取2006年至2016年欧盟GPS位置的气候和土壤信息。

我研究了Worldclim，但是记录似乎要在2010年停止。对于土壤信息，我想使用提供栅格文件的欧洲土壤数据库v2.0。不幸的是，我不知道如何从栅格文件中提取值。对于气候信息，我尝试使用RNCEP，但是找不到有关此错误的信息

  

strsplit（grep（“ missing_value”，x = readLines（scale.offset.missingvals.temp），）中的错误：     下标超出范围

library(RNCEP)
dat <- NCEP.gather(variable='air', level=850, months.minmax=c(1,12), 
years.minmax=c(2006,2016), lat.southnorth=c(51,53), 
lon.westeast=c(4,5),reanalysis2 = FALSE, return.units = TRUE)

答案将简要描述如何从栅格文件中提取点的值以及如何获取2006年至2016年的温度和降雨量信息。

Answer 1

您可以使用raster :: extract

library(raster)
r <- raster(ncol=36, nrow=18, vals=1:(18*36))
xy <- cbind(-50, seq(-80, 80, by=20))

extract(r, xy)
#[1] 626 554 482 410 338 266 194 122  50