Hacer un gráfico con múltiples perfiles superpuestos con ejes personalizados.
El script permite crear múltiples perfiles superpuestos y ajustar los ejes de datos.
0. Pasos Previos: En los entradas previas (11.1 y 11.2) se explicó como extraer datos de grillas, graficarlos en perfiles e incluir las coordenadas de orientación (O y E). Para graficar múltiples perfiles en el mismo gráfico se realizan con los siguientes comandos (la diferencia con el ejemplo de 11.2 es que no incluyen el argumento -Y%H%). Tampoco incluyen los argumentos -B (que escribir los ejes) que se agregarán a continuación (1 y 2).
REM Dibujar Datos
gmt psxy -R0/%KM%/%Topo% -J -O -K "temp_data" -i2,3 -W1,blue >> %OUT%
gmt psxy -R0/%KM%/%Magn% -J -O -K "temp_data" -i2,5 -W1,red >> %OUT%
gmt psxy -R0/%KM%/%Grav% -J -O -K "temp_data" -i2,6 -W1,green >> %OUT%
1. Agregar 2 Ejes: Los siguientes comandos basados en psbasemap permiten agregar el eje de los datos topográficos a la derecha (--MAP_FRAME_AXES=E) y magnetometrícos a la izquierda (--MAP_FRAME_AXES=W). Con el argumento -By se definen intervalos de las anotaciones (a), marco (f) y grillado (g). Con --FONT_ANNOT_PRIMARY se puede ajustar que los valores de los ejes se dibujen de los mismos colores que los datos correspondientes. Los valores del eje topográfico fueron divididos por 1000 para escribirlos como km en lugar de m (-R0/%KM%/-6.2/5.3 y -Bya2.5f0.5g2.5). El título de los eje esta definido por FONT_LABEL.
Nota: --PARAMETRO=Valor permite realizar cambios temporarios de los valores por defecto y sólo afectan el comando en el que se escribe (ver gmtset).
REM Dibujar Escalas
rem gmt psbasemap -R0/%KM%/%Topo% -J -O -K -Bya2500f500g2500+l"Topograf\355a (m)" >> %OUT% --FONT_ANNOT_PRIMARY=8,Helvetica,blue --MAP_FRAME_AXES=E
gmt psbasemap -R0/%KM%/-6.2/5.3 -J -O -K -Bya2.5f0.5g2.5+l"Topograf\355a (km)" >> %OUT% --FONT_ANNOT_PRIMARY=8,Helvetica,blue --MAP_FRAME_AXES=E
gmt psbasemap -R0/%KM%/%Magn% -J -O -K -Bya+l"Magnetometr\355a (nT)" >> %OUT% --FONT_ANNOT_PRIMARY=8,Helvetica,red --MAP_FRAME_AXES=W
El script permite crear múltiples perfiles superpuestos y ajustar los ejes de datos.
0. Pasos Previos: En los entradas previas (11.1 y 11.2) se explicó como extraer datos de grillas, graficarlos en perfiles e incluir las coordenadas de orientación (O y E). Para graficar múltiples perfiles en el mismo gráfico se realizan con los siguientes comandos (la diferencia con el ejemplo de 11.2 es que no incluyen el argumento -Y%H%). Tampoco incluyen los argumentos -B (que escribir los ejes) que se agregarán a continuación (1 y 2).
REM Dibujar Datos
gmt psxy -R0/%KM%/%Topo% -J -O -K "temp_data" -i2,3 -W1,blue >> %OUT%
gmt psxy -R0/%KM%/%Magn% -J -O -K "temp_data" -i2,5 -W1,red >> %OUT%
gmt psxy -R0/%KM%/%Grav% -J -O -K "temp_data" -i2,6 -W1,green >> %OUT%
1. Agregar 2 Ejes: Los siguientes comandos basados en psbasemap permiten agregar el eje de los datos topográficos a la derecha (--MAP_FRAME_AXES=E) y magnetometrícos a la izquierda (--MAP_FRAME_AXES=W). Con el argumento -By se definen intervalos de las anotaciones (a), marco (f) y grillado (g). Con --FONT_ANNOT_PRIMARY se puede ajustar que los valores de los ejes se dibujen de los mismos colores que los datos correspondientes. Los valores del eje topográfico fueron divididos por 1000 para escribirlos como km en lugar de m (-R0/%KM%/-6.2/5.3 y -Bya2.5f0.5g2.5). El título de los eje esta definido por FONT_LABEL.
Nota: --PARAMETRO=Valor permite realizar cambios temporarios de los valores por defecto y sólo afectan el comando en el que se escribe (ver gmtset).
REM Dibujar Escalas
rem gmt psbasemap -R0/%KM%/%Topo% -J -O -K -Bya2500f500g2500+l"Topograf\355a (m)" >> %OUT% --FONT_ANNOT_PRIMARY=8,Helvetica,blue --MAP_FRAME_AXES=E
gmt psbasemap -R0/%KM%/-6.2/5.3 -J -O -K -Bya2.5f0.5g2.5+l"Topograf\355a (km)" >> %OUT% --FONT_ANNOT_PRIMARY=8,Helvetica,blue --MAP_FRAME_AXES=E
gmt psbasemap -R0/%KM%/%Magn% -J -O -K -Bya+l"Magnetometr\355a (nT)" >> %OUT% --FONT_ANNOT_PRIMARY=8,Helvetica,red --MAP_FRAME_AXES=W
2. Agregar 3er Eje: -Xa1.6c permite desplazar el eje gravimétrico 1,6 cm a la derecha del eje topográfico. La "a" indica que el origen del gráfico se desplaza a su posición original después de graficar el eje.
gmt psbasemap -R0/%KM%/%Grav% -J -O -K -Bya+l"Anomal\355a Aire Libre (mGal)" >> %OUT% --FONT_ANNOT_PRIMARY=8,Helvetica,green --MAP_FRAME_AXES=E -Xa1.6c
No hay comentarios.:
Publicar un comentario