Important

Traducerea este un efort al comunității, la care puteți să vă alăturați. În prezent, această pagină este tradusă 66.67%.

17.15. Decuparea și îmbinarea straturilor raster

Notă

În această lecție vom vedea un alt exemplu de pregătire a datelor spațiale, pentru a continua utilizarea geoalgoritmilor în scenarii din lumea reală.

În această lecție, vom calcula un strat de pantă pentru suprafața care înconjoară o zonă dată a orașului, dintr-un strat vectorial cu un singur poligon. DEM-ul de bază este împărțit în două straturi raster care, împreună, acoperă o suprafață mult mai mare decât cea din jurul orașului în care ne dorim să lucrăm. Dacă deschideți proiectul corespunzător acestei lecții, veți vedea ceva de genul următor.

Aceste straturi au două probleme:

• Acestea acoperă o zonă care este prea mare pentru ceea ce dorim (suntem interesați de o regiune mai mică din jurul centrului orașului)

• Ele se află în două fișiere diferite (limitele orașului se încadrează doar într-un singur strat raster, dar, așa cum s-a mai zis, dorim o anumită suprafață adițională în jurul acestuia).

Ambele sunt ușor rezolvabile cu geoalgoritmii corespunzători.

În primul rând, vom crea un dreptunghi care definește zona dorită. Pentru aceasta, vom crea un strat care conține caseta de încadrare a stratului, împreună cu limitele suprafeței orașului, apoi vom crea un tampon, astfel încât să existe un strat raster care îl acoperă un pic mai mult decât este necesar.

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.

Atenționare

Sintaxa s-a schimbat în versiunile recente; setați atât Distanța cât și vertexul Arcului la .25

Aici se află caseta de încadrare obținută, utilizând parametrii de mai sus

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.

Notă

Versiunile recente au o interfață mai complexă. Asigurați-vă că cel puțin o metodă de compresie este selectată.

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.

Notă

Puteți economisi timp prin efectuarea mai întâi a îmbinării, și abia mai apoi a decupării, evitându-se astfel apelarea de două ori a algoritmului de decupare. Totuși, în cazul în care mai multe straturi cu o dimensiune apreciabilă trebuie îmbinate, vă veți alege cu un strat voluminos, care poate fi dificil de prelucrat ulterior. În acest caz, s-ar putea avea să apelați algoritmul de tăiere de mai multe ori, operație consumatoare de timp, însă nu vă îngrijorați pentru că veți vedea că există unele instrumente adiționale, dedicate automatizării acestui proces. Pentru exemplul următor nu este cazul să vă faceți griji, deoarece folosim numai două straturi.

Cu asta, vom obține DEM-ul final pe care ni-l dorim.

Acum este timpul să calculăm stratul pantei.

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.

Panta poate fi de acum calculată, cu ajutorul noului DEM.

Iar aici este stratul pantei rezultate.

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.

Datorită proceselor de reproiectare, stratul final ar putea conține datele din afara casetei de încadrare, pe care am calculat-o într-unul dintre primii pași. Acest lucru poate fi rezolvat prin reluarea decupării, așa cum am procedat la obținerea DEM-ului de bază.