Mapas Geológicos: 13.6. Bloque 3D. Animación.

martes, 28 de mayo de 2019

13.6. Bloque 3D. Animación.

Animación mostrando un bloque 3D de la profundidad del moho en Sudamérica.

Esta animación permite mostrar las zonas profundas del moho ocultas por la superficie. Se arma combinando distintos bloques 3D en un gif.

Script

0. Datos: Se utilizó la grilla de profundidad del moho (Más información en Base de Datos):

Profundidad Moho: Tabla de datos depthtomoho.xyz (aquí comprimido en archivo zip) del Modelo CRUST 1.0. Profundidad de km (positivo hacia arriba).

1. Crear grilla: Se utiliza surface para crear la grilla Moho.nc (ver. 8.5 para más detalles).

gmt surface "depthtomoho.xyz" -I6m -R%REGION% -G"Moho.nc" -fg -r

2. Crear CPT personalizada: con makecpt se crea la paleta de color (CPT) que varía gradualmente entre lightred y lightblue entre los valores de 75 y 0 km de profundidad, respectivamente (ver 5.5 para más detalles).

gmt makecpt -Clightred,lightblue -T-75/0 -Z > %color%

3. Efecto Sombreado: Con grdgradient se crea una grilla para hacer el efecto de sombreado (ver 5.2 para más detalles). La

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

4. Bloques 3D: Para este ejemplo se crearon 35 figuras. La primera corrresponde a un bloque 3D entero, mientras que para el resto fue nivelado a distintas profundidades entre 10 y 74 km (cada 2 km).

A. Bloque 3D entero: Con grdview se crea el bloque 3D (mas detalles en 13.1).

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

B. Bloque 3D Nivelado:
Se define la variable Above que es el valor de profundidad (en positivo) a la cual se quiere nivelar el bloque 3D (también se utiliza para definir el nombre de la imagen de salida). Luego, con grdclip se nivela la grilla con el siguiente argumento:

-Sa-%above%/-%Above%: Setea todos los datos por encima del valor de la variable %Above% a ese mismo valor. (Se agrega el símbolo "-" porque los datos están en negativo).

SET Above=66

gmt grdclip "Moho.nc" -G%CUT% -Sa-%Above%/-%Above%

Nivelado a 32 km de profundidad.

Nivelado a 20 km de profundidad.

Nivelado a 44 km de profundidad.

Nivelado a 56 km de profundidad.

5. Animación: La herramienta convert del programa GraphicsMagik permite combinar los 35 bloques 3D en un gif con las siguientes opciones:

-delay: intervalo entre imágenes en centésimas de segundos.
-loop: cuantas veces se repite (0 = infinitamente).

gm convert -delay 10 -loop 0 3d_*.png 3d_Moho_RedBlue_10.gif

Alternativamente se puede utilizar el sitio gifmaker.