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

Question

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

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

Error in strsplit(grep("missing_value", x=readLines(scale.offset.missingvals.temp), : subscript out of bounds

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

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

Extracting climate and soil data for longitude-latitude points

r

latitude-longitude

r-raster