gmt的投影方式

  其实参见GMT文档就可以了。不过,一般我们画图不需要每个参数都知道吧。下面就列出几个典型的投影方式和例子,免得以后到处找。

  • 蒙卡托投影
    首先是Mercator投影,最开始学的时候就是用这个。GMT主页说它可能是最著名的。它是等角圆柱形投影。

脚本如下:

1
2
3
4
#!/bin/bash
gmt begin mercater png
gmt coast -R90/140/10/50 -JM6i -Bxa10f5 -Bya10f5 -Dl -W0p -A5000 -Ggray
gmt end
  • 米勒投影
    米勒投影把整个地球搞的挺方,比其他的扁的好看一点儿。
1
2
3
4
#!/bin/bash
gmt begin miller png
gmt coast -R-180/180/-90/90 -JJ6i -Ggray -Swhite -Dl -Bxa90 -Bya45 -A30000 -W0p
gmt end
  • 圆柱等积投影
    可以画矩形扁扁的地球。
1
2
3
4
#!/bin/bash
gmt begin cyl1 png
gmt coast -R-180/180/-90/90 -JY0/0/6i -Ggray -Swhite -Dl -Bxa90 -Bya45 -A30000 -W0p
gmt end
  • 圆柱等距投影
    也是扁扁的地球。
1
2
3
4
#!/bin/bash
gmt begin cyl2 png
gmt coast -R-180/180/-90/90 -JQ6i -Ggray -Swhite -Dl -Bxa90 -Bya45 -A30000 -W0p
gmt end
  • 方位等距投影
    一般以某个点为中心的画法,例如接收函数台站和事件分布。
1
2
3
4
#!/bin/bash
gmt begin azi1 png
gmt coast -Rg -JE0/0/6i -Ggray -Swhite -Dl -Bxa90 -Bya45 -A30000 -W0p
gmt end
  • 正投影
    从无限远处看地球。可以简单理解为外星人视角。
1
2
3
4
#!/bin/bash
gmt begin orth png
gmt coast -R-180/180/-90/90 -JG0/0/6i -Ggray -Swhite -Dl -Bxa90 -Bya45 -A30000 -W0p
gmt end
  • 透视投影
    gmt文档说是飞行器视角,我愿意称之为上帝视角,你想怎么看怎么看。

    -JG0/0/6000/90/0/0/60/60/5i

    投影中心的经度和纬度(0ºE/0ºN)。
    观测者海拔高度(以公里为单位,6000 公里)。如果该值小于 10,则假定为观测者到地心的距离(以地球半径为单位)。如果附加“r”,则表示观测者到地心的距离(以公里为单位)。
    方位角以度为单位(90 度为正东)。相当于往东飞。
    倾斜角度(0度)。这是相对于天顶的视角,0度表示垂直向下看,60度倾斜角表示从地平线以上30度角看。
    旋转角度(0度)。这是图像的视轴旋转(顺时针方向)。180度旋转模拟倒飞的情况。
    视场的宽度和高度以度为单位(60)。这个数字取决于你是用肉眼观察(在这种情况下,你的视野大约为 60º 宽),还是用双筒望远镜等工具观察。

1
2
3
4
#!/bin/bash
gmt begin per png
gmt coast -R-180/180/-90/90 -JG0/0/6000/90/0/0/60/60/5i -Ggray -Swhite -Dl -Bxa90 -Bya45 -A30000 -W0p
gmt end
  • 罗宾逊投影
1
2
3
4
#!/bin/bash
gmt begin rob png
gmt coast -R-180/180/-90/90 -JN5i -Ggray -Swhite -Dl -Bxa90 -Bya45 -A30000 -W0p
gmt end