17.15. Decuparea și îmbinarea straturilor raster¶

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

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

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

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

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.

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

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.

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

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.

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

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.

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

Iar aici este stratul pantei rezultate.

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

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

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.