def setup(self):
test_file_path = mm.datasets.get_path("bubenec")
self.df_buildings = gpd.read_file(test_file_path, layer="buildings")
self.df_streets = gpd.read_file(test_file_path, layer="streets")
self.df_tessellation = gpd.read_file(test_file_path,
layer="tessellation")
self.df_streets["nID"] = mm.unique_id(self.df_streets)
self.df_buildings["height"] = np.linspace(10.0, 30.0, 144)
self.df_tessellation["area"] = self.df_tessellation.geometry.area
self.df_buildings["area"] = self.df_buildings.geometry.area
self.df_buildings["fl_area"] = mm.FloorArea(self.df_buildings,
"height").series
self.df_buildings["nID"] = mm.get_network_id(self.df_buildings,
self.df_streets, "nID")
blocks = mm.Blocks(self.df_tessellation, self.df_streets,
self.df_buildings, "bID", "uID")
self.blocks = blocks.blocks
self.df_buildings["bID"] = blocks.buildings_id
self.df_tessellation["bID"] = blocks.tessellation_id
self.swb = Queen.from_dataframe(self.df_buildings)
self.sw5 = mm.sw_high(k=5, gdf=self.df_tessellation, ids="uID")
self.sw3 = mm.sw_high(k=3, gdf=self.df_tessellation, ids="uID")
self.sws = mm.sw_high(k=2, gdf=self.df_streets)
nx = mm.gdf_to_nx(self.df_streets)
nx = mm.node_degree(nx)
self.nodes, self.edges, W = mm.nx_to_gdf(nx, spatial_weights=True)
self.swn = mm.sw_high(k=3, weights=W)
def test_get_node_id(self):
nx = mm.gdf_to_nx(self.df_streets)
nodes, edges = mm.nx_to_gdf(nx)
self.df_buildings["nID"] = mm.get_network_id(self.df_buildings,
self.df_streets, "nID")
ids = mm.get_node_id(self.df_buildings, nodes, edges, "nodeID", "nID")
assert not ids.isna().any()
convex_hull = edges.unary_union.convex_hull
enclosures = mm.enclosures(edges, limit=gpd.GeoSeries([convex_hull]))
enclosed_tess = mm.Tessellation(self.df_buildings,
unique_id="uID",
enclosures=enclosures).tessellation
links = mm.get_network_ratio(enclosed_tess, edges)
enclosed_tess[links.columns] = links
ids = mm.get_node_id(
enclosed_tess,
nodes,
edges,
node_id="nodeID",
edge_keys="edgeID_keys",
edge_values="edgeID_values",
)
assert not ids.isna().any()
def test_get_node_id(self):
nx = mm.gdf_to_nx(self.df_streets)
nodes, edges = mm.nx_to_gdf(nx)
self.df_buildings["nID"] = mm.get_network_id(
self.df_buildings, self.df_streets, "nID"
)
ids = mm.get_node_id(self.df_buildings, nodes, edges, "nodeID", "nID")
assert not ids.isna().any()
def setup_method(self):
test_file_path = mm.datasets.get_path("bubenec")
self.df_buildings = gpd.read_file(test_file_path, layer="buildings")
self.df_streets = gpd.read_file(test_file_path, layer="streets")
self.df_tessellation = gpd.read_file(test_file_path, layer="tessellation")
self.df_buildings["height"] = np.linspace(10.0, 30.0, 144)
self.df_buildings["volume"] = mm.Volume(self.df_buildings, "height").series
self.df_streets["nID"] = mm.unique_id(self.df_streets)
self.df_buildings["nID"] = mm.get_network_id(
self.df_buildings, self.df_streets, "nID"
)
def setup_method(self):
test_file_path = mm.datasets.get_path('bubenec')
self.df_buildings = gpd.read_file(test_file_path, layer='buildings')
self.df_streets = gpd.read_file(test_file_path, layer='streets')
self.df_tessellation = gpd.read_file(test_file_path,
layer='tessellation')
self.df_buildings['height'] = np.linspace(10., 30., 144)
self.df_buildings['volume'] = mm.volume(self.df_buildings, 'height')
self.df_streets['nID'] = mm.unique_id(self.df_streets)
self.df_buildings['nID'] = mm.get_network_id(self.df_buildings,
self.df_streets, 'uID',
'nID')
def setup(self):
test_file_path = mm.datasets.get_path("bubenec")
self.df_buildings = gpd.read_file(test_file_path, layer="buildings")
self.df_tessellation = gpd.read_file(test_file_path, layer="tessellation")
self.df_streets = gpd.read_file(test_file_path, layer="streets")
self.df_streets["nID"] = range(len(self.df_streets))
self.limit = mm.buffered_limit(self.df_buildings, 50)
nx = mm.gdf_to_nx(self.df_streets)
self.nodes, self.edges = mm.nx_to_gdf(nx)
self.df_buildings["nID"] = mm.get_network_id(
self.df_buildings, self.df_streets, "nID"
)
def setup(self):
test_file_path = mm.datasets.get_path("bubenec")
self.df_buildings = gpd.read_file(test_file_path, layer="buildings")
self.df_streets = gpd.read_file(test_file_path, layer="streets")
self.df_tessellation = gpd.read_file(test_file_path, layer="tessellation")
self.df_buildings["height"] = np.linspace(10.0, 30.0, 144)
self.df_buildings["volume"] = mm.Volume(self.df_buildings, "height").series
self.df_streets["nID"] = mm.unique_id(self.df_streets)
self.df_buildings["nID"] = mm.get_network_id(
self.df_buildings, self.df_streets, "nID"
)
self.df_buildings["orient"] = mm.Orientation(self.df_buildings).series
self.df_tessellation["orient"] = mm.Orientation(self.df_tessellation).series
self.sw = Queen.from_dataframe(self.df_tessellation, ids="uID")
self.swh = mm.sw_high(k=3, gdf=self.df_tessellation, ids="uID")
self.swb = Queen.from_dataframe(self.df_buildings, ids="uID")
def setup_method(self):
test_file_path = mm.datasets.get_path('bubenec')
self.df_buildings = gpd.read_file(test_file_path, layer='buildings')
self.df_streets = gpd.read_file(test_file_path, layer='streets')
self.df_tessellation = gpd.read_file(test_file_path,
layer='tessellation')
self.df_streets['nID'] = mm.unique_id(self.df_streets)
self.df_buildings['height'] = np.linspace(10., 30., 144)
self.df_tessellation['area'] = self.df_tessellation.geometry.area
self.df_buildings['area'] = self.df_buildings.geometry.area
self.df_buildings['fl_area'] = mm.floor_area(self.df_buildings,
'height')
self.df_buildings['nID'] = mm.get_network_id(self.df_buildings,
self.df_streets, 'uID',
'nID')
self.df_buildings, self.df_tessellation, self.blocks = mm.blocks(
self.df_tessellation, self.df_streets, self.df_buildings, 'bID',
'uID')
def setup_method(self):
test_file_path = mm.datasets.get_path("bubenec")
self.df_buildings = gpd.read_file(test_file_path, layer="buildings")
self.df_streets = gpd.read_file(test_file_path, layer="streets")
self.df_tessellation = gpd.read_file(test_file_path, layer="tessellation")
self.df_streets["nID"] = mm.unique_id(self.df_streets)
self.df_buildings["height"] = np.linspace(10.0, 30.0, 144)
self.df_tessellation["area"] = self.df_tessellation.geometry.area
self.df_buildings["area"] = self.df_buildings.geometry.area
self.df_buildings["fl_area"] = mm.FloorArea(self.df_buildings, "height").series
self.df_buildings["nID"] = mm.get_network_id(
self.df_buildings, self.df_streets, "nID"
)
blocks = mm.Blocks(
self.df_tessellation, self.df_streets, self.df_buildings, "bID", "uID"
)
self.blocks = blocks.blocks
self.df_buildings["bID"] = blocks.buildings_id
self.df_tessellation["bID"] = blocks.tessellation_id
def test_get_network_id_duplicate(self):
self.df_buildings["nID"] = range(len(self.df_buildings))
buildings_id = mm.get_network_id(self.df_buildings, self.df_streets,
"nID")
assert not buildings_id.isna().any()
def test_get_network_id(self):
buildings_id = mm.get_network_id(self.df_buildings, self.df_streets,
"nID")
assert not buildings_id.isna().any()
def time_get_network_id(self):
mm.get_network_id(self.df_buildings, self.df_streets, "nID")
import geopandas as gpd
import momepy as mm
tess = gpd.read_file('files/prg.gpkg', layer='queen')
blg = gpd.read_file('files/prg.gpkg', layer='buildings')
streets = gpd.read_file('files/streets.gpkg', layer='streets')
streets['nID'] = mm.unique_id(streets)
blg['nID'] = mm.get_network_id(blg, streets, 'uID', 'nID', 150)
tess = tess.merge(blg[['uID', 'nID']], on='uID', how='left')
snapped = mm.snap_street_network_edge(streets, blg, 20, tess, 120)
blocks, blg['bID'], tess['bID'] = mm.blocks(tess, snapped, blg, 'bID', 'uID')
tess.to_file('files/elements.gpkg', layer='tessellation', driver='GPKG')
blg.to_file('files/elements.gpkg', layer='buildings', driver='GPKG')
blocks.to_file('files/elements.gpkg', layer='blocks', driver='GPKG')
streets.to_file('files/elements.gpkg', layer='streets', driver='GPKG')
