15.1. Imagen Satelital + Batimetria.

Cómo recortar una grilla batimétrica para combinarla con una imagen satelital.

El script utiliza la líneas de costa (de la base de datos GSHHG) para recortar (clip) únicamente los datos de las áreas "húmedas" (mares, lagos) de la grilla GMRT y sobreponerlos en una imagen satelital.

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

Imagen satelital SAM_5km.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.

REM    Region Geografica
    SET    REGION=-78/-18/-56/-29

REM   Proyeccion Cilindrica: (M)ercator y Ancho de la figura (W) 
    SET    W=15c
    SET    PROJ=M%W%

REM    Ubicar imagen raster (o eps) en el mapa
    gmt psimage -R -J -O -K >> %OUT% -Dx0/0+w%W% "SAM_5km.png"  

REM    Dibujar frame (-B): Anotaciones (a), frame (f), grilla (g).

    gmt psbasemap -R -J -O -K >> %OUT% -Baf

2. Preparación datos GMRT: A continuación se preparan los datos de la grilla GMRT para sobreponer los datos batimétricos. 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_2.grd" -G%CUT% -R%REGION% -fg
   gmt makecpt -Cbathy -Z > %color%  
   gmt grdgradient %CUT% -A90 -G%SHADOW% -Ne0.5

3. Agregar datos batimétricos: El comando pscoast permite utilizar los datos de línea de costa para recortar la grilla GMRT. El comando pscoast -Sc inicia el recorte de los datos a través de la línea de costa (con resolución máxima, -Df). Luego, se dibuja el mapa batimétrico (grdimage). Finalmente, el comando pscoast -Q cierra el recorte. Todos los comandos escritos entre estos dos comandos son recortados.

REM    Iniciar recorte siguiendo la linea de costa 
    gmt pscoast -R -J -Sc -K -O >> %OUT% -Df
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 pscoast -R -J -Q -K -O >> %OUT%

4. Adornar: Finalmente se puede agregar la línea de costa (pscoast -W) y la escala de colores (psscale).

    gmt psscale -R -J -O -K -C%color% >> %OUT% -Dx15.3/0+w9.3/0.618c -Ba500+l"Profundidad (m)"  -I  -G-8200/0 
    gmt pscoast -R -J -O -K >> %OUT% -Df -W1/thinnest