En esta entrada voy a explicar como modificar el mapa anterior (2. Hemisférico) para crear una mapa polar de la Antártida Argentina con proyección acimutal estereográfica. También se explica como agregar un borde al mapa y anotaciones.
Datos: GSHHG, DCW.
Versión GMT: 5.2.1 (r15220)
Descripción del Script (Link):
Para crear este mapa solamente hay que modificar las variables del sistema de proyección y la región geográfica del mapa:
Variables del Mapas:
1. Proyección: Primero hay que modificar la posición donde se centra el mapa, en este caso el polo sur (-90). También se elige que meridiano va a dibujarse hacia arriba (en este caso -65 para que el sector Antártico Argentino quede hacia arriba).
rem Proyecciones acimutales requieren 3 parametros + 1 opcional (lon0/lat0[/horizon]/width
rem L(a)mbert Equal Area; (S)tereographic; Orto(g)rafica; (E)quidistante
rem SET PROJ=G-65/-30/90/15c
SET PROJ=S-65/-90/90/15c
2. Región: Para los mapas con proyecciones polares los atajos d y g no sirven. Por lo tanto hay que definir la región que queremos que abarque el mapa. Para este ejemplo utilizamos la región comprendida al sur del paralelo 50° Sur (-180/180/-90/-50).
REM Region geografica del mapa (W/E/S/N) d=-180/180/-90/90 g=0/360/-90/90
rem SET REGION=d
SET REGION=-180/180/-90/-50
3. Marco Mapa: Finalmente, modificamos el argumento -B de psbasemap para agregar un borde (f) y anotaciones (a) de las coordenadas al mapa. Con el argumento -Bafg GMT calcula valores por default en base a las características del mapa.
REM Dibujar marco del mapa
rem gmt psbasemap -R -J -O -K -Bg >> %OUT%
gmt psbasemap -R -J -O -K -Bafg >> %OUT%
 |
| Mapa polar con proyección estereográfica centrado en el polo sur. |
Datos: GSHHG, DCW.
Versión GMT: 5.2.1 (r15220)
Descripción del Script (Link):
Para crear este mapa solamente hay que modificar las variables del sistema de proyección y la región geográfica del mapa:
Variables del Mapas:
1. Proyección: Primero hay que modificar la posición donde se centra el mapa, en este caso el polo sur (-90). También se elige que meridiano va a dibujarse hacia arriba (en este caso -65 para que el sector Antártico Argentino quede hacia arriba).
rem Proyecciones acimutales requieren 3 parametros + 1 opcional (lon0/lat0[/horizon]/width
rem L(a)mbert Equal Area; (S)tereographic; Orto(g)rafica; (E)quidistante
rem SET PROJ=G-65/-30/90/15c
SET PROJ=S-65/-90/90/15c
2. Región: Para los mapas con proyecciones polares los atajos d y g no sirven. Por lo tanto hay que definir la región que queremos que abarque el mapa. Para este ejemplo utilizamos la región comprendida al sur del paralelo 50° Sur (-180/180/-90/-50).
REM Region geografica del mapa (W/E/S/N) d=-180/180/-90/90 g=0/360/-90/90
rem SET REGION=d
SET REGION=-180/180/-90/-50
3. Marco Mapa: Finalmente, modificamos el argumento -B de psbasemap para agregar un borde (f) y anotaciones (a) de las coordenadas al mapa. Con el argumento -Bafg GMT calcula valores por default en base a las características del mapa.
REM Dibujar marco del mapa
rem gmt psbasemap -R -J -O -K -Bg >> %OUT%
gmt psbasemap -R -J -O -K -Bafg >> %OUT%
No hay comentarios.:
Publicar un comentario