Posted

De popup informatie van de Veiligheidsregio’s en RES-regio’s kaarten is uitgebreid met het aantal inwoners (CBS) en het aantal onderliggende gemeenten. Daarnaast is de regio-naam nu gelijk de hyperlink naar de regio-website.

De popup wordt vormgegeven in html en de feature attributes zijn toegankelijk via feature.properties.XXXX:

var popupContent = "<p><b><a href=" + feature.properties.website + " target=blank>RES " + feature.properties.naam + "</a></b></p><p>oppervlakte: <b>" +  feature.properties.area_sq_km + "</b> km<sup>2</sup><br><b>" + feature.properties.bevolking2020 + "</b> inwoners (CBS, 2020)<br>bestaat uit <b>" + feature.properties.n_gemeenten + "</b> gemeenten</p>";

Author

Posted

Nieuwe kaart online: de 30 regios van het nationaal programma regionale energie strategie, oftewel de res regios.

De regionale energiestransities zijn niet alleen regelmatig in het nieuws, maar het ‘officiele’ kaartje van de regio indeling is ook nog eens vrij onduidelijk. Reden te meer om deze regio’s op de kaart te zetten.

De RES regio’s zijn samengesteld uit een of meer Gemeenten. Net als bij de Metropoolregio’s worden Provinciegrenzen overschreden: RES regio Food Valley ligt in zowel Gelderland als Utrecht. In Friesland, Groningen, Drenthe, Flevoland en Zeeland valen de RES regio’s samen met de Provinciegrenzen. Tegelijkertijd zijn er twee RES regio’s niet groter dan een gemeente: RES Goeree-Overflakkee en RES Hoeksche Waard.

Hoewel de begrenzing van RES Rotterdam Den Haag gelijk is aan die van de gelijknamige Metropoolregio, lijken dit toch twee verschillende organisaties te zijn. RES Arnhem-Nijmegen lijkt identiek aan de “Groene Metropool Arnhem Nijmegen”, maar bevat niet Gemeente Mook en Middelaar.

De kaart hanteert dezelfde Leaflet functionaliteit als de Veiligheidsregio’s kaart.

Author
Categories openbaar bestuur, webmaps

Posted

Some interesting snippets from a larger batch script where several gdal raster utilities are called upon for rasterizing, raster calculations and raster aggregations.

create timestamp

C:\> SET timestamp=%date:~0, 4%%date:~5,2%%date:~8,2%t%time:~0,2%%time:~3,2%%time:~6,2%
C:\> echo %timestamp%
20210623t135528

vector 2 raster

gdal_rasterize -at -l  %input_layer% -of GTiff -te 10000 300000 280000 625000 -tr 2.5 2.5 -ot Byte -a_srs epsg:28992 -co COMPRESS=LZW -burn 1 %input_gdb% %input_layer% %output_name%_%timestamp%.tif

Target resolution and extent are specified with -tr and -te. Vector shapes are set to value 1 with -burn 1.

raster calculations

python "%gdal_dir%/gdal_calc.py" --co COMPRESS=LZW -A %sinput_A% -B %input_B% -C %input_C% --outfile=%output_name%_%timestamp%.tif --calc="numpy.where((A+B+C)>=1, 1, 0)"

Add three 0/1 rasters together and set to 1 where the sum is equal to or larger than 1, else 0. The gdal_calc.py utility is called with python, which of course only works if a proper python environment is activated:

call c:\apps\Miniconda3\Scripts\activate.bat
call conda activate environment_name

raster aggregations

gdalwarp -co COMPRESS=LZW -tr 25 25 -te 10000 300000 280000 625000 -r sum %input%.tif %output_name%_%timestamp%_25m.tif 

Aggregating from 2.5 to 25m, again by setting -tr and -te. Resampling method is sum, which is available in gdalwarp but not in gdal_translate

Author
Categories gdal

Posted

Copy ESRI Filegeodatabase featureclass to Shapefile with Selection from Attributes and for selected Attribute columns only using GDAL ogr2ogr

ogr2ogr -sql "SELECT * FROM IN_LAYER WHERE ATTRIBUTE_X IN ('Value1', 'Value2')" -select ATTRIBUTE_X, ATTRIBUTE_Y output.shp input.gdb IN_LAYER

Author
Categories gdal, GIS