Graficar múltiples perfiles en un mapa.
Introducción:En el ejemplo anterior (11.4) se muestra como graficar un perfil de anomalías sobre un mapa a partir de una linea. En este, se muestra como graficar múltiples perfiles perpendiculares a un perfil guía.
Script
gmt grdcut "GMRTv3_5.grd" -G%CUT% -R
gmt makecpt -Cibcso -T-8400/0/10 -Z -N > %color%
gmt makecpt -Celevation -T0/6050/50 -Z >> %color%
gmt grdgradient %CUT% -A270 -G%SHADOW% -Ne0.5
gmt grdimage -R -J -O -K %CUT% -C%color% -I%SHADOW% >> %OUT%
gmt psscale -R -J -O -K -DJRM+o0.3c/0+w15/0.618c -C%color% -Ba1+l"Elevaciones (km)" -I >> %OUT% -W0.001
gmt psxy -R -J -O -K >> %OUT% "Trench_PlateProject.txt" -W1,red
2. Datos Gravimétricos: Con img2grd se convierten los datos gravimétricos (anomalías de Aire Libre) para su posterior extracción. Ver el punto 4 de 7.1 para más detalles.
gmt img2grd %GRA% -R -G%TGRA% -T1 -I1 -E -S0.1
3. Definir Perfil Guía: Definir una línea N-S según el meridiano 55 ºW entre 10 ºS y 55ºS. Los datos se guardan en el archivo "temp_line" con los siguientes comandos. La unión de los puntos se realiza siguiendo un círculo máximo.
REM Perfil: Crear archivo para dibujar perfil (Long Lat)
echo -55 -10 > "temp_line"
echo -55 -55 >> "temp_line"
4. Extraer datos de los perfiles: con grdtrack se crean los perfiles perpendiculares al perfil guía ("temp_line") con las siguientes opciones. Los datos de los perfiles se guardan en el archivo "temp_data"
Introducción:En el ejemplo anterior (11.4) se muestra como graficar un perfil de anomalías sobre un mapa a partir de una linea. En este, se muestra como graficar múltiples perfiles perpendiculares a un perfil guía.
Script
0. Datos:
1. Mapa Base: Se crea con los siguientes comandos (ver 5.1 y 5.2 para más detalles). - DEM: GMRTv3_5.grd
- Anomalías de Aire Libre: grav.img.27.1 (en formato img). Ver Base de Datos.
- Trinchera: Datos de Plates Project (Trench_PlateProject.txt). Ver Base de Datos.
gmt grdcut "GMRTv3_5.grd" -G%CUT% -R
gmt makecpt -Cibcso -T-8400/0/10 -Z -N > %color%
gmt makecpt -Celevation -T0/6050/50 -Z >> %color%
gmt grdgradient %CUT% -A270 -G%SHADOW% -Ne0.5
gmt grdimage -R -J -O -K %CUT% -C%color% -I%SHADOW% >> %OUT%
gmt psscale -R -J -O -K -DJRM+o0.3c/0+w15/0.618c -C%color% -Ba1+l"Elevaciones (km)" -I >> %OUT% -W0.001
gmt psxy -R -J -O -K >> %OUT% "Trench_PlateProject.txt" -W1,red
2. Datos Gravimétricos: Con img2grd se convierten los datos gravimétricos (anomalías de Aire Libre) para su posterior extracción. Ver el punto 4 de 7.1 para más detalles.
gmt img2grd %GRA% -R -G%TGRA% -T1 -I1 -E -S0.1
3. Definir Perfil Guía: Definir una línea N-S según el meridiano 55 ºW entre 10 ºS y 55ºS. Los datos se guardan en el archivo "temp_line" con los siguientes comandos. La unión de los puntos se realiza siguiendo un círculo máximo.
REM Perfil: Crear archivo para dibujar perfil (Long Lat)
echo -55 -10 > "temp_line"
echo -55 -55 >> "temp_line"
4. Extraer datos de los perfiles: con grdtrack se crean los perfiles perpendiculares al perfil guía ("temp_line") con las siguientes opciones. Los datos de los perfiles se guardan en el archivo "temp_data"
- -G%TGRA%: Grilla con datos de anomalías de Aire Libre.
- -C5000k/50k/200k: Perfiles perpendiculares equiespaciados cada 200 km según el perfil guía, cada uno de 5000 km de longitud y con datos cada 50 km.
5. Graficar perfiles en el mapa: pswiggle permite graficar los perfiles con las siguientes opciones:
6. Argentina-Chile: Ejemplos acotando la región a Argentina-Chile.
SET REGION=-88/-42/-56/-21
- -i0,1,4: Selecciona las columnas con los datos del archivo "temp_data"
- -F: Agrega recuadro blanco debajo de la escala (+gwhite), con un borde (+p) y efecto sombreado (+s).
- -Gred@50+p: Pinta de rojo las anomalías positivas con 50% de transparencia (@).
- -Gblue@50+n: Pinta de azul las anomalías negativas con 50% de transparencia (@).
- -Z500: Escala de las anomalías.
- -D: Dibuja una escala de 100 mGal abajo a la derecha (en este ejemplo).
- -T: Dibuja el perfil.
- -W: Dibuja el borde de las anomalías.
 |
| -C2000/100/100 |
 |
| -C2000k/10k/500k |
 |
| -C2000/10/200 |
 |
| -C5000k/100k/1000k |
 |
| -C5000k/50k/200k |
SET REGION=-88/-42/-56/-21
 |
| Sólo se pinta la anomalía (sin -T ni -W). |
No hay comentarios.:
Publicar un comentario