5.2. Mapa Topográfico II. Sombreado.

Esta entrada explica como adaptar el mapa topográfico anterior para agregar un sombreado (hill-shaded).

Mapa topográfico de la región de cuyo con iluminación (hill-shaded) desde el este.

Mapa (Link):

1. Gradiente: Para crear un mapa topográfico con sombras (shaded relief o hill-shaded) es necesario primero calcular un grilla de gradientes. Para ello se utiliza grdgradient que se guarda con el nombre de la variable "shadow".
 La opción -A especifica un azimut para realizar una derivada direccional de la grilla. En la práctica indica la dirección desde donde se ilumina la grilla.
La opción -N permite normalizar los valores. Por default elijo -Ne0.5 pero uno debería probar distintos valores (de -N y -A) para obtener un mejor resultado. En este ejemplo la iluminación principal es desde el E (90).

SET    SHADOW=temp_%title%_shadow.grd
gmt grdgradient %CUT% -G%SHADOW% -A90 -Ne0.5

2. Imagen: Para que la imagen incluya al archivo de gradiente con los valores de intensidad se indica con la opción -I y el nombre de la grilla previamente creada.

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

3. Escala de Color: Agregando la opción -I a psscale, se agrega el efecto de sombra a la escala vertical.

gmt psscale -O -K -Dx15.5/0+w11/0.618c -C%color% -Ba500+l"Topograf\355a (m)"  -I >> %OUT% 

Ejemplos con otras iluminaciones:

Iluminación desde el Norte (-A0).

Iluminación desde el Sur (-A180).

Iluminación desde el Oeste (-A270).

PD 24-06-2017: A partir de la versión 5.4.1 de GMT (creo) es posible definir el sombreado directamente desde grdimage con -I, evitando el paso 1. Con +a se define el azimut y con +n el método de normalización y el parámetro. El argumento -I+a90+ne0.5 permite obtener el mismo efecto de sombreado del paso 1.


gmt grdimage -R -J -O -K %CUT% -C%color% -I+a90+ne0.5 >> %OUT%