Parallel code refinement

[gfdrr]

Parallel code refinement

[matamadio]

• The zonal step seem to introduce an overestimation error ONLY for total Pop; when tested on DOM against original code and QGIS zonal, the error is between 1% and 50%. The raster creation is checked to be correct according to exported byproducts.

     with mp.Pool(cores) as p:
        # Get total exposure for each ADM region
        func = partial(zonal_stats_partial, raster=exp_ras, stats="sum")
        stats_parallel = p.map(func, np.array_split(adm_data.geometry, cores))

    exp_per_ADM = list(it.chain(*stats_parallel))

[matamadio]

Thanks Arthur for refining the code.