Contents

Regional aggregates for Cambodia

Regional aggregates for Cambodia

import pandas, geopandas
import contextily
import xarray, rioxarray
from rasterstats import zonal_stats

Boundaries

Important

Please see here for more details

url = (
    "https://geonode.wfp.org/geoserver/wfs"\
    "?format_options=charset:UTF-8&"\
    "typename=geonode:khm_adm2_un&"\
    "outputFormat=SHAPE-ZIP&"\
    "version=1.0.0&"\
    "service=WFS&"\
    "request=GetFeature"
)
url

'https://geonode.wfp.org/geoserver/wfs?format_options=charset:UTF-8&typename=geonode:khm_adm2_un&outputFormat=SHAPE-ZIP&version=1.0.0&service=WFS&request=GetFeature'

cam = geopandas.read_file(url)

bb = cam.total_bounds

Regional friction

  • Motorised

agg_m = pandas.DataFrame(
    zonal_stats(
        cam,
        "cambodia_2020_motorized_friction_surface.tif"
    ),
    index = cam.index
)

  • Walking

agg_w = pandas.DataFrame(
    zonal_stats(
        cam,
        "cambodia_2020_walking_friction_surface.tif"
    ),
    index = cam.index
)

Regional pollution

agg_p = pandas.DataFrame(
    zonal_stats(
        cam,
        "cambodia_s5_no2.tif",
        all_touched = True
    ),
    index = cam.index
)

Join together

db = cam[[
    "adm2_name", "adm2_altnm", "geometry"
]].join(
    agg_m[["mean"]].rename(columns={"mean": "motor_mean"})
).join(
    agg_w[["mean"]].rename(columns={"mean": "walk_mean"})
).join(
    agg_p[["mean"]].rename(columns={"mean": "no2_mean"})
)
db.info()

<class 'geopandas.geodataframe.GeoDataFrame'>
RangeIndex: 198 entries, 0 to 197
Data columns (total 6 columns):
 #   Column      Non-Null Count  Dtype   
---  ------      --------------  -----   
 0   adm2_name   198 non-null    object  
 1   adm2_altnm  122 non-null    object  
 2   geometry    198 non-null    geometry
 3   motor_mean  198 non-null    float64 
 4   walk_mean   198 non-null    float64 
 5   no2_mean    198 non-null    float64 
dtypes: float64(3), geometry(1), object(2)
memory usage: 9.4+ KB

Write out

out_gj = "cambodia_regional.gpkg"
! rm -f $out_gj
db.to_file(out_gj, driver="GPKG")
! du -h $out_gj

2.8M	cambodia_regional.gpkg