读取特定的栅格文件并在R中创建平均栅格

Question

我很绝望，因为我的问题看起来很简单，但是我找不到解决方法。

目标：

1. 我想从一个文件夹中读取1到4个光栅文件。我需要的一个名称在列表中以字符形式列出。
2. 打开文件后，我想创建一个对应于文件均值的新栅格。

我可以在QGIS上对其进行管理，但是由于我个人很多，我需要使流程自动化！

1）它应该与blake2_256一起使用，但是由于名称在列表中，因此我不知道该怎么做。

例如：对于第一个人，我必须读取2个分别名为list.files(pattern = )和2018-12-27_sic.tif的文件

我尝试用2018-12-27_sic_con.tif阅读，readGDAL无效

非常感谢您的宝贵帮助

Answer 1

我将使用stack包中的calc和raster函数。函数stack创建了所有具有相同分辨率和范围的栅格堆栈，并使操作变得容易，例如对每个像元取均值。所以：

library(raster)

fs <- list.files(pattern='tif$')

rasterstack <- stack(fs)

rastermean <- calc(rasterstack, fun=mean)

请注意，如果栅格的分辨率不同，则必须使用resample函数；如果栅格的分辨率不相同，则必须使用crop。在RStudio中输入?resample和?crop会向您显示有关使用这些功能的说明。