绘制地震 M-T 图
- 示例贡献者:
戴宗辉(作者)、王亮(作者)、陈箫翰(修订)
绘制地震 M-T 图的前提条件是要学会利用 -B 选项 选项绘制时间轴底图。关于时间轴底图的绘制方法请参考 时间轴。此外,用户还必须按照如下格式准备好绘制地震 M-T 图所用的数据:
# 日期T时间 震级
2017-02-14T04:19:49 1.3
2017-02-14T06:22:05 1.9
2017-02-14T17:31:24 2.9
2017-02-14T17:42:41 1.8
2017-02-14T19:46:05 1.9
2017-02-15T01:38:02 1.8
2017-02-15T16:58:43 1.7
...
数据文件为两列,第一列是以 “YYYY-mm-ddThh:mm:ss” 格式给出的地震发震时刻,第二列是地震震级。使用 plot 模块的 -Sb 选项可以将以上格式的数据绘制为地震 M-T 图。
示例数据:mt.txt
#!/usr/bin/env bash
gmt begin mt
# 设置 SECONDARY TICK 长度为 0
gmt set MAP_TICK_LENGTH_SECONDARY 0p
# 设置地图上的时间只显示年份和月份
gmt set FORMAT_DATE_MAP yyyy-mm
# -BSW 只显示左边和下边坐标
# -Bsxa1O 设置X轴次级标注的间隔为1个月,注意最后一个字符是大写字母O不是数字0
# -Bpxa7Rf1d 设置X轴主级标注的间隔为7天,刻度间隔为1天
# -B+l"Date" 设置X轴标签
# -Bya1+l 设置Y轴标注间隔为1,并设置Y轴标签,标签使用斜体(2号字体),L 为下标
gmt basemap -R2017-02-10T00:00:00/2017-05-04T00:00:00/0/5 -JX15c/5c -BSW -Bsxa1O -Bpxa7Rf1d -Bx+l"Date" -Bya1+l@%2%M@-L@-@%%
# -Sb0.005c 中 0.005c 代表每条线的宽度
gmt plot mt.txt -Sb0.005c -Gblack
gmt end show