Esempio n. 1
0
def make_irisrequest(station, begin, end, channel, loc):

    print 'IRISREQUEST FOR STATION ', station

    net, sta = station.split('_')
    url = 'http://service.iris.edu/fdsnws/dataselect/1/query?'
    st = Stream()

    for i in channel:

        parameter = urllib.urlencode({
            'net': net,
            'sta': sta,
            'loc': loc,
            #'cha': i,                               #hs v2
            'starttime': begin,
            'endtime': end,
            'nodata': '404',
        })

        u = ('%s%s') % (url, parameter)
        saveUrl(station, u)

        #data = urllib.urlopen(u).read()                          #hs
        try:
            data = urllib.urlopen(u).read()  #hs
        except:
            continue

        #tf  = NamedTemporaryFile()                               #hs
        tf = NamedTemporaryFile(suffix='xtmp')  #hs
        tf.write(data)
        tf.seek(0)

        t = isMSEED(tf.name)

        if t == True:
            st += read(tf.name, 'MSEED')

        tf.close()
        os.remove(tf.name)
        break  #hs v2

#       if t == True:     size = proof_file_v3(st,channel)      #hs
    if len(st) > 0: size = proof_file_v3(st, channel)  #hs
    else: size = 0

    print 'SIZE: ----> ', size
    return size
Esempio n. 2
0
def make_irisrequest(station, begin, end, channel, loc):

    logger.info('\033[31m IRISREQUEST \033[0m\n' % ())

    net, sta = station.split('_')
    url = 'http://service.iris.edu/fdsnws/dataselect/1/query?'
    st = Stream()
    for i in channel:
        parameter = urllib.urlencode({
            'net': net,
            'sta': sta,
            'loc': loc,
            'cha': i,
            'starttime': begin,
            'endtime': end,
            'nodata': '404',
        })
        u = ('%s%s') % (url, parameter)
        data = urllib.urlopen(u).read()
        tf = NamedTemporaryFile()
        tf.write(data)
        tf.seek(0)
        st += read(tf.name, 'MSEED')
        tf.close()

    print st
    output = ('%s.%s-IRIS.mseed') % (net, sta)
    st.write(output, format='MSEED')

    size = proof_file_v2(output, channel)
    #print 'SIZE: ----> ',size
    try:
        #os.remove(fname)
        os.remove(output)
    except:
        print ''

    return size