crowdedness_create_network

crowdedness_create_network.py

network = trips.drop_duplicates('route_id').dropna()
network = network[~(network.trip_long_name == 'Stopbus i.p.v. trein')]

features=[]
for si in network.shape_id.unique():
    section = shapes[shapes.shape_id == si]
    ls =LineString(section[['shape_pt_lon','shape_pt_lat']].to_numpy().tolist())
    features.append(
        Feature(geometry=ls, properties=dict(label=str(si)))
    )

with open('network.geojson', 'w') as outfile:
        dump(FeatureCollection(features), outfile)