[R-sig-Geo] Plot NUTS regions with color

[Phillips Rogfield]

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 using r-project.org
>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo using r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo