Geocube

from functools import partial

from geocube.api.core import make_geocube

from geocube.rasterize import rasterize_image, rasterize_points_griddata

out_grid = make_geocube(

    vector_data=gdf,

    measurements=["COL"],

    resolution=(-1, 1), # -1 will invert lat

    rasterize_function=partial(rasterize_image, all_touched=True),

)

out_grid = make_geocube(

    vector_data=gdf,

    measurements=["COL", "datetime"],

    resolution=(-1, 1), # -1 will invert lat

    group_by="datetime",

)