Importante

La traduzione è uno sforzo comunitario you can join. Questa pagina è attualmente tradotta al 66.67%.

17.15. Tagliare e unire raster

Nota

In questa lezione vedremo un altro esempio di preparazione di dati spaziali, per continuare a utilizzare geoalgoritmi in scenari reali.

Per questa lezione, andremo a calcolare un raster delle pendenze per una zona intorno a un’area cittadina, che è fornita in un vettore costituito da un singolo poligono. Il DEM di base è diviso in due raster che, insieme, coprono un’area molto più estesa di quella attorno alla città con cui vogliamo lavorare. Se apri il progetto corrispondente a questa lezione, vedrai qualcosa del genere.

Questi raster hanno due problemi:

• Coprono un’area che troppo estesa per i nostri scopi (ci interessa solo una zona più piccola intorno al centro cittadino)

• Si trovano in due file diversi (i confini cittadini si trovano all’interno di un solo raster, ma, come si è già detto, vogliamo avere dello spazio aggiuntivo attorno ad esso).

Entrambi i problemi sono facilmente risolvibili con i geoalgoritmi appropriati.

Innanzitutto, creiamo un rettangolo che definisca l’area di cui abbiamo bisogno. Per fare ciò, creiamo un vettore contenente il perimetro di delimitazione del vettore con i confini dell’area cittadina, e in seguito eseguiamo un buffer su di esso, in modo da avere un raster che si estenda un po” oltre lo spazio minimo necessario.

To calculate the bounding box , we can use the Extract layer extent algorithm

To buffer it, we use the Buffer algorithm, with the following parameter values.

Avvertimento

La sintassi è cambiata nelle ultime versioni; imposta .25 sia per Distanza che per Vertice dell’arco.

Ecco il perimetro di delimitazione risultante ottenuto utilizzando i parametri sopra mostrati.

It is a rounded box, but we can easily get the equivalent box with square angles, by running the Extract layer extent algorithm on it. We could have buffered the city limits first, and then calculate the extent rectangle, saving one step.

You will notice that the rasters has a different projection from the vector. We should therefore reproject them before proceeding further, using the Warp (reproject) tool.

Nota

Le ultime versioni hanno un’interfaccia più complessa. Assicurati che sia selezionato almeno un metodo di compressione.

With this layer that contains the bounding box of the raster layer that we want to obtain, we can crop both of the raster layers, using the Clip raster by mask layer algorithm.

Once the layers have been cropped, they can be merged using the SAGA Mosaic raster layers algorithm.

Nota

Puoi risparmiare del tempo prima unendo e poi tagliando, evitando di richiamare due volte l’algoritmo di ritaglio. Comunque, se ci sono diversi raster da unire ed essi hanno una dimensione piuttosto elevata, otterrai un raster esteso che sarà successivamente difficile da procesare. In tal caso, potrebbe essere necessario lanciare l’algoritmo di ritaglio diverse volte, il che potrebbe essere dispendioso in termini di tempo, ma non preoccuparti, presto vedremo che esistono strumenti aggiuntivi per automatizzare tale operazione. In questo esempio abbiamo solo due raster, per cui al momento non dovresti preoccuparti di ciò.

Così facendo, otteniamo il DEM finale desiderato.

Adesso è il momento di creare il raster delle pendenze.

A slope layer can be computed with the Slope algorithm, but the DEM obtained in the last step is not suitable as input, since elevation values are in meters but cellsize is not expressed in meters (the layer uses a CRS with geographic coordinates). A reprojection is needed. To reproject a raster layer, the Warp (reproject) algorithm can be used again. We reproject into a CRS with meters as units (e.g. 3857), so we can then correctly calculate the slope, with either SAGA or GDAL.

Con il nuovo DEM, la pendenza può essere ora calcolata.

E questo è il raster delle pendenze risultante.

The slope produced by the Slope algorithm can be expressed in degrees or radians; degrees are a more practical and common unit. In case you calculated it in radians, the Metric conversions algorithm will help us to do the conversion (but in case you didn’t know that algorithm existed, you could use the raster calculator that we have already used).

Reprojecting the converted slope layer back with the Reproject raster layer, we get the final layer we wanted.

I processi di riproiezione potrebbero aver fatto in modo che il raster finale contenga alcuni dati al di fuori del perimetro di delimitazione che avevamo calcolato in uno dei passaggi iniziali. Questo può essere risolto tagliandolo di nuovo, così come abbiamo fatto per ottenere il DEM di base.