GMT5的-B选项提供了很多种标注日期的方式。但要实现我们平时常用的标注需要通过自定义轴实现。
参考了seisman的例子
http://seisman.info/gmt5-custom-axes.html
http://docs.gmt-china.org/option/B.html
(一)关键
(1)-Bpxc和-Bsxc
设置x轴的primary和secondary属性,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:windowsFonts del gmt.* |
(二)出现的错误
- gawk不起作用,是因为没将其放在gmt的bin目录下,或者放到系统环境变量中。
- 没办法画出x轴,是因为文件名不对,多加了txt后缀,或者是编码格式问题,应将其utf8格式转为ansi格式。