Extract key instrument response information from a RESP file, which can be extracted from a dataless SEED volume by, for example, using the script obspy-dataless2resp or the rdseed program. At the moment, you have to determine yourself if the given response is for velocity or displacement and if the values are given in rad or Hz. This is still experimental code (see also documentation for obspy.sac.sacio.attach_paz()). Attaches to a trace a paz AttribDict containing poles, zeros, and gain.
Parameters: |
|
---|
>>> from obspy import Trace
>>> tr = Trace()
>>> respfile = os.path.join(os.path.dirname(__file__), 'tests', 'data',
... 'RESP.NZ.CRLZ.10.HHZ')
>>> attach_resp(tr, respfile, torad=True, todisp=False)
>>> print tr.stats.paz.keys()
['sensitivity', 'digitizer_gain', 'seismometer_gain', 'zeros', 'gain',
't_shift', 'poles']
>>> print tr.stats.paz.poles
[(-0.15931644664884559+0.15931644664884559j),
(-0.15931644664884559-0.15931644664884559j),
(-314.15926535897933+202.31856689118268j),
(-314.15926535897933-202.31856689118268j)]