blazing216的个人博客分享 http://blog.sciencenet.cn/u/blazing216

博文

Python绘制SAC文件

已有 773 次阅读 2025-10-27 15:00 |个人分类:地球物理相关基础|系统分类:科研笔记

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://blog.sciencenet.cn/blog-725052-1507660.html


下一篇:Python绘制以距离排序的地震图
收藏 IP: 113.55.110.*| 热度|

1 王涛

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...
扫一扫,分享此博文

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2025-12-6 04:09

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部