为何无法重命名rnaturalearth的sf类数据列？

错误原因 #

这个报错是因为sf类型对象自带agr属性（全称attribute-geometry relationship，用于标记非几何列与空间几何的关联类型）。当用dplyr::rename()直接重命名列时，没有同步更新agr属性的名称映射，导致sf内部无法找到对应列的agr记录，从而触发报错。

解决方案 #

提供三种可行的解决方法：

方法一：使用sf包原生重命名方法 #

sf包的rename()方法会自动同步更新agr属性，避免报错：

library(rnaturalearth)
library(dplyr)
library(sf)

baseBorderMap <- ne_countries(returnclass = 'sf')
# 用sf的rename方法重命名
baseBorderMap <- sf::rename(baseBorderMap, countryCode = iso_a2)

方法二：先移除agr属性再重命名 #

如果不需要保留agr属性，可以先移除再执行重命名：

baseBorderMap <- ne_countries(returnclass = 'sf')
# 移除agr属性
attr(baseBorderMap, "agr") <- NULL
# 执行重命名
baseBorderMap <- rename(baseBorderMap, countryCode = iso_a2)
# 可选：若需重新设置agr属性，可使用st_set_agr
baseBorderMap <- st_set_agr(baseBorderMap, "countryCode", "constant")

方法三：使用dplyr的rename_with函数 #

rename_with同样能适配sf对象的重命名操作：

baseBorderMap <- ne_countries(returnclass = 'sf')
baseBorderMap <- rename_with(baseBorderMap, ~"countryCode", .cols = iso_a2)

额外提示 #

注意rnaturalearth返回的iso_a2列存在部分缺失值（如南极洲），后续关联数据集时建议提前处理这些缺失值，避免关联失败。

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