R语言sf对象裁剪时出现TopologyException拓扑异常问题求助

成因 #

1. 拓扑错误触发报错：maps::map提供的湖泊数据中，在接近北纬50°的位置存在**嵌套面（Nested shells）**的拓扑缺陷。当st_crop裁剪范围设置为ymax=50时，刚好截取到这个有问题的区域，导致生成的几何对象无效，后续st_transform处理时触发TopologyException报错。
2. 平面几何处理的局限性：你关闭了S2球面几何（sf_use_s2(FALSE)），此时sf使用平面几何引擎（GEOS）处理经纬度数据，不仅会弹出"assumes planar"的警告，还对拓扑错误的容忍度极低，直接触发报错。而将ymax改为45时，裁剪范围避开了这个有拓扑缺陷的区域，因此代码能正常运行。

解决方法 #

方法1：修复几何有效性 #

在裁剪后、投影转换前，用st_make_valid()修复拓扑错误，这是最直接的解决方案：

library(ggplot2)
library(sf)
library(tidyverse)
library(showtext)
library(showtextdb)
library(ggtext)

sf_use_s2(FALSE)

lakes <- maps::map("lakes", fill=TRUE, plot =FALSE) %>%
  st_as_sf() %>%
  st_crop(xmin = -140, xmax = -55, ymin = 15, ymax = 50) %>%
  st_make_valid() %>%  # 修复拓扑错误
  st_transform(crs = st_crs("+proj=laea +lat_0=45 +lon_0=-100"))

方法2：启用S2球面几何处理 #

开启S2球面几何，它对球面几何的拓扑错误容忍度更高，且更适配经纬度数据的处理：

library(ggplot2)
library(sf)
library(tidyverse)
library(showtext)
library(showtextdb)
library(ggtext)

sf_use_s2(TRUE)  # 启用S2球面几何

lakes <- maps::map("lakes", fill=TRUE, plot =FALSE) %>%
  st_as_sf() %>%
  st_crop(xmin = -140, xmax = -55, ymin = 15, ymax = 50) %>%
  st_transform(crs = st_crs("+proj=laea +lat_0=45 +lon_0=-100"))

方法3：先转换投影再裁剪 #

先将湖泊数据转换到目标投影（LAEA），再进行裁剪，避免在经纬度平面处理时的拓扑问题：

library(ggplot2)
library(sf)
library(tidyverse)
library(showtext)
library(showtextdb)
library(ggtext)

sf_use_s2(FALSE)

target_crs <- st_crs("+proj=laea +lat_0=45 +lon_0=-100")
# 先转换投影，再裁剪
lakes <- maps::map("lakes", fill=TRUE, plot =FALSE) %>%
  st_as_sf() %>%
  st_transform(target_crs) %>%
  st_crop(st_bbox(c(xmin = -140, xmax = -55, ymin = 15, ymax = 50) %>% 
                    st_sfc(crs = 4326) %>% 
                    st_transform(target_crs)))

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