def test_evalresp(self, plot=False): resp_fpath = common.test_data_file('test2.resp') freqs = num.logspace(num.log10(0.001), num.log10(10.), num=1000) transfer = evalresp.evalresp( sta_list='BSEG', cha_list='BHZ', net_code='GR', locid='', instant=util.str_to_time('2012-01-01 00:00:00'), freqs=freqs, units='DIS', file=resp_fpath, rtype='CS')[0][4] pz_fpath = common.test_data_file('test2.sacpz') zeros, poles, constant = pz.read_sac_zpk(pz_fpath) resp = trace.PoleZeroResponse(zeros, poles, constant) transfer2 = resp.evaluate(freqs) if plot: plot_tfs(freqs, [transfer, transfer2]) assert numeq(transfer, transfer2, 1e-4)
def test_evalresp(self, plot=False): testdir = os.path.dirname(__file__) freqs = num.logspace(num.log10(0.001), num.log10(10.), num=1000) transfer = evalresp.evalresp(sta_list='BSEG', cha_list='BHZ', net_code='GR', locid='', instant=util.str_to_time('2012-01-01 00:00:00'), freqs=freqs, units='DIS', file=os.path.join(testdir, 'response', 'RESP.GR.BSEG..BHZ'), rtype='CS')[0][4] pzfn = 'SAC_PZs_GR_BSEG_BHZ__2008.254.00.00.00.0000_2599.365.23.59.59.99999' zeros, poles, constant = pz.read_sac_zpk(filename=os.path.join( testdir, 'response', pzfn)) resp = trace.PoleZeroResponse(zeros, poles, constant) transfer2 = resp.evaluate(freqs) if plot: import pylab as lab lab.plot(freqs, num.imag(transfer)) lab.plot(freqs, num.imag(transfer2)) lab.gca().loglog() lab.show() assert numeq(transfer, transfer2, 1e-4)