Question

官方文档

sac

简介

绘制 SAC 格式的地震波形数据

sac 模块可以读取SAC文件并绘制波形数据。

备注

sac 模块修改自原 pssac 与 pssac2，其功能类似，但语法不同。

sac 模块实现波形绘制的步骤是：

1. 读入 SAC 文件列表

2. 根据

   gmt sac [ saclist|SACfiles ]

   SACfiles

   要绘制在图上的一系列SAC文件名，目前仅支持等间隔SAC数据

   saclist

   SAC 文件列表，每行包含一个SAC文件名及其对应的控制参数。

   文件格式为:

   filename [X  Y [pen]]
   

   • filename 是要绘制的SAC文件名

   • X 和 Y 控制SAC波形的第一个数据点在地图上的位置。若省略 X 和 Y ，则使用其默认值，否则此处指定的 X 和 Y 将具有最高优先级

     • 对于线性投影而言，X 默认是SAC文件的开始时间。使用

       -Bparameters (more …)

       设置底图边框和轴属性

       -C[t0/t1]

       只读取并绘制时间窗 t0 到 t1 范围内的波形。

       t0 和 t1 均是相对于参考时间的秒数，参考时间由

       利用 SAC 的命令 funcgen seismogram 生成了波形，想要绘制单个波形，并分别为正负部分涂色:

       gmt sac seis.SAC -JX10c/5c -R9/20/-2/2 -Baf -Fr -Gp+gblack -Gn+gred -png single
       

       利用 SAC 命令 datagen sub tel *.z 生成多个波形，将其绘制在距离剖面上:

       gmt sac *.z -R200/1600/12/45 -JX15c/5c -Bx200+l"T(s)" -By5+lDegree -BWSen \
            -Ed -M1.5c -W0.5p,red -png distance_profile
       

       利用 SAC 命令 datagen sub tel *.z 生成多个波形，将其绘制在地图上:

       gmt begin map pdf
       gmt sac *.z -JM15c -R-120/-40/35/65 -Baf -M1i -S300c
       saclst stlo stla f *.z | gmt plot -St0.4c -Gblack -i1,2
       gmt end show
       

       相关模块

       meca, polar, coupe, basemap, plot