GMT5 自定义坐标轴

GMT5-B选项提供了很多种标注日期的方式。但要实现我们平时常用的标注需要通过自定义轴实现。

参考了seisman的例子

http://seisman.info/gmt5-custom-axes.html

http://docs.gmt-china.org/option/B.html

(一)关键

(1)-Bpxc-Bsxc

设置x轴的primarysecondary属性,c表示自定义轴

(2)跟在-Bpxc-Bsxc后面的文件格式

代码

set input=chuandianeq.txt

set title="Southwest China Seismicity"

set PS=chuandian_seismicity.ps

set J=X6i/2i

set R=1/2192/3/9

gmt psbasemap -J%J% -R%R% -BWSen -Bya1f0.5+lMagnitude -BpxcDateTick_px.txt -P -K > %PS%

gmt psbasemap -J%J% -R%R% -BsxcDateTick_sx.txt -K -O >> %PS%

gawk "{print $1,$3}" %input% | gmt psxy -J%J% -R%R% -Sc0.1c -W1p,black -O -K >>%PS%

gmt psxy -J -R -T -O >> %PS%

gmt psconvert %PS% -Tj -E400 -P -C-sFONTPATH=c:\windows\Fonts

del gmt.*

(二)出现的错误

  1. gawk不起作用,是因为没将其放在gmt的bin目录下,或者放到系统环境变量中。
  2. 没办法画出x轴,是因为文件名不对,多加了txt后缀,或者是编码格式问题,应将其utf8格式转为ansi格式。