SEISAMUSE

  python脚本练习-fk。

  你是不是也想让自己的程序像obspy一样，叫全世界的人pip install就可以用？

  以下是python脚本练习2，功能为读取选定时间段内的hdf5文件，文件名形如20130512_cross_spec.hdf5。该文件总共有24段，每一段是一个小时某台阵平均互相关谱。然后画出这个互相关谱，保存为png格式图片。

  今天学习匹配场处理（Matched Field Processing, MFP），见这里。

  收集了一些有趣的jupyter notebook插件。

  哪里去了解有趣全球新闻？

  $\LaTeX$和markdown的数学算符几乎是一样的，那怎么编辑呢？

  以下是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，并画出统计结果。
  阅读全文 »

  那么Hexo下如何画流程图呢？

  今天学习Sequencer。