R语言求取大量遥感影像的平均值、标准差：raster库

  本文介绍基于R语言中的raster包，批量读取多张栅格图像，对多个栅格图像计算平均值、标准差，并将所得新的栅格结果图像保存的方法。
  在文章基于R语言的raster包读取遥感影像中，我们介绍了基于R语言raster包，对单张或多张栅格图像加以平均值、标准差计算的方法；但这一篇文章中的标准差计算方法仅仅可以对一张栅格图像的全部像元加以计算，即标准差计算结果是一个具体的数值，而不是一景结果影像；无法对多张、多时相的栅格图像进行计算。本文就介绍另一种方法，可以对多个时相的大量栅格影像加以逐像元平均值、标准差的计算，从而使得最终的结果是一景表示各个像元在全部时相的图像中的平均值或标准差的图像。
  首先，我们按照文章基于R语言的raster包读取遥感影像中提到的方法，配置、加载raster包，并通过stack()函数读取同一文件夹下的全部栅格图像，具体代码如下所示。其中，代码的含义我们在上述这一篇文章中已经加以介绍，这里就不再赘述。
[code]library(raster)tif_file_path

分享、互助 让互联网精神温暖你我