11.4. Perfil Geofísico en mapa

Dibujar perfil gravimétrico de anomalías de aire libre en un mapa.

Introducción: En ocasiones puede ser útil graficar datos de un perfil en un mapa. Se explica como extraer los datos de anomalías de Aire Libre a lo largo de un perfil y graficarlos sobre un mapa batimétrico-topográfico.

Nota: Es el primer ejemplo con la versión GMT6. Puede haber ligeros cambios en la sintaxis en algunas herramientas.

Versión GMT: 6.0.0rc4

Script

0. Datos:

• DEM: GEBCO_2014
• Anomalías de Aire Libre: grav.img.23.1 (en formato img). Ver Base de Datos.

1. Mapa Base: Se crea con los siguientes comandos (ver 5.1 y 5.2 para más detalles).

grdcut %DEM% -G%CUT% -R%REGION% 

makecpt -Crelief -G-5850/2500 -T-5850/2500/50 -Z -D >%color%
gmt grdgradient %CUT% -A90 -G%SHADOW% -Ne0.5h
gmt grdimage -R -J -O -K %CUT% -C%color% -I%SHADOW% >> %OUT%

 

2. Definir Perfil: Crear un archivo (temp_line) con las coordenadas del perfil. La unión de los puntos se realiza siguiendo un círculo máximo.

REM    Perfil: Crear archivo para dibujar perfil (Long Lat)
echo -76 -32 > "temp_line"
echo -46 -32 >> "temp_line"

3. Datos Gravimétricos: Se crea el archivo "temp_data" con los datos de las anomalias de aire libre. Ver 11.1 para la explicación de los comandos y 7.1 para el comando img2grd.

gmt img2grd %GRA% -R -G%TGRA% -T1 -I1 -E -S0.1
gmt sample1d "temp_line" -I0.2k > "temp_sample1d"
gmt mapproject "temp_sample1d" -G+k > "temp_track"
gmt grdtrack "temp_track" -G%TGRA% > "temp_data"
gmtinfo "temp_data"

3. Graficar perfil en el mapa: El comando pswiggle permite graficar el perfil con las siguientes opciones:

• -Gred+p: Pinta de rojo las anomalías positivas.
• -Gblue+n: Pinta de azul las anomalías negativas.
• -Z300: 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
• -t50: transparencia (50%).
• -i0,1,3: Selecciona las columnas donde están los datos en el archivo "temp_data" de longitud, latitud y anomalía de Aire Libre.

gmt pswiggle -R -J -O -K "temp_data" -i0,1,3 >> %OUT% -Gred+p -Gblue+n -Z300 -DjRB+o0.1/0.1+w100+lmGal -t50 -T -W

pswiggle sin -T ni -W.

pswiggle con -T y sin -W.

pswiggle con -T y -W.

Con proyección UTM 20 sur (U-20/15c).