Gebieden in een raster vullen door interpolatie
Inleiding
Sommige modelberekeningen vereisen een hoogteraster zonder gaten onder het rekendomein. Dit geldt bijvoorbeeld voor 2D-simulaties met HEC-RAS. Dit hoofdstuk beschrijft hoe HydroToolbox kan helpen om gaten in rasters te dichten middels ruimtelijke interpolatie.
Maar ook voor bijvoorbeeld SOBEK 1D/2D geldt dat soms watergangen in 2D moeten worden ‘dichtgesmeerd’ om te voorkomen dat aan het model de berging en afvoercapciteit ten onrechte dubbel wordt toegekend. Dit kan gebeuren als een watergang zowel in 1D als in 2D is geschematiseerd én deze watergang in 2D meerdere cellen breed is. De berging en afvoercapaciteit zitten dan zowel in het 1D-dwarsprofiel als in de naburige 2D-cellen.
- Kies in het startscherm de optie ‘GIS’
- Kies ‘grid operations’ - ‘adjust’ - ‘fill up gaps by spatial interpolation’
Hier krijgt u twee opties:
- from grid
- from polygon outlines
from grid
Interpolation from grid houdt in dat alleen gaten (nodata-values) worden dichtgeïnterpoleerd en dat dit gebeurt door in de directe omgeving te zoeken naar actieve cellen en hiertussen ruimtelijk te interpoleren.
- Stel de opties naar wens in. Zie onderstaand voor een toelichting
- Druk op ‘Execute’ om het proces te starten. Kies een locatie en naam voor het resultaatbestand
- Grid
- Klik op de button <..> rechts van het veld om te bladeren naar een raster waarin gaten gedicht moeten worden.
- Interpolation
- Op dit moment ondersteunt HydroToolbox uitsluitend ruimtelijke interpolatie volgens Inverse Distance Weighting
- Options:
- Activeer de checkbox als uitsluitend gaten mogen worden gevuld die in hun geheel omringd zijn door geldige waarden (geen nodata_value). Stel een maximale zoekstraal in naar naburige cellen die mee mogen wegen in de interpolatie, uitgedrukt in het aantal cellen afstand.
From polygon outlines
Interpolatie vanuit de omtrek van polygonen. In dit geval berekent HydroToolbox voor ieder polygoon uit een opgegeven shapefile nieuwe opvulwaarden. Het programma doet dit door ruimtelijk te interpoleren tussen de vectorpunten die de omtrek van het polygoon vormen.
- Stel de opties naar wens in. Zie onderstaand voor een toelichting
- Druk op ‘Execute’ om het proces te starten. Kies een locatie en naam voor het resultaatbestand
- Grid
- Klik op de button <..> om te bladeren naar een rasterbestand waarin gebieden moeten worden aangepast.
- Polygon shapefile:
- Klik op de button <..> om te bladeren naar een shapefile met polygonen
Het programma geeft twee opties:
interpolation from elevation grid
In dit geval haalt het programma onder ieder vectorpunt van het onderhavige polygoon de Z-waarde uit het gekozen onderliggende raster. Deze hoogtewaardes worden ruimtelijk geïnterpoleerd en zo wordt het polygoon gevuld.
interpolation from Z-value from shape
In dit geval haalt het programma onder ieder vectorpunt van het onderhavige polygoon de Z-waarde uit het vectorpunt zelf.
Interpolatie uit de Z-waarden van een polygoonbestand is alleen mogelijk wanneer het polygoonbestand van het type PolygonZ is en waarin de vectorpunten dus Z-waarden bevatten.