Combinar 2 o mas grillas en una misma figura 3-D.
Introducción: Se explica como combinar en un bloque 3-D distintas grillas que representan distintas superficies. Las superficies se grafican agregando una separación vertical (despiece) de manera que las grillas no se superpongan y se puedan ver los detalles de cada una.
Graficar las superficies tiene la ventaja (con respecto a mapas 2-D) que se aprecian cada grilla de forma mas intuitiva. Por ejemplo, es más fácil correlacionar máximos entre superficies.
Script
0. Datos: Se utilizan las siguientes 3 grillas (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 medidas de km.
- Relieve: Grilla GMRTv3_5.grd (ver Anexo 4. Descargar grilla GMRT).
- Anomalías Gravimétricas de Aire Libre: Versión 24.1 de la Grilla de Sandwell y Smith (link ftp grav.img.24.1).
1. Variables principales: Es útil definir las siguientes variables que se utilizarán para dibujar los bloques 3D.
- REGION: área geográfica (longitud y latitud) del bloque 3D.
- PROJ: Proyección geográfica del mapa (Mercator) y ancho de la figura (14 cm).
- PROZ: Escala vertical del bloque (3 cm).
- persp: acimut y elevación de la vista 3D.
- A: dirección de acimut de la fuente de luz para el efecto de sombreado.
SET REGION=-88/-18/-60/20
SET PROJ=M14c
SET PROZ=3c
SET persp=210/25
SET A=160
2. Bloque 3-D Moho: El primer es la profundidad a la discontinuidad de Mohorovičić (moho). El bloque 3D incluye los ejes de longitud y latitud.
A. Grilla: surface permite crear una grilla a partir de la tabla de datos depthtomoho.xyz usando una superficie continua.
B. Información Grilla: grdinfo permite saber el rango de valores de la grilla que luego se agrega a la variable REGION3D. El valor mínimo corresponde a la mínima profundidad del moho (-70 km en este caso) y el máximo a 10 km (para que incluya al relieve).
grdinfo %CUT% -T1
SET REGION3D=%REGION%/-70/10
C. Paleta de Color: grd2cpt crea la paleta de color (CPT) %color% que se ajusta a los valores de la grilla (%CUT%).
gmt grd2cpt -Crainbow -Z %CUT% >%color%
D. Gradiente: grdgradient crea una grilla (%SHADOW%) que se utilizará para el efecto sombreado.
gmt grdgradient %CUT% -A160 -G%SHADOW% -Ne0.5
E. Bloque 3D: grdview crea el bloque 3D (ver 13.1 para más detalles).
gmt grdview %CUT% -R -J -JZ -O -K -p -Qi300 -C%color% -I%SHADOW% >> %OUT% -Wf0.5 -BnSwEZ+b -Baf -Bzaf+l"Altura (Km)"
SET PROJ=M14c
SET PROZ=3c
SET persp=210/25
SET A=160
2. Bloque 3-D Moho: El primer es la profundidad a la discontinuidad de Mohorovičić (moho). El bloque 3D incluye los ejes de longitud y latitud.
A. Grilla: surface permite crear una grilla a partir de la tabla de datos depthtomoho.xyz usando una superficie continua.
- -I6m: resolución de 6 minutos de arco. Se usa una resolución mayor a la de los datos (1 grado) para que la imagen no se vea pixelada.
- -R%REGION%: Región geográfica de la grilla creada.
- -G%CUT%: grilla de salida guardada con en la variable %CUT%.
- -fg: establece que los datos son coordenadas geográficas.
- -r: estable pixel node registration (en lugar de gridline, más info).
B. Información Grilla: grdinfo permite saber el rango de valores de la grilla que luego se agrega a la variable REGION3D. El valor mínimo corresponde a la mínima profundidad del moho (-70 km en este caso) y el máximo a 10 km (para que incluya al relieve).
grdinfo %CUT% -T1
SET REGION3D=%REGION%/-70/10
C. Paleta de Color: grd2cpt crea la paleta de color (CPT) %color% que se ajusta a los valores de la grilla (%CUT%).
gmt grd2cpt -Crainbow -Z %CUT% >%color%
D. Gradiente: grdgradient crea una grilla (%SHADOW%) que se utilizará para el efecto sombreado.
gmt grdgradient %CUT% -A160 -G%SHADOW% -Ne0.5
E. Bloque 3D: grdview crea el bloque 3D (ver 13.1 para más detalles).
gmt grdview %CUT% -R -J -JZ -O -K -p -Qi300 -C%color% -I%SHADOW% >> %OUT% -Wf0.5 -BnSwEZ+b -Baf -Bzaf+l"Altura (Km)"
2. Bloque 3-D Relieve: Se utilizan los mismos comandos con las siguientes modificaciones:
- Con grdmath se convierten los valores de altura en m a km (divido 1000).
- En grdview sólo se dibuja la escala vertical (se omiten los ejes de latitud y longitud) y se agrega un desplazamiento vertical de 7 cm (-Y7c).
- Con pscoast se dibuja la línea de costa (ver ).
gmt grdcut "GMRTv3_5.grd" -G%CUT% -R%REGION% -fg
gmt grdmath %CUT% 1000 DIV = %CUT%
gmt grd2cpt -Crainbow -Z %CUT% >%color% gmt grdgradient %CUT% -A160 -G%SHADOW% -Ne0.5
gmt grdview %CUT% -R%REGION3D% -J -JZ -O -K -p -Qi300 -C%color% -I%SHADOW% >> %OUT% -BZ -Bzaf+l"Altura (km)" -Y7c
gmt pscoast -R -J -JZ -O -K -p%persp%/0 -Df >> %OUT% -W1/0.3,black
 |
| Bloques 3D SIN desplazamiento vertical de las superficies. |
 |
| Bloques 3D CON desplazamiento vertical de las superficies de 7 cm. |
3a. Bloque 3-D Anomalías de Aire Libre: Se utilizan los mismos comandos con las siguientes ligeras modificaciones:
- Se agrega grdclip para recortar los máximos y mínimos de la superficie. En este ejemplo se limitan los máximos y mínimos a 100 y -100 mGal respectivamente.
- Para más detalles de Anomalías de Aire Libre ver 7.1.
gmt img2grd "E:\Facultad\Datos_Geofisicos\Gravimetria\Sandwell-Smith\24.1\grav.img.24.1" -R%REGION% -G%CUT% -T1 -I1 -E -S0.1
gmt grd2cpt %CUT% -Z -L-100/100 -D > %color%
gmt grdclip %CUT% -G%CUT% -Sa100/100 -Sb-100/-100
gmt grdgradient %CUT% -A%A% -G%SHADOW% -Ne0.5
gmt grdview %CUT% -R%REGION%/-370/950 -J -JZ -O -K -p -Qi300 -C%color% -I%SHADOW% >> %OUT% -BZ -Bzaf+l"mGal"-Y7c
gmt grd2cpt %CUT% -Z -L-100/100 -D > %color%
gmt grdclip %CUT% -G%CUT% -Sa100/100 -Sb-100/-100
gmt grdgradient %CUT% -A%A% -G%SHADOW% -Ne0.5
gmt grdview %CUT% -R%REGION%/-370/950 -J -JZ -O -K -p -Qi300 -C%color% -I%SHADOW% >> %OUT% -BZ -Bzaf+l"mGal"-Y7c
 |
| Bloque 3D de anomalías aire libre SIN grdclip. |
 |
| Bloque 3D de anomalías aire libre CON grdclip. |
No hay comentarios.:
Publicar un comentario