SEISAMUSE

Jun Xie的博客

发表于 更新于 分类于 本文字数： 4.7k 阅读时长 ≈ 4 分钟

  以下是python脚本练习1，功能包括：

  • 遍历目录events_20250619下所有子目录中以bhz.SAC_rm结尾的SAC文件；
  • 对这些数据进行窄带滤波，宽度为中心频率（周期分之一）的$\pm$5mHz，滤波器为4个极点0相位的Butterworth，滤波周期为arange(25,145,10)；
  • 计算窄带滤波后的每个周期的信噪比。信噪比定义为信号窗口内，波形包络的最大值比上噪声窗口的均方根。信号窗口定义为2.5-5km/s的到时。噪声窗定义为信号末端之后1000秒开始的1000秒长度的窗口。计算的信噪比写入到user0；
  • 将处理后的数据写到新的文件夹bp_sac中，文件名命名为z.year.jday.00.STA.bhz.period，仅保留信噪比大于3的数据。
  • 采用并行处理(8个cpu)。
  • 统计每个周期信噪比大于3的波形数据。
  • 统计每个周期信噪比大于3的波形的平均信噪比。
  • 将统计结构写入csv，并画出统计结果。
    阅读全文 »

发表于 更新于 分类于 本文字数： 2.7k 阅读时长 ≈ 2 分钟

  安装了seisman的HinetPy，运行脚本数据下载脚本的时候就遇到这个问题：

1
2
3
4
5
6
from HinetPy import Client, win32
  File "/home/junxie/work/Snet/HinetPy-main/HinetPy/__init__.py", line 23, in <module>
    from .client import Client
  File "/home/junxie/work/Snet/HinetPy-main/HinetPy/client.py", line 22, in <module>
    from urllib3.util import create_urllib3_context
ImportError: cannot import name 'create_urllib3_context' from 'urllib3.util' (/home/junxie/.conda/envs/py3/lib/python3.9/site-packages/urllib3/util/__init__.py)

怎么整？

阅读全文 »
0%