6.2. Mapa Aspecto.

Crear mapa de aspecto (orientación de pendientes) a partir de un Modelo de Elevación Digital (DEM), crear una paleta de colores (cpt) y agregar una escala de color.

Mapa de orientación de pendientes (aspecto) de la región de cuyo elaborado a partir de la grilla GMRT. Incluye datos del IGN: ejidos urbanos (áreas verdes), centros poblados (puntos negros), vías férreas (líneas blancas), límites políticos (líneas negras).

Datos:
Los mapas de aspecto se realizan a partir de Modelos de Elevación Digital (DEM) batimétricos-topográficos. En Base de Datos se mencionan distintas fuentes y como descargar los DEM. Para este ejemplo, utilizo la grilla GMRTv3_1.grd descargada del portal GMRT, que es la misma utilizada en ejemplos 5.

Script

Definir Variables

1. Región: Mapa comprendido entre los meridianos 72° y 64° W y los paralelos 35° y 30°S.

rem    Region: Argentina
    SET    REGION=-72/-64/-35/-30

2. Proyección: Mapa con proyección Mercator de 15 cm de ancho.

rem    Proyeccion Mercator (M)
    SET    PROJ=M15c

Hacer Mapa:

3. Mapa Base: Usamos comandos previos para hacer un mapa que incluya datos del IGN, líneas de costas, escala gráfica, pintar áreas océanos y hacer un marco al mapa explicado en entradas previas. El mapa creado es igual al del ejemplo 6.

REM    Datos Instituto Geografico Nacional (IGN)
    gmt psxy -R -J -O -K "E:\Facultad\Datos_Geofisicos\IGN\1_GMT\limites_politico_administrativos.gmt" -Wthin >> %OUT% 
    gmt psxy -R -J -O -K "E:\Facultad\Datos_Geofisicos\IGN\1_GMT\003_Red_Ferroviaria.gmt" -Wthin,white -Am >> %OUT%
    gmt psxy -R -J -O -K "E:\Facultad\Datos_Geofisicos\IGN\1_GMT\005_Centros_Poblados.gmt" -Sc0.04 -Gblack >> %OUT%
    gmt psxy -R -J -O -K "E:\Facultad\Datos_Geofisicos\IGN\1_GMT\010_Ejidos_Urbanos.gmt" -Wfaint -Ggreen >> %OUT%  

REM    Dibujar frame
    gmt psbasemap -R -J -O -K -Baf >> %OUT%

rem    Pintar areas húmedas: Oceanos (-S)
    gmt pscoast -R -J -O -K -Df -Sdodgerblue2  >> %OUT%

REM    Dibujar Linea de Costa (W1)
    gmt pscoast -R -J -O -K -Df -W1/faint >> %OUT%

REM    Dibujar Escala en el mapa centrado en -Lg Lon0/Lat0, calculado en meridiano (+c), ancho (+w), elegante(+f), unidad arriba de escala (+l), unidad con los valores (+u). n: Centrado en el 88% del eje X y 7.5% del eje Y.
    gmt psbasemap -R -J -O -K -Ln0.88/0.075+c-32:00+w100k+f+l   >> %OUT%

4. Gradiente: Usamos grdgradient para calcular el gradiente de la superficie del terreno a partir de la grilla topográfica (GMRTv3_1.grd). Con la opción -Da creamos una nueva grilla (%CUT%) con el acimut (medido en sentido horario desde el norte) de la pendiente. La opción -Da indica que los acimut corresponden a la dirección pendiente-abajo. (Utilizar -D para la dirección pendiente-arriba). Con -fg indicamos que la grilla está en coordenadas geográficas.

REM    Calcular Grilla de orientacion de Pendientes (aspecto). -D(a): direccion con pendiente-arriba (abajo).
    gmt grdgradient GMRTv3_1.grd -Da -G%CUT% -R -fg

5. Paleta de Color: Antes de dibujar el mapa, es necesario crear una paleta de colores (cpt) para las distintas orientaciones de la pendiente del terreno.
Para este ejemplo utilizamos makecpt para crear una nueva cpt (%color%) a partir del cpt maestro cyclic (-Ccyclic) para el rango de 0 a 360°con intervalos cada 5° (-T0/360/5) y valores continuos (-Z).

REM    Crear Paleta de Colores. Paleta Maestra (-C), Definir rango (-Tmin/max/intervalo), CPT continuo (-Z)
    gmt makecpt -Ccyclic -T0/360/5 -Z >%color%

6. Imagen: Con grdimage se crea una imagen a partir de la grilla a partir de la grilla %CUT% (punto 4) y la cpt %color% (punto 5) con el siguiente comando:

REM    Crear Imagen a partir de grilla con sombreado
    gmt grdimage -R -J -O -K %CUT% -C%color% >> %OUT%

7. Escala de Color: Con psscale se agrega una escala de color a partir del cpt previo (-C%color%). Con la opción -Dx indicamos su posición en coordenadas del gráfico (-Dx0/-1.0) y sus dimensiones (15 cm de largo y 0,618 de ancho; +w15/0.618c). Para dibujar una escala horizontal, en lugar de vertical, se agrega una "h" al final de la opción -D. Con la opción -B se agrega anotaciones cada 30°, con marcas cada 15° y con la leyenda "Orientación pendiente (°)".

REM    Agregar escala vertical a partir de CPT (-C). Posición (x,y) +wlargo/ancho. Anotaciones (-Ba). Leyenda (+l). 
    gmt psscale -O -K -Dx0/-1.0+w15/0.618ch -C%color% -Ba30f15:"Orientaci\363n pendiente(\232)": >> %OUT%