如果想从地震数据中截取一个时间段，可使用如下脚本

当截取时间段是相对于这段波形数据本身的，比如波形数据的前10秒，或者第15-30秒，可用

t = tr.times()  # tr是存储波形数据的obspy.Trace对象

# tr = tr.times() + tr.stats.sac.b # 如果有stats.sac.b这个属性

sel = (t >= 15) & (t <= 30)

t_cut, data_cut = t[sel], tr.data[sel]

# t_cut, data_cut即为截取后的波形数据

当截取时间段是以UTC时间定义的，比如从2025年10月20日13:30:30 到之后的100秒，可用

tb = obspy.UTCDateTime('2025-10-20T13:30:30')

te = tb + 100

tr_cut = tr.slice(tb, te)

t_cut = tr_cut.times()

data_cut = tr_cut.data

# t_cut, data_cut即为截取后的波形数据

转载本文请联系原作者获取授权，同时请注明本文来自徐逸鹤科学网博客。

链接地址：https://wap.sciencenet.cn/blog-725052-1507671.html?mobile=1

收藏

分享