Outdated version of the documentation. Find the latest one here.

17.15. Decuparea și îmbinarea straturilor raster

Note

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

../../../_images/medfordarea.png

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.

Pentru a calcula caseta de încadrare, putem folosi algoritmul Poligonului din extinderea stratului

../../../_images/bbox.png

Pentru a-l tampona, vom folosi algorimul Fixed distance buffer, cu următoarele valori pentru parametri.

../../../_images/buffer_dialog.png

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

../../../_images/buffer1.png

Este o casetă rotundă, dar putem obține cu ușurință o casetă echivalentă, cu unghiuri drepte, prin rularea algoritmului Poligon din extinderea stratului asupra ei. Am putea să tamponăm mai întâi limitele orașului, iar apoi să calculăm extinderea dreptunghiului, economisind un pas.

../../../_images/buffer_squared.png

Veți observa că rasterele au o proiecție diferită față de vector. Prin urmare, ar trebui să le reproiectăm înainte de a trece mai departe, folosind instrumentul Warp (reproiectare).

../../../_images/warp1.png

Note

Recent versions have a more complex interface. Select at least one compression method.

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 grid with polygon algorithm.

../../../_images/clip1.png

După straturile au fost decupate, ele se pot uni cu ajutorul algoritmului GDAL Merge.

../../../_images/merge1.png

Warning

The following paragraph is obsolete.

A cellsize is needed for the merged layer. We will use the same one of the input ones. You do not need to know it in advance before calling the algorithm. Just click on the button in the right–hand size of the text field and you will have a dialog to enter small mathematical formulas, and a list of frequently used values, among them the cellsizes and bounding coordinates of all available layers.

Note

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.

../../../_images/finaldem.png

Acum este timpul să calculăm stratul pantei.

Stratul pantei poate fi calculat cu ajutorul algoritmului Pantă, Aspect, Curbură, însă DEM-ul obținut în ultima etapă nu este potrivit ca intrare, deoarece valorile altitudinii sunt exprimate în metri, iar mărimea celulei nu este exprimată în metri (stratul folosește un CRS cu coordonate geografice). De aceea, este nevoie de o reproiectare. Pentru a reproiecta stratul raster, se poate utiliza iarăși algoritmul Warp (reproiectare). Vom efectua reproiectarea într-un CRS având metrul ca unitate (cum ar fi 3857), astfel încât vom putea calcula corect panta, fie cu SAGA, ori cu GDAL.

Here is the reprojected DEM.

../../../_images/warpeddem.png

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

../../../_images/slope.png

Iar aici este stratul pantei rezultate.

../../../_images/slopereproj.png

The slope produced by the Slope, Aspect, Curvature algorithm is expressed in radians, but degrees are a more practical and common unit. 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).

../../../_images/metricconversions.png

Reproiectând stratului pantă convertit, cu ajutorul Reproiectării stratului raster, obținem stratul final pe care l-am dorit.

Warning

todo: De adăugat imaginea

The reprojection processes have caused the final layer to contain data outside the bounding box that we calculated in one of the first steps. This can be solved by clipping it again, as we did to obtain the base DEM.