5.1. Mapa Topográfico

Hacer mapa topográfico a partir de un Modelo de Elevación Digital (DEM). Se explica como recortar una grilla, crear una paleta de colores (CPT) y agregar una escala.

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

0. Datos:
Los mapas batimétricos-topográficos se realizan a partir de Modelos de Elevación Digital (DEM). 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.

Script

Definir Variables

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

REM    Region: Cuyo
    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, pintar áreas océanos y hacer un marco al mapa explicado en entradas previas.


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,red -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%

4. Recortar Grilla: La herramienta grdcut nos permite extraer una subregión de una grilla. Para extraer la región (definida en 1) de la grilla GMRTv3_1.grd y guardarla con el nombre de la variable "CUT" se utiliza el siguiente comando. La opción -fg especifica que los datos de la grilla están en coordenadas geográficas.

SET    CUT=temp_%title%.grd
gmt grdcut GMRTv3_1.grd -G%CUT% -R -fg

5. Paleta de Color I: Para colorear el mapa topográfico es necesario crear una paleta de color (CPT), cuyo nombre definimos con la variable "color". Con makecpt podemos editar una CPT maestro (en este caso dem4) para el rango de 0 a 7000 m con intervalos cada 250 m. La nueva CPT se guarda con el nombre de la variable "color".

gmt makecpt -Cdem4 -T0/7000/250 >%color%


6. Imagen: Con grdimage se crea una imagen a partir de la grilla. Para crear una imagen a partir de la grilla "CUT" (ver punto 5) y con la paleta de colores creada (punto 6) se usa el siguiente comando.

gmt grdimage -R -J -O -K %CUT% -C%color% >> %OUT%

7. Escala de Color: Con psscale se agrega una escala de color a partir de la CPT previo (punto 6) a la derecha del mapa (15.5/0) de 11 cm de largo y 0,618 de ancho. La escala tiene anotaciones cada 500 m y con la leyenda "Topografía (m)".


gmt psscale -O -K -Dx15.5/0+w11/0.618c -C%color% -Ba500+l"Topograf\355a (m)" >> %OUT%

8. Paleta de Color II: Editamos la CPT (creado en el punto 6) para obtener una escala de colores continua (-Z).

gmt makecpt -Cdem4 -T0/7000/250 -Z >%color%