R语言Eckert投影绘制问题：跨180°（0-360范围）出现大量水平线

使用R语言绘制Eckert IV投影地图时，将经度范围从180°到180°转换为0到360°后，地图上出现大量水平线，代码及问题截图如下：

install.packages("sf")
install.packages("ggplot2")
install.packages("rnaturalearth")
install.packages("rnaturalearthdata")

library(sf)
library(ggplot2)
library(rnaturalearth)
library(rnaturalearthdata)

world <- ne_countries(scale = "medium", returnclass = "sf",)
proj_params <- "+proj=eck4 +lon_0=181 +datum=WGS84"
world_wrapped <- st_wrap_dateline(world, options = c("WRAPDATELINE=YES"))

world_eckert <- st_transform(world_wrapped, crs = proj_params)


p <- ggplot() +
  geom_sf(data = world_eckert, fill = "lightyellow", color = "black") +
  labs(title = "Eckert IV Projection World Map (Centered at 180°)") +
  theme_minimal() +
  theme(plot.title = element_text(hjust = 0.5))

出现水平线的核心原因是跨日界线的多边形在投影转换时出现错误连接，以下是针对性修复方法：

方法1：优化日界线包裹参数 #

给st_wrap_dateline添加DATELINEOFFSET=180参数，明确日界线位置，避免多边形错误拆分：

world_wrapped <- st_wrap_dateline(world, options = c("WRAPDATELINE=YES", "DATELINEOFFSET=180"))

方法2：修正投影中心 #

代码中投影中心设为181°，改为标准180°可避免边界计算异常：

proj_params <- "+proj=eck4 +lon_0=180 +datum=WGS84"

方法3：预处理修复无效多边形 #

部分原始多边形可能存在拓扑错误，转换前先修复有效性：

world_clean <- st_make_valid(world)
world_wrapped <- st_wrap_dateline(world_clean, options = c("WRAPDATELINE=YES", "DATELINEOFFSET=180"))

完整修复代码 #

install.packages("sf")
install.packages("ggplot2")
install.packages("rnaturalearth")
install.packages("rnaturalearthdata")

library(sf)
library(ggplot2)
library(rnaturalearth)
library(rnaturalearthdata)

world <- ne_countries(scale = "medium", returnclass = "sf")
# 修复多边形拓扑有效性
world_clean <- st_make_valid(world)
# 设置标准180°中心的Eckert IV投影
proj_params <- "+proj=eck4 +lon_0=180 +datum=WGS84"
# 正确处理日界线包裹
world_wrapped <- st_wrap_dateline(world_clean, options = c("WRAPDATELINE=YES", "DATELINEOFFSET=180"))

world_eckert <- st_transform(world_wrapped, crs = proj_params)

p <- ggplot() +
  geom_sf(data = world_eckert, fill = "lightyellow", color = "black") +
  labs(title = "Eckert IV Projection World Map (Centered at 180°)") +
  theme_minimal() +
  theme(plot.title = element_text(hjust = 0.5))

print(p)

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