script-climate-abfiremap.py

import requests
import geopandas
import json
import datawrappergraphics

# Test chart: "HqkeQ"
# Live chart: "L45df"
CHART_ID = "L45df"

r = requests.get("https://services.arcgis.com/Eb8P5h4CJk8utIBz/ArcGIS/rest/services/Active_Wildfire_Locations/FeatureServer/0/query?where=1%3D1&objectIds=&time=&geometry=&geometryType=esriGeometryEnvelope&inSR=&spatialRel=esriSpatialRelIntersects&resultType=standard&distance=0.0&units=esriSRUnit_Meter&returnGeodetic=false&outFields=*&returnGeometry=true&featureEncoding=esriDefault&multipatchOption=xyFootprint&maxAllowableOffset=&geometryPrecision=&outSR=&defaultSR=&datumTransformation=&applyVCSProjection=false&returnIdsOnly=false&returnUniqueIdsOnly=false&returnCountOnly=false&returnExtentOnly=false&returnQueryGeometry=false&returnDistinctValues=false&cacheHint=false&orderByFields=&groupByFieldsForStatistics=&outStatistics=&having=&resultOffset=&resultRecordCount=&returnZ=false&returnM=false&returnExceededLimitFeatures=true&quantizationParameters=&sqlFormat=none&f=pjson&token=")

data = geopandas.read_file(json.dumps(r.json()))

data["opacity"] = 0.5
data = data.drop(columns="ID")

data["markerColor"] = data["FIRE_STATUS"].replace({"Under Control": "#436170",
                                                   "New": "#F8C325",
                                                   "Out of Control": "#c42127",
                                                   "Being Held": "#000000",
                                                   "Assistance Started": "#1F78B4"
                                                   })

data["type"] = "point"
data["icon"] = "fire"

avg = data["AREA_ESTIMATE"].min()
std = data["AREA_ESTIMATE"].std()

data["scale"] = ((data["AREA_ESTIMATE"] - avg) / (std)) + 1
data["scale"] = data["scale"].apply(lambda x: 2.2 if x > 2.2 else x)

data["tooltip"] = "<big>" + data['RESP_AREA'] + "</big><br><b>Status</b>: " + data['FIRE_STATUS'] + "</span>" + "<br><b>Cause</b>: " + data["GENERAL_CAUSE"] 
data = data.sort_values("scale")

percent_under_control = round(len(data[data['FIRE_STATUS'] == 'Under Control'])/len(data)*100, 0)

chart = (datawrappergraphics.Map(chart_id=CHART_ID)
            .data(data, "./assets/shapes/shapes-abfiremap.json")
            .head(f"There are <b>{len(data)} wildfires</b> burning across Alberta")
            .deck(f"As of today, {percent_under_control}% are listed as under control.")
            .footer(source="Government of Alberta")
            .publish()
            )