Read in CSV data and rasterize it

rasterize.py

from functools import partial
import json

from geocube.api.core import make_geocube
from geocube.rasterize import rasterize_image
from rasterio.enums import MergeAlg
from shapely.wkt import loads
from shapely.geometry import box, mapping, MultiPoint
import pandas as pd
import geopandas as gpd


df = pd.read_csv(
    "testdata.csv",
    index_col=False,
    usecols=[3, 8],
    converters={8: loads},
)
gdf = gpd.GeoDataFrame(
    df,
    crs="epsg:4326",
    geometry=df.geometry,
)
lcc = gdf.to_crs("epsg:8763")
grid = make_geocube(
    vector_data=lcc,
    measurements=["speed"],
    resolution=(-12000, 12000),
    rasterize_function=partial(rasterize_image, merge_alg=MergeAlg.add),
    fill=0,
    align=(9000, 8000),
)

testdata.csv