1.5. Mapas Físico - Político. Argentina.

Mapa del territorio continental americano con datos físico-políticos-culturales del Instituto Geográfico Nacional (IGN). 

Mapa de Argentina con proyección Mercator Transversal con datos del IGN de cursos de Agua (Azul), cuerpos de Agua (azul claro), límites políticos (negro), red ferroviaria (rojo), ejidos urbanos (verde) y centros poblados (círculos negro).

Este ejemplo muestra como realizar un mapa rectangular del territorio continental americano de Argentina con una proyección no rectilínea (paralelos y meridianos), en este caso la proyección cilíndrica Mercator Transversal (o Cassini). Además se muestra como utilizar datos físico-políticos del Instituto Geográfico Nacional (IGN).

 

Datos:

• Global Self-consistent, Hierarchical, High-resolution Geography Database (GSHHG; incluido en GMT, Link): Conjunto de datos de líneas de costa de alta resolución.
• Digital Chart of the World (DCW): Base de datos vectorial de los países del mundo y de las provincias/estados de los 8 más grande (en superficie, incluye Argentina).
• Instituto Geográfico Nacional (IGN). Archivos SHP disponibles. Link para descargar los datos en formato gmt.

0. Preparación de los Datos:
GMT puede leer archivos vectoriales Shapefile (SHP) mediante GDAL. Sin embargo, a veces puede fallar. En el Anexo 1 se explica como adecuar los datos vectoriales (SHP) al formato compatible con GMT.

Script

Variables:
1. Región: En GMT hay 2 métodos para describir la región geográfica de un mapa (-R).

En los ejemplos anteriores, utilizamos la forma de la izquierda. En este ejemplo utilizamos la forma de la derecha, que nos permite obtener mapa rectangulares para proyección no rectilíneas. En este ejemplo utilizaremos las coordenadas 81° W y 55° S para la esquina inferior izquierda, y 53°W y 21° S para la esquina superior derecha. Se debe agregar una "r" al final para indicarle al GMT que se utiliza el segundo método.

REM    Titulo del mapa
    SET    title=EJ1.5_MapaPolitico_Argentina
    echo %title%

rem    Region: Argentina
    SET    REGION=-81/-55/-53/-21r

2. Proyección: Entre las proyecciones cilíndricas de GMT, hay 2 que no son rectilíneas: Mercator Transversal (T) y Cassini (C). Para este ejemplo utilizamos la primera centrada en el meridiano 60° W y en el paralelo 30°S (es opcional), y con un ancho de 13 cm.
 
rem    Proyeccion Cilindrica: (T)ransverse Mercator, (C)assini
    SET    PROJ=T-60/-30/13c

REM     Nombre archivo de salida
    SET    OUT=%title%.ps


Dibujar el Mapa:

3. Comandos previos: Agregamos las líneas de comandos de las entradas previas para pintar las áreas húmedas, resaltar a Argentina, dibujar el marco del mapa, los límites de países, línea de costa.

rem    Pintar areas húmedas: Oceanos (-S) y Lagos (-C+ll/)f
    Set color=dodgerblue2
    gmt pscoast -R -J -O -K -D%D% >> %OUT% -S%color% -Cwhite+l

REM    Resaltar paises DCW
    gmt pscoast -R -J -O -K >> %OUT% -EAR,FK,GS+grosybrown2+p

REM    Dibujar frame
    gmt psbasemap -R -J -O -K >> %OUT% -Baf

REM    Dibujar Bordes Administrativos. N1: paises.(Nborder[/pen])
    gmt pscoast -R -J -O -K -D%D% >> %OUT% -N1/0.75

REM    Dibujar Linea de Costa
    gmt pscoast -R -J -O -K -D%D% >> %OUT% -W1/faint

REM    Convertir ps en otros formatos: EPS (e), PDF (f), JPEG (j), PNG (g), TIFF (t)
    psconvert %OUT% -Tt -A

4. Datos IGN: con psxy se grafican los datos del IGN los cuales son vectoriales de tres tipos: puntos (pueblos), líneas (límites administrativos, departamentos, cursos de agua, red ferroviaria, límites marinos (COPLA)) y áreas (cuerpos de agua, ejidos urbanos). Para los puntos usamos el argumentos -S para dibujar círculos (c) de 0.04 cm (-Sc0.04). Para las líneas usamos -W y -G para definir el color de las áreas. Estas líneas de comando deben ubicarse antes de los datos de GSHHG.
 

Nota: El nombre de los archivos indican la ruta absoluta de la ubicación de los datos del IGN en mi computadora.


REM    Datos Instituto Geografico Nacional (IGN)
REM    -----------------------------------------------------------------------------------------------------------
rem    Cursos y Cuerpos de Agua
    gmt psxy -R -J -O -K "E:\Facultad\Datos_Geofisicos\IGN\1_GMT\004_Cuerpos_De_Agua.gmt" -G%color% >> %OUT%
    gmt psxy -R -J -O -K "E:\Facultad\Datos_Geofisicos\IGN\1_GMT\001_Cursos_De_Agua.gmt" -Wfaint,blue >> %OUT%

rem    Departamentos
    gmt psxy -R -J -O -K "E:\Facultad\Datos_Geofisicos\IGN\1_GMT\Departamentos.gmt" -Wthinnest,- >> %OUT% 

rem    limites administrativos
    gmt psxy -R -J -O -K "E:\Facultad\Datos_Geofisicos\IGN\1_GMT\limites_politico_administrativos.gmt" -Wthinner >> %OUT% 

rem    Limites Marinos (COPLA)
    gmt psxy -R -J -O -K "E:\Facultad\Datos_Geofisicos\IGN\1_GMT\LIMITE_COPLA.gmt" -Wthinner,.- -Am >> %OUT% 

rem    Dibujar red ferroviaria
    gmt psxy -R -J -O -K "E:\Facultad\Datos_Geofisicos\IGN\1_GMT\003_Red_Ferroviaria.gmt" -Wthin,red -Am >> %OUT%

rem    Pueblos y Ejidos Urbanos
    gmt psxy -R -J -O -K "E:\Facultad\Datos_Geofisicos\IGN\1_GMT\005_Centros_Poblados.gmt" -Sc0.04 -Gblack >> %OUT%
    gmt psxy -R -J -O -K "E:\Facultad\Datos_Geofisicos\IGN\1_GMT\010_Ejidos_Urbanos.gmt" -Wfaint -Ggreen >> %OUT%   

Edición de Script
5. Región: Se puede editar la región definida en el punto 1 para obtener un mapa de la provincia de Buenos Aires.

rem    Region: Argentina     
    SET    REGION=-81/-55/-53/-21r
rem    Region: Buenos Aires
    SET    REGION=-64/-42/-56/-33r