giscoR (version 0.2.4)

gisco_get: Get geospatial data from GISCO API

Description

Loads a simple feature (sf) object from GISCO API entry point or your local library.

Usage

gisco_get_coastallines(
  year = "2016",
  epsg = "4326",
  cache = TRUE,
  update_cache = FALSE,
  cache_dir = NULL,
  verbose = FALSE,
  resolution = "20"
)
gisco_get_communes(
  year = "2016",
  epsg = "4326",
  cache = TRUE,
  update_cache = FALSE,
  cache_dir = NULL,
  verbose = FALSE,
  spatialtype = "RG",
  country = NULL
)
gisco_get_countries(
  year = "2016",
  epsg = "4326",
  cache = TRUE,
  update_cache = FALSE,
  cache_dir = NULL,
  verbose = FALSE,
  resolution = "20",
  spatialtype = "RG",
  country = NULL,
  region = NULL
)
gisco_get_lau(
  year = "2016",
  epsg = "4326",
  cache = TRUE,
  update_cache = FALSE,
  cache_dir = NULL,
  verbose = FALSE,
  country = NULL,
  gisco_id = NULL
)
gisco_get_nuts(
  year = "2016",
  epsg = "4326",
  cache = TRUE,
  update_cache = FALSE,
  cache_dir = NULL,
  verbose = FALSE,
  resolution = "20",
  spatialtype = "RG",
  country = NULL,
  nuts_id = NULL,
  nuts_level = "all"
)
gisco_get_urban_audit(
  year = "2020",
  epsg = "4326",
  cache = TRUE,
  update_cache = FALSE,
  cache_dir = NULL,
  verbose = FALSE,
  spatialtype = "RG",
  country = NULL,
  level = NULL
)

Arguments

year

Release year. See Release years available on gisco_get.

epsg

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

• "4258": ETRS89

• "4326": WGS84

• "3035": ETRS89 / ETRS-LAEA

• "3857": Pseudo-Mercator

cache

A logical whether to do caching. Default is TRUE.

update_cache

A logical whether to update cache. Default is FALSE. When set to TRUE it would force a fresh download of the source .geojson file.

cache_dir

A path to a cache directory. See Details on gisco_get.

verbose

Display information. Useful for debugging, default is FALSE.

resolution

Resolution of the geospatial data. One of

• "60": 1:60million

• "20": 1:20million

• "10": 1:10million

• "03": 1:3million

• "01": 1:1million

spatialtype

Type of geometry to be returned:

• "RG": Regions - MULTIPOLYGON/POLYGON object.

• "LB": Labels - POINT object.

• "BN": Boundaries - LINESTRING object.

• "COASTL": coastlines - LINESTRING object.

• "INLAND": inland boundaries - LINESTRING object.

country

Optional. A character vector of country codes. See Details on gisco_get.

region

Optional. A character vector of UN M49 region codes or European Union membership. Possible values are "Africa", "Americas", "Asia", "Europe", "Oceania" or "EU" for countries belonging to the European Union (as per 2021). See About world regions and gisco_countrycode

gisco_id

Optional. A character vector of GISCO_ID LAU values.

nuts_id

Optional. A character vector of NUTS IDs.

nuts_level

NUTS level. One of "0" (Country-level), "1", "2" or "3". See https://ec.europa.eu/eurostat/web/nuts/background.

level

Level of Urban Audit. Possible values are "CITIES", "FUA", "GREATER_CITIES" or NULL, that would download the full dataset.

Value

A sf object specified by spatialtype.

About world regions

Regions are defined as per the geographic regions defined by the UN (see https://unstats.un.org/unsd/methodology/m49/. Under this scheme Cyprus is assigned to Asia. You may use region = "EU" to get the EU members (reference date: 2021).

Release years available

• gisco_get_coastallines: one of "2006", "2010", "2013" or "2016".

• gisco_get_communes: one of "2001", "2004", "2006", "2008", "2010", "2013" or "2016".

• gisco_get_countries: one of "2001", "2006", "2010", "2013", "2016" or "2020".

• gisco_get_lau: one of "2016", "2017", "2018" or "2019"

• gisco_get_nuts: one of "2003", "2006", "2010", "2013", "2016" or "2021".

• gisco_get_urban_audit: one of "2001", "2004", "2014", "2018" or "2020".

Details

country only available on specific datasets. Some spatialtype options (as "BN", "COASTL", "INLAND") may not include country-level identifiers.

country could be either a vector of country names, a vector of ISO3 country codes or a vector of Eurostat country codes. Mixed types (as c("Turkey","US","FRA")) would not work.

Sometimes cached files may be corrupt. On that case, try re-downloading the data setting update_cache = TRUE.

Set cache_dir = "path/to/dir" or options(gisco_cache_dir = "path/to/dir). If you experience any problem on download, try to download the corresponding .geojson file by any other method and save it on your cache_dir.

For a complete list of files available check gisco_db.

See Also

gisco_db, gisco_attributions(), gisco_coastallines

gisco_countrycode(), gisco_countries

gisco_nuts

Examples

# NOT RUN {
library(sf)

##################################
# Example - gisco_get_coastallines
##################################

coastlines <- gisco_get_coastallines()
plot(st_geometry(coastlines), col = "palegreen", border = "lightblue3")
title(
  main = "Coastal Lines",
  sub = gisco_attributions(),
  line = 1
)

###############################
# Example - gisco_get_countries
###############################

sf_world <- gisco_get_countries()
plot(st_geometry(sf_world), col = "seagreen2")
title(sub = gisco_attributions(), line = 1)


sf_africa <- gisco_get_countries(region = "Africa")
plot(st_geometry(sf_africa),
  col = c("springgreen4", "darkgoldenrod1", "red2")
)
title(sub = gisco_attributions(), line = 1)

sf_benelux <-
  gisco_get_countries(country = c("Belgium", "Netherlands", "Luxembourg"))
plot(st_geometry(sf_benelux),
  col = c("grey10", "orange", "deepskyblue2")
)
title(sub = gisco_attributions(), line = 1)

##########################
# Example - gisco_get_nuts
##########################

nuts1 <- gisco_get_nuts(
  resolution = "20",
  year = "2016",
  epsg = "4326",
  nuts_level = "1",
  country = "ITA"
)
nuts2 <- gisco_get_nuts(
  resolution = "20",
  year = "2016",
  epsg = "4326",
  nuts_level = "2",
  country = "ITA"
)
nuts3 <- gisco_get_nuts(
  resolution = "20",
  year = "2016",
  epsg = "4326",
  nuts_level = "3",
  country = "ITA"
)

plot(st_geometry(nuts3),
  border = "grey60",
  lty = 3
)

plot(st_geometry(nuts2),
  lwd = 2,
  border = "red2",
  add = TRUE
)

plot(st_geometry(nuts1),
  lwd = 3,
  border = "springgreen4",
  add = TRUE
)

box()
title(
  main = "NUTS Levels on Italy",
  sub = gisco_attributions(),
  cex.sub = 0.7,
  line = 1
)
legend(
  "topright",
  legend = c("NUTS 1", "NUTS 2", "NUTS 3"),
  col = c("springgreen4", "red2", "grey60"),
  lty = c(1, 1, 3),
  lwd = c(3, 2, 1),
  bty = "n",
  y.intersp = 2
)
# }