gmt begin site_map png,eps

  gmt set FORMAT_GEO_MAP=ddd:mm:ssF

  gmt set MAP_FRAME_TYPE=fancy MAP_FRAME_WIDTH=0.5p

  gmt set MAP_TICK_LENGTH=0.0c

  gmt set FONT_ANNOT_PRIMARY=12p

  set range=-180/180/-90/90

  set map=N6i

  set gridX=a60f30

  set gridY=a30f30

  gmt coast -R%range% -J%map% -Bx%girdX% -By%gridY% -Bwsne -Gwhite -S135/206/235 -X2.5i -Y3.6i

  awk "{print $3, $2}" site_RT1.txt | gmt plot -Sc0.2 -G0/0/255

  awk "{print $3, $2, $5, $1}" site_RT1.txt | gmt text -F+f7p,0+j -Dj0.15c/0.15c

gmt end show

附件：

site_map.txt

site_RT1.txt

注意：GMT 6脚本site_map.txt为windows下bat文件，运行前需将后缀名由txt改为bat。

效果图：

转载本文请联系原作者获取授权，同时请注明本文来自周锋科学网博客。

链接地址：https://wap.sciencenet.cn/blog-350278-1289830.html?mobile=1

收藏

分享