R语言sf包判断跨180°日期变更线几何范围内的点位归属方法

问题本质 #

• 判断失效的核心原因并非sf球面运算版本不足，而是NOAA导出的夏威夷EEZ几何被±180°子午线硬切割为两个独立的多边形碎块；st_shift_longitude()仅完成坐标值的线性平移（将[-180,180]区间转换为[0,360]区间），不会自动缝合碎块的拓扑关系，运算时仍会把两个碎块判定为完全分离的独立几何，自然无法匹配点位位置。
• 此前尝试的EPSG:2163投影方案失效，是因为该投影的有效覆盖范围仅为北美及邻近大西洋区域，研究区所在的中太平洋位置已超出投影适用边界，几何形变会完全扭曲原始边界的真实形状。

正确操作流程 #

1. 统一坐标基准并合并碎块几何 #

先对多边形和点位做完全一致的经度范围平移，再合并被切开的独立碎块为连续整体：

library(sf)
library(tidyverse)

# 读取并筛选目标EEZ要素
eez_raw <- read_sf("USMaritimeLimitsAndBoundariesSHP/USMaritimeLimitsNBoundaries.shp") %>%
  filter(OBJECTID == 1206) %>%
  st_transform(4326) %>%
  st_shift_longitude() # 统一平移到[0,360]经度区间，与world2地图基准对齐

# 关键步骤：合并被日期变更线劈开的碎块，修复拓扑错误
eez_merged <- eez_raw %>%
  st_union() %>% # 合并独立几何碎块为整体
  st_make_valid() %>% # 修复shapefile自带的自相交、微小间隙问题
  st_sf() # 转换为标准sf对象适配后续运算

# 采样点位做和多边形完全相同的坐标转换
pnts_processed <- pnts %>%
  st_shift_longitude()

2. 执行空间归属判断 #

使用预处理完成的几何开展计算，即可得到正确结果：

# 返回每个点位是否落在EEZ边界内的逻辑向量
point_in_eez <- st_within(pnts_processed, eez_merged, sparse = FALSE, prepared = TRUE)[,1]

常见操作误区 #

• 不要跳过碎块合并步骤：仅靠st_shift_longitude()得到的视觉效果对齐，不代表几何拓扑关系连通；此时两个碎块在计算逻辑上仍相距360个经度，完全无法触发相交判断。
• 不要选择适配范围不匹配的投影：若需要使用平面投影计算，可选用覆盖中太平洋区域的投影（如EPSG:3832），但本次研究涉及的EEZ边界形变平缓，合并拓扑后的WGS84几何完全可以直接开展平面运算，无需额外投影转换。
• 遇到偶发判断偏差时，可对几何做快速校验：用st_is_valid()检查几何是否存在自相交、缝隙等错误，修复后再重新运算。

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