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')

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

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

收藏

分享