Ik heb het script aangepast zodat er alleen losse polygonen in het CSV bestand komen. Hiervoor heb ik een extra ogr2ogr stap toegevoegd, om de MULTIPOLYGON geoms te converteren naar POLYGON en ook de geometrie te versimpelen. Door het versimpelen van geometrieen kunnen juist MULTIPOLYGON geometrieen ontstaan, vandaar dat het niet in de (Spatialite) SQL query kan plaatsvinden.
#!/usr/bin/env bash
wget https://service.pdok.nl/cbs/pc4/atom/v1_0/downloads/cbs_pc4_2020.gpkg.zip # zie https://service.pdok.nl/cbs/pc4/atom/v1_0/cbs_pc4_2020.xml
unzip cbs_pc4_2020.gpkg.zip
ogr2ogr -f GPKG cbs_pc4_2020_ex.gpkg cbs_pc4_2020.gpkg -explodecollections -unsetfid -simplify 10 -nlt POLYGON
sqlite3 cbs_pc4_2020_ex.gpkg <<< "
select load_extension('mod_spatialite');
select EnableGpkgMode();
.mode csv
.output test.csv
select postcode,ST_AsText(geom) from cbs_pc4_2020 LIMIT 2;
"
De CSV ziet er nu zo uit:
1184,"POLYGON((121294.892023 478995.29298, 121372.789003 478948.788616, 121429.059332 478936.241411, 121442.786583 478952.842007, 121504.12951 478961.674802, 121600.413749 479005.099134, 121604.488114 479046.855796, 121556.749083 479091.014994, 121488.688699 479124.183289, 121338.962136 479064.383471, 121294.904969 478995.401674, 121147.249231 479085.147127, 121376.099415 479475.123328, 121422.41304 479441.938245, 121462.057772 479478.607639, 121467.026234 479489.303799, 121399.379222 479513.445478, 121446.07658 479594.473451, 121592.710211 479469.068402, 121642.827442 479527.886115, 121618.75783 479557.066357, 121643.170338 479588.420175, 121678.811635 479599.386593, 121664.612157 479616.229207, 121682.217383 479644.848649, 121615.291707 479672.07172, 121600.321737 479709.968354, 121558.860898 479717.418301, 121533.794129 479743.511153, 121651.695065 479948.402046, 121646.838856 479971.110255, 121613.186786 479984.831417, 121674.687111 480167.228264, 121815.580004 480016.999514, 122021.102095 480011.259798, 122225.478064 479839.76493, 121999.53344 479582.35541, 121740.865896 479334.804586, 121717.055031 479355.136851, 121651.090635 479259.919965, 121635.382971 479142.739043, 121677.37997 479047.940422, 121707.36187 479070.686478, 121838.957114 478901.570447, 121850.917067 478846.423806, 121833.496013 478788.31885, 121766.783389 478744.905627, 121276.649687 478695.899025, 121222.900725 478665.210634, 121165.700245 478595.042692, 121141.488505 478512.330092, 121153.735983 478343.955605, 121140.689809 478283.246167, 120937.002243 478170.314984, 120887.190088 478112.766844, 120874.466652 478071.463048, 120878.507609 478033.109585, 120910.780915 477997.902518, 121049.935684 477910.083705, 121056.990486 477837.698413, 120984.25368 477773.006907, 120813.134007 477744.96475, 120766.780623 477705.238914, 120734.851846 477645.57256, 120718.383552 477525.904931, 120688.121569 477445.068219, 120340.430775 476984.035877, 120303.315343 476888.021611, 120310.193839 476804.208688, 120282.472089 476796.845147, 120269.271505 476900.258206, 120306.663788 476989.784326, 120399.010229 477127.331835, 120617.176311 477391.979755, 120677.482445 477499.756884, 120703.363522 477628.069474, 120611.639396 477669.634785, 120628.726231 477711.404736, 120438.778389 477952.766983, 120493.355085 478012.075554, 120615.60531 478064.610414, 120682.024661 478215.297752, 120765.685431 478341.137889, 120741.365823 478431.406044, 120666.98193 478501.99576, 120651.871934 478538.849271, 120970.802939 478490.599916, 120980.75907 478549.880164, 120926.216754 478537.910126, 120905.605778 478594.762462, 120761.645468 478610.343592, 120726.264875 478637.994914, 120434.676019 478680.664912, 120386.674976 478672.632837, 120333.927543 478628.473398, 120267.30136 478642.183867, 120331.577072 478811.346393, 120359.939051 478871.486418, 120409.680923 478868.900481, 120950.039299 478731.378112, 121039.735328 478888.702916, 121064.583659 478878.178247, 121103.551208 478930.206446, 121196.293106 478926.73368, 121234.148353 478977.315399, 121294.892023 478995.29298))"
1695,"POLYGON((136384.812273 518690.957131, 136439.010836 518716.504991, 136488.424541 518637.30669, 136564.04007 518512.936822, 136656.05398 518571.392948, 136723.271433 518426.536312, 136853.944039 518292.580411, 136860.804181 518283.61578, 136865.748176 518273.35799, 136872.962229 518248.755731, 136877.93428 518220.28694, 136886.760843 518219.13979, 136676.546813 518160.984398, 136719.354814 518089.9132, 135713.501156 517791.494577, 135694.275682 517771.926228, 135426.38098 517683.916773, 135420.019738 517667.651474, 135353.5589 517640.650868, 135303.354991 517662.872527, 135310.191193 517691.701398, 135257.093951 517651.426681, 135245.742184 517652.444794, 135230.800455 517679.512383, 135250.395206 517691.192203, 135256.756687 517652.019772, 135288.198946 517685.460549, 135335.467576 517700.169932, 135375.189764 517694.860703, 135352.419156 517738.366502, 135395.551489 517764.449055, 135431.32429 517711.638702, 135638.065769 517767.289231, 135664.503409 517809.032347, 135686.411235 517802.531207, 135784.658635 517832.204375, 135777.523977 517852.55658, 135791.478893 517862.216647, 135810.785561 517838.7453, 135898.156149 517861.531052, 135890.827409 517890.018408, 135956.902934 517905.236847, 136044.858765 517959.818189, 136059.577579 517949.649922, 136149.711704 518007.677697, 136196.7019 517940.939736, 136312.982559 517975.167015, 136288.769746 518023.895655, 136382.856533 518082.46813, 136437.451752 518010.96717, 136470.207643 518021.189147, 136448.175906 518059.919655, 136498.530785 518073.941818, 136566.001599 518120.407174, 136604.849249 518063.430364, 136648.705346 518074.227559, 136623.938498 518125.796885, 136630.307255 518141.255162, 136697.202127 518177.553797, 136875.986735 518220.958756, 136860.327076 518283.330421, 136734.199166 518409.516077, 136723.495227 518425.542092, 136715.343955 518439.346308, 136696.105384 518475.839537, 136655.796798 518568.927787, 136563.913226 518512.356557, 136440.451644 518713.264409, 136387.877497 518685.357228, 136289.643862 518839.152329, 136384.812273 518690.957131))"
De SQL query is overigens ook vanuit QGIS uit te voeren. In het Browser panel rechtermuisklik op de GeoPackage en kies dan Execute SQL… .
Vervolgens SQL query in input veld invoeren op Execute klikken en dan op Load layer klikken:
Vervolgens rechtermuisklik op de laag en dan Export > Save Features As…, om de laag (met enkel WKT geometrieen) op te slaan als CSV:
