R语言sf包st_wrap_dateline处理跨日界线多边形未按预期拆分问题

我尝试生成0级H3六边形网格（共122个单元格），但跨国际日界线的多边形出现反转、远离日界线的异常情况。原本期望用st_wrap_dateline()沿日界线拆分这些多边形，但未达到预期效果，代码和异常效果如下：

install.packages("h3jsr")
install.packages("sf")

library("h3jsr")
library("sf")

# Make number sequence
pts  = as.data.frame(merge(seq (-90,90,1),seq (0,360,1),all=TRUE))
colnames(pts) <- c("LAT","LON")

pts_ <- st_as_sf(x = pts, coords = c("LON","LAT"),crs = 4326)

# resolve H3 index at Resolution
pts_H3 <- point_to_h3(pts_, res = 0, simple = FALSE)

pts_H3 <- as.data.frame(table(pts_H3$h3_resolution_0))

getH3_Poly <- h3_to_polygon(pts_H3$Var1,simple=FALSE)

# timeline wrap
getH3_Poly <- st_wrap_dateline(getH3_Poly, options = c("WRAPDATELINE=YES"))
plot(getH3_Poly)

问题根源 #

1. H3多边形方向冲突：H3返回的多边形坐标为逆时针顺序，而WGS84（EPSG:4326）坐标系中，sf默认将逆时针环解析为“反多边形”（覆盖除六边形本身外的整个地球），视觉上表现为多边形反转、远离日界线。
2. 采样流程冗余误差：通过点采样统计H3索引的方式可能遗漏部分单元格，不如直接生成所有0级H3单元格准确。
3. st_wrap_dateline参数不足：仅指定WRAPDATELINE=YES未明确日界线位置，无法触发正确的拆分逻辑。

修复步骤与代码 #

1. 直接生成所有0级H3单元格

跳过点采样，直接获取完整的0级H3索引列表：

install.packages("h3jsr")
install.packages("sf")

library(h3jsr)
library(sf)

# 获取所有0级H3索引（共122个）
all_h3_res0 <- get_res0_indexes()

2. 转换多边形并修复方向

将H3索引转为多边形后，用st_make_valid()自动修正多边形环方向，解决反转问题：

getH3_Poly <- h3_to_polygon(all_h3_res0, simple = FALSE)

# 修复多边形有效性，自动调整环方向
getH3_Poly <- st_make_valid(getH3_Poly)

3. 配置st_wrap_dateline参数拆分日界线

添加DATELINEOFFSET=180明确指定日界线位置，触发跨边界多边形的拆分：

# 沿180°经线拆分跨日界线的多边形
getH3_Poly_wrapped <- st_wrap_dateline(
  getH3_Poly,
  options = c("WRAPDATELINE=YES", "DATELINEOFFSET=180")
)

# 查看修复后的网格
plot(getH3_Poly_wrapped)

关键说明 #

• get_res0_indexes()是h3jsr包提供的原生函数，可直接获取完整的0级H3单元格列表，避免采样误差。
• st_make_valid()会自动检测并修正多边形的环方向，确保sf将其识别为正确的六边形区域，而非反区域。
• DATELINEOFFSET=180参数明确告知GDAL以180°经线为日界线，确保跨边界多边形被正确拆分为两个部分。

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