martes, 17 de septiembre de 2013

Generic Mapping Tools (GMT): Proyecciones

Proyección- lineal

Vamos a comenzar haciendo un mapa básico para los ejes x-y. Queremos que este vaya de 10 a 70 en el eje x, anotando cada 10 espacios,y de -3 a 8 en y, anotando cada 1. El mapa final debería ser de 4 X 3 pies de tamaño. Una forma de hacerlo es:

psbasemap -R10/70/-3/8 -JX4i/3i -B10/1:."My first plot": -P > plot.ps


Proyección- logarítmica

Vamos a construir un basemap para un gráfico log-log. Vamos a asumir que los datos raw en "x" van desde 3 a 9613 y que "y" va desde 3.2 10^20  a 6.8 10^24. Una manera de hacerlo es:

psbasemap -R1/10000/1e20/1e25 -JX9il/6il \-B2:"Wavelength (m)":/a1pf3:"Power (W)":WS > plot.ps



Proyección- Mercator

A pesar de los problemas en las altas latitudes, la proyección Mercator (-JM) sigue siendo la más usada por los científicos.Esta es una de más múltiples proyecciones cilíndricas ofrecidas por GMT; Nos enfocaremos en esta proyección. La sintaxis es simplemente

-JM(ancho)

Para hacer mapas de costas usamos pscoast el cual automáticamente nos da acceso a las costas, ríos y base de datos de fronteras derivadas de GSHHG. Además de los cambios comunes, vamos a necesitar algunas especificaciones adicionales para pscoast 



Uno de -W,-G,-S debe ser elegido. Nuestro primer ejemplo de Coastline es de Centro américa

pscoast -R-90/-70/0/20 -JM6i -P -B5g5 -Gchocolate > map.ps


Proyección-Eckert VI y Proyección VI


Proyección-Ortográfica







No hay comentarios:

Publicar un comentario