导入栅格堆栈并返回新图层

时间:2013-05-07 17:53:05

标签: r raster

我正在尝试使用rasters包将大约50个.tif文件导入堆栈。由于文件数量很多,我不想全部列出来。我尝试了一些方法,但都没有成功,包括这段代码:

  

试验< -stack(choose.files(),值= T)   .local(x,...)出错:     参数应为Raster *对象或文件名

此外,我想基于50 I导入对新栅格进行分类。例如,我有一组具有相同范围的多边形(可以转换为栅格),并且想要计算每个多边形中出现值“1”的次数。我认为“覆盖”功能最适合这个,但我很难想到如何编码相关的功能。任何帮助将不胜感激。

1 个答案:

答案 0 :(得分:1)

许多文件的简单方法应该是使用(经过一些修改)

f <- list.files(path='???', pattern='.tif$', full.names=TRUE)
s <- stack(f)

但如果你想以交互方式进行:

f <- choose.files()
s <- stack(f)

发生错误是因为您使用“values = TRUE”(这不是堆栈的参数并被解释为文件名)