✍️ 陈箫翰 • 📅 2025-01-25

# 根据发震时刻绘制不同颜色的圆点

本例根据地震发震时刻绘制不同颜色的圆，这个例子展示了GMT拥有强大的时间格式处理能力。

示例数据：`data.txt`

数据格式如下:

    2022-04-03T01:25:11.530 99.459 37.082  9.0 1.4 0.0 0.0 0.0 1
    2022-04-02T18:09:7.260 103.417 36.829  10.0 1.6 0.0 0.0 0.0 2
    2022-03-29T15:37:0.570 103.373 36.882  10.0 1.7 0.0 0.0 0.0 3
    ...

前三列分别为发震年月日时分秒、经纬度

    #!/usr/bin/env bash
    gmt begin ex010
        # 获取发震时刻最大最小值
        # 以 -Tzmin/zmax/dz 的形式报告第一列，每四周一个间隔dz
        T_cpt=$(gmt info data.txt -T4w+c0)

        # 根据发震时刻绘制不同颜色的圆点
        gmt makecpt -Chot ${T_cpt} -Z -Di
        gmt plot data.txt -Sc0.2c -W -C -JM15c -Baf -i1,2,0
        gmt colorbar -Bxaf -C
    gmt end show