「GMT」常用命令选项笔记

Posted on | In | 0 Comments | Visitors:
前一阵子学习GMT,在此对一些常用命令做下备忘

-B 指定绘图边界

1
-Bxinfo[/yinfo][:."title string":][W|w][E|e][S|s][N|n]
  • 其中 xinfo 格式为：
  • [a]tick[m|c][ftick[m|c]][gtick[m|c]][l|p][:"axis label":][:,"unit label":]
  • 其中a,f,g分别为刻度标记，刻度和网格间距

-C 指定绘图拷贝数

缺省为1

-H 告知头记录数

缺省为1，若数据文件头大于1，如4，必须写-H4

-J 选择地图投影

  • 线性投影 -JX width/height
    • -JXwidth[d] 当x,y用地理坐标的度表示时。
    • -JXwidthl 缩放前数值的log10
    • -JXwidthp power 缩放前数值的指数
  • 莫卡托投影 -JMwidth
  • 正形投影 -JGwidth 指定投影中心
  • 圆锥投影 -JBwidth 指定投影中心和两条标准平行线
  • Eckert投影 -JK[f|s]lon0/width 指定中央子午线

-R 指定绘图范围

-Rxmin/xmax/ymin/ymax

-Rxlleft/ylleft/xuright/yuright[r] 加r则是左下角和右上角坐标

-U 标注绘图日期

-U"optinal string"

-K 和 -O

规则

  • 如果一张图只要一个命令就能画出来，则不需要-K -O
  • 如果一张图要两个命令画出，则第一个-K 第二个 -O
  • 如果一张图要多个命令画出，则第一个用-K 中间用-K-O，最后一个用-O
  • 第一个命令使用重定向符号 >创建新文件并写入PS命令
  • 余下所有绘图命令使用重定向符号>>向PS文件中追加新的PS命令。

避免犯错技巧

例

1
2
3
4
5
6
7
8
9
10
11
12
13
14
#!/bin/bash
#
# GMT template
#
PS=example.ps
J=X6i #表明是6英寸
R=0/1/0/1
# Begin GMT Plot
psxy -J$J -R$R -T -K > $PS
# now begin your real plot commands
# every command should be with "-K" "-O" and ">>"
# End GMT Plot
psxy -J$J -R$R -T -O >> $PS
# Plot finished. Now you can convert PS to other format