在 R 中提取 WorldClim 数据但缺少几个月

Question

我有一个关于在 R 中提取 WorldClim 数据的问题并得到了很好的帮助（ 使用 R 为单个国家提取 WORLDCLIM 数据）

该解决方案工作正常，但我现在面临一个单独的问题。

我执行以下代码以 .cvs 格式提取并保存 WorldClim 每月最低温度数据

library(sp)
library(rgeos)
library(rgdal)
r<- getData('CMIP5', var='tmin', res=10, rcp=45, model='HE', year=70)
r <- r[[c(1,12)]]
poly <- shapefile("C:/Users/pc5/Desktop/r1//2011_Dist.shp")
plot(poly)
pr <- extract(r, poly, fun='mean', na.rm=TRUE, df=TRUE, weights = TRUE) 
write.csv(cbind(poly$DISTRICT,ex),"Worldclim.csv", row.names = F)

但是，当我提取并保存数据时，我只得到三个变量：DISTRICT 名称、第 1 个月的 TMIN 和第 12 个月的 TMIN。我在运行代码时看到了

ex <- extract(r, poly, fun='mean', na.rm=TRUE, df=TRUE, weights = TRUE) 

它只提取了 3 个变量，而不是 13 个（12 个月 + I 区名），尽管我可以看到所有 12 个月的数据都以 Geotiff 格式下载到默认文件夹中。 我究竟做错了什么？

Answer 1

运行没有r <- r[[c(1,12)]] ，你会得到想要的结果。 r <- r[[c(1,12)]]部分代码仅从 12 个月数据的栅格堆栈中删除 1 个月和 12 个月。 您可以通过以下代码确定

r<- getData('CMIP5', var='tmin', res=10, rcp=45, model='HE', year=70)
r #will show you dimensions  : 900, 2160, 1944000, 12  (nrow, ncol, ncell, nlayers)
r1 <- r[[c(1,12)]]
r1 #will show you dimensions  : 900, 2160, 1944000, 2  (nrow, ncol, ncell, nlayers)
plot(r) #Plot of all the 12 months
plot(r1) #Plot of only 1 and 12 month

我也编辑了我以前的答案。