使用 geom_sf 跨越日期变更线时地图上的封闭边界

Question

我正在尝试制作一个罗宾逊投影图，它不是以经度 0 为中心，并且穿过 180° 线的两侧。

我的代码如下：

library("rnaturalearthdata")
world <- ne_countries(scale = "medium", returnclass = "sf")
robinson = "+proj=robin +lon_0=-90 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs"
transformed = st_transform(world, robinson)
ggplot(data=transformed) + geom_sf(fill = "gray", colour = "black")+theme_bw()

生成此地图，该地图的方向正确，但所有分割边界都是“由内而外”并环绕整个地图。

9年前就有这个相关问题，但我希望通过最近的软件包更新，有一个不太复杂的解决方案，这将允许我们使用

geom_point()

绘制纬度/经度数据点，而无需额外的转换.

编辑＃2：L. Tyrone 的解决方案对我来说非常有用。要使用此方法绘制点，您必须将它们从数据框转换为具有相同坐标参考系统的 sf 对象，否则纬度/经度将被解释为米，最终与 0,0 无法区分：

# Transform with fake data...
fakedata=data.frame(lat=c(-10,-20,80),lon=c(45,-100,120),tag=c("indian","tropical","arctic"))
transpoint = st_as_sf(fakedata,coords=c("lon","lat"),crs=4326)
dtran = st_transform(transpoint,robinson)

然后您可以通过将特殊的

geometry=

和

stat=

声明添加到

geom_point

甚至

geom_text_repel

来绘制它们：

geom_point(data=dtran,aes(geometry=geometry,color=tag),stat="sf_coordinates") +
geom_text_repel(data=dtran,aes(geometry=geometry,label=tag),stat="sf_coordinates")

最终结果：

Answer 1

在投影之前，您需要“破坏”任何交叉 180 度的多边形。

st_break_antimeridian()

包中的

sf

函数是实现此目的的最简单方法。它会发出警告，您可以安全地忽略该警告。

library(rnaturalearth)
library(sf)
library(dplyr)
library(ggplot2)

world <- ne_countries(scale = "medium", returnclass = "sf") %>%
  st_set_crs(4326)

robinson <- "+proj=robin +lon_0=-90 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs"

world_robinson_90w <- world %>%
  st_break_antimeridian(lon_0 = -90) %>%
  st_transform(crs = robinson)

ggplot() +
  geom_sf(data = world_robinson_90w)