R语言遍历文件夹求取其中全部栅格文件的平均值

  本文介绍基于R语言中的raster包，遍历读取多个文件夹下的多张栅格遥感影像，分别批量对每一个文件夹中的多个栅格图像盘算平均值，并将所得各个结果栅格分别加以生存的方法。
  其中，本文是用R语言来进行操作的；假如希望基于Python语言实现类似的平均值求取操作，大家可以参考Python ArcPy批量盘算多时相遥感影像的各像元平均值与Python忽略NoData盘算多张遥感影像的像元平均值：whitebox库这两篇文章。
  起首，来看一下本文所需实现的需求。如下图所示，现有多个文件夹，其中每一个文件夹内部都含有大量的栅格遥感影像。

  其中，上图中的每一个文件夹的定名都是以遥感影像的分幅条带号为依据的。例如，打开第一个名为47RMN的文件夹，其中均为条带号为47RMN（即同一空间范围）、不同成像时间的遥感影像，如下图所示；其中，紫色框内的遥感影像文件名即可看出，这些图像是同一条带号、不同时间的遥感影像数据。

  我们要做的，就是分别对每一个文件夹中的全部遥感影像盘算平均值，从而得到不同条带号遥感影像的平均值；最终我们将得到多张结果图像，每一景结果图像就是这一条带号、不同成像时间对应的遥感影像的平均值。同时为了方便区分，我们需要将每一景结果图像文件的文件名设置为与条带号有关的内容。
  明确了需求，我们即可开始代码的撰写。本文所用到的代码如下所示。
[code]library(raster)result_path

继续阅读请点击广告