Извлечение климатических и почвенных данных для долготно-широтных точек

Question

Я хотел бы получить информацию о климате и почве для местоположений GPS в ЕС за 2006–2016 годы.

Я заглянул в 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

Вы можете использовать растр :: 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