Querying amenities in Greater Victoria using OSM

fetch_amenities_osm.py

import osmnx as ox
import matplotlib.pyplot as plt

from pathlib import Path

AMENITY_QUERIES = [
    ("bank", {"amenity": "bank"}),
    ("bar", {"amenity": ["bar", "pub"]}),
    ("coffee shop", {"amenity": "cafe", "cuisine": "coffee_shop"}),
    (
        "greater victoria public library",
        {"operator": "Greater Victoria Public Library"},
    ),
    ("medical clinic", {"amenity": "clinic"}),
    ("school", {"amenity": "school"}),
    (
        "store",
        {
            "shop": ["supermarket", "greengrocer", "convenience"],
            "name": ["Walmart", "Wholesale Club"],
        },
    ),
]
COL_REGEX = r"^(name|brand|geometry|(addr:\w+))$"

if __name__ == "__main__":
    data_dir = Path("amenity_data_gva")
    data_dir.mkdir(exist_ok=True)

    for amenity, query in AMENITY_QUERIES:
        places = ox.geometries_from_place("Greater Victoria", query)
        places = places.filter(regex=COL_REGEX)
        places.to_file(data_dir / f"{amenity}.geojson", driver="GeoJSON")