Geofabrik OpenStreetMap Extracts

python -m pip install requests osmium

from pathlib import Path
import osmium
import requests

url = "https://download.geofabrik.de/north-america/canada/prince-edward-island-latest.osm.pbf"
path = Path("prince-edward-island-latest.osm.pbf")
with requests.get(url, stream=True, timeout=60) as response:
    response.raise_for_status()
    with path.open("wb") as output:
        for chunk in response.iter_content(chunk_size=1024 * 1024):
            output.write(chunk)

class RoadCounter(osmium.SimpleHandler):
    def __init__(self):
        super().__init__()
        self.count = 0

    def way(self, way):
        if "highway" in way.tags:
            self.count += 1

roads = RoadCounter()
roads.apply_file(str(path))
print(f"{roads.count} mapped road ways in this dated extract")