如何在R中自动化创建19个data.frame并合并为一个数据框？

首先注意：你原手动代码里存在一个小问题——直接赋值文件路径字符串给raster1是无效的，getValues()需要的是Raster对象，必须用raster()函数加载栅格文件，否则会报错。以下是自动化的完整解决方案：

步骤1：加载依赖包 #

library(raster)

步骤2：自动化处理所有栅格 #

# 生成19个栅格的文件路径列表
file_paths <- paste0("~/Documents/QGIS/Raster/exemple_", 1:19, ".tif")

# 批量处理每个栅格，生成去NA的数据框集合
raster_dfs <- lapply(file_paths, function(path) {
  # 读取栅格文件
  raster_obj <- raster(path)
  # 提取栅格值并转为数据框
  df <- data.frame(raster_val = getValues(raster_obj))
  # 移除含NA的行
  na.omit(df)
})

# 合并所有数据框
raster_19 <- do.call(cbind, raster_dfs)
# 添加Source列
raster_19$Source <- "Native_area"

可选优化：自定义列名 #

如果希望每个数据框的列名和你手动代码的命名（raster1_val、raster2_val...）保持一致，可以修改批量处理的代码：

raster_dfs <- lapply(1:19, function(i) {
  path <- paste0("~/Documents/QGIS/Raster/exemple_", i, ".tif")
  raster_obj <- raster(path)
  # 用序号自定义列名
  df <- data.frame(!!paste0("raster", i, "_val") := getValues(raster_obj))
  na.omit(df)
})

代码说明 #

• paste0自动生成1到19的栅格文件路径，避免手动重复输入
• lapply遍历所有路径，通过匿名函数一次性完成「读取栅格→提取值→转数据框→去NA」的流程，返回包含19个数据框的列表
• do.call(cbind, raster_dfs)替代手动写19个数据框参数，自动完成所有数据框的列合并

内容的提问来源于stack exchange，提问作者Tristan