Changeset 1182 for SHX


Ignore:
Timestamp:
16.02.2016 22:56:38 (4 years ago)
Author:
klaus
Message:

first working version of metadata copy from FDSN into local SHX database

Location:
SHX/trunk/SeismicHandler
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • SHX/trunk/SeismicHandler/commands/meta.py

    r1171 r1182  
    205205            except: 
    206206                print "illegal channel format i '%s'" % chan 
     207                continue 
    207208            if subcmd == 'complete_location': 
    208209                meta = ChannelMeta( 
  • SHX/trunk/SeismicHandler/tools/FdsnMeta2Db.py

    r1181 r1182  
    1515from obspy.fdsn.client import Client 
    1616from obspy.core import UTCDateTime 
     17from SeismicHandler.modules.stations import ChannelMeta, Stations 
    1718 
    1819 
    19 if __name__ == '__main__': 
     20def copyFdsnMeta( net, server ): 
    2021 
    2122    stime = UTCDateTime(1975, 1, 1 ) 
    2223    etime = UTCDateTime() 
    2324     
    24     #client = Client("BGR") 
    25     #inv = client.get_stations( network='GR', starttime=stime, endtime=etime, 
    26     #    level="response" ) 
    27     fp = open( 'metasave.pickle' ) 
    28     inv = pickle.load(fp) 
    29     fp.close() 
     25    stations = Stations() 
     26     
     27    client = Client(server) 
     28    inv = client.get_stations( network=net, starttime=stime, endtime=etime, 
     29        level="response" ) 
     30    #fp = open( 'metasave.pickle' ) 
     31    #inv = pickle.load(fp) 
     32    #fp.close() 
    3033     
    3134    for network in inv: 
    32         print network 
     35        #print dir(network), network.code 
    3336        for sta in network.stations: 
    3437            print sta.code, sta.latitude, sta.longitude 
    3538            print sta.start_date, sta.end_date 
    3639            for chan in sta.channels: 
    37                 print '    ', chan.code, chan.location_code, chan.start_date, chan.end_date 
    38                 #print chan.response 
     40                print '    ', sta.code, chan.location_code, chan.code, chan.start_date, chan.end_date 
     41                #print dir(chan.response) 
     42                print '    ', chan.response.instrument_sensitivity.value 
    3943                paz = chan.response.get_paz() 
    4044                #print dir(paz) 
    41                 print '      ', paz.normalization_factor 
     45                print '      ', paz.normalization_factor, 
    4246                print '      ', paz.normalization_frequency 
    4347                print '      ', paz.poles 
    4448                print '      ', paz.zeros 
     49                print '      ', paz.stage_gain, 
     50                print '      ', paz.stage_gain_frequency, 
     51                print '      ', paz.stage_sequence_number 
     52                meta = ChannelMeta( 
     53                    network=network.code, 
     54                    station=sta.code, 
     55                    location=chan.location_code, 
     56                    stream=chan.code[:-1], 
     57                    component=chan.code[-1], 
     58                    start=chan.start_date, 
     59                    end=chan.end_date, 
     60                    latitude=sta.latitude, 
     61                    longitude=sta.longitude, 
     62                    elevation=sta.elevation, 
     63                    gain=1.e9/chan.response.instrument_sensitivity.value, 
     64                    poles=str(paz.poles), 
     65                    zeros=str(paz.zeros), 
     66                ) 
     67                stations.add(meta, replace=True, local=True) 
    4568     
    46     #fp = open( 'metasave.pickle', 'w' ) 
    47     #pickle.dump( inv, fp ) 
    48     #fp.close() 
     69    fp = open( 'metasave%s.pickle' % net, 'w' ) 
     70    pickle.dump( inv, fp ) 
     71    fp.close() 
    4972     
     73if __name__ == '__main__': 
     74 
     75    for net in ('SX','TH','GE', 'BW'): #'GR' 
     76        copyFdsnMeta( net, 'BGR' ) 
     77 
Note: See TracChangeset for help on using the changeset viewer.