16.3. Máscara Bloque 3-D Cuenca del Paraná.

Como utilizar una máscara en un bloque 3-D para resaltar la Cuenca del Paraná.

Introducción: En 16.1 y 16.2 se vieron ejemplos para resaltar un área irregular. En este ejemplo se explica como utilizar un bloque 3-D para resaltar la Cuenca del Paraná. Para hacer la figura creamos 2 bloques 3-D y luego agregamos datos culturales.

Datos:

1. Grilla GMRT_v3.3low.grd con valores de altura en metros (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:

0. Variables: Se definen las siguientes variables para hacer la figura (ver 0 en 12.1 para más detalles). A diferencia del 12.1, en este ejemplo se utilizan valores en km para el eje Z (-5 a 7 km).

SET    REGION=-70/-42/-36/-12
SET    REGION3D=%REGION%/-5/7
SET    PROJ=M14c
SET    PROZ=1.5c
SET    persp=160/30
SET    persp3d=%persp%/0

1. Conversión m a km: El texto "=/0.001" agregado al nombre de la grilla (%CUT%, %CUT1%) multiplica la grilla por esa escala. En este ejemplo permite convertir los datos de altura en metros de la grilla original a kilómetros.


2. Bloque inferior:
A. Paleta de Color: Se crea una paleta de color (CPT) azul (dodgerblue2) para las profundidades (-5 a 0 km) y con una transición blanco a negro para las alturas (0 a 6,4 km; ver 5.5 para más detalles).

echo    -5    dodgerblue2     0     dodgerblue2 >  %color1%
echo    0     white         6.4     black       >> %color1%

B. Dibujar bloque: Los siguientes comandos crean una grilla para el efecto de sombreado (grdgradient) y el bloque 3-D (grdview; ver 12.1 para más detalles).

gmt grdgradient %CUT1% -A160 -G%SHADOW1% -Ne0.5
gmt grdview %CUT1%=nf/0.001 -R -J -JZ -O -K -p -Qs -C%color1% -I%SHADOW1% >> %OUT% -N-5+g- -Wf0.5 -BnSwEZ+b -Baf -Bzaf+l"Altura (km)"

 

3. Bloque Superior:
A. Crear máscara y recortar grilla: Se utilizan los siguientes comandos para recortar la grilla según la Cuenca del Paraná aplicando una máscara (ver puntos 2 a 4 de 16.1 para más detalles).

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


B. Crear CPT y Dibujar bloque: De la misma manera que en 1, creamos una CPT de colores para la cuenca del Paraná, una grilla para el efecto de sombreado (%SHADOW%) y un bloque 3-D únicamente con la superficie de la cuenca (sin ejes ni fachada).

gmt makecpt -Cdem4 -T0/7.000/0.250 -Z > %color%
gmt grdgradient %CUT%  -A160 -G%SHADOW%  -Ne0.5
gmt grdview %CUT%=nf/0.001  -R -J -JZ -O -K -p -Qs -C%color%  -I%SHADOW%  >> %OUT%

4. Datos Culturales: Los siguientes comandos permiten dibujar los límites de países (-N1) y provincias/estados/departamentos/etc (-N2), la línea de costa (-W1), el borde de la Cuenca del Paraná (%CLIP%) y agregar la escala de colores debajo (ver 12.2 para más detalles). 

gmt pscoast -R%REGION% -Df -M -N1/ | grdtrack -G%CUT1%=nf/0.001 -sa | psxyz -R%REGION3D% -J -JZ -O -K -p >> %OUT% -W0.3,black
gmt pscoast -R%REGION% -Df -M -N2/ | grdtrack -G%CUT1%=nf/0.001 -sa | psxyz -R%REGION3D% -J -JZ -O -K -p >> %OUT% -W0.25,black,-
gmt pscoast -R -J -JZ -O -K -p%persp3d% -Df >> %OUT% -W1/0.3,black
gmt grdtrack -R%REGION% %CLIP% -G%CUT1%=nf/0.001 -sa | psxyz -R%REGION3D% -J -JZ -O -K -p >> %OUT% -Wthin,blue
gmt psscale -R%REGION% -J -O -K -p%persp% >> %OUT% -C%color% -I -DJCB+o0/0.7c+w14/0.618ch -Ba+l"Altura (km)"