Raster package

library(raster)
## Loading required package: sp
fr <- getData(name = "GADM", country = "FRA", level = 1)
plot(fr)

names(fr@data)
##  [1] "PID"       "ID_0"      "ISO"       "NAME_0"    "ID_1"     
##  [6] "NAME_1"    "NL_NAME_1" "VARNAME_1" "TYPE_1"    "ENGTYPE_1"
fr$NAME_1
##  [1] "Île-de-France"              "Champagne-Ardenne"         
##  [3] "Picardie"                   "Haute-Normandie"           
##  [5] "Centre"                     "Basse-Normandie"           
##  [7] "Bourgogne"                  "Nord-Pas-de-Calais"        
##  [9] "Lorraine"                   "Alsace"                    
## [11] "Franche-Comté"              "Pays de la Loire"          
## [13] "Bretagne"                   "Poitou-Charentes"          
## [15] "Aquitaine"                  "Midi-Pyrénées"             
## [17] "Limousin"                   "Rhône-Alpes"               
## [19] "Auvergne"                   "Languedoc-Roussillon"      
## [21] "Provence-Alpes-Côte d'Azur" "Corse"
object.size(fr) / 1000000
## 4.017712 bytes
frs <- rgeos::gSimplify(fr, tol = 0.001)
frs <- SpatialPolygonsDataFrame(frs, fr@data)
object.size(frs) / 1000000
## 0.869484 bytes
# library(devtools)
# install_github("rstudio/leaflet")
library(leaflet)
leaflet() %>% 
  addTiles() %>% 
  addPolygons(data = frs, popup = fr$NAME_1)

Challenge: download data for 2 different countries and 2 different levels.

Challenge: Try to merge some data you’ve found from the internet with this geographical data and plot the result using your prefered method out of tmap, base graphics, leaflet or ggmap.