✍️ 戴宗辉 , 王亮 , 陈箫翰 • 📅 2022-04-22

# 绘制地震 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