-
Notifications
You must be signed in to change notification settings - Fork 2
/
create_graph.py
51 lines (37 loc) · 1.59 KB
/
create_graph.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
import geopandas
import logging
from graphmaker.graph import Graph
from clean_up import correct_islands
from main import plans
log = logging.getLogger(__name__)
log.addHandler(logging.StreamHandler())
log.setLevel(logging.INFO)
columns = ['rounded11', '2011', 'T16PRESR', 'T16PRESD',
'T16SENR', 'T16SEND', 'population', 'GOV_4_1',
'TS_4_1', 'Remedial']
def create_graph_from_shapefile(shapefile_path='./data/wes_with_districtings/wes_with_districtings.shp'):
log.info(f"Creating the queen adjacency graph")
pa_queen = Graph.from_shapefile(
shapefile_path, adjacency_type='queen', data_columns=None, id_column='wes_id')
log.info(f"Saving the queen adjacency graph")
pa_queen.add_columns_from_shapefile(
shapefile_path, columns, id_column='wes_id')
# Fix nodes whose assignment is different from all their neighbors, and
# whose neighbors all have the same assignment.
# In practice, this affects one node's assignment in two plans. (Node '1')
for plan in plans:
correct_islands(pa_queen.graph, plan)
pa_queen.save('./PA_queen.json')
log.info(f"Creating the rook adjacency graph")
pa_rook = Graph.from_shapefile(
shapefile_path, adjacency_type='rook', data_columns=None, id_column='wes_id')
pa_rook.add_columns_from_shapefile(
shapefile_path, columns, id_column='wes_id')
log.info(f"Saving the queen adjacency graph")
for plan in plans:
correct_islands(pa_rook.graph, plan)
pa_rook.save('./PA_rook.json')
def main():
create_graph_from_shapefile()
if __name__ == '__main__':
main()