在地图周围绘制各州省份，并用德国邮政编码着色

我正在尝试创建一张地图，在其中按德国邮政编码进行着色，并且我想添加德国各州作为边界。

# for loading our data
library(raster)
library(readr)
library(readxl)
library(sf)
library(dplyr)

# for datasets
library(maps)
library(spData)

# for plotting
library(grid)
library(tmap)
library(viridis)

获取德国的形状文件（link https://opendata-esri-de.opendata.arcgis.com/datasets/5b203df4357844c8a6715d7d411a8341_0/explore?location=51.099785%2C10.454062%2C7.47）。在德国，邮政编码称为 Postleitzahlen (PLZ)。

germany <- read_sf("data/OSM_PLZ.shp")

将 PLZ 分为任意组以进行绘图。

germany <- germany %>% 
  mutate(plz_groups = case_when(
    substr(plz, 1, 1) == "1" ~ "Group A",
    substr(plz, 2, 2) == "2" ~ "Group B",    
    TRUE ~ "Group X" # rest
  ))

由PLZ进行绘图填充：

tm_shape(germany) +
  tm_fill(col = "plz_groups") 

尝试在顶部绘制德国各州（“Bundesland”）的边界：

tm_shape(germany) +
  tm_fill(col = "plz_groups") +
  tm_borders(col = "bundesland")

col2rgb(col, alpha = TRUE) 中的错误：“bundesland”中的 ungültiger Farbname

要在顶部绘制德国各州的边界，您首先必须获得德国各州的形状文件，例如从ESRI https://opendata-esri-de.opendata.arcgis.com/datasets/esri-de-content::bundesl%C3%A4ndergrenzen-2018/explore?location=50.910218%2C10.454033%2C7.00。之后您可以通过添加该形状文件tm_shape并通过绘制边界tm_borders 哪里的col参数用于设置边框线的颜色。来自文档：

对于tm_borders，它是指定边框线颜色的单一颜色值。

library(sf)
library(dplyr)
library(tmap)

germany <- read_sf("Postleitzahlengebiete_-_OSM/OSM_PLZ.shp")

german_states <- read_sf("Bundesl%C3%A4ndergrenzen_2018/Bundesländergrenzen_2018.shp")

germany <- germany %>%
  mutate(plz_groups = case_when(
    substr(plz, 1, 1) == "1" ~ "Group A",
    substr(plz, 2, 2) == "2" ~ "Group B",
    TRUE ~ "Group X" # rest
  ))

tm_shape(germany) +
  tm_fill(col = "plz_groups") +
  tm_shape(german_states) +
  tm_borders(col = "black")