地震数据进行预处理：滤波和添加地震事件信息

在处理地震数据时，对地震数据进行预处理：滤波以及添加地震事件信息。

    for i in range(520):
        st = obspy.read(filepath + str(i) + '.sac')
        tr = st[0]
        tr.filter('bandpass', freqmin=8, freqmax=15, corners=4, zerophase=True)
        # tr.filter('highpass', freq=8, corners=4, zerophase=True)
        # print(st[0].data)
        tr = (tr.data) / np.max(abs(tr.data))
        sacfile = Trace()
        sacfile.data = tr[:]
        
        sac = SACTrace.from_obspy_trace(sacfile)
        sac_data = sac.data
        sac.stla=35
        sac.stlo=110+(80+j*500+i*25)/111000
       
        sac.delta=0.0006
        sac.evla=35
        sac.evlo=110+(6568+500*j)/111000
        sac.evdp=0.05
        sac.write(overpath + str(j*520+i) + ".sac")

来源：CSDN

作者：朋哥的求职之路

链接：https://blog.csdn.net/weixin_46518295/article/details/104885568