【绘图】GMT绘制世界地图的底图（不同投影方式）+ GMT站点分布图

版权声明：本文为博主原创文章，未经博主允许不得转载。保留追究法律责任，转载联系博主。交流欢迎加QQ群463339221。 https://blog.csdn.net/m0_37362454/article/details/82111368

主要是直观的记录不同的投影方式其投影效果，便于以后再进行投影时可以快速的选择投影类型

pscoast -Rg -JN0/5i -A10000 -Gwhite -Slightgray -Wthinnest -P -Dc -Bg30/g15:.Robinson: >robinson.ps

pscoast -Rg -JH180/4.5i -A10000 -Gwhite -Slightgray -Wthinnest -P -Dc -Bg30/g15:.Hammer: >Hammer.ps

pscoast -R-180/180/-85/85 -JM4.5i -A10000 -Gwhite -Slightgray -Wthinnest -P -Dc -Bg30/g15:.Mercator: >Mercator.ps

pscoast -R-180/180/-90/90 -JQ80/4.5i -A10000 -Gwhite -Slightgray -Wthinnest -P -Dc -B60/30:.Cylin-Equidistant: >cylin_equidistant.ps

pscoast -R-180/180/-90/90 -JG85/40/4.5i -A10000 -Gwhite -Slightgray -Wthinnest -P -Dc -Bg30/g15:.Orthographic: >Orthographic.ps

 

代码

准备文件：

点文件：example_point.txt

-53.854 -1.999  0.2c
116.594 40.065  0.2c
-96.111 39.672  0.2c

文字文件：example_text.txt，文字的坐标由点坐标的纬度加2度得到。

-53.854 0.001   10p 0   BC  Amazon
116.594 42.065  10p 0   BC  BeiJing
-96.111 41.672  10p 0   BC  America

•  

# 在全球地图上标一些点
ps=example.ps
R=-180/180/-80/80
J=X16d/9d #d : degree of arc
point=example_point.txt
text=example_text.txt

# ps文件头
gmt psxy -J$J -R$R -T -K > $ps

# 海岸线
gmt pscoast -R$R -J$J -Ba60  -K -Ggrey -Swhite -O >> $ps

# 文字
gmt pstext $text -R$R -J$J -F+f+a+j -K -O >> $ps
# 画点
gmt psxy $point -R$R -J$J -Sa -Gblack -K -O >> $ps

# ps文件尾
gmt psxy -R$R -J$J -T -O >> $ps

# 保存
gmt psconvert $ps -A -Tg -P

结果