求助：用ggplot复刻墨卡托投影地图，解决顶部裁剪与比例问题

问题根源 #

你遇到的顶部裁剪，本质是coord_map()处理墨卡托投影时，对高纬度区域的坐标计算存在局限，手动设置ylim后绘图区域未预留足够空间；而墨卡托投影的宽高比需严格匹配其数学特性，不能随意拉伸。

修复方案 #

方案1：改用coord_sf（推荐，更稳定） #

ggplot2目前更推荐用sf包处理空间数据，投影计算更精准，能从根源避免裁剪问题：

library(ggplot2)
library(sf)
# 先安装依赖包：install.packages("rnaturalearth")
world_sf <- ne_countries(scale = "medium", returnclass = "sf")

ggplot(world_sf) +
  geom_sf(fill = "gray80", color = "white") +
  # 指定墨卡托投影（EPSG:3857），匹配目标坐标范围
  coord_sf(
    crs = st_crs(3857),
    xlim = c(-20037508, 20037508), # 墨卡托全球标准x范围
    ylim = c(-15000000, 15000000) # 对应你设置的lat范围(-52,83.6)的投影后y值
  ) +
  # 强制匹配墨卡托投影的宽高比
  theme(aspect.ratio = 1)

方案2：修复原有coord_map代码 #

若坚持使用maps包的map_data，可调整参数并强制宽高比：

library(maps)
library(ggplot2)
world <- map_data("world")

ggplot(world, aes(long, lat, group = group)) +
  geom_polygon(fill = "gray80", color = "white") +
  # 指定墨卡托投影，调整ylim避免接近投影极限，开启expand预留空间
  coord_map(
    projection = "mercator",
    xlim = c(-180, 180),
    ylim = c(-52, 83), # 缩小顶部范围，避开墨卡托投影极限（~85°）
    expand = TRUE # 关键：自动预留边缘空间防止裁剪
  ) +
  # 根据纬度范围匹配宽高比，此值适配你设置的lat区间
  theme(aspect.ratio = 1.3)

关键说明 #

• 墨卡托投影在纬度接近85°时会无限拉伸，ylim不要设置太接近该值，否则会出现裁剪或异常变形。
• expand = TRUE会自动在绘图区域边缘预留缓冲空间，直接解决顶部裁剪问题。
• 墨卡托投影的宽高比并非固定1:1：纬度范围越小，宽高比越接近1；纬度范围越靠近两极，宽高比需调大。

内容的提问来源于stack exchange，提问作者Steve'sConnect