为什么某些 xlims 和 ylims 在 ggplot 和 sf 中产生此错误？

我正在学习使用 ggplot 和 sf 来处理空间数据。

当我尝试绘制以下图时，出现错误：

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

world <- ne_countries(scale = "medium", returnclass = "sf")

ggplot(data = world) +
   geom_sf() +
   coord_sf(xlim = c(-102.15, -74.12), ylim = c(0, 33.97), expand = FALSE)

# Error in st_cast.POINT(x[[1]], to, ...) : 
#  cannot create MULTILINESTRING from POINT

然而，如果我只是在任一方向上稍微调整 ylims，它就会起作用！：

ggplot(data = world) +
   geom_sf() +
   coord_sf(xlim = c(-102.15, -74.12), ylim = c(0.01, 33.97), expand = FALSE)

or

ggplot(data = world) +
   geom_sf() +
   coord_sf(xlim = c(-102.15, -74.12), ylim = c(-0.01, 33.97), expand = FALSE)

所以这是一个很好的解决方案，但我想知道问题是什么？在最后的情节中，我注意到一些小岛屿悄悄进入了地图（加拉帕戈斯群岛？） 这个多边形的尖端看起来像一个点而不是一个多边形，是否以某种方式干扰了第一段代码？或者发生了其他事情，我怎样才能更优雅地解决呢？

我不相信你描述的问题是由加拉帕戈斯群岛引起的 - 因为它存在于scale = "small"其中省略了这些岛屿。

我尝试了世界国家数据集的另一个来源（来自自然地球的数据集不是球形有效的，所以我想排除这个问题）。错误仍然存​​在，现在我怀疑它与中的代码有关{ggplot2}- 特别是在默认情况下它如何处理赤道（即零）作为限制expand = TRUE被用户覆盖（可能值得在他们的 github 上提交问题）。

同时我建议两种可能的方法：

• 在 ggplot 调用之外裁剪您感兴趣的区域（出于多种原因，这可能更可取 - 请参阅ggplot2 和 sf：geom_sf_text 在 coord_sf 设置的限制内 https://stackoverflow.com/questions/66031935/ggplot2-and-sf-geom-sf-text-within-limits-set-by-coord-sf/66032492#66032492)
• 避免expand = FALSE在你的coord_sf() call

对于在数据级别（而不是在 ggplot 调用中的表示级别）裁剪世界对象，请考虑以下代码：

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

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

bounds <- matrix(c(-102.15, 0, 
                   -74.12, 0, 
                   -74.12, 33.97, 
                   -102.15, 33.97,
                   -102.15, 0), 
                 byrow = TRUE, ncol = 2) %>%
  list() %>% 
  st_polygon() %>% 
  st_sfc(crs = 4326) 

small_world <- st_intersection(world, bounds) 

ggplot(data = small_world) +
  geom_sf()