def __init__(self, shapefile=Path("assets/wrs/WRS2_descending.shp")): speedups.disable() self.shapefile = ogr.Open(str(shapefile)) self.layer = self.shapefile.GetLayer(0) self.polygons = [] for i in range(self.layer.GetFeatureCount()): feature = self.layer.GetFeature(i) path = feature['PATH'] row = feature['ROW'] geom = feature.GetGeometryRef() shape = shapely.wkt.loads(geom.ExportToWkt()) self.polygons.append((shape, path, row))
# Here I apply a global model on bus OD data in Singapore. The bus OD data is a small subset of bus trip data available via LTA Data Mall: https://datamall.lta.gov.sg/content/datamall/en.html. Preprocessing has been performed by matching O/D to bus stop locations, and line features connecting each pair of OD are created. Weight column carries the expected volume of passengers on the pair of OD nodes. # # %% [markdown] # ## Global Model # %% import geopandas as gpd import networkx as nx import pandas as pd import numpy as np import matplotlib.pyplot as plt import os, sys # %% from shapely import speedups speedups.disable() # %% # import line shapefile (representing OD flows between each pair of bus stops) orig_flow = gpd.read_file(r"data\SampleODFlow.shp") orig_flow.head() # %% orig_flow.info() # %% orig_flow.plot(linewidth=0.5, alpha=0.01) # %% [markdown] # Some of the nodes are in erroneous locations (X/Y of 0,0). We will remove these OD links from the dataset.
def tearDown(self): self.assertTrue(speedups._orig) speedups.disable() self.assertFalse(speedups._orig)
def tearDown(self): self.assertTrue(speedups._orig) speedups.disable() self.assertFalse(speedups._orig)