pssac2 使用教程

pssac2pssac 的升级版。相比于 pssac 而言增加了更多的功能,选项也更加符合 GMT 的风格。

pssac2pssac 的选项之间有很多相似之处,所以本文只介绍二者之间的不同之处,读者应首先阅读《pssac 使用教程》。

在地图上绘制波形

将 SAC 波形数据绘制在地图上也算是常见的需求之一。最常见的情况是,在地图上有一些台站,需要将波形画在台站的附近。

此问题的难点在于,在画地图时,-R 指定的是区域范围,即横轴是经度,纵轴是纬度,而 pssac 在绘制波形时,横轴是时间,纵轴是振幅。因而想要将波形放在特定的位置,不是一件容易的事情。

pssac 使用教程

注意:本文仅针对 Prof. Lupei Zhu 的 pssac 程序,其仅可在 GMT 4下使用。 对于GMT5用户,GMT 在 5.3.1 之后自带了 pssac,可以使用 gmt pssac 直接调用。

相关文档见:

pssac 可以读取 SAC 波形数据,并利用 GMT 的绘图功能将波形数据绘制在 PS 文件中。 pssac 绘制的波形图,大致可以分为三类:

  1. 波形振幅图:横轴是时间,纵轴是振幅;
  2. 剖面图:横轴是时间,纵轴是其他物理量(比如震中距、方位角、反方位角以及文件号);
  3. 任意位置图:将波形文件绘制在图上的任意位置;

本文首先简单介绍 pssac 的各个选项,然后依次介绍如何绘制波形振幅图、剖面图以及任意位置图,在这个过程中会详细介绍较复杂的选项。

GMT 绘制无刻度轴

提出问题

要绘制一张底图,要求

  1. 左边和下边有标注和刻度;
  2. 右边和上边无标注且无刻度;

最直接的想法是使用 -BWSen,绘图得到的效果图如左图所示,右边和上边虽然没有了标注,但是刻度线依然存在。实际想要实现的效果是右图。

SEED 格式转 SAC 格式

本文会介绍将 SEED、miniSEED 等格式转换为 SAC 格式所使用的工具和用法。需要注意, 本文并不会详细介绍工具的具体语法,只介绍常用的功能,其他功能及选项,读者应自行阅读相关文档。