r.slope.aspectgrass — online w chmurze

PROGRAM:

IMIĘ

r.slope.aspekt - Generuje mapy rastrowe nachylenia, aspektu, krzywizny i częściowe
pochodne z rastrowej mapy wysokościowej.
Aspekt jest obliczany w kierunku przeciwnym do ruchu wskazówek zegara od wschodu.

SŁOWA KLUCZOWE

raster, teren, aspekt, nachylenie, krzywizna

STRESZCZENIE

r.slope.aspekt
r.slope.aspekt --help
r.slope.aspekt [-a] podniesienie=Nazwa [nachylenie=Nazwa] [aspekt=Nazwa] [format=ciąg]
[precyzja=ciąg] [krzywizna=Nazwa] [krzywizna=Nazwa] [dx=Nazwa] [dy=Nazwa]
[dxx=Nazwa] [barwiony=Nazwa] [dxy=Nazwa] [zskala=unosić się] [min_nachylenie=unosić się] [--przepisać]
[--pomoc] [--gadatliwy] [--cichy] [--ui]

Flagi:
-a
Nie wyrównuj bieżącego regionu do rastrowej mapy wysokości

--przepisać
Zezwalaj plikom wyjściowym na zastępowanie istniejących plików

--help
Wydrukuj podsumowanie wykorzystania

--gadatliwy
Pełne wyjście modułu

--cichy
Cichy moduł wyjściowy

--UI
Wymuś uruchomienie okna GUI

Parametry:
podniesienie=Nazwa [wymagany]
Nazwa wejściowej mapy rastrowej elewacji

nachylenie=Nazwa
Nazwa wyjściowej mapy rastrowej nachylenia

aspekt=Nazwa
Nazwa wyjściowej mapy rastrowej proporcji

format=ciąg
Format zgłaszania nachylenia
Opcje: stopnie, procent
Zaniedbanie: stopni

precyzja=ciąg
Rodzaj wyjściowych map aspektów i nachyleń
Opcje: KOMÓRKA, FCELL, DCELL
Zaniedbanie: FCELL

krzywizna=Nazwa
Nazwa mapy rastrowej krzywizny profilu wyjściowego

krzywizna=Nazwa
Nazwa wyjściowej mapy rastrowej krzywizny stycznej

dx=Nazwa
Nazwa wyjściowej mapy rastrowej pochodnej cząstkowej pierwszego rzędu dx (nachylenie EW).

dy=Nazwa
Nazwa wyjściowej mapy rastrowej pochodnej cząstkowej pierwszego rzędu dy (nachylenie NS).

dxx=Nazwa
Nazwa wyjściowej mapy rastrowej dxx pochodnej cząstkowej drugiego rzędu

barwiony=Nazwa
Nazwa wyjściowej mapy rastrowej dyy pochodnej cząstkowej drugiego rzędu

dxy=Nazwa
Nazwa wyjściowej mapy rastrowej dxy pochodnej cząstkowej drugiego rzędu

zskala=unosić się
Mnożnik do konwersji jednostek wysokości na jednostki poziome
Zaniedbanie: 1.0

min_nachylenie=unosić się
Minimalna wartość nachylenia (w procentach), dla której obliczany jest aspekt
Zaniedbanie: 0.0

OPIS

r.slope.aspekt generuje mapy rastrowe nachylenia, aspektu, krzywizn oraz pierwszej i drugiej
uporządkować pochodne cząstkowe z mapy rastrowej rzeczywistych wartości wysokości. Użytkownik musi
określ wejście podniesienie mapa rastrowa i co najmniej jedna wyjściowa mapa rastrowa. Użytkownik może
również określ format dla nachylenia (stopnie, procent; domyślnie = stopnie) i zskala:
mnożnik do konwersji jednostek wysokości na jednostki poziome; (domyślnie 1.0).

podniesienie wejściowa mapa rastrowa określona przez użytkownika musi zawierać rzeczywiste wartości wysokości,
nie dane przeskalowane lub skategoryzowane. Jeśli wartości wysokości są w innych jednostkach niż w
jednostek poziomych, należy je przekonwertować na jednostki poziome za pomocą parametru zskala.
In TRAWA GIS 7, pionowy jednostek jest nie przypuszczalny do be metrów każdy więcej. Dla litu szacuje się przykład, if
obie Twój pionowy oraz poziomy jednostek jest stóp, parametr zskala nie wolno używać.

aspekt wyjściowa mapa rastrowa wskazuje kierunek, w którym skierowane są zbocza. Aspekt
kategorie reprezentują liczbę stopni na wschód. Pliki kategorii i tabeli kolorów są również
generowane dla mapy rastrowej proporcji. Kategorie aspektu reprezentują liczbę stopni
wschód i zwiększają się w kierunku przeciwnym do ruchu wskazówek zegara: 90 stopni to północ, 180 to zachód, 270 to południe
360 to wschód.
Uwaga: Wartości te można przekształcić na wartości azymutu (0 to północ, 90 to wschód itd.).
używając r.mapcalc:
# przekonwertuj kąty z CCW na północ w górę
r.mapcalc "aspekt_azymutu = (450 - ccw_aspekt) % 360"

Aspekt nie jest zdefiniowany dla nachylenia równego zeru. Tak więc większość komórek o bardzo małych
nachylenie ma kategorię 0, 45, ..., 360 cali aspekt wyjście. Istnieje możliwość zmniejszenia
odchylenie w tych kierunkach poprzez odfiltrowanie aspektu w obszarach, w których znajduje się teren
prawie płaski. Opcja min_nachylenie można użyć do określenia minimalnego nachylenia dla którego aspektu
jest obliczany. Aspekt dla wszystkich komórek z nachyleniem min_nachylenie jest ustawione na zero (brak danych).

nachylenie wyjściowa mapa rastrowa zawiera wartości nachylenia, wyrażone w stopniach nachylenia od
poziomy jeśli formatWybrana jest opcja =stopnie (domyślna) i procentowy wzrost jeśli
format= wybrano opcję procentową. Generowane są pliki kategorii i tabeli kolorów.

Krzywizny profilowe i styczne to krzywizny w kierunku najbardziej stromego zbocza
i odpowiednio w kierunku stycznej konturu. Krzywizny są wyrażone jako
1/metr, np. krzywizna 0.05 odpowiada promieniowi krzywizny 20m. Wypukły
wartości kształtu są dodatnie, a wartości kształtu wklęsłego są ujemne.

Przykład DEM

Nachylenie (stopień) z przykładu DEM Aspekt (stopień) z przykładu DEM

Krzywizna styczna (m-1) z przykładu DEM Krzywizna profilu (m-1) z przykładu DEM

W przypadku niektórych zastosowań użytkownik będzie chciał użyć przeklasyfikowanej rastrowej mapy nachylenia
grupuje wartości nachylenia w zakresy nachylenia. Można to zrobić za pomocą r. przeklasyfikuj. Przykład
przydatna reklasyfikacja jest podana poniżej:
zakres kategorii etykiety kategorii
(w stopniach) (w procentach)
1 0- 1 0- 2%
2 2- 3 3- 5%
3 4- 5 6- 10%
4 6- 8 11- 15%
5 9- 11 16- 20%
6 12- 14 21- 25%
7 15-90 26% i więcej
Poniższa tabela kolorów działa dobrze z powyższym
ponowna klasyfikacja.
kategoria czerwony zielony niebieski
0 179 179 179
1 0 102 0
2 0 153 0
3 128 153 0
4 204 179 0
5 128 51 51
6 255 0 0
7 0 0 0

UWAGI

Aby upewnić się, że rastrowa mapa wysokości nie zostanie nieprawidłowo ponownie próbkowana, ustawienia dla
bieżący region są nieznacznie modyfikowane (tylko w celu wykonania programu): the
rozdzielczość jest ustawiona tak, aby odpowiadała rozdzielczości mapy rastrowej elewacji i krawędzi mapy
regiony (tj. północ, południe, wschód i zachód) są przesuwane, jeśli to konieczne, w celu wyrównania
krawędzie najbliższych komórek na mapie wysokościowej. Jeśli użytkownik naprawdę chce raster
mapa wysokości ponownie próbkowana do bieżącej rozdzielczości regionu, -a należy określić flagę.

Bieżąca maska ​​jest ignorowana.

Algorytm używany do określania nachylenia i proporcji wykorzystuje sąsiedztwo 3x3 wokół każdej komórki
na rastrowej mapie wysokościowej. W związku z tym nie jest możliwe określenie nachylenia i aspektu dla
komórki sąsiadujące z krawędziami w warstwie mapy wysokościowej. Komórkom tym przypisano a
wartość „zerowego nachylenia” (kategoria 0) zarówno na mapach rastrowych nachylenia, jak i proporcji.

Formuła Horna służy do znajdowania pochodnych pierwszego rzędu w kierunkach x i y.

Tylko w przypadku korzystania z całkowitych modeli elewacji aspekt jest obciążony w 0, 45, 90, 180, 225,
kierunki 270, 315 i 360; tzn. rozkład kategorii aspektów jest bardzo nierówny,
ze szczytami w 0, 45,..., 360 kategoriach. Podczas pracy z elewacją zmiennoprzecinkową
modeli, nie występuje taka stronniczość aspektu.

PRZYKŁADY

Obliczenie of nachylenie, aspekt, profil oraz styczny krzywizna
W tym przykładzie mapa nachylenia, kształtu, profilu i stycznej krzywizny jest obliczana na podstawie
mapa rastrowa wysokości (przykładowy zestaw danych Karoliny Północnej):
g.region raster=elewacja
r.slope.aspect elewacja=elewacja nachylenie=nachylenie aspekt=aspekt pcurvature=pcurvure tcurvature=tcurv
# ustaw ładne tabele kolorów dla wyjściowych map rastrowych
r.colors -n map=nachylenie kolor=sepia
r.colors map=aspekt kolor=aspektkolor
r.colors map=pcurv color=krzywizna
r.colors map=tcurv color=krzywizna

Rysunek: Mapa rastrowa nachylenia, kształtu, profilu i stycznej krzywizny (Karolina Północna
zbiór danych)

Klasyfikacja of poważny aspekt kierunki in kompas orientacja
W poniższym przykładzie (na podstawie przykładowego zestawu danych z Karoliny Północnej) najpierw generujemy
standardową mapę proporcji (ukierunkowaną na CCW od wschodu), a następnie przekonwertować ją na orientację kompasu,
i ostatecznie sklasyfikować cztery główne kierunki aspektu (N, E, S, W):
g.region raster=elewacja -p
# generuj mapę aspektów z orientacją CCW
r.slope.aspect elewacja=elewacja aspekt=mójaspekt
# wygeneruj orientację kompasu i sklasyfikowaj cztery główne kierunki (N, E, S, W)
r.mapcalc "aspect_4_directions = eval( \\
kompas = (450 - mój aspekt ) % 360, \\
if(kompas >=0. && kompas < 45., 1) \\
+ if(kompas >=45. && kompas < 135., 2) \\
+ if(kompas >=135. && kompas < 225., 3) \\
+ if(kompas >=225. && kompas < 315., 4) \\
+ if(kompas >=315., 1) \\
)"
# przypisz etykiety tekstowe
r.category aspekt_4_kierunki separator=przecinek reguły=- << EOF
1, północ
2, wschód
3, południe
4, zachód
EOF
# przypisz tabelę kolorów
r.colors Aspect_4_directions Rules=- << EOF
1 253,184,99
2 178,171,210
3 230,97,1
4 94,60,153
EOF
Mapa aspektowa sklasyfikowana według czterech głównych kierunków kompasu (pokazany podzbiór w powiększeniu)

LITERATURA

· Róg, BKP (1981). Wzgórze Zacienienie oraz dotychczasowy Odblask Mapa, Postępowanie
IEEE, 69(1): 14-47.

· Mitasowa, H. (1985). Kartograficzny aspekty of komputer powierzchnia modelowanie. Dr
Praca dyplomowa. Słowacki Uniwersytet Techniczny w Bratysławie

· J. Hofierka, H. Mitasowa, M. Neteler, 2009. Geomorfometria in TRAWA GIS. W:
Hengl, T. i Reuter, HI (red.), Geomorfometria: Koncepcje, Oprogramowanie,
Aplikacje. Rozwój w gleboznawstwie, tom. 33, Elsevier, 387-410 s.,
http://www.geomorphometry.org

Korzystaj z r.slope.aspectgrass online, korzystając z usług onworks.net