R语言循环为sf空间数据框新增原始地图标识列方法

原有循环代码的问题在于全程仅通过paste()拼接赋值语句的文本字符串，没有将文本转化为实际的赋值操作：

• 没有通过传入的字符串定位到全局环境中存储的真实sf对象
• 没有执行新增列的实际逻辑
• 没有将修改后的结果回写到原对象中，最终生成的内容只是展示用的文本，不会对实际数据产生任何影响。

核心逻辑是通过get()从全局环境按字符串名取出真实sf对象，完成列修改后通过assign()将结果回写，替换原对象，完全匹配手动赋值的效果：

# 待处理的sf对象名向量
list_df <- c("map_g1_r1_airport", "map_g1_r1_hotel", "map_g1_r2_bank", "map_g1_r2_market")

for (df_name in list_df) {
  # 从全局环境取出当前对应的sf对象
  current_sf <- get(df_name, envir = .GlobalEnv)
  # 从对象名中提取所属原始地图名称（截断最后一个下划线及之后的要素类型标识）
  belong_map_id <- gsub("(.*)_.*", "\\1", df_name)
  # 为sf对象新增mapid列
  current_sf$mapid <- belong_map_id
  # 将修改完成的对象回写到全局环境，覆盖原有内容
  assign(x = df_name, value = current_sf, envir = .GlobalEnv)
}

执行代码后，可任意抽取一个对象检查结果，例如：

# 查看机场sf对象的mapid列值
map_g1_r1_airport$mapid

输出结果为与对象行数等长的、值为"map_g1_r1"的字符向量，和给出的预期sf结构完全一致。

如果后续需要批量处理大量同命名规则的sf对象，更推荐将零散对象整理为命名列表统一管理，避免全局环境存储过多零散变量：

# 从全局环境批量匹配所有符合命名规则的sf对象，整理为列表
sf_object_list <- mget(
  x = ls(pattern = "^map_g\\d_r\\d_"),
  envir = .GlobalEnv
)
# 批量为列表内所有sf对象新增mapid列
sf_object_list <- lapply(names(sf_object_list), function(each_name) {
  sf_object_list[[each_name]]$mapid <- gsub("(.*)_.*", "\\1", each_name)
  return(sf_object_list[[each_name]])
})

这种方式不需要循环单独回写全局对象，后续做空间要素合并、批量导出操作时效率更高。

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