Извлечение данных WORLDCLIM с использованием R для отдельной страны - PullRequest
Новая
       71

Извлечение данных WORLDCLIM с использованием R для отдельной страны

0 голосов
/ 20 февраля 2020

Я хочу извлечь данные о мировом климате для минимальной и максимальной температуры только для одной страны Индии, используя R, и сохранить его как набор данных (для использования с моим собственным набором данных, который содержит урожайность на районном уровне). Я просмотрел несколько постов и вижу, что это легко сделать в R, однако посты, которым я пытался следовать, немного отличаются с точки зрения команд или последовательностей, и я запутался. (https://gis.stackexchange.com/questions/259478/worldclim-data-na-for-my-coordinates, https://gis.stackexchange.com/questions/227585/using-r-to-extract-data-from-worldclim

Я попытался использовать следующее. 

library(raster)
library(sp)
r<- getData('CMIP5', var='tmin', res=10, rcp=45, model='HE', year=70)
r <- r[[c(1,12)]]
values <- extract(r,points)
df <- cbind.data.frame(coordinates(points),values)
head(df)

Однако я могу запустить только первые две строки и значения строк 

<- extract(r,points) gives the error Error in (function (classes, fdef, mtable)  : 
  unable to find an inherited method for function ‘extract’ for signature ‘"RasterStack", "function"’

Есть предложения?

1 Ответ

0 голосов
/ 20 февраля 2020

Вот решение для этого 

library(raster)
library(sp)
library(rgeos)
library(rgdal)
library(sf)

r<- getData('CMIP5', var='tmin', res=10, rcp=45, model='HE', year=70)

#Using Zonal statistics
poly <- shapefile("Provide_your_drive_name" e.g. "F:\\Kriging in R\\India Shape files\\2011_Dist.shp")
plot(poly)

#This will take considerable time
ex <- extract(r, poly, fun='mean', na.rm=TRUE, df=TRUE, weights = TRUE) 

write.csv(cbind(poly$DISTRICT,ex),"Worldclim.csv", row.names = F)

# using centroids
nc <- st_read(dsn="Provide_your_drive_name" e.g. "F:\\Kriging in R\\India Shape files", layer="2011_Dist")
# just view the attributes & first 6 attribute values of the data
head(nc)
sp_cent <- gCentroid(as(nc, "Spatial"), byid = TRUE)
values <- extract(r,sp_cent)

write.csv(cbind(as.data.frame(nc$DISTRICT),as.data.frame(values)),"Worldclim_centroids.csv", row.names = F)
...