8.1. Mapa Magnetométrico. EMAG2 v2.

Elaborar mapa magnetométrico a partir de grilla EMAG2 (Earth Magnetic Anomaly Grid, V2). Se explican opciones avanzadas para crear y graficar paletas de colores (CPT) con áreas sin datos (NaN; Not a Number).

Mapa de las Anomalías Magnetométricas Corticales (EMAG2 v2) con datos del IGN (provincias, administrativos), COPLA (límites marítimos), del GSHHG (países, línea de costa).

Datos: Se utilizan la grilla EMAG2_V2.grd en formato NETCDF (compatible con GMT). Ver Base de Datos.

Script: Link

Mapa Base: Usamos comandos previos para hacer un mapa que incluya datos del IGN (límite administrativos), límites marinos (COPLA), límites países, líneas de costas, escala gráfica, y marco al mapa explicado en entradas previas.

1. Preparación Datos: Primero utilizo grdcut para recortar la grilla según la variable %REGION% y lo guardo con el nombre %CUT%_2. 

REM    Recortar Grilla EMAG2     
gmt grdcut -R "EMAG2_V2.grd" -G%CUT%_2 -fg

Luego, con grdconvert asigno los valores -99999 como áreas Sin Datos (=nf/1/0/-99999) y los guardo con la variable %CUT%.

REM    Asignar Datos -99999 a NaN. nf=    GMT netCDF format (32-bit float, COARDS, CF-1.5. nd=
    gmt grdconvert %CUT%_2=nf %CUT%=nf/1/0/-99999

2. Paleta de Color: Se presentan 3 cpt elaborados con las herramientas makecpt y grd2cpt. Ver entradas anteriores.

REM    Crear Paleta de Colores. Paleta Maestra (-C), Definir rango (-Tmin/max/intervalo), CPT continuo (-Z)
rem    gmt grd2cpt -Crainbow -Z %CUT% >%color%
rem    gmt makecpt -Crainbow -Z -T-200/200/10 -D >%color%
    gmt grd2cpt -Crainbow -Z %CUT% -Di -L-200/200 >%color%

3. Escala de Color: Se presentan 4 escala de colores verticales con la leyenda (+l) "Magnetrometría (nT)" elaborados con psscale. La única particularidad que se agrega a las explicadas anteriormente es la opción -D+n que agrega un recuadro al inicio de la escala para los áreas sin datos.

REM    Agrega escala de colores. (-E triangles). Posición (-D) (horizontal = h). Posición (x,y) +wlargo/ancho. Anotaciones (-Ba). Leyenda (+l). 
rem    gmt psscale -O -K -C%color% -Baf+l"Magnetometr\355a (nT)" -I >> %OUT% -Dx15.7/0.0+w13.7/0.618c 
rem    gmt psscale -O -K -C%color% -Baf+l"Magnetometr\355a (nT)" -I >> %OUT% -Dx15.7/1.2+w12.5/0.618c+n"Sin datos"       
    gmt psscale -O -K -C%color% -Baf+l"Magnetometr\355a (nT)" -I >> %OUT% -Dx15.7/1.5+w11.9/0.618c+n"Sin datos"+e 

grd2cpt.

grd2cpt y psscale +n"Sin Datos".

makecpt -200 a 200, psscale +d +n"Sin Datos".

grd2cpt -200 a 200, psscale +d +n"Sin Datos".

4. Color Sin Datos: Por default, las áreas sin datos (NaN) son pintadas de gris (gray 127.5). Con gmtset se puede modificar los parámetros por default. En este ejemplo, para asignar el color blanco a NaN se debe anteponer a grd2cpt/makecpt el siguiente comando:

    gmtset COLOR_NAN white

P.D. 19-04-2017:
5. Sin Datos Transparente: El argumento -Q de grdimage hace que las áreas sin datos (NaN) sean transparentes.

gmt grdimage -R -J -O -K %CUT% -C%color% -I%SHADOW% >> %OUT% -Q