6.3. Mapa combinado Aspecto + Pendiente

Mapa combinado de Pendiente y Aspecto.

Introducción: En los ejemplos anteriores se mostró como elaborar mapas de pendientes (6.1) y aspecto (6.2) por separado. En este ejemplo se explica como elaborar un mapa que refleje al mismo tiempo ambos parámetros según lo descripto por Brewer y Marlow (1993). En la figura resultante, los colores indican la dirección (aspecto). Si el color es más fuerte indica pendientes más fuertes, y se vuelve gris a medida que la pendiente es más suave.

 

Los pasos a seguir son: 1) crear ambas grillas por separado, 2) clasificarlas según los tipos y orientaciones de las pendientes, 3) combinarlas las grillas clasificadas, 4) dibujar la grilla.

Script 1 (pendiente en porcentaje)
Script 2 (pendiente en grados)

0. Datos:

• Grilla GMRTv3_1.grd descargada del portal GMRT.
• Paleta de Color Brewer_Aspect-Slope.cpt. Paleta creada con los colores definidos por Brewer y Marlow (1993).

SET    DEM="GMRTv3_1.grd"

1. Grillas de pendiente y aspecto: grdgradient crea las grillas de pendiente (%CUT%) y aspecto (%CUT2%). Las pendientes se calculan en porcentaje (%). Ver 6.1 y 6.2 para más detalles.

gmt grdgradient %DEM% -D  -S%CUT%  -R -fg

gmt grdmath %CUT% 100 MUL = %CUT%
gmt grdgradient %DEM% -Da -G%CUT2% -R -fg

2. Clasificar grillas: grdclip permite clasificar las grillas, es decir modificar los valores de z de las grillas según los siguientes argumentos:

• -Sblow/below: reemplaza valores menores a los low a below.
• -Silow/high/between:: reemplaza valores intermedios a between.
• -Sahigh/above: reemplaza valores mayores a high a above.

Nota: Es importante remarcar que al clasificar las grillas, éstas se convierte en grillas categóricas. Es decir que sus valores representan una categoría.

A. Clasificar pendientes en 4 clases:

Intervalo de inclinación	Clase	Pendiente
< 5%	10	Plana
5 - 20%	20	Suave
20 - 40%	30	Moderada
> 40%	40	Fuerte

gmt grdclip %CUT% -G%CUT% -Sb5/10 -Si5/20/20  -Si20/40/30  -Sa40/40 -V

Nota: Los valores de pendientes de Brewer y Marlow están expresados en porcentaje y los convertí a grados.

B. Clasificar orientaciones en 8 clases:

Intervalo de Dirección	Clase	Orientación
< 22,5°	1	N
22,5 -  67,5°	2	NE
67,5 - 112,5 °	3	E
112,5 - 157,5°	4	SE
157,5 - 202,5°	5	S
202,5 - 247,5 °	6	SO
247,5 - 292,5 °	7	O
292,5 - 337,5 °	8	NO
> 337,5° (hasta 360°)	1	N

gmt grdclip %CUT2% -G%CUT2% -Sb22.5/1 -Si22.5/67.5/2 -Si67.5/112.5/3 -Si112.5/157.5/4 -Si157.5/202.5/5 -Si202.5/247.5/6 -Si247.5/292.5/7 -Si292.5/337.5/8 -Sa337.5/1 -V

3. Combinar grillas:  Las 2 grillas se suman y guardan con el nombre %CUT%. Al combinar las grillas se obtiene una nueva con valores donde la decena (10, 20, 30, 40) indica el tipo de pendiente y la unidad (1 a 8) la orientación.

gmt grdmath %CUT% %CUT2% ADD = %CUT%

4. Paleta de Colores: La CPT Brewer_Aspect-Slope.cpt tiene un color definido para cada uno de los valores de la grilla combinada.

5. Dibujar Mapa: con grdimage se dibuja el mapa. El argumento -nn indica que se utilice como método de interpolación el vecino más cercano (nearest-neighbor) que se utiliza para datos categóricos (como en este caso).

gmt grdimage -R -J -O -K %CUT% -C"Brewer_Aspect-Slope.cpt" >> %OUT% -nn


6. Mapas con inclinación en grados: En estos ejemplos muestro los mapas rehechos con las pendientes calculadas en grados (°).
A. Clasificación pendientes:

Intervalo de inclinación	Clase	Pendiente
< 3°	10	Plana
3 - 12°	20	Suave
12 - 22°	30	Moderada
> 22°	40	Fuerte

gmt grdclip %CUT% -G%CUT% -Sb3/10 -Si6/12/20 -Si12/20/30 -Sa20/40 -V

B. Clasificación de pendiente:

Intervalo de inclinación	Clase	Pendiente
< 1°	10	Plana
1 - 5°	20	Suave
5 - 10°	30	Moderada
> 10°	40	Fuerte

gmt grdclip %CUT% -G%CUT% -Sb1/10 -Si1/5/20  -Si5/10/30  -Sa10/40 -V

 

Referencias:

• Brewer, C.A. & Marlow, K.A. (1993): Color Representation of Aspect and Slope simultaneously. Proceedings, Eleventh International Symposium on Computer-Assisted Cartography (Auto-Carto-11), Minneapolis, October/November 1993, pp. 328-337. Link.
• Mapas combinados de pendiente y orientación de laderas. Gisandbeers. Explicación paso a paso en español. Link.