13.1. Bloque 3D.

Crear un Bloque 3D a partir de DEM y dibujar el mar y la línea de costa en perspectiva.

Se explica como crear un bloque 3-D, pintada con una paleta de colores y con iluminación a partir de un DEM, agregar un frente (o fachada) y ejes. Además se dibujar el mar y la línea de costa en perspectiva. La entrada 5.1 muestra un mapa 2D.

Script

0. Definir Variables: En los scripts anteriores se definían 2 variables (REGION y PROJ). Para hacer un bloque 3D es necesario definir 3 variables nuevas:

• REGION3D: incluye a REGION y agrega el rango de valores de altura (-10000/10000 en este  ejemplo).
• persp: indica la perspectiva o punto de vista desde la cual se observará el bloque 3D. Se define con el acimut (grados contados en sentido horario desde el norte) y la elevación sobre el horizonte. Se puede agregar un tercer valor que indica a que altura el material en 2D se dibujará. Para este ejemplo se eligió eligió 160/30.
• PROZ: indica la escala vertical. En este ejemplo corresponde a 4 cm.

REM    Region: Cuyo
    SET    REGION=-72/-64/-35/-30
    SET    REGION3D=%REGION%/-10000/10000

REM    Proyeccion Mercator (M)
    SET    PROJ=M14c
    SET    PROZ=4c
    SET    persp=160/30

1. Procesamiento Datos: Uso los mismos comandos (y grilla GMRTv3.1.grd) que en la entrada 5.1 para recortar la grilla, crear una paleta de colores y una grilla para el sombreado. Para el acimut de la posición del sol se eligió la misma que %persp% (-A160).

gmt psxy -R%REGION3D% -J%PROJ% -JZ%PROZ% -p%persp% -T -K -P > %OUT%
gmt grdcut GMRTv3_1.grd -G%CUT% -R%REGION% -fg
gmt makecpt -Cdem4 -T0/7000/250 -Z >%color%
gmt grdgradient %CUT% -A160 -G%SHADOW% -Ne0.5

2. Bloque 3D - Superficie: El módulo grdview permite crear una vista 3D elaborado a partir de las variables (-R, -J, -JZ y -p). El argumento -Qi300 indica que dibuje una imagen (con una resolución de 300 dpi) a partir del DEM y con la paleta de colores (-C%color%) y la grilla para el sombreado (-I%SHADOW%). El área en en negro corresponde al Pacífico (la paleta de colores fue definida para áreas emergidas).

gmt grdview -R -J -JZ -O -K %CUT% -p -Qi300 -C%color% -I%SHADOW% >> %OUT%

3. Bloque 3D - Fachada: El argumento -N-10000 dibuja un plano en ese nivel (en metros). Opcionalmente se puede pintar la fachada (o frente) del bloque (+glightgray) y definir el estilo de la pluma (pen attributes) utilizada para las líneas (-Wf0.5).

gmt grdview -R -J -JZ -O -K %CUT% -I%SHADOW% -p%persp% -C%color% -Qi300 >> %OUT% -N-10000+glightgray -Wf0.5

4. Bloque 3D - Ejes: Con el argumento -B se agregan los ejes al bloque 3D. -Baf dibuja los ejes X y Y en la base del bloque (definido por la variable REGION3D). Para definir el eje Z se sigue la misma sintaxis (-Bzaf) y se puede agregar un título (+l"Altura (m)"). El argumento -BnSeWZ1 indica que los ejes (en mayúscula) incluyan las anotaciones. El 1 indica en que esquina se dibuja el eje vertical (1: inferior izquierda, 2: inf derecha, 3: superior der, 4; sup izq). Si no se agrega, el programa define en que esquina ubicarlo. "+b" dibuja un recuadro al bloque con la dimensiones definidas por REGION3D.

gmt grdview -R -J -JZ -O -K %CUT% -I%SHADOW% -p%persp% -C%color% -Qi300 >> %OUT% -N-10000+glightgray -Wf0.5 -BnSwEZ+b -Baf -Bzaf+l"Altura (m)"

5. Mar y Linea de Costa: pscoast permite dibujar visto en perspectiva. Es necesario indicar la altura del plano del mar (-p%perp%/0). -A indica se dibujen los rasgos entre nivel 0 y 1 (excluyendo lagos, islas en lagos, y lagos en islas en lagos que pueden estar a otra altura; ver GSHHS).

gmt pscoast -R -J -JZ -O -K -p%perp%/0 -Df -Sdodgerblue2 -A0/0/1 >> %OUT%
gmt pscoast -R -J -JZ -O -K -p%perp%/0 -Df -W1/0.3,black >> %OUT%