带有 tmap 的构面图答案

【问题标题】：Facet Map with tmap带有 tmap 的构面图
【发布时间】：2020-01-31 22:56:37
【问题描述】：

我正在尝试使用“tmap”包创建多面地图。我的数据包含很多缺失的国家。当我运行下面的代码时，输出显示了很多缺失的国家。怎样才能得到显示失踪国家的情节？

library(tmap)

    world<-data("World")

    MAP<- read_excel("SSPMAP.xlsx")
    MAP2<-MAP %>% group_by(iso_a3,scenario ) %>%
      summarise(cibase=median(gr_C_2030)*100)
    table(MAP2$cibase)
    cut<-c(-4,-2,0,2,4,6,8)

    World_n<-left_join(World, MAP2, by="iso_a3") 

    ###########################
    Base<-tm_shape(World_n, projection = "eck4") +
      tm_polygons("cibase", breaks = cut,
                  palette = "RdYlGn", 
                  border.col = "black", 
                  border.alpha = 1,midpoint = NA ) +
      tm_legend(legend.position = c("left", "bottom")) +
      tm_facets(by=c("scenario"), ncol  =1, showNA = FALSE)+
      tm_layout(panel.labels=c("Baseline : SSP1","Baseline : SSP2","Baseline : SSP4","Baseline : SSP5"),panel.label.color = "black"
                ,panel.label.size=1.3,panel.label.fontface="bold")
    Base
    tmap_save(Base,"Base.png")

【问题讨论】：

标签： dictionary plot global tmap

【解决方案1】：

tmap 的工作原理类似于 GIS，您可以添加图层。因此，如果您想显示缺少数据的国家/地区，您可以将整个世界添加为一个图层并将您的数据放在上面。类似的东西：

Base <- tm_shape(world) + tm_polygons() + 
      tm_shape(World_n, projection = "eck4") +
      tm_polygons("cibase", breaks = cut,
              palette = "RdYlGn", 
              border.col = "black", 
              border.alpha = 1,midpoint = NA ) +
  tm_legend(legend.position = c("left", "bottom")) +
  tm_facets(by=c("scenario"), ncol  =1, showNA = FALSE)+
  tm_layout(panel.labels=c("Baseline : SSP1","Baseline : SSP2","Baseline : SSP4","Baseline : SSP5"),panel.label.color = "black"
            ,panel.label.size=1.3,panel.label.fontface="bold")
Base

【讨论】：

我将使用您对 tmap 的描述为“像 GIS”，因为它将帮助我更好地向人们解释为什么我更喜欢 tmap 而不是 ArcGIS。谢谢！

【解决方案2】：

我认为如果您可以包含一些数据来使用会有所帮助吗？眼下： MAP 正在请求本地文件。所以很难看出发生了什么。我猜想当您创建 MAP2 时，您正在查看人口增长预测的中位数？也许 world_n 中的某些行在 cibase 中没有值？

【讨论】：