Downloads either a simple features (sf), SpatialPolygonDataFrame or a data_frame preprocessed using broom::tidy().

get_eurostat_geospatial(
  output_class = "sf",
  resolution = "60",
  nuts_level = "all",
  year = "2016",
  cache = TRUE,
  update_cache = FALSE,
  cache_dir = NULL,
  crs = "4326"
)

Arguments

output_class

A string. Class of object returned, either sf simple features, df (data_frame) or spdf (SpatialPolygonDataFrame)

resolution

Resolution of the geospatial data. One of "60" (1:60million), "20" (1:20million) "10" (1:10million) "03" (1:3million) or "01" (1:1million).

nuts_level

Level of NUTS classification of the geospatial data. One of "0", "1", "2", "3" or "all" (mimics the original behaviour)

year

NUTS release year. One of "2003", "2006", "2010", "2013", "2016" or "2021"

cache

a logical whether to do caching. Default is TRUE. Affects only queries from the bulk download facility.

update_cache

a logical whether to update cache. Can be set also with options(eurostat_update = TRUE)

cache_dir

a path to a cache directory. The directory have to exist. The NULL (default) uses and creates 'eurostat' directory in the temporary directory from tempdir. Directory can also be set with option eurostat_cache_dir.

crs

projection of the map: 4-digit EPSG code. One of:

  • "4326" - WGS84

  • "3035" - ETRS89 / ETRS-LAEA

  • "3857" - Pseudo-Mercator

Value

a sf, data_frame or SpatialPolygonDataFrame.

Details

The data source URL is http://ec.europa.eu/eurostat/web/gisco/geodata/reference-data/administrative-units-statistical-units.

Author

Markus Kainu <markuskainu@gmail.com>

Examples

sf <- get_eurostat_geospatial(output_class = "sf", resolution = "60", nuts_level = "all")
#> sf at resolution 1:60 read from local file
df <- get_eurostat_geospatial(output_class = "df", resolution = "20", nuts_level = "0")
#> No encoding supplied: defaulting to UTF-8.
#> Regions defined for each Polygons
#> Joining, by = "id"
#> data_frame at resolution 1: 20 cached at: /tmp/Rtmp61Gkth/eurostat/df20020164326.RData
if (FALSE) { spdf <- get_eurostat_geospatial(output_class = "spdf", resolution = "10", nuts_level = "3") }