如何将含坐标列表列的DataFrame转换为MultiPolygon？

要让st_as_sf()识别你的数据，得先把列表格式的经纬度转换成sf包认可的几何对象结构，核心是把每行的经纬度列表组合成坐标矩阵，再构建MultiPolygon，最后生成标准的几何列。

1. 先加载依赖包 #

确保装了dplyr和sf，没装的话先执行install.packages(c("dplyr", "sf"))，然后加载：

library(dplyr)
library(sf)

2. 逐行转换坐标为MultiPolygon #

假设你的数据框叫df，用rowwise()逐行处理：

df_sf <- df %>%
  rowwise() %>%
  # 把lng和lat列表合并成坐标矩阵（sf默认是经度在前、纬度在后）
  mutate(coords = list(cbind(lng, lat))) %>%
  # 先构建Polygon，再包装成MultiPolygon（统一格式，方便后续操作）
  mutate(geometry = list(st_multipolygon(list(list(coords))))) %>%
  # 把geometry列转成sf认可的sfc类型，同时转成sf对象
  st_sf(sfc_column_name = "geometry") %>%
  # 删掉中间生成的coords列（可选）
  select(-coords)

关键注意点 #

• st_multipolygon()要求的输入是三层嵌套列表：最外层对应MultiPolygon，中间层对应里面的单个Polygon，最内层是Polygon的坐标环（这里只用外环就行）。
• 如果你的坐标没闭合（第一个点和最后一个点不一样），得手动在坐标矩阵末尾复制第一个点，不然可能报错。
• 转完之后的df_sf就是标准sf对象，直接能用来做交互式地图（比如搭配leaflet）。

验证转换是否正确 #

可以用下面的代码检查几何类型：

st_geometry_type(df_sf)

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