Terug van vakantie! Een van de activiteiten was een 9-daagse wandeltocht door de Franse Alpen met als hoogtepunt (hihi) de beklimming van Mont Thabor (kaart), zie foto hierboven (foto door FvL).

Met dank aan de organisatie van Spring reizen!

Nieuwe kaart online: de 29 Regionale Uitvoeringsdiensten.

De regionale uitvoeringsdiensten zijn verantwoordelijk voor vergunningverlening, toezicht en handhaving op het gebied van natuurwetgeving, bouwtoezicht, veiligheid, luchtkwaliteit, geluid, energie, afval, asbest en bodem. De omgevingsdiensten werken in opdracht van Gemeenten en Provincies. Omgevingsdiensten zijn niet gelijk aan de afvalverwijderingsdiensten, zoals AVRI en ACV.

Twee omgevingsdiensten bedienen een niet-aaneengesloten gebied. RUD Utrecht is werkzaam in 3 clusters van Gemeenten. De resterende gemeenten in Utrecht worden bediend door de Omgevingsdienst Regio Utrecht, wiens werkgebied maar net aaneengesloten is. In Noord-Holland bediend Omgevingsdienst IJmond twee clusters van gemeenten.

De gemeentelijke indelingen van de omgevingsdiensten wordt gespecificeerd op de site van omgevingsdienst NL. De regio atlas bevat geen informatie over de regionale diensten zelf, maar laat wel zien hoe de 29 diensten zijn geclusterd tot zes Besluit Risico’s Zware Ongevallen regios.

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>";

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.

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