Import the ESRI shapefile of German districts and the district attribute table. Join the two data frames, transform the CRS to “EPSG:3035” and check your changes.
You need to rename one of the id variables or adjust your join accordingly (“id = district_id”).
We want a first descriptive visual of the distribution of Covid-19 cases in Cologne and the surrounding districts. Calculate the number of Covid-19 cases in the last seven days (cases_7days) by population (population) and multiply with 100,000 (in Germany usually called “7 Tages Inzidenzzahl”: number of persons who were infected with Covid-19 in the last seven days per 100k).
Select Cologne (district_id == 5315), find the surrounding districts, and plot Cologne and its surrounding districts.
You can use the dplyr function bind_rows to combine the two spatial objects, “Cologne” and “Cologne Surroundings”.
Save your data set of Cologne and its surrounding districts as an ESRI Shapefile.