Zonal statistics using dask_geomodeling

dask_geomodeling_zonal_stats.py

# Import modules 
from dask_geomodeling.raster import RasterBlock
from dask_geomodeling.raster import base
from dask_geomodeling.geometry import base

# Define paths 
arg_rast_path = 'raster.tif'
arg_vec_path = 'vector.shp'

"""
# How do we define the blocks? 
arg_ras = base.RasterBlock(arg_rast_path)
arg_vec = base.GeometryBlock(arg_vec_path)
"""

# Apply zonal statistics 
geometry.aggregate.AggregateRaster(source=arg_vec,raster=arg_rast, column_name='sum', statistic='sum')