将数据与SpatialPolygonsDataFrame对象列表进行匹配，以在传单中运行等值区域映射

Question

我正在教自己如何在R中使用映射工具，在这种情况下，我希望根据下面mydata数据框的结果在各自的区域级别设置各个国家的等值线图。

请注意，我实际上是在30个国家/地区运行此功能。 我只是为了重现目的而限制在2个国家。 另请注意， getData函数将从GADM.org网站下载。

library(sp)
library(raster)
library(tidyverse)
library(leaflet)

France_map <- getData(name = "GADM", country = "FRA", level = 1)
Germany_map <- getData(name = "GADM", country = "DEU", level = 1)

all_maps <- list(France_map, Germany_map)

mydata <- dput(structure(list(cntry = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 
  1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
  2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), .Label = c("Germany", "France"), 
  class = "factor"), region_name = structure(1:37, .Label = c("Baden-Württemberg", 
  "Bayern", "Berlin", "Brandenburg", "Bremen", "Hamburg", "Hessen", "Mecklenburg-Vorpommern", 
  "Niedersachsen", "Nordrhein-Westfalen", "Rheinland-Pfalz", "Saarland", "Sachsen", 
  "Sachsen-Anhalt", "Schleswig-Holstein", "Thüringen", "Île-de-France", "Champagne-Ardenne", 
  "Picardie", "Haute-Normandie", "Centre", "Basse-Normandie", "Bourgogne", 
  "Nord-Pas-de-Calais", "Lorraine", "Alsace", "Franche-Comté", "Pays de la Loire", 
  "Bretagne", "Poitou-Charentes", "Aquitaine", "Midi-Pyrénées", "Limousin", 
  "Rhône-Alpes", "Auvergne", "Languedoc-Roussillon", "Provence-Alpes-Côte d'Azur"), 
  class = "factor"), prop = c(0.0771812080536913, 0.0808625336927224, 0.154200131820421, 
  0.143712574850299, 0.125, 0.1875, 0.160220994475138, 0.201612903225806, 
  0.149122807017544, 0.140562248995984, 0.111111111111111, 0.130434782608696, 
  0.162454873646209, 0.141891891891892, 0.164835164835165, 0.0810810810810811, 
  0.320361006883289, 0.173944002182893, 0.267984998166555, 0.145696294975402, 
  0.307927600677265, 0.218154962499145, 0.249630808322935, 0.276741747762545, 
  0.0959323496890737, 0.144213110146746, 0.0929990840802643, 0.190088309954592, 
  0.159029556704208, 0.204186267910012, 0.26199124557682, 0.171309738332821, 
  0.125707503945437, 0.241559205433323, 0.170811734710593, 0.275876128675216, 
  0.205483749159249)), .Names = c("cntry", "region_name", "prop"), class = c("tbl_df", 
  "tbl", "data.frame"), row.names = c(NA, -37L)))

我现在要做的是将mydata的prop变量的结果显示在等值区域图上，使用类似于：

leaflet() %>%
  addPolygons()

我的绊脚石是，如何让R匹配prop变量中的结果到SPDF文件中的相应区域来运行等值区域图？

Answer 1

我结合了Germany_map和France_map这样我就可以一次绘制多边形。 然后，我为prop创造了颜色。 当你使用addPolygons() ，你指定fillColor = ~mypal(mydata$prop) 。 prop值不必保留在多边形数据中。 同样，您在addLegend()指定values = mydata$prop 。 您将下面的结果视为屏幕截图。

library(sp)
library(raster)
library(maptools)
library(leaflet)

# Combine German and French polygons
mymap <- rbind(Germany_map, France_map, makeUniqueIDs = TRUE)

### Create colors for prop
mypal <- colorNumeric(palette = "Reds", domain = mydata$prop)


leaflet() %>% 
addProviderTiles("OpenStreetMap.Mapnik") %>%
setView(lat = 49.2402, lng = 6.9969, zoom = 6) %>%
addPolygons(data = mymap,
            stroke = FALSE, smoothFactor = 0.2, fillOpacity = 0.5,
            fillColor = ~mypal(mydata$prop),
            popup = paste("Region: ", mymap$NAME_1, "<br>",
                          "Prop: ", mydata$prop, "<br>")) %>%
addLegend(position = "bottomright", pal = mypal, values = mydata$prop,
          title = "Prop",
          opacity = 0.5)