||
SAC文件是常见的地震波形数据文件,存储的是单道的地震波形数据,亦可在文件里存储该波形相关的信息,如台站名,台网名,台站位置,地震事件位置等。
Python可使用obspy读取SAC文件。以下是相关代码示例。
import obspy
import matplotlib.pyplot as plt
# 读取SAC文件
# sac_file = 'xxxx.sac' # 替换'xxxx.sac'为你想绘制的SAC文件
st = obspy.read(sac_file)
tr = st[0]
t = tr.times() + tr.stats.sac.b
# t是时间, tr.data是波形数据
# 绘制波形数据
plt.plot(t, tr.data)
plt.savefig('seis.png', dpi=300, bbox_inches='tight')
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2025-12-6 04:09
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社