R语言sf包处理shapefile将几何字段转为经纬度列的问题

你现有代码仅完成了空间数据读取、投影转换、子集筛选和几何对象打印操作，未执行任何几何坐标提取逻辑，因此不会自动生成独立的经纬度列。
此外需要注意：MULTIPOLYGON（多面）是面状矢量要素，单个要素可包含多个多边形、每个多边形由多组边界顶点坐标构成，无法直接自动拆分为单值的经纬度列，需要根据实际业务需求选择对应的提取方式。你调用的as.data.frame()仅会将几何列以嵌套列表形式存入数据框，不会自动拆分坐标值。

根据两类常见需求，提供对应可直接运行的代码：

方案1：提取每个多面要素的质心经纬度（单要素对应单对经纬度） #

这是最常用的场景，适合需要给每个Reef斑块匹配唯一代表坐标的需求，提取后不会破坏原始多面几何结构，每个要素仍保留一行：

library(sf)
library(dplyr)

# 读取原始shp数据
shp <- st_read("/home/rdfleay/Desktop/R/WestAustralia/GeographeBay/data/gb_gensub50a.shp")
# 转换为WGS84经纬度坐标系（EPSG:4326）
shp_4326 <- st_transform(shp, "EPSG:4326")
# 筛选底物类型为Reef的要素
reef_GBay <- shp_4326 %>% 
  filter(Substrate == "Reef")

# 提取质心经纬度，生成独立lon、lat列
reef_GBay <- reef_GBay %>%
  # 计算每个多面的质心点
  mutate(geom_centroid = st_centroid(geometry)) %>%
  # 从质心点中提取经纬度：X对应经度lon，Y对应纬度lat
  mutate(
    lon = st_coordinates(geom_centroid)[, 1],
    lat = st_coordinates(geom_centroid)[, 2]
  ) %>%
  # 可选：移除临时生成的质心几何列，保留原始多面geometry
  select(-geom_centroid)

运行后调用head(reef_GBay)即可查看新增的经纬度列，如果需要转为不带空间几何属性的普通数据框，在代码末尾加%>% st_drop_geometry()即可。

方案2：提取多面所有边界顶点的经纬度（单顶点对应一行） #

如果需要获取构成每个多面斑块的所有边界点坐标，可先将多面拆解为单个点要素再提取坐标，注意该操作会扩张数据行数，每个边界顶点单独占一行：

reef_GBay_vertices <- reef_GBay %>%
  # 依次将多面拆解为多点、单点要素
  st_cast("MULTIPOINT") %>%
  st_cast("POINT") %>%
  # 提取每个顶点的经纬度
  mutate(
    lon = st_coordinates(geometry)[, 1],
    lat = st_coordinates(geometry)[, 2]
  )

• 你之前运行的st_transform(shp, crs = st_crs(shp))属于无效操作：投影转换需要指定和原始坐标系不同的目标CRS，转换为要素自身坐标系不会对坐标值产生任何修改，自然无法从米单位转为经纬度单位。
• EPSG:4326坐标系下，st_coordinates()返回的第一列为X轴对应经度，第二列为Y轴对应纬度，注意不要颠倒经纬度顺序。

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