def METHOD_NAME(
    base,
    population,
    method,
    output_demand,
    vector_base_map=None,
    vector_methods=None,
    vector_column_prefix=None,
):
    """ """
    r.stats_zonal(
        base=base,
        flags="r",
        cover=population,
        method=method,
        output=output_demand,
        overwrite=True,
        quiet=True,
    )
    # copy 'reclassed' as 'normal' map (r.mapcalc)
    # so as to enable removal of it and its 'base' map
    demand_copy = output_demand + "_copy"
    copy_expression = "{input_raster}"
    copy_expression = copy_expression.format(input_raster=output_demand)
    copy_equation = EQUATION.format(result=demand_copy, expression=copy_expression)
    r.mapcalc(copy_equation, overwrite=True)
    # remove reclassed map 'output_demand'
    # and rename 'copy' back to 'output_demand'
    g.remove(flags="f", type="raster", name=output_demand, quiet=True)
    g.rename(raster=(demand_copy, output_demand), quiet=True)
    if output_demand and vector_base_map:
        update_vector(
            vector=vector_base_map,
            raster=output_demand,
            methods=vector_methods,
            column_prefix=vector_column_prefix,
        )