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

示例贡献者:

陈箫翰

本例根据地震发震时刻绘制不同颜色的圆，这个例子展示了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 -Ra -JM15c -Baf -i1,2,0
    gmt colorbar -Bxaf -C
gmt end show
../../_images/8c5f348b0529ec48c9309de8fe53a1f6.png