从R循环中的多个文件夹中提取文件

时间:2016-02-23 18:41:10

标签: r r-raster

我尝试使用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的文件夹的目录,但是然后无法同时从每个文件夹中提取单个裁剪文件。

我确定这很简单,但我不能在生活中考虑如何同时从单独的文件夹中提取文件。

0 个答案:

没有答案