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

Question

我想获得2006至2016年欧盟全球定位系统地点的气候和土壤信息。

然而，我查了一下Worldclim，这些记录似乎在2010年就停止了。关于土壤信息，我想使用欧洲土壤数据库v2.0，它提供光栅文件。不幸的是，我不知道如何从光栅文件中提取值。对于气候信息，我尝试使用RNCEP，但是无法找到有关此错误的信息。

字符串拆分中的错误(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

你可以使用光栅：：提取

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