Instrument correction for Iris data

Question

Hi

I was wondering if it might be possible to see a script where someone has downloaded and applied an instrument correction for an iris seismometer.

I am able to download the data, and the RESP file using the client but when I try and use this file with invsim.evalresp I get an syntax error. My invsim.evalresp line is as follows:

invsim.evalresp(0.025, 512, IU_RAR_response.txt, t1, station='RAR', channel='BHZ', network='IU', locid='10', units='ACC', freq=False)

Is it possible to use the RESP file directly within seisSim, or is it always neccessary to pull out the relevant numbers and poles using evalresp.

Any help would be much appreciated.

Thanks' Dan

该提问来源于开源项目：obspy/obspy

Answer 1

[yanmail] Hi Dan,

I apologize for the lack of documentation on this function. I've been meaning to write a short tutorial for a while but haven't gotten around to it. To answer your question, here is a code snippet on how to use the evalresp option: {{{

!python

from obspy.core import read, UTCDateTime from obspy.signal import seisSim rawf = 'CRLZ.HHZ.SAC' respf = 'RESP.NZ.CRLZ.10.HHZ'

read in your raw data

tr = read(rawf)[0]

define a filter band to

prevent amplifying noise

during the deconvolution

fl1=0.00588 fl2=0.00625 fl3=30. fl4=35.

this can be the date of your raw data or

any date for which the SEED Resp-file is valid

date=UTCDateTime(2003,11,1,0,0,0) seedresp = {'filename':respf,'date':date,'units':'VEL'} tr.data = seisSim(tr.data,tr.stats.sampling_rate,paz_remove=None, remove_sensitivity=False,pre_filt=(fl1,fl2,fl3,fl4), seedresp=seedresp)

}}}

The evalresp function is called from seisSim if the seedresp dictionary is not None. This functionality is still in the testing phase, so you should check the results. Hope this helps.

Yannik