能否在R中直接重命名GeoPackage图层（无需读写删操作）

sf包本身没有提供像你伪代码里那样的st_rename函数来直接重命名GeoPackage图层，但我们可以利用GeoPackage本质是SQLite数据库的特性，直接操作底层数据库来实现，完全不需要执行读入图层、复制重命名、删除原图层这三步。

单图层重命名步骤 #

library(DBI)
library(RSQLite)

# 连接到你的GeoPackage文件
geopackage <- "example.gpkg"
con <- dbConnect(RSQLite::SQLite(), dbname = geopackage)

# 定义新旧图层名
old_layer <- "Layer1_old"
new_layer <- "Layer1_new"

# 1. 重命名图层对应的SQLite表
dbExecute(con, sprintf("ALTER TABLE %s RENAME TO %s", old_layer, new_layer))

# 2. 更新GeoPackage的元数据表（否则st_layers无法识别新名称）
dbExecute(con, sprintf("UPDATE gpkg_contents SET table_name = '%s' WHERE table_name = '%s'", new_layer, old_layer))

# 关闭数据库连接
dbDisconnect(con)

# 验证结果
st_layers(geopackage)

批量循环重命名的实现 #

如果要批量处理多个图层，可以按照以下方式循环执行：

library(sf)
library(DBI)
library(RSQLite)

geopackage <- "example.gpkg"
# 获取所有图层名称
lyrs <- st_layers(geopackage)$name

# 定义新旧名称的映射关系，这里示例是给原名称加_new后缀，你可以自定义规则
name_mapping <- setNames(paste0(lyrs, "_new"), lyrs)

# 连接数据库
con <- dbConnect(RSQLite::SQLite(), dbname = geopackage)

# 循环处理每个图层
for (old_name in names(name_mapping)) {
  new_name <- name_mapping[old_name]
  
  # 重命名表
  dbExecute(con, sprintf("ALTER TABLE %s RENAME TO %s", old_name, new_name))
  
  # 更新元数据
  dbExecute(con, sprintf("UPDATE gpkg_contents SET table_name = '%s' WHERE table_name = '%s'", new_name, old_name))
}

# 关闭连接
dbDisconnect(con)

# 查看重命名后的图层列表
st_layers(geopackage)

注意事项 #

• 操作前一定要备份GeoPackage文件，避免误操作导致数据丢失
• 如果图层有对应的空间索引（名称以rtree_开头的表），重命名主表后索引会失效，你可以手动重建索引：

# 循环重建所有新图层的空间索引
for (new_name in name_mapping) {
  st_create_index(geopackage, layer = new_name)
}

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