读取特定光栅文件并在R中创建平均光栅
我很沮丧,因为我的问题看起来很简单,但我不知道如何处理它 目标:读取特定光栅文件并在R中创建平均光栅,r,design-patterns,raster,qgis,R,Design Patterns,Raster,Qgis,我很沮丧,因为我的问题看起来很简单,但我不知道如何处理它 目标: 我想从文件夹中读取1到4个光栅文件。我需要的人的名字列在一个字符列表中 打开文件后,我想创建一个与文件平均值相对应的新光栅 我可以在QGIS上管理它,但我需要自动化这个过程,因为我有很多人 1) 它应该可以处理list.files(pattern=)但是由于名称在列表中,我不知道怎么做 例如:对于第一个人,我必须阅读名为2018-12-27_-sic.tif和2018-12-27_-sic.tif 我试着用readGDAL,ope
list.files(pattern=)
但是由于名称在列表中,我不知道怎么做
例如:对于第一个人,我必须阅读名为2018-12-27_-sic.tif
和2018-12-27_-sic.tif
我试着用readGDAL
,open阅读。GDAL
不起作用
非常感谢您的宝贵帮助我将使用
光栅
软件包中的堆栈
和计算
功能。函数stack
创建一个光栅堆栈,所有光栅都具有相同的分辨率和范围,并且可以轻松执行诸如获取每个单元的平均值之类的操作。因此:
library(raster)
fs <- list.files(pattern='tif$')
rasterstack <- stack(fs)
rastermean <- calc(rasterstack, fun=mean)
库(光栅)
非常感谢。我终于设法用另一种方法做了。有时候,你只需要扭转问题,找到一个新的解决方案!在本例中,我没有直接读取目标文件,而是读取所有文件,并使用函数[whitch]选择所需的文件