2.3. Mapa Satelital. Blue Marble

Mapas continentales a partir de imagen satelital (Blue Marble).

Se explica como hacer mapamundis o hemisféricos a partir de la imagen satelital Blue Marble: Next Generation de la NASA.

Script

Datos:
Imagen Satelital: Utilizo la siguiente imagen de la NASA (link) en formato geotiff de 0,1º de resolución. Esta imagen contiene datos en los polos y asigna colores azules a los océanos según la profundidad.

• BlueMarbleNG-TB_2004-06-01_rgb_3600x1800.tiff

1. Imagen a RGB: Con mirrone se descompone la imagen en grillas RGB (ver 2.2 para más detalles).

• BlueMarble_r.grd
• BlueMarble_g.grd
• BlueMarble_b.grd

2. Mapas: grdimage permite recombinar la grillas RGB para distintas proyecciones (-J) y regiones geográficas (-R).

SET    PROJ=W-65/15c
SET    REGION=-110/30/-90/20

gmt grdimage -R -J -O -K >> %OUT% "BlueMarble_r.grd" "BlueMarble_g.grd" "BlueMarble_b.grd"

3. Límites y marco: se agregan los límites de los países, línea de costa y marco del mapa

gmt pscoast -R -J -O -K -N1/0.2,- >> %OUT%
gmt pscoast -R -J -O -K -Df -W1/ >> %OUT%
gmt psbasemap -R -J -O -K -B0 >> %OUT%

Mollweide mostrando Sudamérica, Atlántico y parte de África).

4A. Mapamundis: Otros ejemplos de mapas con distintas proyecciones (-J) y regiones geográficas (-R).

Mollweide centrado en Greenwich (W).

Mollweide centrado en Sudamérica (W-65).

Robinson centrado en Sudamérica (N-65).

Eckert VI centrado en Sudamérica (Ks-65).

4B. Mapas hemisféricos:

Proyección estereográfica (-S-65/0/90)

Proyección ortográfica (-G-65/0/90)

PD 08-08-2019: Con la versión de GMT 6.0 es posible hacer mapas con proyección hemisférica NO centrada en el ecuador.

Proyección ortográfica (-G-65/-30/90)