Doorlatendheid baseren op GeoTOP

[leendertvanwolfswinkel]

Dag collega's,

Voor de gemeente Zevenaar en Venlo bieden we een afkoppelkansenkaart aan. Om de infiltratiekansen naar de ondergrond in beeld te brengen willen we de GeoTOP gebruiken. Daarom hebben Ivar en Leendert het één en ander technisch uitgezocht (voor details kan je dus bij hun zijn). N.a.v. dat uitzoekwerk heb ik nog contact gehad met de Geologische Dienst Nederland. Hieronder de mail met daarin het overzicht.

Omdat het voor meerdere mensen interessant kan zijn, gaat dit over de schutting van mijn 'cubicle':
GeoTOP
De ondergrond van bijna heel Nederland is bij Geologische Dienst Nederland (GDN, onderdeel van TNO) via het DINO-loket beschikbaar als open data. Deze GeoTOP geeft een gedetailleerd driedimensionaal beeld van de ondergrond tot een diepte van 30 meter onder maaiveld, het door de mens meest intensief benutte deel van de ondergrond (Figuur 1). Het is een verfijning van de bovenkant van DGM en REGIS-II. In dit model is de ondergrond in miljoenen voxels (blokjes) van 100 x 100 x 0,5 m opgedeeld. Iedere voxel bevat informatie over de grondsoort en daarmee samenhangende fysische en chemische eigenschappen (Bron: GDN).

Helaas is de GeoTOP nog niet voor heel Nederland beschikbaar. Momenteel zijn ze bezig met Brabant en Limburg. Daarna gaan ze weer verder met het noorden/oosten. Er is geen tijdspad bekend. Voor het ontbrekende deel zouden we gebruik kunnen maken van Nederland 3D.

Groeten,
Jordie
---------- Forwarded message ---------
From: Stafleu, J. (Jan) <[email protected]>
Date: ma 26 nov. 2018 om 16:16
Subject: GeoTOP
To: [email protected] <[email protected]>

Beste Jordie,

Naar aanleiding van ons telefoongesprek stuur ik je een aantal links met informatie over GeoTOP en Nederland 3D.

DINOloket Ondergrondmodellen: https://www.dinoloket.nl/ondergrondmodellen

Als je GeoTOP aanklikt zie je het huidige dekkingsgebied en in oranje de locatie van de boringen die we bij het maken van het model gebruikt hebben.

Vervolgens kun je links op de knop met het vierkantje klikken waarna je op de kaart een interessegebied kan afbakenen. Voor de GIS-bestanden van GeoTOP is dat niet heel spannend omdat je altijd de kaarten van het hele dekkingsgebied ontvangt. Kies na het aangeven van een gebied voor “ArcGIS – landsdekkend” en volg de aanwijzingen op het scherm.

Als alternatief is het ook mogelijk om kaartbladen en software van de “Subsurface viewer” te downloaden. Hiermee kun je het model in 3D bekijken. GeoTOP is hiervoor in kaartbladen van 20 x 25 km opgedeeld en nu is het wel van belang om het juiste interessegebied te kiezen. Elk kaartblad bevat ook een 3D-dataset van het model; hoe je die uit het kaartblad haalt staat beschreven in de handleiding die met de bestelling meegeleverd wordt.

Toelichting – op de pagina waar het huidige dekkingsgebied van GeoTOP staat afgebeeld is ook een tabblad “Toelichting”. Navigeer naar “Detaillering van de bovenste lagen met GeoTOP” en vervolgens “Bekijken en aanvragen GeoTOP”. Je vindt hier uitleg over de GIS-bestanden. Bij “Meer weten?” vind je meer uitgebreide informatie over GeoTOP.

Directe links:

https://www.dinoloket.nl/bekijken-en-aanvragen-geotop

https://www.dinoloket.nl/meer-weten

Nederland 3D:

Dit model staat niet op DINOloket maar op https://www.grondwatertools.nl/nl3d

Een toelichting staat hier: https://www.grondwatertools.nl/toelichting-nl3d-viewer

En een viewer en downloadpagina is hier: https://www.grondwatertools.nl/nl3d-viewer

Legenda’s

In bijgaande twee bestanden vind je de legenda’s van de GIS-bestanden van GeoTOP, een voor de lithoklasse en een voor de geologische eenheden. De kolom voxel-nummer geeft de getalswaarde van de gridcellen weer. De legenda’s voor de Subsurface-viewer staan in de eerder genoemde handleiding.

Mocht je aanvullende vragen hebben, laat het me gerust weten.

Vriendelijke groeten,

Jan Stafle
REF_GTP_LITHO_CLASS.XLSX
REF_GTP_STR_UNIT.XLSX
u

[leendertvanwolfswinkel]

Voorbereidingen voor dit ticket staan op branch 'geotop'

• Het bepalen van de doorlatendheid gaan we niet meer doen met het algoritme dat daarvoor is gemaakt. Dat bepaalt nu alleen nog het bodemtype en heb ik hernoemd.

• geotop data downloaden: staat op utr-qgis-01:
  D:\leendert.vanwolfswin\geotop\GTP14\Model_GTM000000000041\GeoTOP_v01r4s1\GeoTOP_v01r4s1\lagenmodel\grids

• geotop rasters samplen met percelen (statistic "majority"). resultaat moet één vectorlaag zijn waarin van elke laag de majority is weggeschreven. De kolomnamen worden:

geotop_0cm_majority,
geotop_50cm_majority,
geotop_100cm_majority,
geotop_150cm_majority,
geotop_200cm_majority,
geotop_250cm_majority,
geotop_300cm_majority,
geotop_350cm_majority,
geotop_400cm_majority,
geotop_450cm_majority,
geotop_500cm_majority

• De resulterende laag als "geotop" naar postgis database schrijven
• Script 'geotop_infiltration.sql' runnen