r - 使用 dismo 包创建生物气候变量的问题

我想获得物种分布建模的 future 情景的所有生物气候变量。所以我使用 worldclim 数据库中的三个变量在“dismo”包中运行“biovars”函数，我得到了一个 12 层的 RasterBrick:

prec<-stack(paste(getwd(),"/prec_2080/wc_2_5m_HADCM3_B2a_2080_prec_",1:12,".bil",sep=""))
tmin<-stack(paste(getwd(),"/tmin_2080/wc_2_5m_HADCM3_B2a_2080_tmin_",1:12,".bil",sep=""))
tmax<-stack(paste(getwd(),"/tmax_2080/wc_2_5m_HADCM3_B2a_2080_tmax_",1:12,".bil",sep=""))
x<-biovars(prec=prec,tmin=tmin,tmax=tmax)
 x
class       : RasterBrick
dimensions  : 3600, 8640, 12  (nrow, ncol, nlayers)
resolution  : 0.04166667, 0.04166667  (x, y)
extent      : -180, 180, -60, 90  (xmin, xmax, ymin, ymax)
projection  : NA
values      : C:/DOCUME~1/Marco/LOCALS~1/TMP/R_raster_tmp/raster_tmp_8984740455.grd
min values  :     42 -65458  -1017      0     71      0 -65439     22     23     56 ...
max values  :  65456    213      1  34159  65534  65513  65534  65507  65503  65518 ...

但是，我认为应该有 19 个生物气候变量。正如你提到的，除了那里之外，生物变种还有更多的争论，但我不知道它们是什么。你能帮我吗？
另一个问题是我在编写这些变量时出错:

writeRaster(x,paste(getwd(),"/wc_2_5m_HADCM3_B2a_2080_1.grd",sep=""))

而且，当我尝试逐个编写它们时，出现以下错误:

for (i in 10:12) {
writeRaster(x[[i]],paste(getwd(),"/wc_2_5m_HADCM3_B2a_2080_",i,".grd",sep=""),overwrite=TRUE)
}

三个输入变量具有相同的维度，例如:

prec
class       : RasterStack
dimensions  : 3600, 8640, 12  (nrow, ncol, nlayers)
resolution  : 0.04166667, 0.04166667  (x, y)
extent      : -180, 180, -60, 90  (xmin, xmax, ymin, ymax)
projection  : NA
min values  : 0 0 0 0 0 0 0 0 0 0 ...
max values  : 65535 65535 65535 65535 65535 65535 65535 65535 65535 65535 ...

有人能解释一下为什么吗？先谢谢了~

最佳答案

这确实是一个错误。它已在 0.5-19 版本中得到修复，该版本应该会在 24 小时内从 R-Forge 和很快从 CRAN 获得。相对湿度

关于r - 使用 dismo 包创建生物气候变量的问题，我们在Stack Overflow上找到一个类似的问题：https://stackoverflow.com/questions/5116448/