SAC 中将位移记录转换成速度记录

将位移记录转换成速度记录,本质上就是微分,在 SAC 中有三种方法可以实现。

先生成测试用的数据:

SAC> fg seis
SAC> rmean; rtr; taper
SAC> w seis.sac
  1. 直接对 diff 命令对波形做时间域微分:

    SAC> r ./seis.sac
    SAC> dif
    SAC> w seis.sac.diff
  2. transfer 命令去仪器响应的功能实现微分,该操作本质上在频率域完成:

    SAC> r ./seis.sac
    SAC> transfer from none to vel
       Station (CDV), Channel (-12345)
    SAC> w seis.sac.transfer
  3. 时间域的微分在频率域等效于乘以叫频率:

    SAC> r ./seis.sac
    SAC> fft
     DC level after DFT is -0.0051373
    SAC> mulomega
    SAC> ifft
    SAC> w seis.sac.mulomega

测试结果表明,三种算法的结果差别很小,可以忽略。不过需要注意的是,使用 diff 做时间域微分时数据点数会少一个。