15.2. Imagen Satelital + Topografia País.

Recortar los datos topográficos de un país y agregarlos sobre un mapa satelital.

En la entrada 15.1 el recorte de datos solo se podía realizar a través de la línea de costa. Este script, en cambio, permite utilizar cualquier polígono para  realizar el recorte. En este caso es particular, utilizo el límite de México (tomado de la base de datos Digital Chart of World; DCW), para recortar (clip) los datos de la grilla GMRT y sobreponerlos en una imagen satelital. Los pasos 0, 1, 2 y 5 son similares a los explicados en 15.1.

0. Datos: Se utiliza la imagen satelital (Mexico_2000m.png) y la grilla GMRT (GMRTv3_3_High.grd). Los anexos 2 y 4 explican como descargar la imagen satelital y la grilla GMRT, respectivamente.

Imagen satelital Mexico_2000m.png

Mapa batimétrico-topográfico creado con la grilla GMRT.

1. Variables y Mapa Satelital: Definimos la región geográfica y la  proyección (Mercator) que corresponden a la imagen satelital descargada. También se define el ancho de la figura (W=15c) y se define el marco del mapa (-Baf). En 2. Mapa Satelital se explican los comandos en detalle.

    SET    REGION=-118/-86/14/33
    SET    W=15c
    SET    PROJ=M%W%

2. Preparación datos GMRT: A continuación se preparan los datos de la grilla GMRT para sobreponer los datos topográficos. Se recorta la grilla GMRT y se crean una paleta de color (CPT) y una grilla para el sombreado. En 5.1 Mapa Topográfico se explican los comandos en detalle, incluyendo grdimage (ver punto 3).


   gmt grdcut "GMRTv3_3_High.grd" -G%CUT% -R%REGION% -fg 

   gmt makecpt -Cdem4 -T0/5400/100 -Z > %color%
   gmt grdgradient %CUT% -A90 -G%SHADOW% -Ne0.5


3. Crear trayectoria de recorte (clip path): pscoast permite extraer el límite de un país para que luego sean utilizados para recortar la grilla GMRT (punto 4). El argumento pscoast -E selecciona los datos de México (MX) y mediante la opción -M los guarda en el archivo "temp_clip". Este archivo define la forma en la que se hará el recorte.


REM    Extraer datos pais
    gmt pscoast -EMX -M > "temp_clip"

Datos parciales del archivo temp_clip.

4. Recortar/Agregar datos topográficos: psclip permite iniciar el recorte de los datos según el archivo "temp_clip". Luego, se dibuja el mapa topográfico (grdimage). Finalmente, el comando psclip -C cierra el recorte. Todos los comandos escritos entre estos dos comandos son recortados.

REM    Iniciar recorte siguiendo la linea de costa 
    gmt psclip -R -J -O -K "temp_clip" >> %OUT%
REM    ***********************************************************************

REM    Dibujar imagen a partir de la grilla
    gmt grdimage -R -J -O -K %CUT% -C%color% -I%SHADOW% >> %OUT%

REM    ***********************************************************************
REM    Finalizar recorte siguiendo la linea de costa 
    gmt psclip -R -J -O -K >> %OUT% -C

5. Adornar: Finalmente se puede agregar la línea de costa (pscoast -W), límite de los países (-N), límite de México (-EMX+pred) y la escala de colores (psscale).


    gmt psscale -R -J -O -K -C%color% >> %OUT% -Dx15.3/0+w9.7/0.618c -Ba+l"Alturas (m)"  -I

    gmt pscoast -R -J -O -K >> %OUT% -Df -W1/thinnest
    gmt pscoast -R -J -O -K >> %OUT% -Df -N1/thinnest
    gmt pscoast -R -J -O -K >> %OUT% -EMX+pred