Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

Add 2021 data and add new epsg options #190

Merged
merged 2 commits into from
Sep 19, 2020
Merged

Add 2021 data and add new epsg options #190

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
85d4ca0
Add 2021 data and add new epsg options
dieghernan Sep 17, 2020
f2209e7
Fix typo on man
dieghernan Sep 18, 2020
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
34 changes: 23 additions & 11 deletions R/get_eurostat_geospatial.R
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,4 +1,4 @@
#' @title Download Geospatial Data from GISGO
#' @title Download Geospatial Data from GISCO
#' @description Downloads either a simple features (sf), SpatialPolygonDataFrame or a
#' data_frame preprocessed using
#' \code{broom::tidy()}.
Expand All @@ -14,7 +14,7 @@
#' @param nuts_level Level of NUTS classification of the geospatial data. One of
#' "0", "1", "2", "3" or "all" (mimics the original behaviour)
#' @param year NUTS release year. One of
#' "2003", "2006", "2010", "2013" or "2016"
#' "2003", "2006", "2010", "2013", "2016" or "2021"
#' @param cache a logical whether to do caching. Default is \code{TRUE}. Affects
#' only queries from the bulk download facility.
#' @param update_cache a locigal whether to update cache. Can be set also with
Expand All @@ -24,6 +24,12 @@
#' 'eurostat' directory in the temporary directory from
#' \code{\link{tempdir}}. Directory can also be set with
#' \code{option} eurostat_cache_dir.
#' @param crs projection of the map: 4-digit \href{http://spatialreference.org/ref/epsg/}{EPSG code}. One of:
#' \itemize{
#' \item "4326" - WGS84
#' \item "3035" - ETRS89 / ETRS-LAEA
#' \item "3857" - Pseudo-Mercator
#' }
#' @export
#' @details The data source URL is \url{http://ec.europa.eu/eurostat/web/gisco/geodata/reference-data/administrative-units-statistical-units}.
#' @author Markus Kainu <markuskainu@gmail.com>
Expand All @@ -37,7 +43,7 @@ get_eurostat_geospatial <- function(output_class="sf",
resolution="60",
nuts_level = "all", year = "2016",
cache = TRUE, update_cache = FALSE,
cache_dir = NULL){
cache_dir = NULL, crs = "4326"){
# Check if you have access to ec.europe.eu.
if (!check_access_to_data()){
message("You have no access to ec.europe.eu.
Expand Down Expand Up @@ -66,13 +72,19 @@ Please check your connection and/or review your proxy settings")

# Check year is of correct format
year <- as.character(year)
if (!as.numeric(year) %in% c(2003, 2006, 2010, 2013, 2016)) {
stop("Year should be one of 2003, 2006, 2010, 2013 or 2016")
if (!as.numeric(year) %in% c(2003, 2006, 2010, 2013, 2016, 2021)) {
stop("Year should be one of 2003, 2006, 2010, 2013, 2016 or 2021")
}

if (as.numeric(year) == 2003 & as.numeric(resolution) == 60) {
stop("NUTS 2003 is not provided at 1:60 million resolution. Try 1:1 million, 1:3 million, 1:10 million or 1:20 million")
}

# Check crs is of correct format
crs <- as.character(crs)
if (!as.numeric(crs) %in% c(4326, 3035, 3857)) {
stop("crs should be one of 4326, 3035 or 3857")
}

# message("
# COPYRIGHT NOTICE
Expand Down Expand Up @@ -101,7 +113,7 @@ Please check your connection and/or review your proxy settings")
# information regarding their licence agreements.
# ")

if (resolution == "60" && year == 2016){
if (resolution == "60" && year == 2016 && crs == "4326"){

if (nuts_level %in% c("all")){
shp <- eurostat_geodata_60_2016
Expand Down Expand Up @@ -157,7 +169,7 @@ Please check your connection and/or review your proxy settings")
if (!cache || update_cache || !file.exists(cache_file)){

if (nuts_level %in% c("0","all")){
url <- paste0("http://ec.europa.eu/eurostat/cache/GISCO/distribution/v2/nuts/geojson/NUTS_RG_",resolution,"M_",year,"_4326_LEVL_0.geojson")
url <- paste0("http://ec.europa.eu/eurostat/cache/GISCO/distribution/v2/nuts/geojson/NUTS_RG_",resolution,"M_",year,"_",crs,"_LEVL_0.geojson")
resp <- RETRY("GET", url, terminate_on = c(404))
if (httr::http_error(resp)) {
stop(paste("The requested url cannot be found within the get_eurostat_geospatial function:", url))
Expand All @@ -167,7 +179,7 @@ Please check your connection and/or review your proxy settings")
}
}
if (nuts_level %in% c("1","all")){
url <- paste0("http://ec.europa.eu/eurostat/cache/GISCO/distribution/v2/nuts/geojson/NUTS_RG_",resolution,"M_",year,"_4326_LEVL_1.geojson")
url <- paste0("http://ec.europa.eu/eurostat/cache/GISCO/distribution/v2/nuts/geojson/NUTS_RG_",resolution,"M_",year,"_",crs,"_LEVL_1.geojson")
resp <- RETRY("GET", url, terminate_on = c(404))
if (httr::http_error(resp)) {
stop(paste("The requested url cannot be found within the get_eurostat_geospatial function:", url))
Expand All @@ -177,7 +189,7 @@ Please check your connection and/or review your proxy settings")
}
}
if (nuts_level %in% c("2","all")){
resp <- RETRY("GET", paste0("http://ec.europa.eu/eurostat/cache/GISCO/distribution/v2/nuts/geojson/NUTS_RG_",resolution,"M_",year,"_4326_LEVL_2.geojson"), terminate_on = c(404))
resp <- RETRY("GET", paste0("http://ec.europa.eu/eurostat/cache/GISCO/distribution/v2/nuts/geojson/NUTS_RG_",resolution,"M_",year,"_",crs,"_LEVL_2.geojson"), terminate_on = c(404))
if (httr::http_error(resp)) {
stop(paste("The requested url cannot be found within the get_eurostat_geospatial function:", url))
} else {
Expand All @@ -186,7 +198,7 @@ Please check your connection and/or review your proxy settings")
}
}
if (nuts_level %in% c("3","all")){
resp <- RETRY("GET", paste0("http://ec.europa.eu/eurostat/cache/GISCO/distribution/v2/nuts/geojson/NUTS_RG_",resolution,"M_",year,"_4326_LEVL_3.geojson"), terminate_on = c(404))
resp <- RETRY("GET", paste0("http://ec.europa.eu/eurostat/cache/GISCO/distribution/v2/nuts/geojson/NUTS_RG_",resolution,"M_",year,"_",crs,"_LEVL_3.geojson"), terminate_on = c(404))
if (httr::http_error(resp)) {
stop(paste("The requested url cannot be found within the get_eurostat_geospatial function:", url))
} else {
Expand Down Expand Up @@ -251,7 +263,7 @@ Please check your connection and/or review your proxy settings")
}
}

if (resolution == "60" & year == 2016){
if (resolution == "60" & year == 2016 & crs == "4326"){
if (output_class == "sf") message("sf at resolution 1:60 read from local file")
if (output_class == "df") message("data_frame at resolution 1:60 read from local file")
if (output_class == "spdf") message("SpatialPolygonDataFrame at resolution 1:60 read from local file")
Expand Down
14 changes: 11 additions & 3 deletions man/get_eurostat_geospatial.Rd
View file
Open in desktop

Some generated files are not rendered by default. Learn more about how customized files appear on GitHub.