Simplification code

[martinfleis]

This is a snippet that can be used to simplify the network. It requires sgeop development version available from https://github.com/uscuni/sgeop and this pysal/momepy#652.

import geopandas
import sgeop

roads = geopandas.read_parquet("/data/uscuni-ulce/processed_data/streets/streets_69300.parquet")
buildings = geopandas.read_parquet(
    "/data/uscuni-ulce/processed_data/buildings/buildings_69300.parquet", columns=["geometry"]
)

new_roads = sgeop.simplify_network(
    roads,
    exclusion_mask=buildings.geometry,
    artifact_threshold_fallback=7,
)