Python:向IRIS发送邮件申请事件或者连续数据
此脚本可以自动向IRIS发送邮件申请事件或者连续数据(本人暂时只是测试了126邮箱,且脚本只在ubuntu环境下测试过)。运行脚本前需要开起邮箱的SMTP服务,详细可参考csdn中的名为:利用Python代码自动发邮件/出现550,535错误??一文。脚本包的下载链接为:https://download.csdn.net/download/u011563036/18469128
一. 申请地震事件数据
- 准备台站信息目录,如station_info_public所示,第一列是台网名,第二列是台站名,第三、四列是台站的纬度和经度,第四列是要下载的分量名。
BI DHAK 23.730 90.400 SHE,SHN,SHZ
RM SLV 21.325 103.907 HHE,HHN,HHZ
IN SHL 25.567 91.856 HHE,HHN,HHZ
IM CMAR 18.458 98.943 SHZ
- 准备事件目录,如event_catalog_public所示,第一列到第六列分别是年,儒略日,小时,分,秒,毫秒,第七、八列是事件的纬度和经度,第九列是震源深度(单位为m),第九列是震级。
2019 004 06 22 45 0630 -6.518 130.407 127900 5.6
2019 005 18 47 11 0740 51.335 -178.118 30000 5.9
2019 005 19 25 38 0520 -8.144 -71.587 570410 6.8
2019 005 22 54 14 0440 25.823 144.569 14280 5.7
2019 006 13 41 59 0600 34.125 45.679 14000 5.6
2019 006 17 27 18 0980 2.258 126.758 43210 6.6
2019 007 03 11 54 0280 7.866 121.950 19430 5.5
2019 008 12 39 30 0950 30.587 131.044 35000 6.3
2019 010 17 00 37 0630 -26.904 -175.545 10000 5.8
- 然后在main.py脚本中填写台站目录和事件目录的路径,以及一些相关的信息。随即运行main.py即可。
二. 申请连续数据
- 同上,准备台站信息目录。
- 运行create_fate_catalog_for_download_continue_data.py,生成虚拟的事件目录信息,如fate_catalog所示。
2019 001 00 00 00 0000
2019 002 00 00 00 0000
2019 003 00 00 00 0000
2019 004 00 00 00 0000
2019 005 00 00 00 0000
2019 006 00 00 00 0000
- 然后同上,运行main.py脚本,但这里注意time_length需要设置为一天的时间长度,脚本中给的是86399.99,即下载的是一天天的连续数据。
不管申请事件或者连续数据,如果申请成功最后IRIS会提供一个下载数据的链接。
注:脚本默认截取的事件数据的起始事件为事件目录中的发震时刻,可以通过设置main.py中的time_length变量的大小设定截取事件数据的长度。output文件夹中生成的每个文本就是发送给IRIS申请数据用的文本内容。如果不想重复下载数据,在每一次运行main.py之前先删除output文件夹。还有一点必须注意,main.py中填写的邮箱密码是开启smtp的授权码,不是邮箱的密码。切记!切记!切记!
参考:
1、https://blog.seisman.info/python-send-mail(其中auto_send_mail.py脚本修改自seisman的sendmail.py脚本)