0

我正在尝试使用 rasters 包将大约 50 个 .tif 文件导入堆栈。由于文件太多,我不想一一列举。我尝试了几种方法,但都没有成功,包括以下代码:

test<-stack(choose.files(),values=T) .local(x, ...) 中的错误:参数应该是 Raster* 对象或文件名

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

4

1 回答 1

1

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

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

但是,如果您想以交互方式进行操作:

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

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

于 2013-05-11T19:12:18.117 回答