18.16. Hydrologische analyse¶

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.

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.

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.

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).

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).

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.

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.

En dit is wat u zult krijgen.

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

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,

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

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

U zult iets krijgen zoals dit.

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

De resulterende statistieken zijn de volgende.

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.