如何将GeoJSON格式坐标转换为R sp所需的坐标格式？

Hey there! 很高兴你觉得《多边形重叠百分比》里的方法好用～ 针对你提到的几百个GeoJSON多边形转R sp坐标格式的问题，我整理了两种靠谱的方法，帮你高效搞定：

方法一：用sf包（推荐，最省心） #

sf是R里处理空间数据的现代工具，对GeoJSON支持非常友好，还能和sp包无缝转换，完全不用手动折腾嵌套坐标的问题。步骤超简单：

1. 先安装并加载sf包
2. 直接读取GeoJSON文件为sf对象
3. 一键转成sp格式

代码示例：

# 安装sf包（第一次用需要安装）
install.packages("sf")
library(sf)

# 读取单个GeoJSON文件
sf_polygon <- st_read("your_single_polygon.geojson")

# 转换为sp对象
sp_polygon <- as(sf_polygon, "Spatial")

# 批量处理文件夹里的所有GeoJSON
# 先获取所有GeoJSON文件路径
geo_file_paths <- list.files(
  path = "your_geojson_directory", 
  pattern = "\\.geojson$", 
  full.names = TRUE
)

# 循环读取并转换
sp_polygon_list <- lapply(geo_file_paths, function(file_path) {
  # quiet=TRUE避免打印过多信息
  sf_obj <- st_read(file_path, quiet = TRUE)
  as(sf_obj, "Spatial")
})

方法二：手动处理坐标（适合不想依赖sf的场景） #

如果不想用sf包，也可以手动解析GeoJSON的嵌套坐标列表，转换成sp要求的矩阵格式（sp的Polygon对象需要每行是(x,y)的矩阵，且首尾点必须相同，保证多边形闭合）。

代码示例：

# 加载jsonlite用来读取GeoJSON
install.packages("jsonlite")
library(jsonlite)
library(sp)

# 读取单个GeoJSON文件
geo_json_data <- fromJSON("your_polygon.geojson")

# 提取坐标：注意GeoJSON的坐标是多层嵌套，索引根据你的数据结构调整
# 一般结构是 features$geometry$coordinates，MultiPolygon会多一层
coords_nested <- geo_json_data$features$geometry$coordinates[[1]][[1]]

# 转换为矩阵，并确保多边形闭合（首尾点一致）
if (!all(coords_nested[1,] == coords_nested[nrow(coords_nested),])) {
  coords_nested <- rbind(coords_nested, coords_nested[1,])
}
coords_matrix <- as.matrix(coords_nested)

# 创建sp的Polygon和SpatialPolygons对象
single_poly <- Polygon(coords_matrix)
poly_collection <- Polygons(list(single_poly), ID = "poly_1") # 给个唯一ID
sp_polygon <- SpatialPolygons(list(poly_collection))

# 批量处理的话，循环每个文件，注意给每个多边形设置不同的ID即可

小提示 #

• 如果你处理的是MultiPolygon类型的GeoJSON，需要循环提取每个子多边形的坐标，分别创建Polygon对象，再组合到Polygons里
• 始终确保转换后的多边形是闭合的，否则sp会报错

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