iris_pull.py 文件源码

def plot_responses(options, config, event_names):
    import pylab as lab

    conf = config['iris_pull_config']
    if not event_names:
        sys.exit('need event name')

    for event_name in event_names:
        conf.event_name = event_name
        stations = _get_stations(conf)
        event = _get_event_infos(conf)
        fband = conf.inv_response_frequencyband


        combi_get_responses(stations, event.time, conf.path('resp_path'))

        for station in stations:

            for cha in station.get_channels():
                resp = cha.inv_response

                fmin, fmax = fband[0], fband[3]
                freqs = num.exp(num.linspace(num.log(fmin), num.log(fmax), 400))
                amps = num.abs(resp.evaluate(freqs))

                lab.loglog(freqs, amps)


    lab.show()