16.2. Mapa Topográfico de la Cuenca del Paraná.

Utilizar una máscara para crear un mapa topográfico con colores dentro de la Cuenca del Paraná.

Introducción: En 5.1 a 5.2 se explica como crear un mapa topográfico con colores según una paleta de color (CPT). Sin embargo en ocasiones uno puede querer resaltar un área irregular como una cuenca hidrográfica o un país. Una forma (presentada aquí) es utilizar una máscara para restringir el color a esa área dejando el resto en blanco y con efecto de sombreado.

Datos:

1. Grilla GMRT_v3.3low.grd (ver Anexo 4 para descargar otra grilla).
2. Archivo Cuenca_Parana.txt datos de HydroSHEDS (archivo sa_bas_15s_beta.shp; ver Base de Datos). 

Script: Para crear la figuras vamos a superponer 2 mapas: el mapa inferior en blanco y el superior en color.
Es importante que para ambos mapas se utilicen los mismos parámetros para crear el efecto de sombreado (-A y -Ne0.5). Ver punto 1 en 5.2 para más detalles.

gmt grdgradient %CUT%  -A270 -G%SHADOW%  -Ne0.5

gmt grdgradient %CUT1% -A270 -G%SHADOW1% -Ne0.5

Ejemplo Cuenca de Paraná:

0. Mapa Base: Dibujar marco del mapa, línea de costa, límites de países y divisiones internar (provincias/estados/departamentos/territorios) y pintar los océanos. Se utilizan los siguientes comandos explicados en 1.1 y 1.2.


gmt psbasemap -R -J -O -K -Baf >> %OUT%
gmt pscoast -R -J -O -K >> %OUT% -Df -W1/0.5
gmt pscoast -R -J -O -K >> %OUT% -Df -N1/0.30 
gmt pscoast -R -J -O -K >> %OUT% -Df -N2/0.25,-.
gmt pscoast -R -J -O -K >> %OUT% -Sdodgerblue2 -C-

1. Mapa inferior: Es un mapa topográfico pintado con una paleta de color blanca y con efecto de sombreado. Se utilizaron los siguientes comandos (explicados en 5.5 1A).

gmt grdcut "GMRTv3_3_low.grd" -R -G%CUT1%
gmt grdinfo %CUT1% -T50 
echo    -5000     white    6400     white    > %color1%
gmt grdgradient %CUT1%  -A270 -G%SHADOW1% -Ne0.5
gmt grdimage -R -J -O -K %CUT1% -C%color1% -I%SHADOW1% >> %OUT%

2. Recortar Grilla Irregular: Antes de dibujar el mapa superior hay que recortar la grilla según la forma de la cuenca del Paraná utilizando una máscara. Se utilizaron los siguientes comandos explicados en 16.1.

SET    CLIP="Cuenca_Parana.txt"
gmt grdmask -R%CUT1% %CLIP% -G%MASK% -NNaN/NaN/1
gmt grdmath %CUT1% %MASK% MUL = %CUT%


3. Mapa superior: Una vez creada la grilla irregular, se dibuja con grdimage. Hay que agregar -Q para que queden transparentes las áreas NaN permitiendo ver el mapa inferior.

gmt makecpt > %color% -Cdem4 -Z -T0/5700/50 

gmt grdgradient %CUT%  -A270 -G%SHADOW%  -Ne0.5

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

gmt psscale -R -J -O -K -C%color% >> %OUT% -I -DJRM+o0.3c/0+w14/0.618c   -Ba+l"Alturas (km)" -W0.001

gmt psxy -R -J -O -K >> %OUT% %CLIP% -W1,red