R语言tmap多图层地图实现单个统一图例的方法问询

合并重复图例的实现方法 #

你遇到的多图层重复图例问题可以通过以下两步解决：

第一步：固定全局统一的颜色映射规则 #

首先基于所有年份的完整数据计算统一的数值范围/分箱，保证同一个颜色在所有年份对应相同的数值，避免不同图层颜色映射逻辑不一致：

# 计算所有数据的全局数值范围，可根据需求自定义分箱规则
val_min <- min(pesticides$values, na.rm = TRUE)
val_max <- max(pesticides$values, na.rm = TRUE)
# 示例：自定义5个等距分箱，也可以使用分位数分箱等其他规则
val_breaks <- seq(val_min, val_max, length.out = 6)

第二步：调整tmap绘图参数合并图例 #

对于你使用as.layers = TRUE生成交互式切换图层的场景，修改绘图代码如下即可合并为单个图例：

mx <- tm_shape(World)+
  tm_view(set.view = c(9.7, 45.5, 8), view.legend.position = c("left", "top"))+
  tm_basemap("Esri.WorldTopoMap")+
  tm_borders("white", lwd = .5) +
  tm_shape(pesticides) +
  tm_dots(
    col="values",
    # 加入固定的分箱规则，保证所有图层颜色映射一致
    breaks = val_breaks,
    legend.hist = TRUE,
    size=0.3,
    legend.show = TRUE,
    popup.vars=c("Id Station"="station",
                 "Valore"= "values",
                 "Substance" = "PARAM")
  )+
  tm_facets(by= "year",
            nrow = 9,
            sync = TRUE,
            as.layers = T,
            drop.units = T,
            # 新增参数：合并所有分面/图层的相同图例
            legend.combine = TRUE)

如果你的tmap版本不支持legend.combine参数，也可以只保留最后一个图层的图例显示，其余图层隐藏图例：将tm_dots中的legend.show替换为条件判断legend.show = ~year == max(year)即可。 #

不同年份数据分布差异下的图例选择建议 #

是否需要调整图例取决于你的分析目的：

• 如果你需要跨年份比较农药残留的绝对变化：必须使用统一图例，否则同一颜色在不同年份代表的数值不同，会完全误导对比结论，这也是绝大多数环境监测空间展示场景的选择。
• 如果你仅需要观察单一年份内部的残留高低空间分布，不需要跨年份对比：可以为每个年份设置独立图例，但需要在图层说明处明确标注每个年份的数值范围，避免使用者混淆。

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