r.slope.aspectgrass: en línea en la nube

PROGRAMA:

NOMBRE

r.pendiente.aspecto - Genera mapas ráster de pendiente, aspecto, curvaturas y parciales.
derivados de un mapa ráster de elevación.
El aspecto se calcula en sentido antihorario desde el este.

PALABRAS CLAVE

ráster, terreno, aspecto, pendiente, curvatura

SINOPSIS

r.pendiente.aspecto
r.pendiente.aspecto --ayuda
r.pendiente.aspecto [-a] elevación=nombre [pendiente=nombre ] [aspecto=nombre ] [formato=cadena]
[precisión=cadena] [curvatura=nombre ] [tcurvatura=nombre ] [dx=nombre ] [dy=nombre ]
[dxx=nombre ] [dyy=nombre ] [dxy=nombre ] [escala z=flotar] [pendiente_min=flotar] [-exagerar]
[-ayuda] [-verboso] [-tranquilo] [-ui]

Banderas
-a
No alinee la región actual con el mapa de elevación ráster

--Sobrescribir
Permitir que los archivos de salida sobrescriban los archivos existentes

--ayuda
Resumen de uso de impresión

--verboso
Salida del módulo detallado

--tranquilo
Salida de módulo silencioso

--ui
Forzar el inicio del cuadro de diálogo GUI

parámetros:
elevación=nombre [requerido]
Nombre del mapa ráster de elevación de entrada

pendiente=nombre
Nombre del mapa ráster de pendiente de salida

aspecto=nombre
Nombre del mapa ráster de aspecto de salida

formato=cadena
Formato para informar la pendiente
Opciones: grados, por ciento
Por defecto: grados

precisión=cadena
Tipo de mapas de aspecto y pendiente de salida
Opciones: CELDA, FCEL, DCELL
Por defecto: FCCEL

curvatura=nombre
Nombre del mapa ráster de curvatura del perfil de salida

tcurvatura=nombre
Nombre del mapa ráster de curvatura tangencial de salida

dx=nombre
Nombre del mapa ráster de derivada parcial de primer orden dx (pendiente EW) de salida

dy=nombre
Nombre del mapa ráster de derivada parcial de primer orden dy (pendiente NS) de salida

dxx=nombre
Nombre del mapa ráster dxx de derivada parcial de segundo orden de salida

dyy=nombre
Nombre del mapa ráster dyy de derivada parcial de segundo orden de salida

dxy=nombre
Nombre del mapa ráster dxy derivado parcial de segundo orden de salida

escala z=flotar
Factor multiplicativo para convertir unidades de elevación en unidades horizontales.
Por defecto: 1.0

pendiente_min=flotar
Valor mínimo de pendiente (en porcentaje) para el que se calcula el aspecto
Por defecto: 0.0

DESCRIPCIÓN

r.pendiente.aspecto genera mapas ráster de pendiente, aspecto, curvaturas y primero y segundo
ordenar derivadas parciales de un mapa ráster de valores de elevación reales. El usuario debe
especificar la entrada elevación mapa ráster y al menos un mapa ráster de salida. El usuario puede
también especifique el formato para pendiente (grados, porcentaje; predeterminado = grados), y el escala z:
factor multiplicativo para convertir unidades de elevación en unidades horizontales; (predeterminado 1.0).

Los elevación El mapa ráster de entrada especificado por el usuario debe contener valores de elevación reales,
no datos reescalados o categorizados. Si los valores de elevación están en otras unidades que en el
unidades horizontales, deben convertirse a unidades horizontales mediante el parámetro escala z.
In HIERBA GIS 7, vertical unidades que están no ficticio a be metros any más. Para ejemplo, if
ambas a tu manera vertical y horizontal unidades que están pies, parámetro escala z no debe utilizarse.

Los aspecto El mapa ráster de salida indica la dirección en la que miran las pendientes. El aspecto
las categorías representan el número de grados del este. Los archivos de tablas de categorías y colores también
generado para el mapa ráster de aspecto. Las categorías de aspecto representan el número de grados de
al este y aumentan en sentido antihorario: 90 grados es el norte, 180 es el oeste, 270 es el sur
360 es el este.
Nota: Estos valores se pueden transformar en valores de acimut (0 es el norte, 90 es el este, etc.)
usando r.mapcalc:
# convertir ángulos de CCW a norte hacia arriba
r.mapcalc "azimuth_aspect = (450 - ccw_aspect)% 360"

El aspecto no está definido para una pendiente igual a cero. Por lo tanto, la mayoría de las células con un tamaño muy pequeño
pendiente terminan teniendo categoría 0, 45, ..., 360 en aspecto producción. Es posible reducir
el sesgo en estas direcciones al filtrar el aspecto en áreas donde el terreno es
casi plano. Una opcion pendiente_min se puede utilizar para especificar la pendiente mínima para qué aspecto
se calcula. El aspecto de todas las celdas con pendiente. pendiente_min se establece a nulo (sin datos).

Los pendiente El mapa ráster de salida contiene valores de pendiente, expresados ​​en grados de inclinación desde
la horizontal si formato= Se elige la opción de grados (la predeterminada), y en porcentaje de aumento si
formato= se elige la opción de porcentaje. Se generan archivos de tablas de categorías y colores.

Las curvaturas de perfil y tangenciales son las curvaturas en la dirección de la pendiente más pronunciada.
y en la dirección de la tangente del contorno respectivamente. Las curvaturas se expresan como
1 / metro, por ejemplo, una curvatura de 0.05 corresponde a un radio de curvatura de 20 m. Convexo
los valores de forma son positivos y los valores de forma cóncava son negativos.

Ejemplo DEM

Pendiente (grados) del ejemplo DEM Aspecto (grados) del ejemplo DEM

Curvatura tangencial (m-1) del ejemplo DEM Curvatura del perfil (m-1) del ejemplo DEM

Para algunas aplicaciones, el usuario deseará utilizar un mapa ráster de pendiente reclasificado que
agrupa los valores de pendiente en rangos de pendiente. Esto se puede hacer usando r.reclasificar. Un ejemplo de
a continuación se ofrece una reclasificación útil:
etiquetas de categoría de rango de categoría
(en grados) (en porcentaje)
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% y más
La siguiente tabla de colores funciona bien con la anterior.
reclasificación.
categoría rojo verde azul
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

NOTAS

Para asegurarse de que el mapa de elevación ráster no sea muestreado de manera inapropiada, la configuración de
la región actual se modifican ligeramente (solo para la ejecución del programa):
La resolución está configurada para que coincida con la resolución del mapa ráster de elevación y los bordes de la
región (es decir, el norte, sur, este y oeste) se desplazan, si es necesario, para alinearse a lo largo de
bordes de las celdas más cercanas en el mapa de elevación. Si el usuario realmente quiere el ráster
mapa de elevación remuestreado a la resolución de la región actual, el -a se debe especificar la bandera.

La máscara actual se ignora.

El algoritmo utilizado para determinar la pendiente y el aspecto utiliza una vecindad de 3x3 alrededor de cada celda.
en el mapa de elevación ráster. Por lo tanto, no es posible determinar la pendiente y el aspecto para
las celdas adyacentes a los bordes en la capa del mapa de elevación. A estas celdas se les asigna un
valor de "pendiente cero" (categoría 0) en los mapas ráster de pendiente y aspecto.

La fórmula de Horn se usa para encontrar las derivadas de primer orden en las direcciones xey.

Solo cuando se utilizan modelos de elevación de enteros, el aspecto está sesgado en 0, 45, 90, 180, 225,
Direcciones 270, 315 y 360; es decir, la distribución de las categorías de aspectos es muy desigual,
con picos en 0, 45, ..., 360 categorías. Cuando se trabaja con elevación de punto flotante
modelos, no se produce tal sesgo de aspecto.

EJEMPLOS

Cálculo of Pendiente, aspecto, perfiles y tangencial curvatura
En este ejemplo, un mapa de pendiente, aspecto, perfil y curvatura tangencial se calcula a partir de un
mapa ráster de elevación (conjunto de datos de muestra de Carolina del Norte):
g.region raster = elevación
r.lope.aspect elevación = elevación pendiente = pendiente aspecto = aspecto pcurvatura = pcurv tcurvature = tcurv
# establecer tablas de colores agradables para los mapas ráster de salida
r.colors -n map = color de pendiente = sepia
mapa de r.colors = color de aspecto = color de aspecto
r.colors map = pcurv color = curvatura
r.colors map = tcurv color = curvatura

Figura: Mapa ráster de pendiente, aspecto, perfil y curvatura tangencial (Carolina del Norte
conjunto de datos)

Clasificación of gran aspecto direcciones in brújula La orientación
En el siguiente ejemplo (basado en el conjunto de datos de muestra de Carolina del Norte), primero generamos
el mapa de aspecto estándar (orientado en sentido antihorario desde el este), luego conviértalo a la orientación de la brújula,
y finalmente clasifique cuatro direcciones de aspecto principales (N, E, S, W):
g.region ráster = elevación -p
# generar mapa de aspecto con orientación CCW
r.slope.aspect elevación = elevación aspecto = myaspect
# generar la orientación de la brújula y clasificar cuatro direcciones principales (N, E, S, W)
r.mapcalc "aspect_4_directions = eval (\\
brújula = (450 - myaspect)% 360, \\
if (brújula> = 0. && brújula <45., 1) \\
+ if (brújula> = 45. && brújula <135., 2) \\
+ if (brújula> = 135. && brújula <225., 3) \\
+ if (brújula> = 225. && brújula <315., 4) \\
+ si (brújula> = 315., 1) \\
)"
# asignar etiquetas de texto
r.category aspect_4_directions separator = comma rules = - << EOF
1, norte
2, este
3, sur
4, oeste
EOF
# asignar tabla de colores
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 de aspecto clasificado en cuatro direcciones principales de la brújula (se muestra un subconjunto ampliado)

Referencias

· Horn, BKP (1981). Hill Sombreado y los reflectancia Mapa, Actas de la
IEEE, 69(1): 14-47.

· Mitasova, H. (1985). Cartográfico aspectos of computadora superficie modelado. Doctorado
tesis. Universidad Técnica Eslovaca, Bratislava

· Hofierka, J., Mitasova, H., Neteler, M., 2009. Geomorfometria in HIERBA SIG. En:
Hengl, T. y Reuter, HI (Eds), Geomorfometría: Conceptos, software
Aplicaciones. Desarrollos en la ciencia del suelo, vol. 33, Elsevier, 387-410 págs.
http://www.geomorphometry.org

Utilice r.slope.aspectgrass en línea utilizando los servicios de onworks.net