Mapa del Gradiente Vertical de la Gravedad (VGG) para el sur de Sudamérica y el Atlántico Sur. Además se explica como optimizar una paleta de colores a partir de los datos de una grilla.
Introducción: Este mapa es similar al anterior (7.1), con la única diferencia
que se grafican los datos del Gradiente Vertical de la Gravedad (VGG) en
lugar de las Anomalías Gravimétricas de Aire Libre. Es por eso que se
utiliza el mismo script, y solamente se adaptaron los comandos de convertir datos, Paleta de Color y Escala de Color. Se presentan 4 paletas de colores (cpt), combinados con 2 Escalas de Color.
Datos:
Script: Link
4. Convertir datos: img2grd nos permite extraer un subconjunto definidos por %REGION% (punto 1) del archivo curv.img.23.1. Los datos se guardan en la grilla %CUT% con una resolución de 1 minuto (-I1). Los valores se multiplican por 0.02 (-S0.02) para obtener Eötvös. -T1 permite obtener datos en todos los puntos. -E fuerza a que la grilla tenga exactamente la región solicitada (%REGION%).
gmt img2grd "E:\Facultad\Datos_Geofisicos\Gravimetria\Sandwell-Smith\23.1\curv.img.23.1" -R%REGION% -G%CUT% -D -T1 -I1 -E -S0.02
5. Paleta de Color: Se presentan 4 cpt elaborados con las herramientas makecpt y grd2cpt que permitan mostrar mejor los datos (ver entrada anterior para más detalles). Todas las cpt tienen valores continuos (-Z) y se guardan con el nombre %color%. Los 2 primeros abarcan todo el rango de valores (-140 a 140 Eötvös). Para los otros 2 cpt, los valores se restringen a -20 a 20 (Eötvös) e incluyen el parámetro -D que produce que los valores que queden fuera del rango sean pintados con los colores de los extremos del cpt.
REM Crear Paleta de Colores. Paleta Maestra (-C), Definir rango (-Tmin/max/intervalo), CPT continuo (-Z)
gmt makecpt -T-140/140/10 -Z >%color%
gmt grd2cpt %CUT% -Z >%color%
gmt makecpt -T-20/20/1 -Z -D >%color%
gmt grd2cpt %CUT% -Z -L-20/20 -D >%color%
 |
| Mapa del Gradiente Vertical de la Gravedad (VGG, Sandwell et al. v 23.1) con datos del IGN (provincias, administrativos), COPLA (límites marítimos), del GSHHG (países, línea de costa). |
Datos:
- gradiente Vertical de la Gravedad (VGG, Sandwell et al., 2014) en formato img (curv.img.23.1). Ver Base de Datos.
Script: Link
4. Convertir datos: img2grd nos permite extraer un subconjunto definidos por %REGION% (punto 1) del archivo curv.img.23.1. Los datos se guardan en la grilla %CUT% con una resolución de 1 minuto (-I1). Los valores se multiplican por 0.02 (-S0.02) para obtener Eötvös. -T1 permite obtener datos en todos los puntos. -E fuerza a que la grilla tenga exactamente la región solicitada (%REGION%).
gmt img2grd "E:\Facultad\Datos_Geofisicos\Gravimetria\Sandwell-Smith\23.1\curv.img.23.1" -R%REGION% -G%CUT% -D -T1 -I1 -E -S0.02
5. Paleta de Color: Se presentan 4 cpt elaborados con las herramientas makecpt y grd2cpt que permitan mostrar mejor los datos (ver entrada anterior para más detalles). Todas las cpt tienen valores continuos (-Z) y se guardan con el nombre %color%. Los 2 primeros abarcan todo el rango de valores (-140 a 140 Eötvös). Para los otros 2 cpt, los valores se restringen a -20 a 20 (Eötvös) e incluyen el parámetro -D que produce que los valores que queden fuera del rango sean pintados con los colores de los extremos del cpt.
REM Crear Paleta de Colores. Paleta Maestra (-C), Definir rango (-Tmin/max/intervalo), CPT continuo (-Z)
gmt makecpt -T-140/140/10 -Z >%color%
gmt grd2cpt %CUT% -Z >%color%
gmt makecpt -T-20/20/1 -Z -D >%color%
gmt grd2cpt %CUT% -Z -L-20/20 -D >%color%
7. Escala de Color: Se presentan 2 escala de colores verticales con la leyenda (+l) "Gradiente Vertical de la Gravedad (Eötvös)" elaborados con psscale a partir del cpt %color%. La primera, que se utilizó para las 2 primeras cpt, tiene una longitud de 13,7 cm y anotaciones cada 10 unidades (-Ba10). La segunda escala, utilizada para las otras 2, incluyen triángulos en los extremos (+e), una longitud de 13,1 cm y anotaciones cada 5 unidades (-Ba5). El parámetro -Dx indica la posición de la escala en coordenadas del gráfico (-15,5 en eje X, 0,0 y 0,3 en eje Y).
REM Agrega escala de colores. (-E triangles). Posición (-D) (horizontal = h)
rem gmt psscale -O -K -Dx15.5/0.0+w13.7/0.618c -C%color% -I -Ba10+l"Gradiente Vertical de la Gravedad (E\366tv\366s)" >> %OUT%
gmt psscale -O -K -Dx15.5/0.3+w13.1/0.618c+e -C%color% -I -Ba5+l"Gradiente Vertical de la Gravedad (E\366tv\366s)" >> %OUT%
 |
| makecpt -T-140/140/10 -Z |
 |
| grd2cpt %CUT% -Z |
 |
makecpt -T-20/20/1 -Z -D
|
|
| grd2cpt %CUT% -Z -L-20/20 -D |
No hay comentarios.:
Publicar un comentario