- 博客(5)
- 资源 (2)
- 收藏
- 关注
RSS订阅原创 IRIS数据的初步处理
Python：IRIS数据的初步处理此前一文分享了如何向IRIS申请数据，详细可参考(Python：向IRIS发送邮件申请事件或者连续数据)一文。但下载好数据后是一个个以.tar.mseed为结尾的数据包，该怎么对这个数据处理，得到sac文件？本文将简单介绍一下数据的转换过程，并提供一个简单的例子。一. 将所有下载好的.tar.mseed为结尾的数据包放入mseed_data文件夹，执行01batch_decompress_data.py脚本，针对每一个.tar.mseed的数据包会生成.mseed，.
2022-03-06 11:13:08
1144
原创 HYPODD输入文件.pha读取
Python：HYPODD输入文件.pha读取此脚本意在读取HypoDD的输入.pha文件将结果保存为字典。因为可能会碰到对.pha文件进行修改，比如修改拾取数据的权重，可能需要修改台站名或台网名，所以此脚本可能可以作为辅助帮助修改.pha文件。.pha的文件格式如下：# 2016 10 14 00 00 09.270 42.8085 13.2121 6.63 0 0 0 0 1T1214 1.3196 1 PED03 7.7896 1 SED24 3.0596 1 PED01 5.6196
2022-02-24 17:45:37
988
1
原创 Python：向IRIS发送邮件申请事件或者连续数据
Python：向IRIS发送邮件申请事件或者连续数据此脚本可以自动向IRIS发送邮件申请事件或者连续数据（本人暂时只是测试了126邮箱，且脚本只在ubuntu环境下测试过）。运行脚本前需要开起邮箱的SMTP服务，详细可参考csdn中的名为：利用Python代码自动发邮件/出现550，535错误？？一文。脚本包的下载链接为：xxxxxx一. 申请地震事件数据准备台站信息目录，如station_info_public所示，第一列是台网名，第二列是台站名，第三、四列是台站的纬度和经度，第四列是要下载的分量
2021-05-08 14:27:28
717
原创 在obspy中获得Wood-Anderson仪器振幅
在obspy中获得Wood-Anderson仪器振幅在地震学学习中，有时候需要对地震的震级大小进行确定，这个时候可能需要将原始波形进行转换，得到Wood-Anderson仪器振幅。这里简单举个例子介绍一下如何通过obspy获得Wood-Anderson仪器振幅。导入需要的包：from obspy import readfrom obspy.io.sac import sacpz.attach_paz定义变量paz_wa保存Wood-Anderson仪器响应文件：paz_wa = {'sensi
2021-04-09 16:43:05
714
原创 Python 儒略日和公历互转代码
Python 儒略日和公历互转代码# python3#-------------------------------------------------------------------------------# Purpose: conversion between Julian date and the Gregorian calendar.# Author: plovess# Created: 2020-03-01# update: #------
2021-04-01 22:10:58
677
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人