使用ggplot2绘制北温带研究站点地图时遇问题求助

我太懂你这个困扰了——用coord_map()设置北温带（23.5°N至66.5°N）纬度范围后地图乱成一团，这其实是ggplot处理多边形投影和范围裁剪时的常见问题，根源在于**coord_map()只是在绘图阶段裁剪显示范围，并没有真正修改底层的多边形数据**，那些跨越你设置的纬度边界的国家/地区多边形会被截断，从而出现线条重叠、错乱的情况。

下面给你两种可靠的解决思路，优先推荐第一种：

方法一：先裁剪空间数据，再绘图（最稳妥） #

我们可以用sf包（空间数据处理的标准工具）先把世界地图数据裁剪到北温带范围内，再用ggplot绘制，这样就能从根源上避免多边形截断的问题。

步骤代码： #

# 加载必要的包
library(ggplot2)
library(sf)
library(dplyr)

# 获取世界地图数据并转换为sf格式（空间数据格式）
world_sf <- st_as_sf(map_data("world"), coords = c("long", "lat"), crs = 4326) %>%
  group_by(group, region) %>%
  summarise(do_union = FALSE) %>%
  st_cast("POLYGON")

# 定义北温带的边界范围（经度-180到180，纬度23.5到66.5）
north_temperate_bbox <- st_bbox(
  c(xmin = -180, xmax = 180, ymin = 23.5, ymax = 66.5), 
  crs = 4326
) %>% st_as_sfc()

# 裁剪地图数据到北温带区域
world_north_temperate <- st_intersection(world_sf, north_temperate_bbox)

# 绘制地图
ggplot() +
  geom_sf(data = world_north_temperate, fill = "#f0f0f0", color = "#333333") +
  coord_sf(ylim = c(23.5, 66.5)) +
  theme_minimal() +
  labs(title = "北温带区域地图")

为什么这个方法有效？ #

st_intersection()会真正对空间多边形进行裁剪，只保留落在北温带范围内的部分，这样绘图时就不会出现截断的多边形线条，地图显示完全正常。

方法二：调整coord_map参数（快速临时解决） #

如果你不想引入新的包，可以尝试指定投影类型并补充完整的经度范围，有时候能缓解错乱问题：

library(ggplot2)

world <- map_data("world")
map <- ggplot() + 
  geom_polygon(data = world, aes(x=long, y = lat, group = group), fill = "#f0f0f0", color = "#333333") +
  # 指定墨卡托投影，同时设置完整的经度范围
  coord_map(projection = "mercator", xlim = c(-180, 180), ylim = c(23.5, 66.5)) +
  theme_minimal()
map

注意事项： #

这种方法只是通过调整投影参数减少错乱，但部分跨纬度边界的多边形还是可能存在显示问题，适合快速测试，长期使用还是推荐方法一。

另外，如果你需要更准确的世界地图边界数据，可以用rnaturalearth包获取最新数据，比如：

library(rnaturalearth)
world_sf <- ne_countries(scale = "medium", returnclass = "sf")

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