zoukankan      html  css  js  c++  java
  • 下载EU台网(欧洲台网)的地震波数据

    retrievedata.py

    ### here first to check the existence of the focal mechanism event file in the NDK directory,
    ###  if existence, mostly useful for the waveforms inversion,if not,download the event quakeml from
    ###  iris without focal mechanism. Then download the data from IRIS and GFZ, and before run, should 
    ###  update the username and password for the restricted data.
    ###  http://geofon.gfz-potsdam.de/waveform/archive/auth/auth-example.php
    ###  wget --post-file /data2/yjgao/token.asc https://geofon.gfz-potsdam.de/fdsnws/dataselect/1/auth -O cred.txt
    ###  Yajian Gao, 2017,12,26,GFZ
    
    from obspy import read_events
    from obspy.clients.fdsn import Client
    from obspy import UTCDateTime
    import os
    import sys
    import datetime
    import obspy
    from pathlib import Path
    import os.path
    from obspy.clients.fdsn.mass_downloader import CircularDomain, 
        Restrictions, MassDownloader
    from obspy.clients.seedlink.easyseedlink import create_client
    #def handle_data(trace):
        #print('Received the following trace:')
        #print(trace)
        #print()
    #client = create_client('geofon.gfz-potsdam.de', on_data=handle_data)
    client = Client("IRIS")
    if len(sys.argv) != 8:
        sys.exit("Usage: python %s year month day hour:minute:seconds eventid evla evlo"% sys.argv[0])
    year,month,day,origin,eventid,evla,evlo = sys.argv[1:]
    hour,minute,second=origin.split(':')
    #############################################download quakeml from IRIS without FOCAL MECHANISM
    
    EVENTFILENAME='C'+str(year)+str(month)+str(day)+str(hour)+str(minute)+'A.ndk'
    PATH="../NDK/"+str(EVENTFILENAME)
    print(PATH)
    ts=os.path.exists(PATH)
    if ts== True:
        print(ts)
        cat = obspy.read_events("../NDK/"+str(EVENTFILENAME))
        print(cat)

    cmd,并进入桌面(取决于下载数据的脚本在哪里) ,例:cd C:UsersDesktop

    输入下载命令:python retrievedata.py 2018 09 28 10:02:45 eventid  evla evlo

  • 相关阅读:
    Deltix Round, Summer 2021
    CF1558F Strange Sort 题解
    将文本复制到粘贴板上
    Mac安装2021ps
    腾讯云TDSQL-C云原生数据库技术
    腾讯云分布式数据库TDSQL再获认可
    T-TDSQL的核心技术
    TDSQL原创技术的出发点
    拓路前行-TDSQL追求极致体验的这一路
    TDSQL多集群下的自动化和交付
  • 原文地址:https://www.cnblogs.com/seisjun/p/10061344.html
Copyright © 2011-2022 走看看