5.4. Mapa Topografico. Paleta de Color (CPT).

Para hacer un mapa topográfico se utiliza una paleta de colores (CPT) para pintar las distintas alturas. En general ocurre que el rango de valores de la CPT es mayor o menor que los de la grilla. A continuación se explica como ajustar la CPT en ambos casos y como utilizar una CPT del sitio cpt-city.

Primero se detalla como descargar una nueva CPT del sitio cpt-city (punto 0). En 1 se explica como saber el rango de valores de la grilla. Luego se explican opciones de makecpt para ajustar una CPT con un rango de valores MENORES que la grilla (pasos 2 a 4). Finalmente, se explica para cuando los rangos son MAYORES (5 y 6). El script utilizado es una adaptación del explicado en 5.1. Mapa Topográfico. Sólo se explican los comandos que fueron modificados. 

0. cpt-city: En el sitio cpt-city hay disponible CPT adicionales. Para descargarlas, una vez elegida, hay que hacer clic en la opción "cpt" (GMT colour palette table). El archivo descargado conviene guardarlo con el resto de las CPT maestras (gmt5/share/cpt/). Si se guarda en en otro carpeta, se debe indicar la ruta (relativa o absoluta; ver 2. Scripts) en el argumento -C. En este ejemplo vamos a utilizar la CPT "wiki-schwarzwald-cont" que esta definida para el rango 0-1500.

1. Rango de valores de la grilla: grdinfo -T permite saber cual es el rango de valores de la grilla (ver punto 6 de la entrada 6 para más detalles). Para este ejemplo, el rango de elevaciones es -3450/6900. Los valores aparecen en la consola al ejecutar el script.

grdinfo %CUT% -T50
pause

2. CPT 1 original: -C indica la CPT a utilizar. Si no se utilizan otros argumentos, makecpt abarca el rango original de valores de la CPT (0-1500). Para las áreas con valores mayores y menores a ese rango se utilizan respectivamente los colores COLOR_FOREGROUND (F) y COLOR_BACKGROUND (B) que aparecen en los triángulos de la escala de color.

gmt makecpt -Cwiki-schwarzwald-cont > %color% -Z

3. CPT 1 -D: La opción -D produce que COLOR_BACKGROUND y COLOR_FOREGROUND coincida con los máximos y mínimos de la CPT.

gmt makecpt -Cwiki-schwarzwald-cont > %color% -Z -D

4. CPT 1 -T: -T permite redefinir el rango de valores que abarca la CPT. Ver grdinfo para saber el rango de valores de la grilla. El siguiente comando permite estirar la CPT para valores entre 0 y 6900. Para las zonas bajo el mar se utiliza el BACKGROUND_COLOR.

gmt makecpt -Cwiki-schwarzwald-cont > %color% -Z -T0/6900

5. CPT 2 original: Para el segundo caso, usamos la CPT maestra geo, la cual esta definida para el rango -8000/8000 (abrir el archivo con gedit). El comando psscale dibuja la escala de colores al costado.

gmt makecpt -Cgeo > %color% -Z
gmt psscale -R -J -O -K -C%color% >> %OUT% -DJRM+o1.0c/0.8+w8.5/0.618c+e+n"Sin Datos" -Ba1000+l"Topograf\355a (m)" -I

6. CPT 2 -G: El argumento -G de psscale permite ajustar el rango de valores de la escala de color a los representados en el mapa. En este ejemplo se utiliza -G-3450/6900 (ver punto 1).

gmt psscale -R -J -O -K -C%color% >> %OUT% -DJRM+o1.0c/0.8+w8.5/0.618c+e+n"Sin Datos" -Ba1000+l"Topograf\355a (m)" -I -G-3450/6900

PD: 17-Abril-2017
7. Anotaciones Escala de Color: El argumento -W en psscale (agregado en versiones más recientes) permite multiplicar todos los valores de Z de la CPT por un factor. Por lo tanto, -W0.001 cambia las anotaciones de la escala de m a km.

    gmt psscale -R -J -O -K -C%color% >> %OUT% -DJRM+o1.0c/0.8+w8.5/0.618c+e+n"Sin Datos"  -I -Ba1+l"Topograf\355a (km)" -G-3450/6900 -W0.001