Hello Mathias,
thank you for your support.
I was able to solve it.
I fortified the Polygons and added the frequency to the dataframe.
Best regards.
On 24/06/2021 08:33, Mathias Moser wrote:
You need to create a frequency table from l and merge this to the shapefile.
Using sf/ggplot2 would also simplify things here:
download.file("
https://gisco-services.ec.europa.eu/distribution/v2/nuts/download/ref-nuts-2021-60m.geojson.zip
", "nuts.zip")
unzip("nuts.zip", exdir = tempdir())
library(dplyr)
library(sf)
EU_NUTS <- read_sf(paste0(tempdir(), "/NUTS_RG_60M_2021_4326_LEVL_3.geojson")
l <- c("AT223", "AT212", "AT212", "AT121")
l_freq <- as.data.frame(table(l))
EU_NUTS <- left_join(EU_NUTS, l_freq, by = c("id" = "l"))
library(ggplot2)
ggplot(filter(EU_NUTS, CNTR_CODE == "AT")) + geom_sf(aes(fill = Freq))
//m
On Wed, 2021-06-23 at 21:11 +0200, Phillips Rogfield wrote:
Dear R-SIG-Geo ML,
I have downloaded the NUTS shapefiles from here:
https://ec.europa.eu/eurostat/web/gisco/geodata/reference-data/administrative-units-statistical-units/nuts
I can load them with the following code:
library(rgdal)
shp_bdir <- [PATH TO SHAPE FILE]
layername <- "NUTS_RG_01M_2021_4326"
shp_folder <- file.path(shp_bdir, paste0(layername,".shp"))
EU_NUTS <- readOGR(dsn = shp_folder, layer = layername)
Then I can plot the regions with:
plot(EU_NUTS)
Now I have a list of NUTS - 3 Level, for example:
l <- c("AT223", "AT212", "AT212", "AT121")
I would like a plot where the NUTS regions are colored, and the more a
particular NUTS is present is the list, the darker its color.
For example "AT212" should be darker than "AT223", because the former is
present two times in the list.
How can I achieve that?
Thank you and best regards.
[[alternative HTML version deleted]]
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
