20km UTM tile splitter

[jdries]

Example grid: https://artifactory.vgt.vito.be/ui/repos/tree/General/auxdata-public/grids/utm-tiling-grid-20km.gpkg

Define a function like this:

def split_area(aoi: Union[bounds, shapely.Geometry], projection="AUTO:42001", tile_size:float=20.0)

That simply returns a list of geometries that define grid cells, coordinates should be in projection.

Then we want to make it trivial to construct a job_db based on this, so wrap this code in a function:

df = gpd.GeoDataFrame({
    "geometry": split_area(aoi),
})
job_db = create_job_db("jobs.csv").initialize_from_df(df)