绘制地震 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
../../_images/c2facdce29be7e155362daa6df78101e.png