17.16. Analisi idrologica¶
Nota
In questa lezione eseguiremo alcune analisi idrologiche. Questa analisi verrá usata in alcune delle prossime lezioni, dato che costituisce un buon esempio di flusso di lavoro per l’analisi, e verrá utilizzato per dimostrare alcune funzionalitá avanzate.
In questa lezione, eseguiremo alcune analisi idrologiche. Iniziando con un DEM, estrarremo una rete di canali, delineeremo gli spartiacque e calcoleremo alcune statistiche.
La prima cosa da fare é caricare il progetto con i dati della lezione, che contengono il DEM.
The first module to execute is Catchment area. You can use anyone of the others named Catchment area. They have different algorithms underneath, but the results are basically the same.
Selezionare il DEM nel campo Elevazione, e lasciare i valori di default per i restanti parametri.
Alcuni algoritmi calcolano diversi layer, ma il Bacino di utenza é l’unico che utilizzeremo.
Potete eliminare gli altri se volete.
La visualizzazione del layer non é molto informativa.
Per sapere il perché, si puó controllare l’istogramma e notare che i valori non sono uniformemente distribuiti (alcune celle hanno dei valori molto alti, ovvero quelle che corrispondono alla rete di canali). Il calcolo del logaritmo del bacino di utenza porta ad un layer che contiene molta piú informazione (tale risultato si puó ottenere tramite il calcolatore raster).
Il bacino di utenza (anche noto come accumulazione di flusso), puó essere utilizzato per impostare una soglia per l’inizializzazione dei canali. Ció si puó fare utilizzando l’algoritmo Rete di canali. Qui le istruzioni per la configurazione (si noti che Soglia di inizializzazione Maggiore di 10.000.000).
Utilizzare il layer di bacino di utenza originale, non quello logaritmico. Quest’ultimo serve solo per la visualizzazione.
Se si aumenta il valore della Soglia di inizializzazione, si ha una rete di canali piú sparsa. Se si diminuisce, si ha una rete piú densa. Con il valore proposto, questo é ció che si ottiene.
L’immagine precedente mostra il layer vettoriale risultante ed il DEM, ma ci dovrebbe anche essere un layer raster con la stessa rete di canali. Il raster layer sará effettivamente quello che verrá utilizzato.
Useremo ora l’algoritmo Bacini spartiacque per delineare i sottobacini corrispondenti alla rete di canali, utilizzando come punti di sbocco tutte le giunzioni. Ecco come impostare i parametri.
E questo é il risultato.
Questo é il risultato raster. Si puó vettorizzare utilizzando l’algoritmo Vectorising grid classes.
Si calcolino adesso le statistiche dei valori di elevazione in uno dei sottobacini. L’idea é quella di avere un layer che rappresenti l’elevazione nei sottobacini e quindi passarlo al modulo che calcola tali statistiche.
First, let’s clip the original DEM with the polygon representing a subbasin. We will use the Clip grid with polygon algorithm. If we select a single subbasin polygon and then call the clipping algorithm, we can clip the DEM to the area covered by that polygon, since the algorithm is aware of the selection.
Selezionare un poligono,
ed invocare l’algoritmo con i seguenti parametri:
L’elemento selezionato nel campo di input é, ovviamente, il DEM che si vuole ritagliare.
Si otterrá qualcosa del genere.
Il layer é pronto per essere utilizzato nell’algoritmo Raster layer statistics.
Di seguito le statistiche risultanti.
Utilizzeremo sia la procedura per il calcolo del bacino che il calcolo delle statistiche in altre lezioni, per trovare altri elementi che ci possano aiutare ad autmatizzarli e lavorare in modo piú efficiente.