2

I have a dataset of around 15000 events (protests) with latitude+longitude coordinates, and I would like to know in which US congressional district each event occurred. I've found a shapefile of US congressional districts, and would like to use it to identify the congr. district for each individual point.

I understand I can use the "Join attributes to location" function in QGIS to do this, but I don't understand how to use it.

Do both files (my dataset file with the events and the congr. district file) need to be shapefiles (.shp)?

Is using this basically going to add another column to my dataset file with info about which district each point is located in?

Never used QGIS or any other geographic program before, so please be specific if possible.

The dataset can be found here: https://acleddata.com/download/22846/

District shapefiles here: http://cdmaps.polisci.ucla.edu/shp/districts114.zip

1
  • 2
    What have you tried? Have you read up on how join attributes by location works?
    – Erik
    Commented Nov 9, 2020 at 9:02

1 Answer 1

4

The short answer: a) both files can be in different format (one as csv, one as shapefile) b) as you can read in the description of the algorithm (in the dialogue window when you start "Join attributes to location"), it creates a new layer, based on your points (protests), adding the attributes of the district the point lies in.

See the following picture (even if its in german, it should become clear how to use).

enter image description here

The longer answer: there are some problems with your data (as often when downloading data, you often have to adjust it for proper use).

First, the protests are in an excel file (xlsx). To import it in QGIS, you need it as a CSV file. Open it in Excel. However, there are a lot of semicolons (;) around in the data, you have to remove them as they are also used as text-delimiter character of the csv. So before exporting your file, in Excel search and replace the semicolon with another character, like hyphen (-). Than save as... and choose csv, best use the CSV UTF-8 option. Now you can import this csv in QGIS: again, see screenshot that should illustrate the options, even if it's in german. Be sure to choose the semicolon as text delimiter and make sure the other characters are not checked. Geometry definition: point coordinates, Geometry CRS EPSG:4326, make sure the x and y-fields are correctly interpreted. Than add to your project.

enter image description here

Loading the shapefile with districts should not be a problem.

However, when using the "Join attributes to location", there is an error message, as there are invalid geometries. So first, you have to repair the geometries of your shapefile. Use the "repair geometries" from the processing toolbox and select your district-layer as input. The algorithm will create a new layer with the repaired geometries. Use this one for processing with "Join attributes to location".

This takes quite a while, on my machine (QGIS 3.16, Win 10) it took 6:42 min. Better first create a spatial index for your vector layer: right click on the layer - properties - source tab - under the entry Geometry: create spatial index. With this, "Join attributes to location" took only 26 seconds!

enter image description here

1
  • Thank you sooo much for that answer, it worked like a charm! Commented Nov 9, 2020 at 19:53

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.