Hoe kan ik de exacte postcodes zien gekoppeld aan wijken/gebieden.
Ik heb aanstaan CBS Wijken en Buurten 2022 versie
en CBS postcode 4 statistieken WMS 2020
PDOK heeft geen dataproduct waarin postcodes gekoppeld zijn aan wijken/buurten. Dus als je dit wilt zul je het zelf moeten maken.
Dit is relatief eenvoudig met een SQL query te bereiken in combinatie met de GeoPackage downloads van Wijken en Buurten en CBS Postcode 4. Zie het bash script hieronder:
# downloaden geopackages via atom services met wget
wget "https://service.pdok.nl/cbs/wijkenbuurten/2022/atom/downloads/wijkenbuurten_2022_v1.gpkg"
wget "https://service.pdok.nl/cbs/pc4/atom/v1_0/downloads/cbs_pc4_2020.gpkg.zip"
unzip cbs_pc4_2020.gpkg.zip
tmp_gpkg=$(mktemp -u --suffix .gpkg) # maak tijdelijk gpkg bestand aan
ogr2ogr -f GPKG "$tmp_gpkg" wijkenbuurten_2022_v1.gpkg wijken # sla wijken en cbs_pc4_2020 tabellen op in tijdelijk bestand
ogr2ogr -append -f GPKG "$tmp_gpkg" cbs_pc4_2020.gpkg.zip cbs_pc4_2020
# voeg postcodes toe aan wijken laag dmv spatial intersects join
ogr2ogr -f GPKG "wijkenbuurten-pc4-join.gpkg" "$tmp_gpkg" -nln wijken -sql "SELECT
a.*,
GROUP_CONCAT(b.postcode) AS postcodes
FROM
wijken AS a,
cbs_pc4_2020 AS b
WHERE
ST_INTERSECTS(a.geom, b.geom)
GROUP BY
a.wijkcode"
Resultaat van bovenstaande code kan je hier downloaden: wijkenbuurten-pc4-join.gpkg
Ik heb ook nog geprobeerd bovenstaande query te gebruiken in een QGIS virtual layer, zoals hier beschreven. Maar kreeg dit niet werkend, QGIS leek daar op vast te lopen.
