ggplot2组合geom_map与geom_point的美学映射报错排查

这个问题是ggplot2多图层叠图时非常典型的美学映射继承错误：

1. 你大概率把label = year、甚至经纬度这类只属于地震数据集earthquakes1的映射写在了最外层ggplot()的全局aes()里。
2. 单独运行geom_point和文本标注层时，全局映射的字段都能在earthquakes1里找到，所以能正常出图。
3. 叠加geom_map层时，该图层默认会继承全局所有美学映射，但你传入的边界数据world1里根本没有year字段；加上geom_map内部会调用坐标处理函数，找不到对应列时会误把函数对象当成映射值传入，就会先弹出"Don't know how to automatically pick scale for object of type function"的警告，后续渲染文本标注时拿不到合法的label列值，就直接抛出列无效的报错。

选任意一种即可解决问题：

• 方案1：分层写映射。不要把单个图层才用到的字段写在全局aes里，geom_map、geom_point、文本标注层各自在层内的aes()里只写自己数据集存在的字段，从根源避免错误继承。
• 方案2：给geom_map层加参数inherit.aes = FALSE，强制该图层不继承全局设置的所有美学映射，只用层内自定义的映射规则。

library(ggplot2)
# 可选，加载ggrepel避免地震点的年份标注重叠
# library(ggrepel)

ggplot() +
  # 海地行政边界图层
  geom_map(
    data = world1,
    map = world1,
    # 只配置边界数据自带的字段映射
    aes(map_id = region, fill = subregion),
    color = "gray40",
    linewidth = 0.3,
    # 如果坚持把映射写在全局，打开下面这行注释即可
    # inherit.aes = FALSE
  ) +
  # 地震散点图层
  geom_point(
    data = earthquakes1,
    # 散点专属的经纬度、震级映射
    aes(x = Longitude, y = Latitude, size = Mag),
    color = "#c0392b",
    alpha = 0.8
  ) +
  # 地震年份标注图层
  geom_text(
    # 换成geom_text_repel可自动避让重叠文本
    data = earthquakes1,
    aes(x = Longitude, y = Latitude, label = year),
    size = 3,
    nudge_y = 0.05
  ) +
  # 固定地图坐标比例，避免边界拉伸变形
  coord_map() +
  # 按需裁剪显示范围到海地周边
  xlim(range(earthquakes1$Longitude) + c(-1,1)) +
  ylim(range(earthquakes1$Latitude) + c(-1,1)) +
  theme_bw()

• 多数据源叠图时，最稳妥的写法是每个geom层单独指定自己的data参数和专属aes映射，全局ggplot()里只放所有图层共用的映射，能规避绝大多数跨图层的映射报错。
• 如果后续处理空间数据，更推荐把边界数据转成sf格式后用geom_sf()绘图，比老版本的geom_map()兼容性更好，映射逻辑更清晰，很少出现这类奇怪的继承问题。

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