1

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 (AGS = district_id).

2

We want a first descriptive visual of the distribution of charging stations in Cologne (or any other district of your choice) and the surrounding districts. Filter the district of Cologne (district_id == "05315") and find the surrounding districts. Calculate the number of chargers per district (charger_count) and the number of chargers per 1,000 inhabitants in each district (charger_dens). Plot the two columns for Cologne and its surrounding districts.

You can use the dplyr function sf::bind_rows() to combine the two spatial objects, “Cologne” and “Cologne Surroundings”.

3

Save your data set of Cologne and its surrounding districts as an ESRI Shapefile.