r.slope.aspectgrass - Online în cloud

PROGRAM:

NUME

r.pantă.aspect - Generează hărți raster de pantă, aspect, curburi și parțiale
derivate dintr-o hartă raster de altitudine.
Aspectul este calculat în sens invers acelor de ceasornic dinspre est.

CUVINTE CHEIE

raster, teren, aspect, pantă, curbură

REZUMAT

r.pantă.aspect
r.pantă.aspect --Ajutor
r.pantă.aspect [-a] elevație=nume [pantă=nume] [aspect=nume] [format=şir]
[precizie=şir] [pcurvatura=nume] [tcurvatura=nume] [dx=nume] [dy=nume]
[dxx=nume] [vesel=nume] [dxy=nume] [zscale=pluti] [min_slope=pluti] [-suprascrie]
[--ajutor] [-prolix] [-liniştit] [-ui]

Steaguri:
-a
Nu aliniați regiunea curentă la harta raster al cotei

--sobrescrie
Permiteți fișierelor de ieșire să suprascrie fișierele existente

--Ajutor
Imprimați rezumatul utilizării

--verbos
Ieșire modulară

--Liniște
Ieșire silențioasă a modulului

--ui
Forțați lansarea dialogului GUI

parametri:
elevație=nume [necesar]
Numele hărții raster a cotei introduse

pantă=nume
Nume pentru harta raster a pantei de ieșire

aspect=nume
Nume pentru harta raster a aspectului de ieșire

format=şir
Format pentru raportarea pantei
Opțiuni: grade, la sută
Mod implicit: grade

precizie=şir
Tipul de aspect de ieșire și hărți de pantă
Opțiuni: CELULA, FCELL, DCELL
Mod implicit: FCELL

pcurvatura=nume
Nume pentru harta raster a curburii profilului de ieșire

tcurvatura=nume
Nume pentru harta raster cu curbură tangențială de ieșire

dx=nume
Nume pentru harta raster cu derivată parțială de ordinul întâi dx (pantă EW).

dy=nume
Nume pentru harta raster cu derivate parțiale de ordinul întâi dy (panta NS).

dxx=nume
Nume pentru harta raster dxx derivată parțială de ordinul doi

vesel=nume
Nume pentru ieșire hartă raster cu derivate parțiale de ordinul doi

dxy=nume
Nume pentru harta raster dxy derivată parțială de ordinul doi

zscale=pluti
Factorul multiplicativ pentru a converti unitățile de cotă în unități orizontale
Mod implicit: 1.0

min_slope=pluti
Valoarea minimă a pantei (în procente) pentru care este calculat aspectul
Mod implicit: 0.0

DESCRIERE

r.pantă.aspect generează hărți raster de pantă, aspect, curbură și prima și a doua
ordonați derivatele parțiale dintr-o hartă raster a valorilor adevărate de altitudine. Utilizatorul trebuie
specificați intrarea elevație hartă raster și cel puțin o hărți raster de ieșire. Utilizatorul poate
specificați de asemenea format pentru pantă (grade, procente; implicit=grade) și zscale:
factor multiplicativ pentru a converti unitățile de cotă în unități orizontale; (implicit 1.0).

elevație Harta raster de intrare specificată de utilizator trebuie să conțină valori reale de altitudine,
nu date redimensionate sau clasificate. Dacă valorile cotei sunt în alte unități decât în
unități orizontale, acestea trebuie convertite în unități orizontale folosind parametrul zscale.
In GRASS GIS 7, vertical de unităţi sunt nu asumat la be de metri Orice mai mult. Pentru exemplu, if
atât ta vertical și orizontală de unităţi sunt picioare, parametru zscale nu trebuie folosit.

aspect Harta raster de ieșire indică direcția în care sunt îndreptate pantele. Aspectul
categoriile reprezintă numărul de grade de est. Fișierele de categorii și de tabel de culori sunt, de asemenea
generat pentru harta raster de aspect. Categoriile de aspect reprezintă numărul de grade de
est și cresc în sens invers acelor de ceasornic: 90 de grade este nord, 180 este vest, 270 este sud
360 este Est.
Notă: Aceste valori pot fi transformate în valori azimut (0 este nord, 90 este est etc.)
folosind r.mapcalc:
# convertiți unghiurile de la CCW la nord în sus
r.mapcalc "azimuth_aspect = (450 - ccw_aspect) % 360"

Aspectul nu este definit pentru panta egală cu zero. Astfel, cele mai multe celule cu un foarte mic
panta ajunge sa aiba categoria 0, 45, ..., 360 in aspect ieșire. Este posibil să se reducă
părtinirea în aceste direcții prin filtrarea aspectului în zonele în care se află terenul
aproape plat. O opțiune min_slope poate fi folosit pentru a specifica panta minimă pentru care aspect
este calculat. Aspectul pentru toate celulele cu panta min_slope este setat la zero (nu există date).

pantă Harta raster de ieșire conține valorile pantei, exprimate în grade de înclinare de la
orizontalul dacă formatEste aleasă opțiunea =grade (prestabilită), iar în procente de creștere dacă
format= optiunea procentuala este aleasa. Sunt generate fișiere de categorii și tabel de culori.

Curbururile de profil și tangențiale sunt curburele în direcția pantei celei mai abrupte
şi respectiv în direcţia tangentei conturului. Curbururile sunt exprimate ca
1/metru, de exemplu, o curbură de 0.05 corespunde unei raze de curbură de 20 m. Convex
valorile formei sunt pozitive, iar valorile formei concave sunt negative.

Exemplu DEM

Pantă (grad) din exemplu DEM Aspect (grad) din exemplu DEM

Curbură tangenţială (m-1) din exemplul DEM Curbura profilului (m-1) din exemplul DEM

Pentru unele aplicații, utilizatorul va dori să folosească o hartă raster reclasificată a pantei care
grupează valorile pantei în intervale de pantă. Acest lucru se poate face folosind r.reclasare. Un exemplu de
o reclasificare utilă este prezentată mai jos:
categorie etichete categorie
(în grade) (în procente)
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 mai mult
Următorul tabel de culori funcționează bine cu cele de mai sus
reclasificare.
categoria roșu verde albastru
+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

NOTE

Pentru a vă asigura că harta de cotă raster nu este reeșantionată inadecvat, setările pentru
regiunea curentă sunt ușor modificate (doar pentru execuția programului): the
rezoluția este setată să se potrivească cu rezoluția hărții raster de cotă și cu marginile
regiunea (adică nordul, sudul, estul și vestul) sunt deplasate, dacă este necesar, pentru a se alinia de-a lungul
marginile celor mai apropiate celule din harta cotelor. Dacă utilizatorul dorește cu adevărat rasterul
harta altitudinii reeșantionată la rezoluția actuală a regiunii, the -a ar trebui specificat steag.

Masca curentă este ignorată.

Algoritmul folosit pentru a determina panta și aspectul folosește o vecinătate 3x3 în jurul fiecărei celule
în harta raster al cotei. Astfel, nu este posibil să se determine panta și aspectul pt
celulele adiacente marginilor din stratul hartă de cotă. Aceste celule sunt atribuite a
valoarea „pantă zero” (categoria 0) atât în ​​hărțile raster cu panta, cât și în hărțile de aspect.

Formula lui Horn este folosită pentru a găsi derivatele de ordinul întâi în direcțiile x și y.

Numai când se utilizează modele de cote întregi, aspectul este părtinit în 0, 45, 90, 180, 225,
direcțiile 270, 315 și 360; adică, distribuția categoriilor de aspect este foarte neuniformă,
cu vârfuri la 0, 45,..., 360 categorii. Când lucrați cu cotă în virgulă flotantă
modele, nu există o astfel de părtinire de aspect.

EXEMPLE

Calcul of pantă, aspect, profil și tangenţial curbură
În acest exemplu, o hartă de pantă, aspect, profil și curbură tangențială sunt calculate din an
Hartă raster de altitudine (set de date eșantion din Carolina de Nord):
g.regiune raster=altitudine
r.slope.aspect elevation=pantă de cotă=pantă aspect=aspect pcurvature=pcurv tcurvature=tcurv
# setați tabele de culori frumoase pentru hărți raster de ieșire
r.culori -n hartă=pantă culoare=sepia
r.colors map=aspect color=aspectcolr
r.colors map=pcurv culoare=curbatură
r.colors map=tcurv culoare=curbatură

Figura: Hartă raster de pantă, aspect, profil și curbură tangențială (Carolina de Nord
set de date)

Clasificare of major aspect traseu in busolă orientare
În exemplul următor (pe baza setului de date eșantion din Carolina de Nord) generăm mai întâi
harta aspectului standard (orientată în sens invers din partea estică), apoi convertiți-o în orientarea busolă,
și, în cele din urmă, clasificați patru direcții de aspect majore (N, E, S, V):
g.regiune raster=altitudine -p
# generați o hartă de aspect cu orientare CCW
r.slope.aspect elevation=elevation aspect=myaspect
# generează orientarea busolei și clasifică patru direcții majore (N, E, S, V)
r.mapcalc "aspect_4_directions = eval( \\
compas=(450 - myaspect ) % 360, \\
if(busola >=0. && busola < 45., 1) \\
+ if(busola >=45. && busola < 135., 2) \\
+ if(busola >=135. && busola < 225., 3) \\
+ if(busola >=225. && busola < 315., 4) \\
+ if(busolă >=315., 1) \\
)"
# atribuiți etichete de text
r.category aspect_4_directions separator=virgulă reguli=- << EOF
1, nord
2, est
3, sud
4, vest
EOF
# atribuiți tabelul de culori
r.culori aspect_4_directions reguli=- << EOF
1 253,184,99
2 178,171,210
3 230,97,1
4 94,60,153
EOF
Harta aspect clasificată în patru direcții principale ale busolei (subsetul mărit este afișat)

REFERINȚE

· Horn, BKP (1981). Deal Umbrire și il Reflexia Hartă, Proceedings of the
IEEE, 69(1): 14-47.

· Mitasova, H. (1985). Cartografic aspecte of calculator suprafaţă modelare. Dr.
teză. Universitatea Tehnică Slovacă, Bratislava

· Hofierka, J., Mitasova, H., Neteler, M., 2009. Geomorfometrie in GRASS GIS. În:
Hengl, T. și Reuter, HI (Eds). Geomorfometrie: Concepte, Software-ul,
Aplicații. Evoluții în știința solului, voi. 33, Elsevier, 387-410 p.,
http://www.geomorphometry.org

Utilizați r.slope.aspectgrass online folosind serviciile onworks.net