带有排列国家多边形的图表

Question

我正在寻找使用以下数据集（下）的特殊图表的解决方案。该图必须包含来自数据集的国家的多边形，但它们必须（1）彼此相邻放置，而不是长、纬度坐标； （2）每个国家的多边形的大小必须对应相对大小，这是一个t$rs变量（计算如下）； (3) 每个国家的填充颜色取决于 t$value 变量的值 -- 如果为正则颜色为绿色，如果为负则为红色。

到目前为止我的代码是：

 library(ggmap)
 library(dplyr)
 library(sf) 
 library(tidyverse)
 library(gridExtra)
 library(rgdal)
 library(tmap)

 # The data set
 t <- data.frame(id    = c(136,142,172,567),
            name  = c("Italy","Norway","Finland","Singapore"),
            value = c(-0.921253632,245.6713064,4.049413648,207.5896534))

 # Min value in absolute terms
 min = min(abs(t$value))

 # Relative value w.r.t. min value 
 t$rs <- t$value / min

 # Shape files for countries
 # https://hub.arcgis.com/datasets/esri::world-countries-generalized/about

 # Unzip files from the ArcGIS archive file
 unzip("World_Countries_(Generalized).zip", exdir = ".")
 shp.file <- "World_Countries__Generalized_.shx"
 countries <- readOGR(shp.file)

 qtm(countries[countries$COUNTRY %in% t$name,])
 

我的输出图附后。距离想要的形状有点远。

Answer 1

我使用 sf 包中的 st_read() 稍微调整了您的数据导入：

library(tidyverse)
library(sf)
library(tmap)

# Unzip files from the ArcGIS archive file
unzip("World_Countries_(Generalized).zip", exdir = ".")
shp.file <- "World_Countries__Generalized_.shx"
countries <- st_read(shp.file) 
   
countries %>% 
   left_join(t, by = c("COUNTRY" = "name")) %>% 
   filter(!is.na(id)) %>% 
   st_as_sf() %>% 
   tm_shape() + 
   tm_fill("value") +
   tm_facets(by = "COUNTRY")