[Feature request] Allow to cache by country

[ColinFay]

Expected behavior

Running gisco_get_nuts(country = "Belgium") then gisco_get_nuts(country = "Austria") caches both, not only the latest.

Observed behavior

The package only caches the latest value.

> tp_dir <- fs::path_temp("gisco_cache")
> fs::dir_create(tp_dir)
> gisco_get_nuts(
    year = "2024",
  nuts_level = 2,
  epsg = "3035",
  resolution = "60",
  country = "Austria",
  cache = TRUE,
  cache_dir = tp_dir
)
> fs::dir_ls(tp_dir)
/var/folders/9w/zdlv83ws6csdjnfc5x819gtr0000gn/T/Rtmp8GVGgl/gisco_cache/NUTS_RG_60M_2024_3035_LEVL_2.geojson

> gisco_get_nuts(
  year = "2024",
  nuts_level = 2,
  epsg = "3035",
  resolution = "60",
  country = "France",
  cache = TRUE,
  cache_dir = tp_dir
)
> fs::dir_ls(tp_dir)
/var/folders/9w/zdlv83ws6csdjnfc5x819gtr0000gn/T/Rtmp8GVGgl/gisco_cache/NUTS_RG_60M_2024_3035_LEVL_2.geojson

[dieghernan]

Hi @ColinFay

Please see https://ropengov.github.io/giscoR/reference/gisco_get_units.html

Download individual shapefiles of units. Unlike gisco_get_countries(), gisco_get_nuts() or gisco_get_urban_audit(), that downloads a full dataset and applies filters, gisco_get_units() downloads a single shapefile for each unit.

gisco_get_nuts() downloads the aggregated file "NUTS_RG_60M_2024_3035_LEVL_2.geojson" and filter it by country, that's why the file you see is the same, but results are different.

This is because GISCO does not provide NUTS files at country level (https://gisco-services.ec.europa.eu/distribution/v2/nuts/nuts-2024-files.html), only aggregated or at per unit level. So at this moment is not possible to implement this feature request, as our policy is to perform minimal modifications in the original source (only validity checks are performed).

If you want individual files use gisco_get_units() but this function is going to be replaced in giscoR v1.0.0:

Current CRAN giscoR version (0.6.1):

library(giscoR)


all_nuts2_codes <- gisco_get_nuts(
  year = "2024",
  nuts_level = 2,
  epsg = "3035",
  resolution = "60"
) |> 
  sf::st_drop_geometry()

# Austria data 
tp_dir_a <- fs::path_temp("Austria")

at <- gisco_get_units(
  id_giscoR = "nuts",
  unit = all_nuts2_codes[all_nuts2_codes$CNTR_CODE == "AT", "NUTS_ID"],
  year = "2024",
  epsg = "3035",
  resolution = "60",
  cache_dir = tp_dir_a
)

list.files(tp_dir_a, full.names = TRUE)
#> [1] "C:/Users/diego/AppData/Local/Temp/RtmpqqAzpN/Austria/AT11-region-60m-3035-2024.geojson"
#> [2] "C:/Users/diego/AppData/Local/Temp/RtmpqqAzpN/Austria/AT12-region-60m-3035-2024.geojson"
#> [3] "C:/Users/diego/AppData/Local/Temp/RtmpqqAzpN/Austria/AT13-region-60m-3035-2024.geojson"
#> [4] "C:/Users/diego/AppData/Local/Temp/RtmpqqAzpN/Austria/AT21-region-60m-3035-2024.geojson"
#> [5] "C:/Users/diego/AppData/Local/Temp/RtmpqqAzpN/Austria/AT22-region-60m-3035-2024.geojson"
#> [6] "C:/Users/diego/AppData/Local/Temp/RtmpqqAzpN/Austria/AT31-region-60m-3035-2024.geojson"
#> [7] "C:/Users/diego/AppData/Local/Temp/RtmpqqAzpN/Austria/AT32-region-60m-3035-2024.geojson"
#> [8] "C:/Users/diego/AppData/Local/Temp/RtmpqqAzpN/Austria/AT33-region-60m-3035-2024.geojson"
#> [9] "C:/Users/diego/AppData/Local/Temp/RtmpqqAzpN/Austria/AT34-region-60m-3035-2024.geojson"

# Belgium data
tp_dir_b <- fs::path_temp("Belgium")
# Getting the data is correct
be <- gisco_get_units(
  id_giscoR = "nuts",
  unit = all_nuts2_codes[all_nuts2_codes$CNTR_CODE == "BE", "NUTS_ID"],
  year = "2024",
  epsg = "3035",
  resolution = "60",
  cache_dir = tp_dir_b
)


list.files(tp_dir_b, full.names = TRUE)
#>  [1] "C:/Users/diego/AppData/Local/Temp/RtmpqqAzpN/Belgium/BE10-region-60m-3035-2024.geojson"
#>  [2] "C:/Users/diego/AppData/Local/Temp/RtmpqqAzpN/Belgium/BE21-region-60m-3035-2024.geojson"
#>  [3] "C:/Users/diego/AppData/Local/Temp/RtmpqqAzpN/Belgium/BE22-region-60m-3035-2024.geojson"
#>  [4] "C:/Users/diego/AppData/Local/Temp/RtmpqqAzpN/Belgium/BE23-region-60m-3035-2024.geojson"
#>  [5] "C:/Users/diego/AppData/Local/Temp/RtmpqqAzpN/Belgium/BE24-region-60m-3035-2024.geojson"
#>  [6] "C:/Users/diego/AppData/Local/Temp/RtmpqqAzpN/Belgium/BE25-region-60m-3035-2024.geojson"
#>  [7] "C:/Users/diego/AppData/Local/Temp/RtmpqqAzpN/Belgium/BE31-region-60m-3035-2024.geojson"
#>  [8] "C:/Users/diego/AppData/Local/Temp/RtmpqqAzpN/Belgium/BE32-region-60m-3035-2024.geojson"
#>  [9] "C:/Users/diego/AppData/Local/Temp/RtmpqqAzpN/Belgium/BE33-region-60m-3035-2024.geojson"
#> [10] "C:/Users/diego/AppData/Local/Temp/RtmpqqAzpN/Belgium/BE34-region-60m-3035-2024.geojson"
#> [11] "C:/Users/diego/AppData/Local/Temp/RtmpqqAzpN/Belgium/BE35-region-60m-3035-2024.geojson"

^{Created on 2025-12-08 with reprex v2.1.1}

Session info

sessioninfo::session_info()
#> Warning in system2("quarto", "-V", stdout = TRUE, env = paste0("TMPDIR=", : el
#> comando ejecutado '"quarto"
#> TMPDIR=C:/Users/diego/AppData/Local/Temp/RtmpqqAzpN/fileae05fab2c4f -V' tiene
#> el estatus 1
#> ─ Session info ───────────────────────────────────────────────────────────────
#>  setting  value
#>  version  R version 4.5.2 (2025-10-31 ucrt)
#>  os       Windows 11 x64 (build 26200)
#>  system   x86_64, mingw32
#>  ui       RTerm
#>  language (EN)
#>  collate  Spanish_Spain.utf8
#>  ctype    Spanish_Spain.utf8
#>  tz       Europe/Madrid
#>  date     2025-12-08
#>  pandoc   3.6.3 @ C:/Program Files/RStudio/resources/app/bin/quarto/bin/tools/ (via rmarkdown)
#>  quarto   NA @ C:\\Users\\diego\\AppData\\Local\\Programs\\Quarto\\bin\\quarto.exe
#> 
#> ─ Packages ───────────────────────────────────────────────────────────────────
#>  package     * version date (UTC) lib source
#>  class         7.3-23  2025-01-01 [2] CRAN (R 4.5.2)
#>  classInt      0.4-11  2025-01-08 [1] CRAN (R 4.5.0)
#>  cli           3.6.5   2025-04-23 [1] CRAN (R 4.5.0)
#>  DBI           1.2.3   2024-06-02 [1] CRAN (R 4.5.0)
#>  digest        0.6.39  2025-11-19 [1] CRAN (R 4.5.2)
#>  e1071         1.7-16  2024-09-16 [1] CRAN (R 4.5.0)
#>  evaluate      1.0.5   2025-08-27 [1] CRAN (R 4.5.1)
#>  fastmap       1.2.0   2024-05-15 [1] CRAN (R 4.5.0)
#>  fs            1.6.6   2025-04-12 [1] CRAN (R 4.5.0)
#>  geojsonsf     2.0.5   2025-11-26 [1] CRAN (R 4.5.2)
#>  giscoR      * 0.6.1   2025-01-27 [1] CRAN (R 4.5.2)
#>  glue          1.8.0   2024-09-30 [1] CRAN (R 4.5.0)
#>  htmltools     0.5.8.1 2024-04-04 [1] CRAN (R 4.5.0)
#>  KernSmooth    2.23-26 2025-01-01 [2] CRAN (R 4.5.2)
#>  knitr         1.50    2025-03-16 [1] CRAN (R 4.5.0)
#>  lifecycle     1.0.4   2023-11-07 [1] CRAN (R 4.5.0)
#>  magrittr      2.0.4   2025-09-12 [1] CRAN (R 4.5.1)
#>  proxy         0.4-27  2022-06-09 [1] CRAN (R 4.5.0)
#>  Rcpp          1.1.0   2025-07-02 [1] CRAN (R 4.5.1)
#>  reprex        2.1.1   2024-07-06 [1] CRAN (R 4.5.0)
#>  rlang         1.1.6   2025-04-11 [1] CRAN (R 4.5.0)
#>  rmarkdown     2.30    2025-09-28 [1] CRAN (R 4.5.1)
#>  rstudioapi    0.17.1  2024-10-22 [1] CRAN (R 4.5.0)
#>  sessioninfo   1.2.3   2025-02-05 [1] CRAN (R 4.5.0)
#>  sf            1.0-23  2025-11-28 [1] CRAN (R 4.5.2)
#>  units         1.0-0   2025-10-09 [1] CRAN (R 4.5.1)
#>  withr         3.0.2   2024-10-28 [1] CRAN (R 4.5.0)
#>  xfun          0.54    2025-10-30 [1] CRAN (R 4.5.2)
#>  yaml          2.3.11  2025-11-28 [1] CRAN (R 4.5.2)
#> 
#>  [1] C:/Users/diego/AppData/Local/R/win-library/4.5
#>  [2] C:/Program Files/R/R-4.5.2/library
#>  * ── Packages attached to the search path.
#> 
#> ──────────────────────────────────────────────────────────────────────────────

giscoR development version (upcoming v1.0.0)

Use gisco_get_unit_nuts()

# pak::pak("ropengov/giscoR")


library(giscoR)

all_nuts2_codes <- gisco_get_nuts(
  year = "2024",
  nuts_level = 2,
  epsg = "3035",
  resolution = "60"
) |> 
  sf::st_drop_geometry()

# Austria data 
tp_dir_a <- fs::path_temp("Austria")

at <- gisco_get_unit_nuts(
  unit = all_nuts2_codes[all_nuts2_codes$CNTR_CODE == "AT",]$NUTS_ID,
  year = "2024",
  epsg = "3035",
  resolution = "60",
  cache_dir = tp_dir_a
)

list.files(tp_dir_a, full.names = TRUE, recursive = TRUE)
#> [1] "C:/Users/diego/AppData/Local/Temp/Rtmpo1bQga/Austria/nuts/units/AT11-region-60m-3035-2024.geojson"
#> [2] "C:/Users/diego/AppData/Local/Temp/Rtmpo1bQga/Austria/nuts/units/AT12-region-60m-3035-2024.geojson"
#> [3] "C:/Users/diego/AppData/Local/Temp/Rtmpo1bQga/Austria/nuts/units/AT13-region-60m-3035-2024.geojson"
#> [4] "C:/Users/diego/AppData/Local/Temp/Rtmpo1bQga/Austria/nuts/units/AT21-region-60m-3035-2024.geojson"
#> [5] "C:/Users/diego/AppData/Local/Temp/Rtmpo1bQga/Austria/nuts/units/AT22-region-60m-3035-2024.geojson"
#> [6] "C:/Users/diego/AppData/Local/Temp/Rtmpo1bQga/Austria/nuts/units/AT31-region-60m-3035-2024.geojson"
#> [7] "C:/Users/diego/AppData/Local/Temp/Rtmpo1bQga/Austria/nuts/units/AT32-region-60m-3035-2024.geojson"
#> [8] "C:/Users/diego/AppData/Local/Temp/Rtmpo1bQga/Austria/nuts/units/AT33-region-60m-3035-2024.geojson"
#> [9] "C:/Users/diego/AppData/Local/Temp/Rtmpo1bQga/Austria/nuts/units/AT34-region-60m-3035-2024.geojson"

# Belgium data
tp_dir_b <- fs::path_temp("Belgium")
# Getting the data is correct
be <- gisco_get_unit_nuts(
  unit = all_nuts2_codes[all_nuts2_codes$CNTR_CODE == "BE", ]$NUTS_ID,
  year = "2024",
  epsg = "3035",
  resolution = "60",
  cache_dir = tp_dir_b
)


list.files(tp_dir_b, full.names = TRUE,recursive = TRUE)
#>  [1] "C:/Users/diego/AppData/Local/Temp/Rtmpo1bQga/Belgium/nuts/units/BE10-region-60m-3035-2024.geojson"
#>  [2] "C:/Users/diego/AppData/Local/Temp/Rtmpo1bQga/Belgium/nuts/units/BE21-region-60m-3035-2024.geojson"
#>  [3] "C:/Users/diego/AppData/Local/Temp/Rtmpo1bQga/Belgium/nuts/units/BE22-region-60m-3035-2024.geojson"
#>  [4] "C:/Users/diego/AppData/Local/Temp/Rtmpo1bQga/Belgium/nuts/units/BE23-region-60m-3035-2024.geojson"
#>  [5] "C:/Users/diego/AppData/Local/Temp/Rtmpo1bQga/Belgium/nuts/units/BE24-region-60m-3035-2024.geojson"
#>  [6] "C:/Users/diego/AppData/Local/Temp/Rtmpo1bQga/Belgium/nuts/units/BE25-region-60m-3035-2024.geojson"
#>  [7] "C:/Users/diego/AppData/Local/Temp/Rtmpo1bQga/Belgium/nuts/units/BE31-region-60m-3035-2024.geojson"
#>  [8] "C:/Users/diego/AppData/Local/Temp/Rtmpo1bQga/Belgium/nuts/units/BE32-region-60m-3035-2024.geojson"
#>  [9] "C:/Users/diego/AppData/Local/Temp/Rtmpo1bQga/Belgium/nuts/units/BE33-region-60m-3035-2024.geojson"
#> [10] "C:/Users/diego/AppData/Local/Temp/Rtmpo1bQga/Belgium/nuts/units/BE34-region-60m-3035-2024.geojson"
#> [11] "C:/Users/diego/AppData/Local/Temp/Rtmpo1bQga/Belgium/nuts/units/BE35-region-60m-3035-2024.geojson"

^{Created on 2025-12-08 with reprex v2.1.1}

Session info

sessioninfo::session_info()
#> Warning in system2("quarto", "-V", stdout = TRUE, env = paste0("TMPDIR=", : el
#> comando ejecutado '"quarto"
#> TMPDIR=C:/Users/diego/AppData/Local/Temp/Rtmpo1bQga/file3c7829811cb9 -V' tiene
#> el estatus 1
#> ─ Session info ───────────────────────────────────────────────────────────────
#>  setting  value
#>  version  R version 4.5.2 (2025-10-31 ucrt)
#>  os       Windows 11 x64 (build 26200)
#>  system   x86_64, mingw32
#>  ui       RTerm
#>  language (EN)
#>  collate  Spanish_Spain.utf8
#>  ctype    Spanish_Spain.utf8
#>  tz       Europe/Madrid
#>  date     2025-12-08
#>  pandoc   3.6.3 @ C:/Program Files/RStudio/resources/app/bin/quarto/bin/tools/ (via rmarkdown)
#>  quarto   NA @ C:\\Users\\diego\\AppData\\Local\\Programs\\Quarto\\bin\\quarto.exe
#> 
#> ─ Packages ───────────────────────────────────────────────────────────────────
#>  package     * version    date (UTC) lib source
#>  class         7.3-23     2025-01-01 [2] CRAN (R 4.5.2)
#>  classInt      0.4-11     2025-01-08 [1] CRAN (R 4.5.0)
#>  cli           3.6.5      2025-04-23 [1] CRAN (R 4.5.0)
#>  curl          7.0.0      2025-08-19 [1] CRAN (R 4.5.1)
#>  DBI           1.2.3      2024-06-02 [1] CRAN (R 4.5.0)
#>  digest        0.6.39     2025-11-19 [1] CRAN (R 4.5.2)
#>  e1071         1.7-16     2024-09-16 [1] CRAN (R 4.5.0)
#>  evaluate      1.0.5      2025-08-27 [1] CRAN (R 4.5.1)
#>  fastmap       1.2.0      2024-05-15 [1] CRAN (R 4.5.0)
#>  fs            1.6.6      2025-04-12 [1] CRAN (R 4.5.0)
#>  giscoR      * 0.9.9.9999 2025-12-08 [1] Github (ropengov/giscoR@3032faf)
#>  glue          1.8.0      2024-09-30 [1] CRAN (R 4.5.0)
#>  htmltools     0.5.8.1    2024-04-04 [1] CRAN (R 4.5.0)
#>  httr2         1.2.2      2025-12-08 [1] CRAN (R 4.5.2)
#>  jsonlite      2.0.0      2025-03-27 [1] CRAN (R 4.5.0)
#>  KernSmooth    2.23-26    2025-01-01 [2] CRAN (R 4.5.2)
#>  knitr         1.50       2025-03-16 [1] CRAN (R 4.5.0)
#>  lifecycle     1.0.4      2023-11-07 [1] CRAN (R 4.5.0)
#>  magrittr      2.0.4      2025-09-12 [1] CRAN (R 4.5.1)
#>  otel          0.2.0      2025-08-29 [1] CRAN (R 4.5.1)
#>  pillar        1.11.1     2025-09-17 [1] CRAN (R 4.5.1)
#>  pkgconfig     2.0.3      2019-09-22 [1] CRAN (R 4.5.0)
#>  proxy         0.4-27     2022-06-09 [1] CRAN (R 4.5.0)
#>  R6            2.6.1      2025-02-15 [1] CRAN (R 4.5.0)
#>  rappdirs      0.3.3      2021-01-31 [1] CRAN (R 4.5.0)
#>  Rcpp          1.1.0      2025-07-02 [1] CRAN (R 4.5.1)
#>  reprex        2.1.1      2024-07-06 [1] CRAN (R 4.5.0)
#>  rlang         1.1.6      2025-04-11 [1] CRAN (R 4.5.0)
#>  rmarkdown     2.30       2025-09-28 [1] CRAN (R 4.5.1)
#>  rstudioapi    0.17.1     2024-10-22 [1] CRAN (R 4.5.0)
#>  sessioninfo   1.2.3      2025-02-05 [1] CRAN (R 4.5.0)
#>  sf            1.0-23     2025-11-28 [1] CRAN (R 4.5.2)
#>  tibble        3.3.0      2025-06-08 [1] CRAN (R 4.5.0)
#>  units         1.0-0      2025-10-09 [1] CRAN (R 4.5.1)
#>  vctrs         0.6.5      2023-12-01 [1] CRAN (R 4.5.0)
#>  withr         3.0.2      2024-10-28 [1] CRAN (R 4.5.0)
#>  xfun          0.54       2025-10-30 [1] CRAN (R 4.5.2)
#>  yaml          2.3.11     2025-11-28 [1] CRAN (R 4.5.2)
#> 
#>  [1] C:/Users/diego/AppData/Local/R/win-library/4.5
#>  [2] C:/Program Files/R/R-4.5.2/library
#>  * ── Packages attached to the search path.
#> 
#> ──────────────────────────────────────────────────────────────────────────────

[dieghernan]

FYI gisco 1.0.0 has been submitted to CRAN, see #114