计算一个文件夹下所有tiff的均值,输出为Excel表
上一篇博客中我们将一张tiff按照shp研究区进行了掩膜,如果我们想要计算掩膜后每一个小tiff的值,该如何实现?
下面用R语言实现
library(sf)
library(raster)
library(tidyverse)
rm(list = ls())
folder_path <- ""
tif_files <- list.files(path = folder_path, pattern = "\\.tif$", full.names = TRUE)
result_df <- data.frame(File_Name = character(), Mean_Value = numeric(), stringsAsFactors = FALSE)
for (tif_file in tif_files) {
tif <- raster(tif_file)
valid_values <- tif[!is.na(tif)]
num_valid_values <- length(valid_values)
if (num_valid_values < 5) {
next
}
mean_value <- mean(valid_values, na.rm = TRUE)
file_name <- gsub("\\.tif$", "", basename(tif_file))
result_df <- rbind(result_df, data.frame(File_Name = file_name, Mean_Value = mean_value))
print(file_name)
}
result_df$File_Name <- as.numeric(result_df$File_Name)
result_df$Mean_Value <- as.numeric(result_df$Mean_Value)
sorted_df <- result_df[order(result_df$File_Name), ]
full_sequence <- data.frame(File_Name = seq(min(result_df$File_Name), max(result_df$File_Name)))
complete_data <- full_sequence %>% left_join(sorted_df, by = "File_Name")
write_xlsx(complete_data,".xlsx")
![在这里插入图片描述](https://img-blog.csdnimg.cn/direct/b0e18b54b811462f9f624e3c6cfcbcde.png)
图1 这里可以看到缺少2,我们加入了2,并赋值为了na
![在这里插入图片描述](https://img-blog.csdnimg.cn/direct/4723921ea880472dba6bb627a50ce6fb.png)
图2 这里可以看到加入了2,并赋值为了na