如何将含坐标的数据框转换为Simple Feature数据框？

嘿，这个错误我之前踩过坑！大概率是你在生成Simple Feature（sf）对象时，分组操作和几何列的长度不匹配导致的——尤其是误用了mutate而非summarise来处理分组后的几何特征，或者手动创建几何列时没对齐行与几何对象的数量。

先给你最稳妥的解决方案：直接用st_as_sf() #

这是官方推荐的转换方式，几乎不会出错，不需要手动折腾几何列：

library(sf)
library(dplyr)

# 假设你的原始数据框长这样
df <- tibble(
  location = c("公园", "商场", "地铁站"),
  x = c(116.397, 116.405, 116.412),
  y = c(39.908, 39.912, 39.905),
  crs_code = 4326 # 替换成你实际用的坐标参考系统代码
)

# 一行代码完成转换
sf_df <- df %>%
  st_as_sf(coords = c("x", "y"), crs = df$crs_code[1])

# 查看结果
sf_df

这个函数会自动把x/y列转换成名为geometry的标准几何列，同时帮你设置好坐标参考系统，完全避开手动操作的长度问题。

如果你非要手动创建几何列（比如生成多点/线） #

如果你的需求是把同一地点的多个坐标合并成一个几何对象（比如多点），这时候不能用mutate，得用summarise——因为mutate要求新列的长度和原数据行数一致，而分组后生成的单个几何对象只能对应一行数据：

# 示例：同一地点有多个坐标，要合并成多点
df_multi <- tibble(
  location = c("公园", "公园", "商场"),
  x = c(116.397, 116.398, 116.405),
  y = c(39.908, 39.909, 39.912),
  crs_code = 4326
)

# 用summarise合并分组后的几何对象
sf_multi_df <- df_multi %>%
  group_by(location) %>%
  summarise(
    geometry = st_sfc(st_multipoint(cbind(x, y))),
    .groups = "drop"
  ) %>%
  st_set_crs(4326)

sf_multi_df

为什么你会遇到那个错误？ #

你收到的Column feature must be length 1 (the group size), not 2提示，本质是：

• 你用group_by分组后，每个组的行数是1，但你生成的feature列每个组返回了2个几何对象；
• 或者反过来，每个组有2行，但你只生成了1个几何对象。

mutate要求新列的元素数量和当前数据的行数严格对齐，而如果是要把分组内的多行合并成一个几何对象，必须用summarise来压缩行数，让几何对象和合并后的行一一对应。

最后检查小细节 #

• 确认x/y列是数值型，没有NA值；
• 坐标参考系统（CRS）代码要正确（比如WGS84是4326，Web墨卡托是3857）。

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