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

在R中提取WorldClim数据,但缺少几个月,r,extract,R,Extract,我有一个关于在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 <- shapefi

我有一个关于在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)
库(sp)
图书馆(rgeos)
图书馆(rgdal)

r在不使用
r的情况下运行代码您使用的
na.rm=TRUE
,因此如果缺少数据,将跳过这些数据,请尝试在不使用此选项的情况下查看您是否在几个月内没有数据
,这将意味着您不幸的是,数据在几个月内不存在(但至少解释一下为什么您什么都没有得到)。谢谢您的回答。我试过了,但得到了同样的结果。在上面的帖子中,我没有提到第一次运行上述代码时,我能够以.cvs格式保存整个数据。然而,第二次尝试时,我遇到了上述问题。如果省略这一行会发生什么:
r谢谢你的回答。我删除了该行,并能够保存所有12个月的数据。是否有任何方法可以循环此命令,以便我可以提取最低和最高温度的数据并保存它,而无需重复运行此命令?例如,
rAlso对于替代质心方法,我是否需要更改任何内容。因为当我运行它时,我得到以下错误
#使用质心nc,我更新了我以前的答案。请检查一下。对于质心选项,shapefile应以
nc的形式读入。代码将一直工作到我执行**sp\n然后尝试回答前面的问题。我现在使用了
sf
库来读取shapefile,并且没有出现错误。
ex <- extract(r, poly, fun='mean', na.rm=TRUE, df=TRUE, weights = TRUE) 
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