| Title: | Administrative Boundaries of Spain |
|---|---|
| Description: | Administrative Boundaries of Spain at several levels (Autonomous Communities, Provinces, Municipalities) based on the 'GISCO' 'Eurostat' database <https://ec.europa.eu/eurostat/web/gisco> and 'CartoBase SIANE' from 'Instituto Geografico Nacional' <https://www.ign.es/>. It also provides a 'leaflet' plugin and the ability of downloading and processing static tiles. |
| Authors: | Diego Hernangómez [aut, cre, cph] (<https://orcid.org/0000-0001-8457-4658>, rOpenSpain) |
| Maintainer: | Diego Hernangómez <[email protected]> |
| License: | GPL-3 |
| Version: | 0.8.0 |
| Built: | 2023-09-10 20:09:00 UTC |
| Source: | https://github.com/rOpenSpain/mapSpain |
Include tiles of public Spanish organisms to a
leaflet::leaflet() map.
addProviderEspTiles(
map,
provider,
layerId = NULL,
group = NULL,
options = providerEspTileOptions()
)
providerEspTileOptions(...)map |
A map widget created from |
provider |
Name of the provider, see esp_tiles_providers for values available. |
layerId |
the layer id |
group |
The name of the group the newly created layers should belong to
Human-friendly group names are permitted–they need not be short,
identifier-style names. Any number of layers and even different types of
layers (e.g. markers and polygons) can share the same group name. See
|
options |
a list of extra options for tile layers, popups, paths (circles, rectangles, polygons, ...), or other map elements |
... |
Arguments passed on to |
providerEspTileOptions() is a wrapper of
leaflet::providerTileOptions().
A modified leaflet::leaflet() map object.
https://dieghernan.github.io/leaflet-providersESP/ leaflet plugin, v1.3.3.
leaflet::leaflet(), leaflet::addTiles()
leaflet::providerTileOptions(), leaflet::tileOptions()
Other imagery utilities:
esp_getTiles(),
esp_make_provider(),
esp_tiles_providers
library(leaflet)
PuertadelSol <-
leaflet() %>%
setView(
lat = 40.4166,
lng = -3.7038400,
zoom = 18
) %>%
addProviderEspTiles(provider = "IGNBase.Gris") %>%
addProviderEspTiles(provider = "RedTransporte.Carreteras")
PuertadelSolCheck if R has access to resources at https://github.com/rOpenSpain/mapSpain/tree/sianedata.
esp_check_access()a logical.
esp_check_access()Use this function with caution. This function would clear your cached data and configuration, specifically:
Deletes the mapSpain config directory
(rappdirs::user_config_dir("mapSpain", "R")).
Deletes the cache_dir directory.
Deletes the values on stored on Sys.getenv("MAPSPAIN_CACHE_DIR") and
options(mapSpain_cache_dir).
esp_clear_cache(config = FALSE, cached_data = TRUE, verbose = FALSE)config |
if |
cached_data |
If this is set to |
verbose |
Logical, displays information. Useful for debugging,
default is |
This is an overkill function that is intended to reset your status as it you would never have installed and/or used mapSpain.
Invisible. This function is called for its side effects.
Other cache utilities:
esp_set_cache_dir()
# Don't run this! It would modify your current state
## Not run:
esp_clear_cache(verbose = TRUE)
## End(Not run)
Sys.getenv("MAPSPAIN_CACHE_DIR")A data frame used internally for translating codes and names of the different subdivisions of Spain. The data frame provides the hierarchy of the subdivisions including NUTS1 level, Autonomous Communities (equivalent to NUTS2), Provinces and NUTS3 level. See Note.
A data frame with 59 rows codes and columns:
NUTS 1 code
NUTS 1 name
NUTS 1 alternative name
NUTS1 1 short (common) name (Spanish)
INE code of the autonomous community
ISO2 code of the autonomous community
NUTS 2 Code
INE name of the autonomous community
ISO2 name of the autonomous community (Spanish)
ISO2 name of the autonomous community (Catalan)
ISO2 name of the autonomous community (Galician)
ISO2 name of the autonomous community (Basque)
NUTS 2 name
CLDR name of the autonomous community (English)
CLDR name of the autonomous community (Spanish)
CLDR name of the autonomous community (Catalan)
CLDR name of the autonomous community (Galician)
CLDR name of the autonomous community (Basque)
Short (common) name of the autonomous community (English)
Short (common) name of the autonomous community (Spanish)
Short (common) name of the autonomous community (Catalan)
Short (common) name of the autonomous community (Galician)
Short (common) name of the autonomous community (Basque)
INE code of the province
ISO2 code of the province
NUTS code of the province
INE name of the province
ISO2 name of the province (Spanish)
ISO2 name of the province (Catalan)
ISO2 name of the province (Galician)
ISO2 name of the province (Basque)
CLDR name of the province (English)
CLDR name of the province (Spanish)
CLDR name of the province (Catalan)
CLDR name of the province (Galician)
CLDR name of the province (Basque)
Short (common) name of the province (English)
Short (common) name of the province (Spanish)
Short (common) name of the province (Catalan)
Short (common) name of the province (Galician)
Short (common) name of the province (Basque)
NUTS 3 code
NUTS 3 name
NUTS 3 short (common) name
Although NUTS2 matches the first subdivision level of Spain (CCAA - Autonomous Communities), it should be noted that NUTS3 does not match the second subdivision level of Spain (Provinces). NUTS3 provides a dedicated code for major islands whereas the Provinces doesn't.
Ceuta and Melilla has an specific status (Autonomous Cities) but are considered as communities with a single province (as Madrid, Asturias or Murcia) on this dataset.
INE: Instituto Nacional de Estadistica: https://www.ine.es/
Eurostat (NUTS): https://ec.europa.eu/eurostat/web/nuts/background
CLDR: https://unicode-org.github.io/cldr-staging/charts/38/index.html
Other datasets:
esp_munic.sf,
esp_nuts.sf,
esp_tiles_providers,
pobmun19
Other political:
esp_get_can_box(),
esp_get_capimun(),
esp_get_ccaa(),
esp_get_comarca(),
esp_get_country(),
esp_get_gridmap,
esp_get_munic(),
esp_get_nuts(),
esp_get_prov(),
esp_get_simpl_prov()
Other dictionary:
esp_dict_region_code()
data("esp_codelist")Converts long subdivision names into different coding schemes and languages.
esp_dict_region_code(sourcevar, origin = "text", destination = "text")
esp_dict_translate(sourcevar, lang = "en", all = FALSE)sourcevar |
Vector which contains the subdivision names to be converted. |
origin, destination
|
One of |
lang |
Language of translation. Available languages are:
|
all |
Logical. Should the function return all names or not?
On |
If no match is found for any value, the function displays a warning and
returns NA for those values.
Note that mixing names of different administrative levels (e.g. "Catalonia"
and "Barcelona") may return empty values, depending on the destination
values.
esp_dict_region_code() returns a vector of characters.
esp_dict_translate() returns a character vector or a named
list with each of the possible names of each sourcevar on the required
language lang.
Other dictionary:
esp_codelist
Other dictionary:
esp_codelist
vals <- c("Errioxa", "Coruna", "Gerona", "Madrid")
esp_dict_region_code(vals)
esp_dict_region_code(vals, destination = "nuts")
esp_dict_region_code(vals, destination = "cpro")
esp_dict_region_code(vals, destination = "iso2")
# From ISO2 to another codes
iso2vals <- c("ES-M", "ES-S", "ES-SG")
esp_dict_region_code(iso2vals, origin = "iso2")
esp_dict_region_code(iso2vals,
origin = "iso2",
destination = "nuts"
)
esp_dict_region_code(iso2vals,
origin = "iso2",
destination = "cpro"
)
# Mixing levels
valsmix <- c("Centro", "Andalucia", "Seville", "Menorca")
esp_dict_region_code(valsmix, destination = "nuts")
## Not run:
# Warning
esp_dict_region_code(valsmix, destination = "codauto")
esp_dict_region_code(valsmix, destination = "iso2")
## End(Not run)
vals <- c("La Rioja", "Sevilla", "Madrid", "Jaen", "Orense", "Baleares")
esp_dict_translate(vals)
esp_dict_translate(vals, lang = "es")
esp_dict_translate(vals, lang = "ca")
esp_dict_translate(vals, lang = "eu")
esp_dict_translate(vals, lang = "ga")
esp_dict_translate(vals, lang = "ga", all = TRUE)When plotting Spain, it is usual to represent the Canary Islands as an inset
(see moveCAN on esp_get_nuts()). These functions provides complementary
lines and polygons to be used when the Canary Islands are displayed
as an inset.
esp_get_can_box() is used to draw lines around the displaced Canary
Islands.
esp_get_can_provinces() is used to draw a separator line between the two
provinces of the Canary Islands.
esp_get_can_box(style = "right", moveCAN = TRUE, epsg = "4258")
esp_get_can_provinces(moveCAN = TRUE, epsg = "4258")style |
Style of line around Canary Islands. Four options available:
|
moveCAN |
A logical |
epsg |
projection of the map: 4-digit EPSG code. One of:
|
A sf polygon or line depending of style parameter.
esp_get_can_provinces returns a LINESTRING object.
While moveCAN is useful for visualization, it would alter the actual
geographic position of the Canary Islands. When using the output for
spatial analysis or using tiles (e.g. with esp_getTiles() or
addProviderEspTiles()) this option should be set to FALSE in order to
get the actual coordinates, instead of the modified ones.
esp_get_can_provinces extracted from CartoBase ANE,
se89_mult_admin_provcan_l.shp file.
Other political:
esp_codelist,
esp_get_capimun(),
esp_get_ccaa(),
esp_get_comarca(),
esp_get_country(),
esp_get_gridmap,
esp_get_munic(),
esp_get_nuts(),
esp_get_prov(),
esp_get_simpl_prov()
Provs <- esp_get_prov()
Box <- esp_get_can_box()
Line <- esp_get_can_provinces()
# Plot
library(ggplot2)
ggplot(Provs) +
geom_sf() +
geom_sf(data = Box) +
geom_sf(data = Line) +
theme_linedraw()
# Displacing Canary
# By same factor
displace <- c(15, 0)
Provs_D <- esp_get_prov(moveCAN = displace)
Box_D <- esp_get_can_box(style = "left", moveCAN = displace)
Line_D <- esp_get_can_provinces(moveCAN = displace)
ggplot(Provs_D) +
geom_sf() +
geom_sf(data = Box_D) +
geom_sf(data = Line_D) +
theme_linedraw()
# Example with poly option
# Get countries with giscoR
library(giscoR)
# Low resolution map
res <- "20"
Countries <-
gisco_get_countries(
res = res,
epsg = "4326",
country = c("France", "Portugal", "Andorra", "Morocco", "Argelia")
)
CANbox <-
esp_get_can_box(
style = "poly",
epsg = "4326",
moveCAN = c(12.5, 0)
)
CCAA <- esp_get_ccaa(
res = res,
epsg = "4326",
moveCAN = c(12.5, 0) # Same displacement factor)
)
# Plot
ggplot(Countries) +
geom_sf(fill = "#DFDFDF") +
geom_sf(data = CANbox, fill = "#C7E7FB", linewidth = 1) +
geom_sf(data = CCAA, fill = "#FDFBEA") +
coord_sf(
xlim = c(-10, 4.3),
ylim = c(34.6, 44)
) +
theme(
panel.background = element_rect(fill = "#C7E7FB"),
panel.grid = element_blank()
)Get a sf point with the location of the political powers for each municipality (possibly the center of the municipality).
Note that this differs of the centroid of the boundaries of the
municipality, returned by esp_get_munic().
esp_get_capimun(
year = Sys.Date(),
epsg = "4258",
cache = TRUE,
update_cache = FALSE,
cache_dir = NULL,
verbose = FALSE,
region = NULL,
munic = NULL,
moveCAN = TRUE,
rawcols = FALSE
)year |
Release year. See Details for years available. |
epsg |
projection of the map: 4-digit EPSG code. One of:
|
cache |
A logical whether to do caching. Default is |
update_cache |
A logical whether to update cache. Default is |
cache_dir |
A path to a cache directory. See About caching. |
verbose |
Logical, displays information. Useful for debugging,
default is |
region |
A vector of names and/or codes for provinces
or |
munic |
A name or |
moveCAN |
A logical |
rawcols |
Logical. Setting this to |
year could be passed as a single year (YYYY format, as end of year) or
as a specific date (YYYY-MM-DD format). Historical information starts as
of 2005.
When using region you can use and mix names and NUTS codes (levels 1,
2 or 3), ISO codes (corresponding to level 2 or 3) or cpro. See
esp_codelist
When calling a superior level (Province, Autonomous Community or NUTS1) , all the municipalities of that level would be added.
A sf point object.
You can set your cache_dir with esp_set_cache_dir().
Sometimes cached files may be corrupt. On that case, try re-downloading
the data setting update_cache = TRUE.
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. Use the option verbose = TRUE for debugging the API query.
While moveCAN is useful for visualization, it would alter the actual
geographic position of the Canary Islands. When using the output for
spatial analysis or using tiles (e.g. with esp_getTiles() or
addProviderEspTiles()) this option should be set to FALSE in order to
get the actual coordinates, instead of the modified ones.
IGN data via a custom CDN (see https://github.com/rOpenSpain/mapSpain/tree/sianedata).
Other political:
esp_codelist,
esp_get_can_box(),
esp_get_ccaa(),
esp_get_comarca(),
esp_get_country(),
esp_get_gridmap,
esp_get_munic(),
esp_get_nuts(),
esp_get_prov(),
esp_get_simpl_prov()
Other municipalities:
esp_get_munic(),
esp_munic.sf
## Not run:
# This code compares centroids of municipalities against esp_get_capimun
# It also download tiles, make sure you are online
library(sf)
# Get shape
area <- esp_get_munic_siane(munic = "Valladolid", epsg = 3857)
# Area in km2
print(paste0(round(as.double(sf::st_area(area)) / 1000000, 2), " km2"))
# Extract centroid
centroid <- sf::st_centroid(area)
centroid$type <- "Centroid"
# Compare with capimun
capimun <- esp_get_capimun(munic = "Valladolid", epsg = 3857)
capimun$type <- "Capimun"
# Get a tile to check
tile <- esp_getTiles(area, "IGNBase.Todo", zoommin = 2)
# Join both point geometries
points <- rbind(
centroid[, "type"],
capimun[, "type"]
)
# Check on plot
library(ggplot2)
library(tidyterra)
ggplot(points) +
geom_spatraster_rgb(data = tile) +
geom_sf(data = area, fill = NA, color = "blue") +
geom_sf(data = points, aes(fill = type), size = 5, shape = 21) +
scale_fill_manual(values = c("green", "red")) +
theme_void() +
labs(title = "Centroid vs. capimun")
## End(Not run)Returns Autonomous Communities of Spain as polygons and points at a specified scale.
esp_get_ccaa() uses GISCO (Eurostat) as source. Please use
giscoR::gisco_attributions()
esp_get_ccaa_siane() uses CartoBase ANE as source, provided by
Instituto Geografico Nacional (IGN), http://www.ign.es/web/ign/portal.
Years available are 2005 up to today.
esp_get_ccaa(ccaa = NULL, moveCAN = TRUE, ...)
esp_get_ccaa_siane(
ccaa = NULL,
year = Sys.Date(),
epsg = "4258",
cache = TRUE,
update_cache = FALSE,
cache_dir = NULL,
verbose = FALSE,
resolution = "3",
moveCAN = TRUE,
rawcols = FALSE
)ccaa |
A vector of names and/or codes for autonomous communities
or |
moveCAN |
A logical |
... |
Arguments passed on to
|
year |
Release year. See |
epsg |
projection of the map: 4-digit EPSG code. One of:
|
cache |
A logical whether to do caching. Default is |
update_cache |
A logical whether to update cache. Default is |
cache_dir |
A path to a cache directory. See About caching. |
verbose |
Logical, displays information. Useful for debugging,
default is |
resolution |
Resolution of the polygon. Values available are
|
rawcols |
Logical. Setting this to |
When using ccaa you can use and mix names and NUTS codes (levels 1 or 2),
ISO codes (corresponding to level 2) or codauto (see esp_codelist).
Ceuta and Melilla are considered as Autonomous Communities on this function.
When calling a NUTS1 level, all the Autonomous Communities of that level would be added.
On esp_get_ccaa_siane(), year could be passed as a single year (YYYY
format, as end of year) or as a specific date (YYYY-MM-DD format).
Historical information starts as of 2005.
A sf object specified by spatialtype.
You can set your cache_dir with esp_set_cache_dir().
Sometimes cached files may be corrupt. On that case, try re-downloading
the data setting update_cache = TRUE.
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. Use the option verbose = TRUE for debugging the API query.
While moveCAN is useful for visualization, it would alter the actual
geographic position of the Canary Islands. When using the output for
spatial analysis or using tiles (e.g. with esp_getTiles() or
addProviderEspTiles()) this option should be set to FALSE in order to
get the actual coordinates, instead of the modified ones.
IGN data via a custom CDN (see https://github.com/rOpenSpain/mapSpain/tree/sianedata).
Other political:
esp_codelist,
esp_get_can_box(),
esp_get_capimun(),
esp_get_comarca(),
esp_get_country(),
esp_get_gridmap,
esp_get_munic(),
esp_get_nuts(),
esp_get_prov(),
esp_get_simpl_prov()
ccaa <- esp_get_ccaa()
library(ggplot2)
ggplot(ccaa) +
geom_sf()
# Random CCAA
Random <- esp_get_ccaa(ccaa = c(
"Euskadi",
"Catalunya",
"ES-EX",
"Canarias",
"ES52",
"01"
))
ggplot(Random) +
geom_sf(aes(fill = codauto), show.legend = FALSE) +
geom_sf_label(aes(label = codauto), alpha = 0.3)
# All CCAA of a Zone plus an addition
Mix <-
esp_get_ccaa(ccaa = c("La Rioja", "Noroeste"))
ggplot(Mix) +
geom_sf()
# Combine with giscoR to get countries
library(giscoR)
library(sf)
res <- 20 # Set same resoluion
europe <- gisco_get_countries(resolution = res)
ccaa <- esp_get_ccaa(moveCAN = FALSE, resolution = res)
# Transform to same CRS
europe <- st_transform(europe, 3035)
ccaa <- st_transform(ccaa, 3035)
ggplot(europe) +
geom_sf(fill = "#DFDFDF", color = "#656565") +
geom_sf(data = ccaa, fill = "#FDFBEA", color = "#656565") +
coord_sf(
xlim = c(23, 74) * 10e4,
ylim = c(14, 55) * 10e4
) +
theme(panel.background = element_rect(fill = "#C7E7FB"))Returns 'comarcas' of Spain as polygons,as provided by the INE (Instituto Nacional de Estadistica).
esp_get_comarca(
region = NULL,
comarca = NULL,
moveCAN = TRUE,
epsg = "4258",
update_cache = FALSE,
cache_dir = NULL,
verbose = FALSE
)region |
A vector of names and/or codes for provinces or |
comarca |
A name or |
moveCAN |
A logical |
epsg |
projection of the map: 4-digit EPSG code. One of:
|
update_cache |
A logical whether to update cache. Default is |
cache_dir |
A path to a cache directory. See About caching. |
verbose |
Logical, displays information. Useful for debugging,
default is |
'Comarcas' (English equivalent: district, county, area or zone) does not always have a formal legal status. They correspond mainly to natural areas (valleys, river basins etc.) or even to historical regions or ancient kingdoms.
When using region you can use and mix names and NUTS codes
(levels 1, 2 or 3), ISO codes (corresponding to level 2 or 3) or
"cpro" (see esp_codelist).
When calling a superior level (Province, Autonomous Community or NUTS1) , all the comarcas of that level would be added.
The use of the information contained on the INE website may be carried out by users or re-use agents, at their own risk, and they will be the sole liable parties in the case of having to answer to third parties due to damages arising from such use.
A sf polygon object.
You can set your cache_dir with esp_set_cache_dir().
Sometimes cached files may be corrupt. On that case, try re-downloading
the data setting update_cache = TRUE.
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. Use the option verbose = TRUE for debugging the API query.
While moveCAN is useful for visualization, it would alter the actual
geographic position of the Canary Islands. When using the output for
spatial analysis or using tiles (e.g. with esp_getTiles() or
addProviderEspTiles()) this option should be set to FALSE in order to
get the actual coordinates, instead of the modified ones.
INE: PC_Axis files.
Other political:
esp_codelist,
esp_get_can_box(),
esp_get_capimun(),
esp_get_ccaa(),
esp_get_country(),
esp_get_gridmap,
esp_get_munic(),
esp_get_nuts(),
esp_get_prov(),
esp_get_simpl_prov()
comarcas <- esp_get_comarca(moveCAN = FALSE)
library(ggplot2)
ggplot(comarcas) +
geom_sf()
# Comarcas of Castille and Leon
comarcas_cyl <- esp_get_comarca("Castilla y Leon")
ggplot(comarcas_cyl) +
geom_sf(aes(fill = ine.prov.name)) +
labs(fill = "Province")
# Comarcas with Mountains or Alt(o,a) in the name
comarcas_alto <- esp_get_comarca(
comarca = "Montaña|Monte|Sierra|Alt",
epsg = 3857
)
ggplot(comarcas_alto) +
geom_sf(aes(fill = ine.ccaa.name)) +
geom_sf_text(aes(label = name), check_overlap = TRUE) +
labs(fill = "CCAA")Returns the boundaries of Spain as a single sf polygon at a specified scale.
esp_get_country(moveCAN = TRUE, ...)moveCAN |
A logical |
... |
Arguments passed on to
|
A sf polygon object.
You can set your cache_dir with esp_set_cache_dir().
Sometimes cached files may be corrupt. On that case, try re-downloading
the data setting update_cache = TRUE.
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. Use the option verbose = TRUE for debugging the API query.
While moveCAN is useful for visualization, it would alter the actual
geographic position of the Canary Islands. When using the output for
spatial analysis or using tiles (e.g. with esp_getTiles() or
addProviderEspTiles()) this option should be set to FALSE in order to
get the actual coordinates, instead of the modified ones.
Other political:
esp_codelist,
esp_get_can_box(),
esp_get_capimun(),
esp_get_ccaa(),
esp_get_comarca(),
esp_get_gridmap,
esp_get_munic(),
esp_get_nuts(),
esp_get_prov(),
esp_get_simpl_prov()
OriginalCan <- esp_get_country(moveCAN = FALSE)
# One row only
nrow(OriginalCan)
library(ggplot2)
ggplot(OriginalCan) +
geom_sf(fill = "grey70")
# Less resolution
MovedCan <- esp_get_country(moveCAN = TRUE, resolution = "20")
library(ggplot2)
ggplot(MovedCan) +
geom_sf(fill = "grey70")Loads a sf polygon with the geographic grids of Spain as provided on the Banco de Datos de la Naturaleza (Nature Data Bank), by the Ministry of Environment (MITECO):
esp_get_grid_BDN() extracts country-wide grids with resolutions
5x5 or 10x10 kms.
esp_get_grid_BDN_ccaa() extracts grids by Autonomous Community with
resolution 1x1 km.
esp_get_grid_BDN(
resolution = 10,
type = "main",
update_cache = FALSE,
cache_dir = NULL,
verbose = FALSE
)
esp_get_grid_BDN_ccaa(
ccaa,
update_cache = FALSE,
cache_dir = NULL,
verbose = FALSE
)resolution |
Resolution of the grid in kms. Could be |
type |
The scope of the grid. It could be mainland Spain ( |
update_cache |
A logical whether to update cache. Default is |
cache_dir |
A path to a cache directory. See About caching. |
verbose |
Logical, displays information. Useful for debugging,
default is |
ccaa |
A vector of names and/or codes for autonomous communities.
See Details on |
A sf polygon
You can set your cache_dir with esp_set_cache_dir().
Sometimes cached files may be corrupt. On that case, try re-downloading
the data setting update_cache = TRUE.
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. Use the option verbose = TRUE for debugging the API query.
BDN data via a custom CDN (see https://github.com/rOpenSpain/mapSpain/tree/sianedata/MTN).
See original metadata and source on https://www.miteco.gob.es/es/biodiversidad/servicios/banco-datos-naturaleza/informacion-disponible/bdn-cart-aux-descargas-ccaa.aspx
Other grids:
esp_get_grid_EEA(),
esp_get_grid_ESDAC(),
esp_get_grid_MTN()
grid <- esp_get_grid_BDN(resolution = "10", type = "main")
library(ggplot2)
ggplot(grid) +
geom_sf() +
theme_light() +
labs(title = "BDN Grid for Spain")Loads a sf polygon with the geographic grids of Spain as provided by the European Environment Agency (EEA).
esp_get_grid_EEA(
resolution = 100,
type = "main",
update_cache = FALSE,
cache_dir = NULL,
verbose = FALSE
)resolution |
Resolution of the grid in kms. Could be |
type |
The scope of the grid. It could be mainland Spain ( |
update_cache |
A logical whether to update cache. Default is |
cache_dir |
A path to a cache directory. See About caching. |
verbose |
Logical, displays information. Useful for debugging,
default is |
A sf polygon
You can set your cache_dir with esp_set_cache_dir().
Sometimes cached files may be corrupt. On that case, try re-downloading
the data setting update_cache = TRUE.
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. Use the option verbose = TRUE for debugging the API query.
Other grids:
esp_get_grid_BDN(),
esp_get_grid_ESDAC(),
esp_get_grid_MTN()
## Not run:
grid <- esp_get_grid_EEA(type = "main", resolution = 100)
grid_can <- esp_get_grid_EEA(type = "canary", resolution = 100)
esp <- esp_get_country(moveCAN = FALSE)
library(ggplot2)
ggplot(grid) +
geom_sf() +
geom_sf(data = grid_can) +
geom_sf(data = esp, fill = NA) +
theme_light() +
labs(title = "EEA Grid for Spain")
## End(Not run)Loads a sf polygon with the geographic grids of Spain as provided by the European Soil Data Centre (ESDAC).
esp_get_grid_ESDAC(
resolution = 10,
update_cache = FALSE,
cache_dir = NULL,
verbose = FALSE
)resolution |
Resolution of the grid in kms. Could be |
update_cache |
A logical whether to update cache. Default is |
cache_dir |
A path to a cache directory. See About caching. |
verbose |
Logical, displays information. Useful for debugging,
default is |
A sf polygon
You can set your cache_dir with esp_set_cache_dir().
Sometimes cached files may be corrupt. On that case, try re-downloading
the data setting update_cache = TRUE.
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. Use the option verbose = TRUE for debugging the API query.
Panagos P., Van Liedekerke M., Jones A., Montanarella L., "European Soil Data Centre: Response to European policy support and public data requirements"; (2012) Land Use Policy, 29 (2), pp. 329-338. doi:10.1016/j.landusepol.2011.07.003
European Soil Data Centre (ESDAC), esdac.jrc.ec.europa.eu, European Commission, Joint Research Centre.
Other grids:
esp_get_grid_BDN(),
esp_get_grid_EEA(),
esp_get_grid_MTN()
## Not run:
grid <- esp_get_grid_ESDAC()
esp <- esp_get_country(moveCAN = FALSE)
library(ggplot2)
ggplot(grid) +
geom_sf() +
geom_sf(data = esp, color = "grey50", fill = NA) +
theme_light() +
labs(title = "ESDAC Grid for Spain")
## End(Not run)Loads a sf polygon with the geographic grids of Spain.
esp_get_grid_MTN(
grid = "MTN25_ETRS89_Peninsula_Baleares_Canarias",
update_cache = FALSE,
cache_dir = NULL,
verbose = FALSE
)grid |
Name of the grid to be loaded. See Details. |
update_cache |
A logical whether to update cache. Default is |
cache_dir |
A path to a cache directory. See About caching. |
verbose |
Logical, displays information. Useful for debugging,
default is |
Metadata available on https://github.com/rOpenSpain/mapSpain/tree/sianedata/MTN.
Possible values of grid are:
| grid_name |
| MTN25_ED50_Peninsula_Baleares |
| MTN25_ETRS89_ceuta_melilla_alboran |
| MTN25_ETRS89_Peninsula_Baleares_Canarias |
| MTN25_RegCan95_Canarias |
| MTN50_ED50_Peninsula_Baleares |
| MTN50_ETRS89_Peninsula_Baleares_Canarias |
| MTN50_RegCan95_Canarias |
A description of the MTN (Mapa Topografico Nacional) grids available:
MTN25_ED50_Peninsula_Baleares
MTN25 grid corresponding to the Peninsula and Balearic Islands, in ED50 and geographical coordinates (longitude, latitude) This is the real MTN25 grid, that is, the one that divides the current printed series of the map, taking into account special sheets and irregularities.
MTN50_ED50_Peninsula_Baleares
MTN50 grid corresponding to the Peninsula and Balearic Islands, in ED50 and geographical coordinates (longitude, latitude) This is the real MTN50 grid, that is, the one that divides the current printed series of the map, taking into account special sheets and irregularities.
MTN25_ETRS89_ceuta_melilla_alboran
MTN25 grid corresponding to Ceuta, Melilla, Alboran and Spanish territories in North Africa, adjusted to the new official geodetic reference system ETRS89, in geographical coordinates (longitude, latitude).
MTN25_ETRS89_Peninsula_Baleares_Canarias
MTN25 real grid corresponding to the Peninsula, the Balearic Islands and the Canary Islands, adjusted to the new ETRS89 official reference geodetic system, in geographical coordinates (longitude, latitude).
MTN50_ETRS89_Peninsula_Baleares_Canarias
MTN50 real grid corresponding to the Peninsula, the Balearic Islands and the Canary Islands, adjusted to the new ETRS89 official reference geodetic system, in geographical coordinates (longitude, latitude).
MTN25_RegCan95_Canarias
MTN25 grid corresponding to the Canary Islands, in REGCAN95 (WGS84 compatible) and geographic coordinates (longitude, latitude). It is the real MTN25 grid, that is, the one that divides the current printed series of the map, taking into account the special distribution of the Canary Islands sheets.
MTN50_RegCan95_Canarias
MTN50 grid corresponding to the Canary Islands, in REGCAN95 (WGS84 compatible) and geographic coordinates (longitude, latitude). This is the real grid of the MTN50, that is, the one that divides the current printed series of the map, taking into account the special distribution of the Canary Islands sheets.
A sf polygon
You can set your cache_dir with esp_set_cache_dir().
Sometimes cached files may be corrupt. On that case, try re-downloading
the data setting update_cache = TRUE.
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. Use the option verbose = TRUE for debugging the API query.
IGN data via a custom CDN (see https://github.com/rOpenSpain/mapSpain/tree/sianedata/MTN).
Other grids:
esp_get_grid_BDN(),
esp_get_grid_EEA(),
esp_get_grid_ESDAC()
grid <- esp_get_grid_MTN(grid = "MTN50_ETRS89_Peninsula_Baleares_Canarias")
library(ggplot2)
ggplot(grid) +
geom_sf() +
theme_light() +
labs(title = "MTN50 Grid for Spain")Loads a hexbin map (sf object) or a map of squares with the boundaries of the provinces or autonomous communities of Spain.
esp_get_hex_prov(prov = NULL)
esp_get_hex_ccaa(ccaa = NULL)
esp_get_grid_prov(prov = NULL)
esp_get_grid_ccaa(ccaa = NULL)prov |
A vector of names and/or codes for provinces
or |
ccaa |
A vector of names and/or codes for autonomous communities
or |
Hexbin or grid map has an advantage over usual choropleth maps. In choropleths, a large polygon data looks more emphasized just because of its size, what introduces a bias. Here with hexbin, each region is represented equally dismissing the bias.
You can use and mix names, ISO codes, "codauto"/"cpro" codes (see esp_codelist) and NUTS codes of different levels.
When using a code corresponding of a higher level (e.g.
esp_get_prov("Andalucia")) all the corresponding units of that level are
provided (in this case , all the provinces of Andalucia).
Results are provided in EPSG:4258, use sf::st_transform()
to change the projection.
A sf POLYGON object.
Other political:
esp_codelist,
esp_get_can_box(),
esp_get_capimun(),
esp_get_ccaa(),
esp_get_comarca(),
esp_get_country(),
esp_get_munic(),
esp_get_nuts(),
esp_get_prov(),
esp_get_simpl_prov()
esp <- esp_get_country()
hexccaa <- esp_get_hex_ccaa()
library(ggplot2)
ggplot(hexccaa) +
geom_sf(data = esp) +
geom_sf(aes(fill = codauto),
alpha = 0.3,
show.legend = FALSE
) +
geom_sf_text(aes(label = label), check_overlap = TRUE) +
theme_void() +
labs(title = "Hexbin: CCAA")
hexprov <- esp_get_hex_prov()
ggplot(hexprov) +
geom_sf(data = esp) +
geom_sf(aes(fill = codauto),
alpha = 0.3,
show.legend = FALSE
) +
geom_sf_text(aes(label = label), check_overlap = TRUE) +
theme_void() +
labs(title = "Hexbin: Provinces")
gridccaa <- esp_get_grid_ccaa()
ggplot(gridccaa) +
geom_sf(data = esp) +
geom_sf(aes(fill = codauto),
alpha = 0.3,
show.legend = FALSE
) +
geom_sf_text(aes(label = label), check_overlap = TRUE) +
theme_void() +
labs(title = "Grid: CCAA")
gridprov <- esp_get_grid_prov()
ggplot(gridprov) +
geom_sf(data = esp) +
geom_sf(aes(fill = codauto),
alpha = 0.3,
show.legend = FALSE
) +
geom_sf_text(aes(label = label), check_overlap = TRUE) +
theme_void() +
labs(title = "Grid: Provinces")Loads a sf polygon object containing areas with the required hydrographic elements of Spain.
esp_get_hydrobasin(
epsg = "4258",
cache = TRUE,
update_cache = FALSE,
cache_dir = NULL,
verbose = FALSE,
resolution = "3",
domain = "land"
)epsg |
projection of the map: 4-digit EPSG code. One of:
|
cache |
A logical whether to do caching. Default is |
update_cache |
A logical whether to update cache. Default is |
cache_dir |
A path to a cache directory. See About caching. |
verbose |
Logical, displays information. Useful for debugging,
default is |
resolution |
Resolution of the polygon. Values available are
|
domain |
Possible values are |
Metadata available on https://github.com/rOpenSpain/mapSpain/tree/sianedata/.
A sf polygon object.
You can set your cache_dir with esp_set_cache_dir().
Sometimes cached files may be corrupt. On that case, try re-downloading
the data setting update_cache = TRUE.
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. Use the option verbose = TRUE for debugging the API query.
IGN data via a custom CDN (see https://github.com/rOpenSpain/mapSpain/tree/sianedata).
Other natural:
esp_get_hypsobath(),
esp_get_rivers()
hydroland <- esp_get_hydrobasin(domain = "land")
hydrolandsea <- esp_get_hydrobasin(domain = "landsea")
library(ggplot2)
ggplot(hydroland) +
geom_sf(data = hydrolandsea, fill = "skyblue4", alpha = .4) +
geom_sf(fill = "skyblue", alpha = .5) +
geom_sf_text(aes(label = rotulo),
size = 3, check_overlap = TRUE,
fontface = "bold",
family = "serif"
) +
coord_sf(
xlim = c(-9.5, 4.5),
ylim = c(35, 44)
) +
theme_void()Loads a sf polygon or line object representing the hypsometry and bathymetry of Spain.
Hypsometry represents the the elevation and depth of features of the Earth's surface relative to mean sea level.
Bathymetry is the measurement of the depth of water in oceans, rivers, or lakes.
esp_get_hypsobath(
epsg = "4258",
cache = TRUE,
update_cache = FALSE,
cache_dir = NULL,
verbose = FALSE,
resolution = "3",
spatialtype = "area"
)epsg |
projection of the map: 4-digit EPSG code. One of:
|
cache |
A logical whether to do caching. Default is |
update_cache |
A logical whether to update cache. Default is |
cache_dir |
A path to a cache directory. See About caching. |
verbose |
Logical, displays information. Useful for debugging,
default is |
resolution |
Resolution of the shape. Values available
are |
spatialtype |
Spatial type of the output. Use |
Metadata available on https://github.com/rOpenSpain/mapSpain/tree/sianedata/.
A sf polygon or line object.
You can set your cache_dir with esp_set_cache_dir().
Sometimes cached files may be corrupt. On that case, try re-downloading
the data setting update_cache = TRUE.
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. Use the option verbose = TRUE for debugging the API query.
IGN data via a custom CDN (see https://github.com/rOpenSpain/mapSpain/tree/sianedata).
Other natural:
esp_get_hydrobasin(),
esp_get_rivers()
# This code would produce a nice plot - It will take a few seconds to run
library(ggplot2)
hypsobath <- esp_get_hypsobath()
# Error on the data provided - There is an empty shape
# Remove:
hypsobath <- hypsobath[!sf::st_is_empty(hypsobath), ]
# Tints from Wikipedia
# https://en.wikipedia.org/wiki/Wikipedia:WikiProject_Maps/Conventions/Topographic_maps
bath_tints <- colorRampPalette(
rev(
c(
"#D8F2FE", "#C6ECFF", "#B9E3FF",
"#ACDBFB", "#A1D2F7", "#96C9F0",
"#8DC1EA", "#84B9E3", "#79B2DE",
"#71ABD8"
)
)
)
hyps_tints <- colorRampPalette(
rev(
c(
"#F5F4F2", "#E0DED8", "#CAC3B8", "#BAAE9A",
"#AC9A7C", "#AA8753", "#B9985A", "#C3A76B",
"#CAB982", "#D3CA9D", "#DED6A3", "#E8E1B6",
"#EFEBC0", "#E1E4B5", "#D1D7AB", "#BDCC96",
"#A8C68F", "#94BF8B", "#ACD0A5"
)
)
)
levels <- sort(unique(hypsobath$val_inf))
# Create palette
br_bath <- length(levels[levels < 0])
br_terrain <- length(levels) - br_bath
pal <- c(bath_tints((br_bath)), hyps_tints((br_terrain)))
# Plot Canary Islands
ggplot(hypsobath) +
geom_sf(aes(fill = as.factor(val_inf)),
color = NA
) +
coord_sf(
xlim = c(-18.6, -13),
ylim = c(27, 29.5)
) +
scale_fill_manual(values = pal) +
guides(fill = guide_legend(
title = "Elevation",
direction = "horizontal",
label.position = "bottom",
title.position = "top",
nrow = 1
)) +
theme(legend.position = "bottom")
# Plot Mainland
ggplot(hypsobath) +
geom_sf(aes(fill = as.factor(val_inf)),
color = NA
) +
coord_sf(
xlim = c(-9.5, 4.4),
ylim = c(35.8, 44)
) +
scale_fill_manual(values = pal) +
guides(fill = guide_legend(
title = "Elevation",
reverse = TRUE,
keyheight = .8
))Returns municipalities of Spain as polygons at a specified scale.
esp_get_munic() uses GISCO (Eurostat) as source. Please use
giscoR::gisco_attributions()
esp_get_munic_siane() uses CartoBase ANE as source, provided by
Instituto Geografico Nacional (IGN), http://www.ign.es/web/ign/portal.
Years available are 2005 up to today.
esp_get_munic(
year = "2019",
epsg = "4258",
cache = TRUE,
update_cache = FALSE,
cache_dir = NULL,
verbose = FALSE,
region = NULL,
munic = NULL,
moveCAN = TRUE
)
esp_get_munic_siane(
year = Sys.Date(),
epsg = "4258",
cache = TRUE,
update_cache = FALSE,
cache_dir = NULL,
verbose = FALSE,
resolution = 3,
region = NULL,
munic = NULL,
moveCAN = TRUE,
rawcols = FALSE
)year |
Release year. See Details for years available. |
epsg |
projection of the map: 4-digit EPSG code. One of:
|
cache |
A logical whether to do caching. Default is |
update_cache |
A logical whether to update cache. Default is |
cache_dir |
A path to a cache directory. See About caching. |
verbose |
Logical, displays information. Useful for debugging,
default is |
region |
A vector of names and/or codes for provinces
or |
munic |
A name or |
moveCAN |
A logical |
resolution |
Resolution of the polygon. Values available are
|
rawcols |
Logical. Setting this to |
The years available are:
esp_get_munic(): year could be one of "2001", "2004", "2006", "2008",
"2010", "2013" and any year between 2016 and 2019.
See giscoR::gisco_get_lau(), giscoR::gisco_get_communes().
esp_get_munic_siane(): year could be passed as a single year ("YYYY"
format, as end of year) or as a specific date ("YYYY-MM-DD" format).
Historical information starts as of 2005.
When using region you can use and mix names and NUTS codes
(levels 1, 2 or 3), ISO codes (corresponding to level 2 or 3) or
"cpro" (see esp_codelist).
When calling a superior level (Province, Autonomous Community or NUTS1) , all the municipalities of that level would be added.
A sf polygon
You can set your cache_dir with esp_set_cache_dir().
Sometimes cached files may be corrupt. On that case, try re-downloading
the data setting update_cache = TRUE.
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. Use the option verbose = TRUE for debugging the API query.
While moveCAN is useful for visualization, it would alter the actual
geographic position of the Canary Islands. When using the output for
spatial analysis or using tiles (e.g. with esp_getTiles() or
addProviderEspTiles()) this option should be set to FALSE in order to
get the actual coordinates, instead of the modified ones.
IGN data via a custom CDN (see https://github.com/rOpenSpain/mapSpain/tree/sianedata).
giscoR::gisco_get_lau(), base::regex().
Other political:
esp_codelist,
esp_get_can_box(),
esp_get_capimun(),
esp_get_ccaa(),
esp_get_comarca(),
esp_get_country(),
esp_get_gridmap,
esp_get_nuts(),
esp_get_prov(),
esp_get_simpl_prov()
Other municipalities:
esp_get_capimun(),
esp_munic.sf
# Get munics
Base <- esp_get_munic(year = "2019", region = "Castilla y Leon")
# Provs for delimiting
provs <- esp_get_prov(prov = "Castilla y Leon")
# Load population data
data("pobmun19")
# Arrange and create breaks
Base_pop <- merge(Base, pobmun19,
by = c("cpro", "cmun"),
all.x = TRUE
)
br <- sort(c(
0, 50, 100, 200, 500,
1000, 5000, 50000, 100000,
Inf
))
Base_pop$cuts <- cut(Base_pop$pob19, br, dig.lab = 20)
# Plot
library(ggplot2)
ggplot(Base_pop) +
geom_sf(aes(fill = cuts), color = NA) +
geom_sf(data = provs, fill = NA, color = "grey70") +
scale_fill_manual(values = hcl.colors(length(br), "cividis")) +
labs(
title = "Population in Castilla y Leon",
subtitle = "INE, 2019",
fill = "Persons"
) +
theme_void()Returns NUTS regions of Spain as polygons and points at a specified scale, as provided by GISCO (Geographic Information System of the Commission, depending of Eurostat).
NUTS are provided at three different levels:
"0": Country level
"1": Groups of autonomous communities
"2": Autonomous communities
"3": Roughly matches the provinces, but providing specific individual objects for each major island
esp_get_nuts(
year = "2016",
epsg = "4258",
cache = TRUE,
update_cache = FALSE,
cache_dir = NULL,
verbose = FALSE,
resolution = "01",
spatialtype = "RG",
region = NULL,
nuts_level = "all",
moveCAN = TRUE
)year |
Release year of the file. One of |
epsg |
projection of the map: 4-digit EPSG code. One of:
|
cache |
A logical whether to do caching. Default is |
update_cache |
A logical whether to update cache. Default is |
cache_dir |
A path to a cache directory. See About caching. |
verbose |
Logical, displays information. Useful for debugging,
default is |
resolution |
Resolution of the geospatial data. One of
|
spatialtype |
Type of geometry to be returned:
|
region |
Optional. A vector of region names, NUTS or ISO codes
(see |
nuts_level |
NUTS level. One of "0" (Country-level), "1", "2" or "3". See Description. |
moveCAN |
A logical |
A sf object specified by spatialtype.
You can set your cache_dir with esp_set_cache_dir().
Sometimes cached files may be corrupt. On that case, try re-downloading
the data setting update_cache = TRUE.
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. Use the option verbose = TRUE for debugging the API query.
While moveCAN is useful for visualization, it would alter the actual
geographic position of the Canary Islands. When using the output for
spatial analysis or using tiles (e.g. with esp_getTiles() or
addProviderEspTiles()) this option should be set to FALSE in order to
get the actual coordinates, instead of the modified ones.
Please check the download and usage provisions on
giscoR::gisco_attributions()
giscoR::gisco_get_nuts(), esp_dict_region_code().
Other political:
esp_codelist,
esp_get_can_box(),
esp_get_capimun(),
esp_get_ccaa(),
esp_get_comarca(),
esp_get_country(),
esp_get_gridmap,
esp_get_munic(),
esp_get_prov(),
esp_get_simpl_prov()
Other nuts:
esp_nuts.sf
NUTS1 <- esp_get_nuts(nuts_level = 1, moveCAN = TRUE)
library(ggplot2)
ggplot(NUTS1) +
geom_sf() +
labs(
title = "NUTS1: Displacing Canary Islands",
caption = giscoR::gisco_attributions()
)
NUTS1_alt <- esp_get_nuts(nuts_level = 1, moveCAN = c(15, 0))
ggplot(NUTS1_alt) +
geom_sf() +
labs(
title = "NUTS1: Displacing Canary Islands",
subtitle = "to the right",
caption = giscoR::gisco_attributions()
)
NUTS1_orig <- esp_get_nuts(nuts_level = 1, moveCAN = FALSE)
ggplot(NUTS1_orig) +
geom_sf() +
labs(
title = "NUTS1",
subtitle = "Canary Islands on the true location",
caption = giscoR::gisco_attributions()
)
AndOriental <-
esp_get_nuts(region = c("Almeria", "Granada", "Jaen", "Malaga"))
ggplot(AndOriental) +
geom_sf()
RandomRegions <- esp_get_nuts(region = c("ES1", "ES300", "ES51"))
ggplot(RandomRegions) +
geom_sf() +
labs(title = "Random Regions")
MixingCodes <- esp_get_nuts(region = c("ES4", "ES-PV", "Valencia"))
ggplot(MixingCodes) +
geom_sf() +
labs(title = "Mixing Codes")Returns provinces of Spain as polygons and points at a specified scale.
esp_get_prov() uses GISCO (Eurostat) as source. Please use
giscoR::gisco_attributions()
esp_get_prov_siane() uses CartoBase ANE as source, provided by Instituto
Geografico Nacional (IGN), http://www.ign.es/web/ign/portal. Years
available are 2005 up to today.
esp_get_prov(prov = NULL, moveCAN = TRUE, ...)
esp_get_prov_siane(
prov = NULL,
year = Sys.Date(),
epsg = "4258",
cache = TRUE,
update_cache = FALSE,
cache_dir = NULL,
verbose = FALSE,
resolution = "3",
moveCAN = TRUE,
rawcols = FALSE
)prov |
A vector of names and/or codes for provinces
or |
moveCAN |
A logical |
... |
Arguments passed on to
|
year |
Release year. See |
epsg |
projection of the map: 4-digit EPSG code. One of:
|
cache |
A logical whether to do caching. Default is |
update_cache |
A logical whether to update cache. Default is |
cache_dir |
A path to a cache directory. See About caching. |
verbose |
Logical, displays information. Useful for debugging,
default is |
resolution |
Resolution of the polygon. Values available are
|
rawcols |
Logical. Setting this to |
When using prov you can use and mix names and NUTS codes (levels 1, 2 or
3), ISO codes (corresponding to level 2 or 3) or "cpro" (see
esp_codelist).
Ceuta and Melilla are considered as provinces on this dataset.
When calling a superior level (Autonomous Community or NUTS1) , all the provinces of that level would be added.
On esp_get_prov_siane(), year could be passed as a single year ("YYYY"
format, as end of year) or as a specific date ("YYYY-MM-DD" format).
Historical information starts as of 2005.
A sf object specified by spatialtype.
You can set your cache_dir with esp_set_cache_dir().
Sometimes cached files may be corrupt. On that case, try re-downloading
the data setting update_cache = TRUE.
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. Use the option verbose = TRUE for debugging the API query.
While moveCAN is useful for visualization, it would alter the actual
geographic position of the Canary Islands. When using the output for
spatial analysis or using tiles (e.g. with esp_getTiles() or
addProviderEspTiles()) this option should be set to FALSE in order to
get the actual coordinates, instead of the modified ones.
IGN data via a custom CDN (see https://github.com/rOpenSpain/mapSpain/tree/sianedata).
Other political:
esp_codelist,
esp_get_can_box(),
esp_get_capimun(),
esp_get_ccaa(),
esp_get_comarca(),
esp_get_country(),
esp_get_gridmap,
esp_get_munic(),
esp_get_nuts(),
esp_get_simpl_prov()
prov <- esp_get_prov()
library(ggplot2)
ggplot(prov) +
geom_sf() +
theme_void()
# Random Provinces
Random <-
esp_get_prov(prov = c(
"Zamora",
"Palencia",
"ES-GR",
"ES521",
"01"
))
ggplot(Random) +
geom_sf(aes(fill = codauto), show.legend = FALSE, alpha = 0.5) +
scale_fill_manual(values = hcl.colors(
nrow(Random), "Spectral"
)) +
theme_minimal()
# All Provinces of a Zone plus an addition
Mix <- esp_get_prov(prov = c(
"Noroeste",
"Castilla y Leon", "La Rioja"
))
Mix$CCAA <- esp_dict_region_code(
Mix$codauto,
origin = "codauto"
)
ggplot(Mix) +
geom_sf(aes(fill = CCAA), alpha = 0.5) +
scale_fill_discrete(type = hcl.colors(5, "Temps")) +
theme_classic()
# ISO codes available
allprovs <- esp_get_prov()
ggplot(allprovs) +
geom_sf(fill = NA) +
geom_sf_text(aes(label = iso2.prov.code),
check_overlap = TRUE,
fontface = "bold"
) +
theme_void()Loads a sf lines or point object representing the nodes and railway lines of Spain.
esp_get_railway(
year = Sys.Date(),
epsg = "4258",
cache = TRUE,
update_cache = FALSE,
cache_dir = NULL,
verbose = FALSE,
spatialtype = "line"
)year |
Release year. |
epsg |
projection of the map: 4-digit EPSG code. One of:
|
cache |
A logical whether to do caching. Default is |
update_cache |
A logical whether to update cache. Default is |
cache_dir |
A path to a cache directory. See About caching. |
verbose |
Logical, displays information. Useful for debugging,
default is |
spatialtype |
Spatial type of the output. Use |
A sf line or point object.
You can set your cache_dir with esp_set_cache_dir().
Sometimes cached files may be corrupt. On that case, try re-downloading
the data setting update_cache = TRUE.
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. Use the option verbose = TRUE for debugging the API query.
IGN data via a custom CDN (see https://github.com/rOpenSpain/mapSpain/tree/sianedata).
Other infrastructure:
esp_get_roads()
provs <- esp_get_prov()
ccaa <- esp_get_ccaa()
# Railways
rails <- esp_get_railway()
# Stations
stations <- esp_get_railway(spatialtype = "point")
# Map
library(ggplot2)
ggplot(provs) +
geom_sf(fill = "grey99", color = "grey50") +
geom_sf(data = ccaa, fill = NA) +
geom_sf(
data = rails, aes(color = tipo),
show.legend = FALSE, linewidth = 1.5
) +
geom_sf(
data = stations,
color = "red", alpha = 0.5
) +
coord_sf(
xlim = c(-7.5, -2.5),
ylim = c(38, 41)
) +
scale_color_manual(values = hcl.colors(
length(unique(rails$tipo)), "viridis"
)) +
theme_minimal()Loads a sf polygon or line object representing rivers, channels, reservoirs and other wetlands of Spain
esp_get_rivers(
epsg = "4258",
cache = TRUE,
update_cache = FALSE,
cache_dir = NULL,
verbose = FALSE,
resolution = "3",
spatialtype = "line",
name = NULL
)epsg |
projection of the map: 4-digit EPSG code. One of:
|
cache |
A logical whether to do caching. Default is |
update_cache |
A logical whether to update cache. Default is |
cache_dir |
A path to a cache directory. See About caching. |
verbose |
Logical, displays information. Useful for debugging,
default is |
resolution |
Resolution of the polygon. Values available are
|
spatialtype |
Spatial type of the output. Use |
name |
Optional. A character or |
Metadata available on https://github.com/rOpenSpain/mapSpain/tree/sianedata/.
A sf polygon or line object.
IGN data via a custom CDN (see https://github.com/rOpenSpain/mapSpain/tree/sianedata).
Other natural:
esp_get_hydrobasin(),
esp_get_hypsobath()
# Use of regex
regex1 <- esp_get_rivers(name = "Tajo|Segura")
unique(regex1$rotulo)
regex2 <- esp_get_rivers(name = "Tajo$| Segura")
unique(regex2$rotulo)
# See the diference
# Rivers in Spain
shapeEsp <- esp_get_country(moveCAN = FALSE)
MainRivers <-
esp_get_rivers(name = "Tajo$|Ebro$|Ebre$|Duero|Guadiana$|Guadalquivir")
sf::st_bbox(MainRivers)
library(ggplot2)
ggplot(shapeEsp) +
geom_sf() +
geom_sf(data = MainRivers, color = "skyblue", linewidth = 2) +
coord_sf(
xlim = c(-7.5, 1),
ylim = c(36.8, 43)
) +
theme_void()
# Wetlands in South-West Andalucia
and <- esp_get_prov(c("Huelva", "Sevilla", "Cadiz"))
Wetlands <- esp_get_rivers(spatialtype = "area")
ggplot(and) +
geom_sf() +
geom_sf(
data = Wetlands, fill = "skyblue",
color = "skyblue", alpha = 0.5
) +
coord_sf(
xlim = c(-7.5, -4.5),
ylim = c(36, 38.5)
) +
theme_void()Loads a sf line object representing the main roads of Spain.
esp_get_roads(
year = Sys.Date(),
epsg = "4258",
cache = TRUE,
update_cache = FALSE,
cache_dir = NULL,
verbose = FALSE,
moveCAN = TRUE
)year |
Release year. See Details for years available. |
epsg |
projection of the map: 4-digit EPSG code. One of:
|
cache |
A logical whether to do caching. Default is |
update_cache |
A logical whether to update cache. Default is |
cache_dir |
A path to a cache directory. See About caching. |
verbose |
Logical, displays information. Useful for debugging,
default is |
moveCAN |
A logical |
year could be passed as a single year ("YYYY" format, as end of year) or
as a specific date ("YYYY-MM-DD" format).
A sf line object.
You can set your cache_dir with esp_set_cache_dir().
Sometimes cached files may be corrupt. On that case, try re-downloading
the data setting update_cache = TRUE.
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. Use the option verbose = TRUE for debugging the API query.
While moveCAN is useful for visualization, it would alter the actual
geographic position of the Canary Islands. When using the output for
spatial analysis or using tiles (e.g. with esp_getTiles() or
addProviderEspTiles()) this option should be set to FALSE in order to
get the actual coordinates, instead of the modified ones.
IGN data via a custom CDN (see https://github.com/rOpenSpain/mapSpain/tree/sianedata).
Other infrastructure:
esp_get_railway()
country <- esp_get_country()
Roads <- esp_get_roads()
library(ggplot2)
ggplot(country) +
geom_sf(fill = "grey90") +
geom_sf(data = Roads, aes(color = tipo), show.legend = "line") +
scale_color_manual(
values = c("#003399", "#003399", "#ff0000", "#ffff00")
) +
guides(color = guide_legend(direction = "vertical")) +
theme_minimal() +
labs(color = "Road type") +
theme(legend.position = "bottom")Loads a simplified map (sf object) with the boundaries of the provinces or autonomous communities of Spain, as provided by the INE (Instituto Nacional de Estadistica).
esp_get_simpl_prov(
prov = NULL,
update_cache = FALSE,
cache_dir = NULL,
verbose = FALSE
)
esp_get_simpl_ccaa(
ccaa = NULL,
update_cache = FALSE,
cache_dir = NULL,
verbose = FALSE
)prov |
A vector of names and/or codes for provinces
or |
update_cache |
A logical whether to update cache. Default is |
cache_dir |
A path to a cache directory. See About caching. |
verbose |
Logical, displays information. Useful for debugging,
default is |
ccaa |
A vector of names and/or codes for autonomous communities
or |
Results are provided without CRS, as provided on source.
You can use and mix names, ISO codes, "codauto"/"cpro" codes (see esp_codelist) and NUTS codes of different levels.
When using a code corresponding of a higher level (e.g.
esp_get_simpl_prov("Andalucia")) all the corresponding units of that level
are provided (in this case , all the provinces of Andalucia).
A sf POLYGON object.
You can set your cache_dir with esp_set_cache_dir().
Sometimes cached files may be corrupt. On that case, try re-downloading
the data setting update_cache = TRUE.
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. Use the option verbose = TRUE for debugging the API query.
INE: PC_Axis files
esp_get_hex_prov(), esp_get_hex_ccaa()
Other political:
esp_codelist,
esp_get_can_box(),
esp_get_capimun(),
esp_get_ccaa(),
esp_get_comarca(),
esp_get_country(),
esp_get_gridmap,
esp_get_munic(),
esp_get_nuts(),
esp_get_prov()
prov_simp <- esp_get_simpl_prov()
library(ggplot2)
ggplot(prov_simp) +
geom_sf(aes(fill = ine.ccaa.name)) +
labs(fill = "CCAA")
# Provs of Single CCAA
and_simple <- esp_get_simpl_prov("Andalucia")
ggplot(and_simple) +
geom_sf()
# CCAAs
ccaa_simp <- esp_get_simpl_ccaa()
ggplot(ccaa_simp) +
geom_sf() +
geom_sf_text(aes(label = ine.ccaa.name), check_overlap = TRUE)Get static map tiles based on a spatial object. Maps can be fetched from various open map servers.
This function is a implementation of the javascript plugin leaflet-providersESP v1.3.3.
esp_getTiles(
x,
type = "IDErioja",
zoom = NULL,
zoommin = 0,
crop = TRUE,
res = 512,
bbox_expand = 0.05,
transparent = TRUE,
mask = FALSE,
update_cache = FALSE,
cache_dir = NULL,
verbose = FALSE,
options = NULL
)x |
An sf or |
type |
This parameter could be either:
|
zoom |
Zoom level. If |
zoommin |
Delta on default |
crop |
|
res |
Resolution (in pixels) of the final tile. Only valid for WMS. |
bbox_expand |
A numeric value that indicates the expansion percentage
of the bounding box of |
transparent |
Logical. Provides transparent background, if supported.
Depends on the selected provider on |
mask |
|
update_cache |
A logical whether to update cache. Default is |
cache_dir |
A path to a cache directory. See About caching. |
verbose |
Logical, displays information. Useful for debugging,
default is |
options |
A named list containing additional options to pass to the query. |
Zoom levels are described on the OpenStreetMap wiki:
| zoom | area to represent |
| 0 | whole world |
| 3 | large country |
| 5 | state |
| 8 | county |
| 10 | metropolitan area |
| 11 | city |
| 13 | village or suburb |
| 16 | streets |
| 18 | some buildings, trees |
For a complete list of providers see esp_tiles_providers.
Most WMS/WMTS providers provide tiles on "EPSG:3857". In case that the tile
looks deformed, try projecting first x:
x <- sf::st_transform(x, 3857)
A SpatRaster is returned, with 3 (RGB) or 4 (RGBA) layers, depending on
the provider. See terra::rast().
.
You can set your cache_dir with esp_set_cache_dir().
Sometimes cached files may be corrupt. On that case, try re-downloading
the data setting update_cache = TRUE.
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. Use the option verbose = TRUE for debugging the API query.
https://dieghernan.github.io/leaflet-providersESP/ leaflet plugin, v1.3.3.
Other imagery utilities:
addProviderEspTiles(),
esp_make_provider(),
esp_tiles_providers
## Not run:
# This script downloads tiles to your local machine
# Run only if you are online
segovia <- esp_get_prov_siane("segovia", epsg = 3857)
tile <- esp_getTiles(segovia, "IGNBase.Todo")
library(ggplot2)
library(tidyterra)
ggplot(segovia) +
geom_spatraster_rgb(data = tile) +
geom_sf(fill = NA)
# Another provider
tile2 <- esp_getTiles(segovia, type = "MDT")
ggplot(segovia) +
geom_spatraster_rgb(data = tile2) +
geom_sf(fill = NA)
# A custom WMS provided
custom_wms <- esp_make_provider(
id = "an_id_for_caching",
q = "https://idecyl.jcyl.es/geoserver/ge/wms?",
service = "WMS",
version = "1.3.0",
format = "image/png",
layers = "geolog_cyl_litologia"
)
custom_wms_tile <- esp_getTiles(segovia, custom_wms)
autoplot(custom_wms_tile) +
geom_sf(data = segovia, fill = NA, color = "red")
# A custom WMTS provider
custom_wmts <- esp_make_provider(
id = "cyl_wmts",
q = "https://www.ign.es/wmts/pnoa-ma?",
service = "WMTS",
layer = "OI.OrthoimageCoverage"
)
custom_wmts_tile <- esp_getTiles(segovia, custom_wmts)
autoplot(custom_wmts_tile) +
geom_sf(data = segovia, fill = NA, color = "white", linewidth = 2)
# Example from https://leaflet-extras.github.io/leaflet-providers/preview/
stamen_water <- list(
id = "Stamen_Water",
q = "https://stamen-tiles-b.a.ssl.fastly.net/watercolor/{z}/{x}/{y}.jpg"
)
stamen <- esp_getTiles(segovia, stamen_water, zoommin = 1)
autoplot(stamen) +
geom_sf(data = segovia, fill = NA, color = "white", linewidth = 1)
## End(Not run)Helper function for esp_getTiles() that helps to create a custom provider.
esp_make_provider(id, q, service, layers, ...)id |
An identifier for the user. Would be used also for identifying cached tiles. |
q |
The base url of the service |
service |
The type of tile service, either |
layers |
The name of the layer to retrieve |
... |
Additional parameters to the query, like |
This function is meant to work with services provided as of the OGC Standard.
Note that:
mapSpain would not provide advice on the parameter q to be
provided.
Currently, on WMTS requests only services with
tilematrixset=GoogleMapsCompatible are supported.
A named list with two elements id and q.
For a list of potential providers from Spain check IDEE Directory.
Other imagery utilities:
addProviderEspTiles(),
esp_getTiles(),
esp_tiles_providers
## Not run:
# This script downloads tiles to your local machine
# Run only if you are online
custom_wms <- esp_make_provider(
id = "an_id_for_caching",
q = "https://idecyl.jcyl.es/geoserver/ge/wms?",
service = "WMS",
layers = "geolog_cyl_litologia"
)
x <- esp_get_ccaa("Castilla y León", epsg = 3857)
mytile <- esp_getTiles(x, type = custom_wms)
tidyterra::autoplot(mytile) +
ggplot2::geom_sf(data = x, fill = NA)
## End(Not run)POLYGON object of Spain (2019)A sf object including all municipalities of Spain as provided by GISCO (2019 version).
A sf object (resolution: 1:1million, EPSG:4258) object with
8,131 rows and columns:
INE code of the autonomous community.
INE name of the autonomous community.
INE code of the province.
INE name of the province.
INE code of the municipality.
Name of the municipality.
LAU Code (GISCO) of the municipality. This is a combination of cpro and cmun, aligned with INE coding scheme.
geometry field.
https://ec.europa.eu/eurostat/web/gisco/geodata/reference-data/, LAU 2019 data.
Other datasets:
esp_codelist,
esp_nuts.sf,
esp_tiles_providers,
pobmun19
Other municipalities:
esp_get_capimun(),
esp_get_munic()
data("esp_munic.sf")
teruel_cpro <- esp_dict_region_code("Teruel", destination = "cpro")
teruel_sf <- esp_munic.sf[esp_munic.sf$cpro == teruel_cpro, ]
teruel_city <- teruel_sf[teruel_sf$name == "Teruel", ]
# Plot
library(ggplot2)
ggplot(teruel_sf) +
geom_sf(fill = "#FDFBEA") +
geom_sf(data = teruel_city, aes(fill = name)) +
scale_fill_manual(
values = "#C12838",
labels = "City of Teruel"
) +
labs(
fill = "",
title = "Municipalities of Teruel"
) +
theme_minimal() +
theme(
text = element_text(face = "bold"),
panel.background = element_rect(colour = "black"),
panel.grid = element_blank(),
legend.position = c(.2, .95)
)POLYGON object of SpainA sf object including all NUTS levels of Spain as provided by GISCO (2016 version).
A sf object (resolution: 1:1million, EPSG:4258) with
86 rows and columns:
NUTS level code (0,1,2,3)
NUTS identifier
Urban Type, see Details
Eurostat Country code ES
NUTS name on Latin characters
NUTS name on local alphabet
Mount Type, see Details
Coast Type, see Details
FID
geometry field
MOUNT_TYPE: Mountain typology:
1: More than 50 % of the surface is covered by topographic mountain areas.
2: More than 50 % of the regional population lives in topographic mountain areas.
3: More than 50 % of the surface is covered by topographic mountain areas and where more than 50 % of the regional population lives in these mountain areas.
4: Non-mountain region / other regions.
0: No classification provided
URBN_TYPE: Urban-rural typology:
1: Predominantly urban region.
2: Intermediate region.
3: Predominantly rural region.
0: No classification provided
COAST_TYPE: Coastal typology:
1: Coastal (on coast).
2: Coastal (less than 50% of population living within 50 km. of the coastline).
3: Non-coastal region.
0: No classification provided
https://gisco-services.ec.europa.eu/distribution/v2/nuts/, file
NUTS_RG_20M_2016_4326.geojson.
Other datasets:
esp_codelist,
esp_munic.sf,
esp_tiles_providers,
pobmun19
Other nuts:
esp_get_nuts()
data("esp_nuts.sf")
nuts <- esp_nuts.sf
# Select NUTS 3
nuts3 <- esp_nuts.sf[esp_nuts.sf$LEVL_CODE == 3, ]
# Combine with full shape
spain <- esp_get_country(moveCAN = FALSE)
# Plot Urban Type: See
# https://ec.europa.eu/eurostat/web/rural-development/methodology
library(ggplot2)
nuts3$URBN_TYPE_cat <- as.factor(nuts3$URBN_TYPE)
levels(nuts3$URBN_TYPE_cat)
levels(nuts3$URBN_TYPE_cat) <- c("Urban", "Intermediate", "Rural")
ggplot(nuts3) +
geom_sf(aes(fill = URBN_TYPE_cat), linewidth = .1) +
scale_fill_manual(values = c("grey80", "#FFC183", "#68AC20")) +
labs(
title = "NUTS3 levels of Spain",
fill = "Urban topology"
) +
theme_linedraw() +
theme(
legend.position = c(.8, .2)
)This function will store your cache_dir path on your local machine and
would load it for future sessions. Type Sys.getenv("MAPSPAIN_CACHE_DIR") to
find your cached path.
Alternatively, you can store the cache_dir manually with the following
options:
Run Sys.setenv(MAPSPAIN_CACHE_DIR = "cache_dir"). You would need to
run this command on each session (Similar to install = FALSE).
Set options(mapSpain_cache_dir = "cache_dir"). Similar to the previous
option. This is not recommended any more, and it is provided for
backwards compatibility purposes.
Write this line on your .Renviron file:
MAPSPAIN_CACHE_DIR = "value_for_cache_dir" (same behavior than
install = TRUE). This would store your cache_dir permanently.
esp_set_cache_dir(
cache_dir,
overwrite = FALSE,
install = FALSE,
verbose = TRUE
)cache_dir |
A path to a cache directory. On missing value the function
would store the cached files on a temporary dir (See |
overwrite |
If this is set to |
install |
if |
verbose |
Logical, displays information. Useful for debugging,
default is |
An (invisible) character with the path to your cache_dir.
Other cache utilities:
esp_clear_cache()
# Don't run this! It would modify your current state
## Not run:
esp_set_cache_dir(verbose = TRUE)
## End(Not run)
Sys.getenv("MAPSPAIN_CACHE_DIR")A named list of length 102 containing the parameters of the url information of different public WMS and WMTS providers of Spain.
Implementation of javascript plugin leaflet-providersESP v1.3.3.
A named list of the providers available with the following structure:
Each item of the list is named with the provider alias.
Each element of the list contains two nested named lists:
static with the parameters to get static tiles plus an additional item
named attribution.
leaflet with additional parameters to be passed onto
addProviderEspTiles().
Providers available to be passed to type on esp_getTiles() are:
"IDErioja"
"IDErioja.Base"
"IDErioja.Relieve"
"IDErioja.Claro"
"IDErioja.Oscuro"
"IGNBase"
"IGNBase.Todo"
"IGNBase.Gris"
"IGNBase.TodoNoFondo"
"IGNBase.Orto"
"MDT"
"MDT.Elevaciones"
"MDT.Relieve"
"MDT.CurvasNivel"
"MDT.SpotElevation"
"PNOA"
"PNOA.MaximaActualidad"
"PNOA.Mosaico"
"OcupacionSuelo"
"OcupacionSuelo.Ocupacion"
"OcupacionSuelo.Usos"
"LiDAR"
"MTN"
"Geofisica"
"Geofisica.Terremotos10dias"
"Geofisica.Terremotos30dias"
"Geofisica.Terremotos365dias"
"Geofisica.ObservedEvents"
"Geofisica.HazardArea"
"VigilanciaVolcanica"
"VigilanciaVolcanica.ErupcionesHistoricas"
"CaminoDeSantiago"
"CaminoDeSantiago.CaminoFrances"
"CaminoDeSantiago.CaminosFrancia"
"CaminoDeSantiago.CaminosGalicia"
"CaminoDeSantiago.CaminosDelNorte"
"CaminoDeSantiago.CaminosAndaluces"
"CaminoDeSantiago.CaminosCentro"
"CaminoDeSantiago.CaminosEste"
"CaminoDeSantiago.CaminosCatalanes"
"CaminoDeSantiago.CaminosSureste"
"CaminoDeSantiago.CaminosInsulares"
"CaminoDeSantiago.CaminosPortugueses"
"Catastro"
"Catastro.Catastro"
"Catastro.Parcela"
"Catastro.CadastralParcel"
"Catastro.CadastralZoning"
"Catastro.Address"
"Catastro.Building"
"Catastro.BuildingPart"
"Catastro.AdministrativeBoundary"
"Catastro.AdministrativeUnit"
"RedTransporte"
"RedTransporte.Carreteras"
"RedTransporte.Ferroviario"
"RedTransporte.Aerodromo"
"RedTransporte.AreaServicio"
"RedTransporte.EstacionesFerroviario"
"RedTransporte.Puertos"
"Cartociudad"
"Cartociudad.CodigosPostales"
"Cartociudad.Direcciones"
"NombresGeograficos"
"UnidadesAdm"
"UnidadesAdm.Limites"
"UnidadesAdm.Unidades"
"Hidrografia"
"Hidrografia.MasaAgua"
"Hidrografia.Cuencas"
"Hidrografia.Subcuencas"
"Hidrografia.POI"
"Hidrografia.ManMade"
"Hidrografia.LineaCosta"
"Hidrografia.Rios"
"Hidrografia.Humedales"
"Militar"
"Militar.CEGET1M"
"Militar.CEGETM7814"
"Militar.CEGETM7815"
"Militar.CEGETM682"
"Militar.CECAF1M"
"ADIF"
"ADIF.Vias"
"ADIF.Nodos"
"ADIF.Estaciones"
"LimitesMaritimos"
"LimitesMaritimos.LimitesMaritimos"
"LimitesMaritimos.LineasBase"
"Copernicus"
"Copernicus.Forest"
"Copernicus.ForestLeaf"
"Copernicus.WaterWet"
"Copernicus.SoilSeal"
"Copernicus.GrassLand"
"Copernicus.RiparianGreen"
"Copernicus.RiparianLandCover"
"Copernicus.Natura2k"
"Copernicus.UrbanAtlas"
"ParquesNaturales"
"ParquesNaturales.Limites"
"ParquesNaturales.ZonasPerifericas"
https://dieghernan.github.io/leaflet-providersESP/ leaflet plugin, v1.3.3.
Other datasets:
esp_codelist,
esp_munic.sf,
esp_nuts.sf,
pobmun19
Other imagery utilities:
addProviderEspTiles(),
esp_getTiles(),
esp_make_provider()
data("esp_tiles_providers")
# Get a single provider
single <- esp_tiles_providers[["IGNBase.Todo"]]
single$static
single$leafletPopulation by municipality (2019)
An example data frame with 8,131 rows containing the population data by municipality in Spain (2019).
INE code of the province.
name of the province.
INE code of the municipality.
Name of the municipality.
Overall population (2019)
Men population (2019)
Women population (2019)
INE: Instituto Nacional de Estadistica https://www.ine.es/
Other datasets:
esp_codelist,
esp_munic.sf,
esp_nuts.sf,
esp_tiles_providers
data("pobmun19")