Deniwelacja

Witaj! :)

Wczoraj wieczorem i znalazłam chwilę, by napisać post. Tym razem tematem naszej pracy będzie... deniwelacja w terenie!
"Deni... ?" Deniwelacja. Nie będę się rozpisywać i powiem to najprościej jak się da - różnica między najniższym i najwyższym punktem wysokościowym w terenie.

ZASTRZEGAM, ŻE PONIŻSZY POST JEST OPARTY NA MOJEJ WIEDZY I MOŻE ZAWIERAĆ NIEWIELKIE BŁĘDY W ODNIESIENIU DO TEGO, CO CHCESZ OSIĄGNĄĆ. (Mój post oparty jest o to, że pierwotna warstwa BĘDZIE docinana)

Od czego zaczynamy? Najpierw zdecydujmy na czym pracujemy. W moim projekcie wykorzystam już wcześniej przygotowany DEM. Zanim jednak zacznę właściwe obliczenia, muszę poświęcić trochę czasu jego obróbce.

Do widoku mapy dodałam interpolację, czyli mój cyfrowy model wysokościowy oraz ramkę obejmującą obszar, który nas interesuje, ona jednak przyda nam się dalej.
Pierwszą znaczącą zmianą DEM będzie przetransponowanie go z warstwy rastrowej do wektorowej. Zrobimy to dość łatwo, wybierając z górnego menu Raster --> Konwersja --> Poligonizuj. Jako plik źródłowy podajemy nasz DEM, w plik wynikowy podajemy lokalizację nowego pliku, który nazwę dem oraz klikamy ok.

(Pojawią się zapewne dwa okienka, które potwierdzamy przyciskiem ok, a pierwotne zamykamy). Warstwa poligonowa utworzona z rastra powinna wyglądać tak:

Jeśli spojrzymy na tę warstwę uważnie zauważymy, że zajmuje ona więcej niż ramka, czyli teren, który nas dotyczy. Aby to zmienić musimy ją przyciąć do ramki. Wektor --> Narzędzia geoprocessingu --> Przytnij. Jako wejściową warstwę wektorową przyjmujemy dem, zaś za maskę przycięcia ramkę. Zapisujemy to w nowym pliku, o nazwie dem_cut. Klikamy ok.

Warstwę dem, DEM oraz ramkę możemy usunąć z widoku mapy. Do projektu dodajemy natomiast warstwę z polami heksagonalnymi.

Teraz już niewiele dzieli nas od końca ;) "Tak szybko?!" No a jak! To QGIS a nie ręczne liczydło (no dobra, czasem bywa powoooooolnie).
Wybieramy z górnego menu Wektor --> Narzędzia zarządzania danymi --> Złącz atrybuty według lokalizacji. W dwa pierwsze okienka wpisujemy nasze dwie warstwy, najpierw heksagony, potem poligon. Zaznaczamy opcję Pobierz podsumowanie przecinających się obiektów, gdzie odznaczamy Średnia, zaś zaznaczamy Min i Max. Zapisujemy to w pliku wyjściowym deni. Wybieramy ok.

Po tym pojawiła się nowa warstwa, która powinna wyglądać jak HEKSAGONY! "Ale... ja mam inną. Dlaczego?" Bo się nie posłuchałeś i dodałeś warstwy w innej kolejności. Zrób to jeszcze raz, uważniej :)

W tabeli atrybutów nowej warstwy są zarówno UniqueID, powierzchnię heksagonów jak i maksymalną i minimalną wysokość.

Teraz musimy tylko policzyć szybko deniwelację. Przyda nam się do tego kalkulator pól.

Jak powyżej, należy stworzyć nowe pole, ja nazwałam je deni. Ustawiamy typ pola na liczby dziesiętne (real), długość na 10, zaś dokładność na 2. Wybieramy Funkcje Pola i wartości, najpierw Maksymalną wysokość, minus, minimalną. klikamy ok. Powstało nowe pole z obliczoną deniwelacją terenu dla każdego heksagonu. Na koniec wchodzimy we właściwości warstwy, wybieramy styl, ustawiamy symbol stopniowy, kolumna: deni,  10 klas, równe przedziały i klikamy klasyfikuj.

Następnie wchodzimy w etykiety, etykietuj tę warstwę z... deni i wybieramy ok. Gotowy projekt powinien wyglądać jak poniżej:

PS. Nie jestem pewna, czy powyższy sposób jest poprawny. Opiera się on jedynie na moim małym doświadczeniu.