Guest User

Untitled

a guest
Jan 5th, 2018
66
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 2.01 KB | None | 0 0
  1. library(sp)
  2. library(raster)
  3. library(tidyverse)
  4.  
  5. France_map <- getData(name = "GADM", country = "FRA", level = 1)
  6. Germany_map <- getData(name = "GADM", country = "DEU", level = 1)
  7.  
  8. all_maps <- list(France_map, Germany_map)
  9.  
  10. mydata <- dput(structure(list(cntry = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L,
  11. 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
  12. 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), .Label = c("Germany", "France"),
  13. class = "factor"), region_name = structure(1:37, .Label = c("Baden-Württemberg",
  14. "Bayern", "Berlin", "Brandenburg", "Bremen", "Hamburg", "Hessen", "Mecklenburg-Vorpommern",
  15. "Niedersachsen", "Nordrhein-Westfalen", "Rheinland-Pfalz", "Saarland", "Sachsen",
  16. "Sachsen-Anhalt", "Schleswig-Holstein", "Thüringen", "Île-de-France", "Champagne-Ardenne",
  17. "Picardie", "Haute-Normandie", "Centre", "Basse-Normandie", "Bourgogne",
  18. "Nord-Pas-de-Calais", "Lorraine", "Alsace", "Franche-Comté", "Pays de la Loire",
  19. "Bretagne", "Poitou-Charentes", "Aquitaine", "Midi-Pyrénées", "Limousin",
  20. "Rhône-Alpes", "Auvergne", "Languedoc-Roussillon", "Provence-Alpes-Côte d'Azur"),
  21. class = "factor"), prop = c(0.0771812080536913, 0.0808625336927224, 0.154200131820421,
  22. 0.143712574850299, 0.125, 0.1875, 0.160220994475138, 0.201612903225806,
  23. 0.149122807017544, 0.140562248995984, 0.111111111111111, 0.130434782608696,
  24. 0.162454873646209, 0.141891891891892, 0.164835164835165, 0.0810810810810811,
  25. 0.320361006883289, 0.173944002182893, 0.267984998166555, 0.145696294975402,
  26. 0.307927600677265, 0.218154962499145, 0.249630808322935, 0.276741747762545,
  27. 0.0959323496890737, 0.144213110146746, 0.0929990840802643, 0.190088309954592,
  28. 0.159029556704208, 0.204186267910012, 0.26199124557682, 0.171309738332821,
  29. 0.125707503945437, 0.241559205433323, 0.170811734710593, 0.275876128675216,
  30. 0.205483749159249)), .Names = c("cntry", "region_name", "prop"), class = c("tbl_df",
  31. "tbl", "data.frame"), row.names = c(NA, -37L)))
  32.  
  33. leaflet() %>%
  34. addPolygons()
Add Comment
Please, Sign In to add comment