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.
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"
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
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
No hay comentarios.:
Publicar un comentario