我尝试使用14个GCM模型进行一些简单的作物分析(134种作物),数据存储在14个文件夹中。我想要做的是从14个模型中提取每个作物的fils,堆叠它们,分析它们然后重复下一个作物。
到目前为止,我写的代码是:
library(raster)
setwd("/Users/Desktop/JERM/Results/Future/RCP4.5/")
for(i in 1:length(CropsUnique))
{
Cropsi <- EcoID[which(Crops %in% CropsUnique[i])]
for(j in 1:length(ldc))
{
CropSpecies <- stack(paste(ldc[j], "Eco", Cropsi, ".grd", sep="/"))
if(length(Cropsi)>1){CropSpecies <- calc(CropSpecies, function(x) max(x, na.rm=TRUE))}
writeRaster(CropSpecies, paste(ldc[j], CropsUnique[i], ".grd", sep=""), overwrite=TRUE)
}
}
但是,我一直收到以下错误&#34;文件错误(x,&#34; rb&#34;):无法打开连接&#34;。我认为问题是R只能到达包含14个GCM的文件夹的目录,但是然后无法同时从每个文件夹中提取单个裁剪文件。
我确定这很简单,但我不能在生活中考虑如何同时从单独的文件夹中提取文件。