使用sf包st_intersection无法沿乌拉尔山脉拆分俄罗斯多边形？

问题出在哪？ #

你当前的思路有两个核心问题：

1. 把俄罗斯多边形转成多线串再求交完全是冗余操作，还会丢失多边形的拓扑属性——拆分多边形不需要先转线，直接针对多边形+切割线操作才是正确逻辑。
2. 即使线线求交得到了点，这些点也无法直接用来拆分多边形；同时投影精度、几何拓扑无效等问题，也会导致st_intersection返回的交点不符合预期。

正确操作步骤 #

1. 修复拓扑：用st_make_valid处理所有几何对象，解决自相交、无效环等常见问题——这是sf操作的前置必备步骤，很多计算失败都源于拓扑无效。
2. 适配切割线：把乌拉尔线裁剪到俄罗斯多边形范围内，给多边形加小缓冲避免线刚好卡边界，确保切割线能完全穿过多边形。
3. 拆分+提取：用st_split拆分多边形，再用st_collection_extract取出拆分后的两个子多边形。

修正后的代码 #

rm(list=ls())
library(rworldmap)
library(tidyverse)
library(sf)

# 设置工作路径（根据你的本地路径调整）
setwd("/Users/mr56267/Documents/UT_web/TLAH_Maps_2023/Mapping_textbook")

# 读取并转换世界地图数据
world_sp <- fortify(getMap())
world_sf <- world_sp %>% 
  st_as_sf(coords = c("long", "lat"), crs = 4326, row.names="group") %>%
  group_by(group) %>% 
  summarise(geometry = st_combine(geometry)) %>% 
  st_cast("POLYGON") %>%
  st_make_valid() # 修复拓扑

# 转换为罗宾逊投影
world_robinson <- st_transform(world_sf, 
                               crs = '+proj=robin +lon_0=0 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs')

# 提取俄罗斯多边形并修复拓扑
Russia.sf <- world_robinson %>% 
  filter(group == "Russia.1") %>%
  st_make_valid()

# 读取乌拉尔边界线并转换投影、修复拓扑
Urals.sf <- st_read("https://raw.githubusercontent.com/sashatrubetskoy/asia_europe_border/master/asia_europe_border.geojson") %>%
  st_transform(crs = st_crs(Russia.sf)) %>%
  st_make_valid()

# 裁剪乌拉尔线到俄罗斯多边形范围内（加1000米缓冲确保线完全穿过多边形）
Urals_clipped.sf <- st_intersection(Urals.sf, st_buffer(Russia.sf, 1000))

# 拆分俄罗斯多边形
split_russia <- st_split(Russia.sf, Urals_clipped.sf)
# 提取拆分后的东西两部分多边形
russia_west <- st_collection_extract(split_russia, "POLYGON")[[1]]
russia_east <- st_collection_extract(split_russia, "POLYGON")[[2]]

# 可视化验证结果
ggplot() +
  geom_sf(data = russia_west, fill = "lightblue", alpha = 0.7) +
  geom_sf(data = russia_east, fill = "lightgreen", alpha = 0.7) +
  geom_sf(data = Urals_clipped.sf, color = "red", size = 1)

关键细节说明 #

• st_make_valid：处理几何对象的拓扑错误，比如自相交、无效环，这是避免sf计算报错的核心操作。
• st_buffer(Russia.sf, 1000)：给多边形加1000米缓冲（和投影单位一致），确保裁剪后的乌拉尔线能完全触及多边形内部，避免因线刚好落在边界上导致拆分失败。
• st_collection_extract：st_split返回的是几何集合，需要从中提取多边形元素，得到两个独立的俄罗斯子区域。

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