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

本文介绍基于R 语言中的raster包,批量读取多张 栅格图像,对多个栅格图像计算平均值标准差,并将所得新的栅格结果图像保存的方法。

在文章基于R语言的raster包读取遥感影像中,我们介绍了基于R 语言raster包,对单张或多张栅格图像加以平均值标准差 计算的方法;但这一篇文章中的标准差 计算方法仅仅可以对一张栅格图像 的全部像元加以计算,即标准差计算结果是一个具体的数值 ,而不是一景结果影像 ;无法对多张、多时相的栅格图像进行计算。本文就介绍另一种方法,可以对多个时相 的大量栅格影像加以逐像元平均值、标准差的计算,从而使得最终的结果是一景表示各个像元在全部时相的图像中 的平均值或标准差的图像

首先,我们按照文章基于R语言的raster包读取遥感影像中提到的方法,配置、加载raster包,并通过stack()函数读取同一文件夹下的全部栅格图像,具体代码如下所示。其中,代码的含义我们在上述这一篇文章中已经加以介绍,这里就不再赘述。

r 复制代码
library(raster)
tif_file_path <- list.files(r"(E:\02_Project\01_Chlorophyll\LCC_SC_2020\SD)", pattern = ".tif$", full.names = TRUE, ignore.case = TRUE)
tif_file_all <- stack(tif_file_path)

运行上述代码,可以看到已经得到了RasterStack格式的结果数据,如下图所示。

接下来,我们通过calc()函数,对多时相栅格遥感影像数据加以计算;其中,其第一个参数tif_file_all就是需要加以计算的多个栅格图像,而第二个参数fun = sd表示我们需要计算标准差;如果我们需要计算平均值,那么就将第二个参数修改为fun = mean即可,我们这里就以标准差为例介绍后续的操作。当然,前述提到的文章基于R语言的raster包读取遥感影像中的方法也是可以对多个栅格图像计算平均值的。

r 复制代码
tif_sd <- calc(tif_file_all, fun = sd)
plot(tif_sd)

此外,上述代码在calc()函数运行时,若某一空间位置上的像元多张栅格遥感影像 中,存在至少一个 无效值(NoData 值),则这一像元在最终的结果图像 中同样为无效值;若希望忽略无效值的这一影响,可以将上述第一句代码修改为如下格式。其中,na.rm = TRUE就表示若某一景栅格遥感影像中某像元为无效值,则忽略这一景影像中的这一个像元。

r 复制代码
tif_sd <- calc(tif_file_all, fun = sd, na.rm = TRUE)

运行calc()函数后,我们可以通过plot()函数将结果图像绘制出来,如下图所示。

上图即为多个栅格图像的像元数值时间序列依次计算标准差所得的结果。

此外,由于我这里的栅格像元数据实际表达的数值 之间有一个缩放系数0.01,因此通过下述代码将其像元值恢复为实际含义的数值。

r 复制代码
tif_sd_new <- tif_sd / 100
plot(tif_sd_new)

随后,重新绘制结果图;确认无误后,即可依据writeRaster()函数,通过如下代码保存我们刚刚得到的标准差结果栅格图像。

r 复制代码
rf <- writeRaster(tif_sd_new, filename = r"(E:\02_Project\01_Chlorophyll\LCC_SC_2020\SD\LCC_SD.tif)", overwrite = TRUE)

运行代码后,如下图所示。其中,writeRaster()函数的第一个参数表示我们将要保存的栅格数据,第二个参数表示保存栅格文件的路径与名称,第三个参数表示,如果第二个参数指定的路径与名称已经有文件存在了,是否直接对其加以覆盖。

随后,我们即可在指定的路径下找到我们刚刚计算得到的多个栅格图像的标准差结果。

至此,大功告成。

相关推荐
没有梦想的咸鱼185-1037-16632 天前
基于R语言机器学习方法在生态经济学领域中的实践技术应用
开发语言·机器学习·数据分析·r语言
zhangfeng11332 天前
R 语法高亮为什么没有,是需要安装专用的编辑软件,R语言自带的R-gui 功能还是比较简单
开发语言·r语言
WangYan20222 天前
经济学+机器学习+R语言:十大原理、熵权法、随机森林、神经网络、因果推断全解析
r语言·生态经济学·经济学
高-老师2 天前
基于R语言的物种气候生态位动态量化与分布特征模拟
开发语言·r语言·物种气候
Tiger Z2 天前
R 语言科研绘图 --- 其他绘图-汇总2
r语言·论文·科研·绘图·研究生
WangYan20222 天前
【物种分布模型】R语言物种气候生态位动态量化与分布特征模拟——气候生态位动态检验、质心转移可视化、适生区预测等
r语言·物种分布模型·物种气候生态位
q567315233 天前
告别低效:构建健壮R爬虫的工程思维
开发语言·爬虫·r语言
没有梦想的咸鱼185-1037-16633 天前
【高分论文密码】大尺度空间模拟预测与数字制图
信息可视化·数据分析·r语言
kebeiovo4 天前
项目必备流程图,类图,E-R图实例速通
开发语言·r语言·流程图
zhangfeng11334 天前
BiocManager下载失败 R语言 解决办法
开发语言·r语言