Warning: file_get_contents(/data/phpspider/zhask/data//catemap/1/typo3/2.json): failed to open stream: No such file or directory in /data/phpspider/zhask/libs/function.php on line 167

Warning: Invalid argument supplied for foreach() in /data/phpspider/zhask/libs/tag.function.php on line 1116

Notice: Undefined index: in /data/phpspider/zhask/libs/function.php on line 180

Warning: array_chunk() expects parameter 1 to be array, null given in /data/phpspider/zhask/libs/function.php on line 181
读取特定光栅文件并在R中创建平均光栅_R_Design Patterns_Raster_Qgis - Fatal编程技术网

读取特定光栅文件并在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

我很沮丧,因为我的问题看起来很简单,但我不知道如何处理它

目标:

  • 我想从文件夹中读取1到4个光栅文件。我需要的人的名字列在一个字符列表中
  • 打开文件后,我想创建一个与文件平均值相对应的新光栅
  • 我可以在QGIS上管理它,但我需要自动化这个过程,因为我有很多人

    1) 它应该可以处理
    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]选择所需的文件