Create shapefile from OpenAQ API V2 locations endpoint

main.py

from urllib import request
import requests
import fiona
from fiona.crs import from_epsg


def result_to_dict(result):
    geom = {
        "type": "Point",
        "coordinates": [
            float(result["coordinates"]["latitude"]),
            float(result["coordinates"]["longitude"]),
        ],
    }
    return {"type": "Feature", "geometry": geom, "properties": {"Name": result["name"]}}


def main():
    url = "https://api.openaq.org/v2/locations?limit=5000&page=1&offset=0&sort=desc&coordinates=34.0522%2C-118.2437&radius=40000&order_by=lastUpdated&dumpRaw=false"
    res = requests.get(url)
    data = res.json()
    records = [result_to_dict(result) for result in data["results"]]
    schema = {"geometry": "Point", "properties": {"Name": "str:48"}}
    with fiona.open(
        "output.shp", "w", crs=from_epsg(4326), driver="ESRI Shapefile", schema=schema
    ) as output:
        output.writerecords(records)


if __name__ == "__main__":
    main()