Numeryczny Model Terenu

Interpolacja rastra

Jeśli dane wysokościowe jakimi dysponujemy istnieją tylko w postaci wektorowej, w pierwszej kolejności powinniśmy wykonać interpolację do postaci rastrowej. Dzięki dostępowi do narzędzi SAGA oraz GRASS, mamy możliwość skorzystania z kilku uznanych w środowisku naukowym metod interpolacyjnych, m.in.:

• Kriging
• Thin plate Spline w odmianach zwykłej i TIN
• Multilevel B-Spline
• Metoda odwrotnych odległości
• Metoda najbliższego sąsiedztwa

Technicznie operacja interpolacji w przypadku większości z tych metod jest bardzo podobna. Powinniśmy wskazać warstwę wektorową opisaną w oknie jako Points, pole atrybutu, oraz rozmiar oczka siatki rastra Gridsize.

Preprocessing "hydrologiczny"

Każdorazowo numeryczny model terenu utworzony przy pomocy interpolacji danych wektorowych, lub pochodzący ze zbiorów takich jak SRTM, AsterDEM, EU-DEM powinniśmy poddać procesowi usuwania zagłębień bezodpływowych (sink removal). Tu również z pomocą przychodzi nam pakiet SAGA, który w zakładce Terrain Analysis - Hydrology udostępnia kilka różnych wersji narzędzia, stosujących różne algorytmy. Wybór konkretnego algorytmu pozostawiam własnym doświadczeniom, lub rekomendacjom w literaturze. W naszym przykładzie skorzystamy z narzędzia opisanego jako Sink removal.

Pozwala nam on dostosować metodę usuwania zagłębień na dwa sposoby - poprzez pogłębianie kanałów spływu powierzchniowego, lub poprzez "zasypywanie".

Kanały spływu powierzchniowego, zlewnie

Kolejne operacje które możemy wykonać przy pomocy narzędzi geoprocessingu to wyznaczenie kanałów spływu powierzchniowego, oraz zlewni. Znajdziemy je w zakładce Terrain Analysis - Channels.

W wyniku tej operacji możemy otrzymać warstwy wektorowe dla zlewni i kanałów, czy rastrowe warstwy określające kierunek spływu powierzchniowego.

Morfometria rzeźby terenu

W zakładce Terrain analysis - morphometry znajdziemy zgrupowane ważniejsze narzędzia i algorytmy do obliczeń cech morfometrycznych rzeźby terenu takich jak nachylenie stoku, ekspozycja, krzywizny, wskaźnik TRI.

Tworzenie warstwic

Jeśli potrzebujemy stworzyć dla celów wizualizacji warstwę wektorową z przebiegiem warstwic, korzystamy w tym celu z narzędzia znajdującego się w menu Raster -> Cięcie i wybieramy pozycję Generuj warstwice.

Następnie wskazujemy przy pomocy przycisku Wybierz plik docelowy, jak również zaznaczymy pole Nazwa atrybutu i wprowadzamy naszą nazwę pola. W polu Cięcie warstwicowe ustawiamy parametr cięcia, odstępu wysokościowego pomiędzy kolejnymi warstwicami.