R语言重投影自定义中央经线地图时如何避免边界回绕异常

异常原因 #

原始全球国家矢量数据的默认经度范围为-180°~180°，当你将中央经线设置为130°时，新投影对应的经度分割线为130°-180°=-50°，该线正好穿过南美大陆。未被切割的跨分割线多边形在重投影时，会被错误地从地图一侧拉伸到另一侧，从而产生图形重叠异常。

修复代码 #

只需要在加载完全球数据后，先对矢量数据做适配新中央经线的切割处理即可，完整代码如下：

library("ggplot2")
theme_set(theme_bw())
library("sf")
library("rnaturalearth")
library("rnaturalearthdata")

# 加载全球国家数据
world <- ne_countries(scale = "medium", returnclass = "sf")

# 核心修复：根据目标中央经线，切割跨经度分割线的多边形
# 参数lon_0和你coord_sf里设置的中央经线保持一致即可
world_cut <- st_wrap_dateline(world, options = c("WRAPDATELINE=YES", "DATELINEOFFSET=50"))

# 绘图时使用处理后的数据即可
ggplot(data = world_cut) +
    geom_sf() +
    coord_sf(crs = "+proj=cea +lon_0=130 +x_0=0 +y_0=0 +lat_ts=45 +ellps=WGS84 +datum=WGS84 +units=m +no_defs")

参数说明 #

• st_wrap_dateline是sf包内置的专门处理日期变更线切割的函数
• 参数里的DATELINEOFFSET=50对应你当前的分割线位置-50°（和130°中央经线正好相差180°），如果你后续修改了中央经线，只需要把这个值调整为180 - 你的中央经线绝对值即可：比如中央经线设为150°，offset就填30
• 处理后的数据会自动把跨分割线的多边形拆成两部分，重投影后不会出现跨整个地图的拉伸重叠，同时你需要的南美区域回绕到左侧的效果也会保留

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