征集令

维护者征集 下面列出了我发起的几个开源项目。挖了坑但是没有时间去填,所以需要征集志愿者参与维护。志愿者越多越好,项目开源的目的就是让所有人都可

判断 SAC 数据是否已去除仪器响应

有些时候,波形数据拿到手了,但是却不知道波形数据是否有被处理过,尤其是有没有去除仪器响应。

先说结论,判断一个 SAC 数据是否已经去除仪器响应的准则是:

若一个 SAC 数据中每个数据点的值都是整数,则这个数据的仪器响应没有被去除的概率为 99.9%

在极坐标下绘制深度剖面

在做大区域层析成像时,通常需要绘制深度剖面。一个简单的示例如下:

这样的深度剖面图实际上是在极坐标下绘制的。先直接给出绘制上图所使用的命令:

gmt psbasemap -R60/80/5571/6371 -JPa15c/70z -Bx5 -By200 -BWNse > test.ps

用 SAC 快速拾取震相

震相拾取,是地震学的一个基本工作。SAC 提供了 plotpk 命令,可以用于手动拾取 震相到时,用起来还算是比较简单方便。

绘制颜色渐变的线条

想要绘制一条颜色变化的线段,下面是演示代码:

#!/bin/bash
# GMT v5.2.1

gmt makecpt -Crainbow -T-2/2/1 > lines.cpt
gmt psxy -JX15c/4c -R0/6/0/4 -B1 -Clines.cpt -W2p > test.ps << EOF
> -Z-1.5
1 2
2 2
> -Z-0.5
2 2
3 2
> -Z0.5
3 2
4 2
> -Z1.5
4 2
5 2
EOF

绘图效果如下图: