如何保存包含多个几何列的Simple Feature对象？

其实你遇到的问题核心是Shapefile格式的先天限制——它只允许每个文件包含单个几何列，所以直接用st_write()保存多几何列的sf对象必然会报错。下面是几种按推荐度排序的最佳解决方案：

1. 优先使用GeoPackage格式（首推） #

GeoPackage（GPKG）是OGC制定的现代GIS标准格式，原生支持多个几何列，完美适配你的需求，操作也最简单。只需要将保存路径的后缀设为.gpkg即可：

library(sf)

# 示例：创建包含多边形、质心、缓冲区的多几何列sf对象
nc <- st_read(system.file("shape/nc.shp", package="sf"))
nc$centroid <- st_centroid(nc$geometry)
nc$buffer <- st_buffer(nc$geometry, 1000)

# 保存为GeoPackage，delete_layer参数避免图层已存在的报错
st_write(nc, "multi_geometry_data.gpkg", delete_layer = TRUE)

读取时直接用st_read("multi_geometry_data.gpkg")就能完整恢复所有几何列和属性，完全不需要额外处理。而且GeoPackage还支持更大的文件体积、更多属性类型，是当前GIS领域主流的Shapefile替代格式。

2. 拆分保存为多个Shapefile（兼容旧系统） #

如果因为项目要求必须使用Shapefile，最清晰不易出错的方式是把不同几何列拆分为独立的sf对象，然后分别保存为带统一前缀的Shapefile文件，保持属性一致：

# 提取多边形几何列并保存
st_write(st_set_geometry(nc, "geometry"), "nc_polygons.shp")

# 提取质心几何列并保存
st_write(st_set_geometry(nc, "centroid"), "nc_centroids.shp")

# 提取缓冲区几何列并保存
st_write(st_set_geometry(nc, "buffer"), "nc_buffers.shp")

你可以把这些文件放在同一个文件夹下，通过文件名前缀关联，读取时也能轻松对应属性信息。

3. 将额外几何列转为WKT字符串保存（折中方案） #

如果不想拆分文件，还可以把额外的几何列转换为WKT（Well-Known Text）格式的字符串，作为普通属性列保存到Shapefile中，读取时再还原为几何对象：

# 将几何列转为WKT字符串
nc$centroid_wkt <- st_as_text(nc$centroid)
nc$buffer_wkt <- st_as_text(nc$buffer)

# 移除原始的sfc类型几何列，只保留主几何列
nc_simple <- nc[, !names(nc) %in% c("centroid", "buffer")]

# 保存为Shapefile
st_write(nc_simple, "nc_with_wkt.shp")

# 读取时还原几何列
nc_read <- st_read("nc_with_wkt.shp")
nc_read$centroid <- st_as_sfc(nc_read$centroid_wkt)
nc_read$buffer <- st_as_sfc(nc_read$buffer_wkt)

这种方式适合需要单个Shapefile文件，但能接受读取时额外转换步骤的场景。

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