如何在sf中合并跨日界线多边形，消除太平洋中心地图的分割线？

核心问题是原始GeoJSON的多边形顶点跨180°经线，即便转换投影或偏移经度，绘制时仍会显示连接日界线的垂直线。以下是几种可行的解决方法：

方法1：用st_wrap_dateline处理日期线跨越 #

该函数可自动调整跨日界线的多边形，将其重新包裹到太平洋中心视角的经度范围，同时消除不必要的分割线：

library(tidyverse)
library(sf)

# 读取数据
regions <- st_read("https://gist.github.com/mhoban/ebdee566c4d1963e9fcd88fae818fc83#file-hawaii-geojson")

# 处理日期线，指定以180°为基准的包裹规则
regions_fixed <- regions %>%
  st_wrap_dateline(options = c("WRAPDATELINE=YES", "DATELINEOFFSET=180")) %>%
  st_set_crs(4326) %>%
  st_shift_longitude()

# 绘制测试
ggplot(regions_fixed %>% st_transform(3832)) + 
  geom_sf(fill = "lightgrey", linewidth = 1)

方法2：临时关闭S2几何处理 #

若S2几何引擎的精度逻辑导致分割线残留，可临时关闭S2后修复合并几何：

# 临时关闭S2几何计算
sf_use_s2(FALSE)

regions_fixed <- regions %>%
  st_shift_longitude() %>%
  st_make_valid() %>%
  st_union() %>%
  st_sf()

# 可选：恢复S2设置
sf_use_s2(TRUE)

# 绘制测试
ggplot(regions_fixed %>% st_transform(3832)) + 
  geom_sf(fill = "lightgrey", linewidth = 1)

方法3：手动调整多边形顶点坐标 #

对跨180°的多边形，手动将西经（负数）经度值加360，转换为东经范围后重构几何：

# 定义坐标调整函数
adjust_coords <- function(geom) {
  coords <- st_coordinates(geom)
  # 将西经坐标转换为东经范围
  coords[, "X"] <- ifelse(coords[, "X"] < 0, coords[, "X"] + 360, coords[, "X"])
  # 重新构建MULTIPOLYGON
  st_multipolygon(list(
    lapply(split(coords, coords[, "L1"]), function(p) {
      matrix(p[, c("X", "Y")], ncol = 2)
    })
  ))
}

# 应用坐标调整
regions_fixed <- regions %>%
  mutate(geometry = map(geometry, adjust_coords)) %>%
  st_set_crs(4326) %>%
  st_transform(3832)

# 绘制测试
ggplot(regions_fixed) + 
  geom_sf(fill = "lightgrey", linewidth = 1)

关键说明 #

• st_wrap_dateline是最简洁的专用方案，推荐优先使用。
• 关闭S2适合复杂多边形场景，但可能影响其他空间运算精度，使用后建议恢复设置。
• 手动调整坐标适合需要完全控制几何修改的场景，需注意处理MULTIPOLYGON的层级结构。

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