obspy.io.sac.sacpz.attach_resp

attach_resp(tr, resp_file, todisp=False, tovel=False, torad=False, tohz=False)[source]

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.io.sac.sacio.attach_paz()). Attaches to a trace a paz AttribDict containing poles, zeros, and gain.

Parameters:
  • tr An ObsPy Trace object
  • resp_file path to RESP-file or file pointer
  • todisp change a velocity transfer function to a displacement transfer function by adding another zero
  • tovel change a displacement transfer function to a velocity transfer function by removing one 0,0j zero
  • torad change to radians
  • tohz change to Hertz
>>> from obspy import Trace
>>> import os
>>> 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)
>>> for k in sorted(tr.stats.paz):  
...     print(k)
digitizer_gain
gain
poles
seismometer_gain
sensitivity
t_shift
zeros
>>> print(tr.stats.paz.poles)  
[(-0.15931...+0.15931...j), (-0.15931...-0.15931...j),
 (-314.159...+202.318...j), (-314.159...-202.318...j)]