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

Posted

Nieuwe kaart online: de vier metropoolregio’s van Nederland.

Aanleiding voor deze kaart was het nieuws dat de directeur van de metropoolregio Rotterdam Den Haag overstapt naar een baan bij het ministerie van Defensie; een overstap die discussie opriep.

De vier Metropoolregio’s zijn opvolgers van de voormalige stads- of plusregios. Metropoolregio’s Rotterdam Den Haag en Eindhoven hebben juridische grondslag op basis van de Wet Gemeenschappelijke Regelingen. Zij zijn verantwoordelijkheid voor regionaal verkeers- en vervoersbeleid. De Metropoolregio’s Arnhem-Nijmegen en Amsterdam daartegen zijn vrijwillige samenwerkingen tussen gemeenten.

Metropoolregio’s volgen gemeentegrenzen, maar overschrijden in tegenstelling tot bijvoorbeeld Veiligheidsregio’s, wel Provinciegrenzen. Metropoolregio Amsterdam omvat twee gemeenten in Flevoland en Mook en Middelaar is de enige Limburgse gemeente tussen 17 Gelderse gemeenten in de metropoolregio Arnhem-Nijmegen.

De kaart is wederom opgebouwd met Leaflet en bevat 2 kaartlagen: de metropoolregio’s en de Provincies. De Provincie-laag ligt ‘onderop’, zodat de mouse-clicks de metropoolregio-kaartlaag bereiken en de gewenste pop-up en highlight veroorzaken. Andersom werkt niet, zo is proefondervindelijk vastgesteld. De pop-up en highlight functionaliteit werkt evenmin als de metropoolregio-kaartlaag deels transparant is.

Author
Categories openbaar bestuur, webmaps

Posted

CC Air55

De voorgestelde fusie tussen gemeenten Barneveld en Scherpenzeel houdt de gemoederen aldaar bezig.

Een interessant aspect van gemeentelijke herindelingen zijn de gevolgen voor de Intergemeentelijke samenwerkingsverbanden. Dit zijn coalities van gemeenten waarbinnen bepaalde gemeentelijke taken gezamenlijk worden uitgevoerd en geregeld, zoals de Veiligheidsregios.

Als 2 gemeentes fuseren kan dat betekenen dat een van de fusie-partners overstapt naar een ander samenwerkingsverband. Zal dat het geval zijn voor Barneveld en Scherpenzeel?

Nee. De Regio Atlas laat zien dat deze gemeenten al samenwerken binnen alle 33 mogelijke samenwerkingsverbanden. De eventuele fusie zal dus geen gevolgen hebben voor de Scherpenzelers (?) en Barnevelders op het gebied van veiligheid (blijft Gelderland-Midden, de omgevingsdienst (blijft de Vallei), Politie-eenheid (blijft Oost-Nederland) en de vele andere indelingen.

De kaarten hieronder illustreren de samenwerkingsintensiteit (van 0 tot 33) van de Gemeenten Barneveld en Scherpenzeel.

Author
Categories openbaar bestuur, webmaps