18.16. Hydrologische analyse

Notitie

In deze les zullen we een hydrologische analyse uitvoeren. Deze analyse zal worden gebruikt in enkele van de volgende lessen, omdat het bestaat uit een goed voorbeeld van een werkstroom voor een analyse, en we het zullen gebruiken voor het demonstreren van enkele geavanceerde mogelijkheden.

In deze les zullen we een hydrologische analyse uitvoeren. Beginnend met een DEM, zullen we een netwerk van kanalen gaan uitnemen, waterkeringen uittekenen en enkele statistieken berekenen.

Het eerste is om het project te laden met de gegevens voor de les, wat slechts een DEM bevat.

../../../_images/dem1.png

De eerste module die moet worden uitgevoerd is Catchment area (in sommige versies van SAGA wordt hij Flow accumulation (Top Down) genoemd). U kunt elke van de andere genaamd Catchment area gebruiken. Zij hebben verschillende onderliggende algoritmen maar de resultaten zijn in de basis hetzelfde.

Selecteer de DEM in het veld Elevation, en laat de rest van de parameters staan op de standaard waarden.

../../../_images/catchmentarea.png

Sommige algoritmen berekenen vele lagen, maar de laag Catchment Area is de enige die we zullen gebruiken.

U kunt de andere weggooien als u dat wilt.

Het renderen van de laag is niet erg informatief.

../../../_images/catchmentlayer.png

To know why, you can have a look at the histogram and you will see that values are not evenly distributed (there are a few cells with very high value, those corresponding to the channel network). Calculating the logarithm of the catchment area value yields a layer that conveys much more information (you can do it using the raster calculator).

../../../_images/catchmentlayerlog.png

Het opvanggebied (ook bekend als flow accumulation), kan worden gebruikt om een drempel in te stellen voor het initiëren van kanalen. Dit kan worden gedaan met behulp van het algoritme Channel network. Hier staat hoe u het in moet stellen (let op de Initiation threshold Greater than 10.000.000).

../../../_images/channelnetwork.png

gebruik de originele laag catchment area, niet die van de logaritme. Die was alleen om te renderen.

Als u de waarde Initiation threshold verhoogt, zult u een spaarzamer netwerk van kanalen verkrijgen. Als u het verlaagt, zult u een dichter verkrijgen. Dit is wat u krijgt met de voorgestelde waarde.

../../../_images/channelnetworklayer.png

De afbeelding hierboven geeft slechts de resulterende vectorlaag en de DEM weer, maar er zou ook een rasterlaag moeten zijn met hetzelfde netwerk van kanalen. Die rasterlaag zal in feite degene zijn die we zullen gaan gebruiken.

Nu zullen we het algoritme Watersheds basins gaan gebruiken om de subbassins uit te tekenen die overeenkomen met dat netwerk van kanalen, waarbij we alle kruisingen daarin gebruiken als uitlaatpunt. Hier staat hoe u het overeenkomende dialoogvensters met parameters moet instellen.

../../../_images/watersheds.png

En dit is wat u zult krijgen.

../../../_images/watershedslayer.png

Dit is een rasterresultaat. U kunt het vectoriseren met behulp van het algoritme Vectorising grid classes.

../../../_images/vectorising.png ../../../_images/watershedslayervector.png

Laten we nu eens proberen statistieken te berekenen over de hoogtewaarden in één van de subbassins. Het idee is om een laag te krijgen die slechts de hoogte binnen dat subbassin weergeeft en die dan doorgeeft naar de module die die statistieken berekent.

Laten we eerst de originele DEM clippen met de polygoon die een subbassin weergeeft. We zullen het algoritme Clip raster with polygon gebruiken. Als we één enkele polygoon voor een subbassin selecteren en dan het algoritme voor het clippen aanroepen, kunnen we de DEM clippen tot het gebied dat wordt bedekt door die polygoon, omdat het algoritme zich bewust is van de selectie.

Selecteer een polygoon,

../../../_images/selectedpolygon.png

en roep het algortime om te clippen aan met de volgende parameters:

../../../_images/clipgrid.png

Het geselecteerde element in het invoerveld is, natuurlijk, de DEM die we willen clippen.

U zult iets krijgen zoals dit.

../../../_images/clippeddem.png

Deze laag is gereed om te worden gebruikt in het algoritme Raster layer statistics.

../../../_images/rasterstats.png

De resulterende statistieken zijn de volgende.

../../../_images/stats.png

We zullen zowel de procedure voor het berekenen van het bassin gebruiken als de berekeningen voor statistieken uit andere lessen, om uit te zoeken hoe andere elementen ons kunnen helpen om beide te automatiseren en meer effectief te werken.