seisamuse
·
PYTHON脚本练习（一）
💡
原文中文，约5100字，阅读约需12分钟。
📝
内容提要
该Python脚本遍历目录中的SAC文件，进行窄带滤波并计算信噪比(SNR)。信号窗口为2.5-5km/s，噪声窗口为信号后1000秒。仅保存SNR大于3的数据，并将结果写入CSV文件和生成可视化图表，使用8个CPU进行并行处理。
🎯
关键要点
- 该Python脚本遍历目录events_20250619下所有子目录中的bhz.SAC_rm结尾的SAC文件。
- 对数据进行窄带滤波，滤波器为4个极点0相位的Butterworth，滤波周期为arange(25,145,10)。
- 计算信噪比(SNR)，信号窗口为2.5-5km/s，噪声窗口为信号后1000秒。
- 仅保存信噪比大于3的数据，并将结果写入新的文件夹bp_sac。
- 采用8个CPU进行并行处理以提高效率。
- 统计每个周期信噪比大于3的波形数据及其平均信噪比。
- 将统计结果写入CSV文件，并生成可视化图表。
❓
延伸问答
这个Python脚本的主要功能是什么？
该脚本遍历指定目录中的SAC文件，进行窄带滤波并计算信噪比(SNR)，仅保存SNR大于3的数据。
如何计算信噪比(SNR)？
信噪比通过信号窗口内波形包络的最大值与噪声窗口的均方根值的比值计算得出。
脚本是如何处理数据的？
脚本使用4个极点0相位的Butterworth滤波器对数据进行窄带滤波，并采用8个CPU进行并行处理以提高效率。
处理后的结果如何保存？
处理后的结果被写入新的文件夹bp_sac中，并生成CSV文件和可视化图表。
信号窗口和噪声窗口的定义是什么？
信号窗口为2.5-5km/s，噪声窗口为信号后1000秒的长度窗口。
脚本生成了哪些可视化图表？
脚本生成了信噪比大于3的波形数量和平均信噪比的可视化图表。
➡️
