Description Usage Arguments Details Author(s) Examples
This function gets the geodata from ECDC GIS server into R SpatialpolygonsDataFrame or SpatialpointsDataFrame.
1 2 | get_GEO_data(layer, STAT_LEVL = c(0), FIELDS = "GEO_ID", isValid = 1,
eu_surround_continents = FALSE)
|
layer |
Select the layer for Hybrid Layer service: 0 for Points, 1 for Polygons. |
STAT_LEVL |
Which data level to choose; see details. Defaults to 0 for EU/EEA. Give a single value or vector. |
FIELDS |
Which data fields to retrieve from the hybrid layer; see details. Defaults to "GEO_ID". Give a single value or character vector of relevant field names. |
isValid |
isValid=1 removes the obsolete NUTS regions. |
eu_surround_continents |
TRUE selects the polygons for Europe, Asia and Africa (if adequate STAT_LEVL selected), to allow mapping European mainland with the surrounding background of Asia and Northern Africa. Requires clipping as e.g. per the SurvMapper. |
STAT_LEVL; 0 for EU/EEA and candidate country level, 10 for the rest of world country level. 1,2,3 = NUTS levels; 11,12 = GAUL levels; 21,22,23,24 = GADM levels. FIELDS; for additional fields to retrieve from the hybrid lay, see: https://gis.ecdc.europa.eu/public/rest/services/UtilitiesServices/EMMA_GoMap_GC_HL_V1_3_WebMerc/MapServer/
Tommi Karki, with thanks to Silviu Ionescu.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 | # get EU/EEA country polygons
plgs <- get_GEO_data(layer = 1)
# get EU/EEA country points
points <- get_GEO_data(layer = 0)
# use the geodata
plot(points, col = "red")
plot(plgs, add = TRUE)
# get all the world country polygons and some additional fields
plgs <- get_GEO_data(layer = 1, FIELDS = c("GEO_ID", "GEO_NAME", "CONTINENT"), STAT_LEVL=c(0,10))
#' # get all Europe, Asia and Africa polygons and some additional fields
plgs <- get_GEO_data(layer = 1, FIELDS = c("GEO_ID", "GEO_NAME", "CONTINENT"), STAT_LEVL=c(0,10),
eu_surround_continents = TRUE)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.