domingo, 15 de mayo de 2016

10. Leyenda, Mapa de ubicación y valores por defecto.

Se explica como agregar una leyenda y un mapa de ubicación, y cómo personalizar los mapas.

En este entrada, explico como agregar, en un mapa previo (en este caso el mapa tectónico de la entrada 4), una leyenda con los símbolos dibujados en el mapa, como agregar una mapa de ubicación (o inset) para localizar y como editar los valores por defecto que utiliza GMT para graficar.


Script: (link). Incluye los comandos previos utilizados para crear el mapa tectónico.

1. Valores por defecto: GMT tiene cerca de 100 parámetros con valores por defecto que definen la apariencia de los mapas (ver manual para una explicación general y figuras que explican los parámetros). Los parámetros se agrupan en 10 sub-secciones temáticas (ver gmt.conf), entre las cuales se encuentran Fuentes, Formato y Marco del Mapa. Los parámetros por defecto pueden modificarse con gmtset. A continuación se muestran los comandos para modificar algunos de los parámetros. 

I: Fuentes: Incluye parámetros que controlan los títulos (FONT_TITLE), etiquetas (FONT_LABEL) y anotaciones primarias (FONT_ANNOT_PRIMARY). Esttán definidos por 3 valores: tamaño, tipografía (en GMT hay 35 fuentes predefinidas que pueden indicarse con el número correspondiente) y color.Los siguientes comandos permiten modificarlos:

REM    Sub-seccion FUENTE
    gmtset FONT_ANNOT_PRIMARY 8,Helvetica,black
    gmtset FONT_LABEL 10,Helvetica,black
    gmtset FONT_TITLE 12,4,Black



II: Formato El siguiente comando se especifica que en las coordenadas se utilicen los símbolos WESN en lugar de valores negativos ("-").


REM    Sub-seccion FORMATO
    gmtset FORMAT_GEO_MAP ddd:mm:ssF



III: Marco del Mapa: Los siguientes comandos modifican los parámetros del marco del mapa. El 1° (AXES) indica que en los ejes se agreguen las anotaciones/coordenadas (en mayúscula). El 2° (TYPE) define un borde elegante (fancy) y con esquinas redondeadas (+). El último redefine el ancho del marco (0.1 cm).

REM    Sub-seccion MAPA
    gmtset MAP_FRAME_AXES WesN
    gmtset MAP_FRAME_TYPE fancy+
    gmtset MAP_FRAME_WIDTH 0.1




2. Leyenda: El proceso de agregar una leyenda a una figura se puede dividir en 2 partes. En la primera se crea un archivo específico (specfile) que contiene la información de la leyenda. En la segunda parte el specfile es graficada con pslegend.

2.A. Crear specfile: Con ECHO (de MS-DOS) se crea el archivo specfile (en este ejemplo llamado "temp_legend"), el cual está conformado por distintos registros. Existen 14 tipo de registros (indicados por la letra inicial y cada uno con su propia sintaxis) que permiten armar la leyenda.

En este ejemplo, el primer registro crea un encabezado (H) con el texto "Leyenda del Mapa".

El segundo registro indica la cantidad de columnas (N 3) que se utilizarán para los próximos 6 registros. En la primera línea se describen los 3 tipos de límite de placas con distintos símbolos (S): línea (-), falla inversa (f+l+t), falla de rumbo (f+l+s) con las mismas características usadas en el mapa y el texto. En la segunda línea se dibujan los círculos (c) de los sismos para eventos someros, intermedios y profundos.

Con N 4 se establece que para la siguientes línea se dibujen 4 registros. Los 2 primeros son rectángulos (r) correspondiente a LIPS y a Ofiolitas, y 2 líneas (zonas de fracturas y dorsales extintas).

El registro G definen el espacio vertical entre líneas (0.075 cm).

Finalmente, con el último registro (M) se agrega una escala medida en la longitud -70, latitud -57, de 500 km con las unidades (+u) y elegante (f).

    echo H 10 Times-Roman Leyenda del Mapa >> temp_legend
    echo N 3 >> temp_legend
    echo S 0.25c - 0.5c - 3.0p,red           0.75c Dorsal                 >> temp_legend
    echo S 0.25c f+l+t 0.5c t 1.0p,green     0.75c Subucci\363n            >> temp_legend
    echo S 0.25c f+l+s 0.5c s 1p,black       0.75c L\355mite Transforme        >> temp_legend
    echo G 0.075c >> temp_legend

    echo S 0.25c c 0.25c red   0.40p     0.5c Sismos someros (0-100 km)         >> temp_legend
    echo S 0.25c c 0.25c green 0.40p     0.5c Sismos intermedios (100-300 km)     >> temp_legend
    echo S 0.25c c 0.25c blue  0.40p     0.5c Sismos profundos (300-700 km)     >> temp_legend
    echo G 0.075c >> temp_legend

    echo N 4 >> temp_legend
    echo S 0.25c r 0.5c p300/29  0.25p    0.75c LIPS          >> temp_legend
    echo S 0.25c r 0.5c purple4  0.25p    0.75c Ofiolitas      >> temp_legend
    echo S 0.25c - 0.5c - 1.0p,violet     0.75c Zonas de Fracturas >> temp_legend
    echo S 0.25c - 0.5c - 0.80p,orange    0.75c Dorsales Extintas  >> temp_legend


    echo G 0.075c >> temp_legend

    echo M -70 -57 500+u f                          >> temp_legend


2.B. Dibujar Leyenda: pslegend permite graficar el archivo temp_legend. Para centrar la leyenda de debajo del mapa (-Dx7.5/-0.2) con una ancho de 15 cm (+w15/0). Al establecer la altura en 0, entonces es calculada automáticamente por el programa. Con -F se dibuja un recuadro a la leyenda con un borde externo (+p), otro interno (+i) y esquinas redondeadas (+r). También sería posible incluir un relleno (+gcolor).


     gmt pslegend -R -J -O -K -Dx7.5/-0.2+w15/0+jTC -F+p+i+r temp_legend >> %OUT%

2.C. Leyenda Auxiliar: Dado que para dibujar las dorsales se utilizó dos símbolos superpuestos (línea roja gruesa y una línea blanca fina) es necesario hacer lo mismo en la leyenda. Para hacerlo, se repite la estructura de la leyenda previa y para el tercer registro se dibuja una línea blanca. Luego se vuelve a dibujar la leyenda. Conviene quitar los textos y bordes para que no se repitan.
Nota: En la primera línea se utiliza ">" para que sobreescriba el "temp_legend" anterior.

REM    Leyenda Auxiliar
    echo H 10 Times-Roman  > temp_legend   
    echo N 3 >> temp_legend
    echo S 0.25c - 0.5c - 1.0p,white >> temp_legend
     gmt pslegend -R -J -O -K -Dx7.5/-0.2+w15/0+jTC temp_legend >> %OUT%



3. Mapa de Ubicación: Para agregar un mapa secundario con la ubicación del mapa principal simplemente se superpone otro mapa sobre el regional. Para esto es importante redefinir las variables de -R (region) y -J (proyección). -X y -Y permiten ajustar la posición del mapa de ubicación. 
En este ejemplo, para dibujar un mapa hemisférico (como el de la entrada 1.2) con proyección azimutal ortográfica de 2.5 cm de (-JG) centrado en la región de Scotia con líneas negras de 0,3 puntos (-W), áreas secas en negro (-G) y húmedas en blanco (-S) y desplazado 5 cm hacia arriba (-Y5.0c) y 12,5 cm hacia la derecha (-X12.5c) se utiliza el siguiente comando:

REM    Dibujar Mapa INSET   
    gmt pscoast -Rd -JG-49.5/-57.5/2.5c -O -K -W0.3p,black -B90 -Dc -Gblack -Swhite -X12.5c -Y5.0c >> %OUT%



4. Área de Estudio: El siguiente comando sirve para agregar un recuadro con el área de estudio en el mapa de ubicación. En la primera parte se genera el archivo "temp_Recuadro" con las coordenadas geográficas de los extremos del mapa principal. En la segunda parte, se grafica dicho archivo con el comando psxy.

REM    Crear y Dibujar Recuadro de zona de estudio
rem    SET    REGION=-79/-20/-63/-50
    echo -79 -63 >> temp_Recuadro
    echo -79 -50 >> temp_Recuadro
    echo -20 -50 >> temp_Recuadro
    echo -20 -63 >> temp_Recuadro

    gmt psxy temp_Recuadro -R -J -Ap -W0.45p,red -O -L -K >> %OUT%



No hay comentarios.:

Publicar un comentario