[obspy-users] Instrument Response Removal Error

ayodeji.kuponiyi at gmail.com ayodeji.kuponiyi at gmail.com
Thu Jun 18 03:46:37 CEST 2015


> Hello,
> 
> I started using obspy very recently and found it a lot easy to use and beneficial - thanks for the program!
> 
> I have a problem with instrument response removal and was not sure what went wrong. I tried the same routine on another station file and it was successful. Please find below the commands I used. I also attach the instrument response file for the station. 
> 
> >>> tr = read('YOSQ.CN.00.CHN.2010.059')
> >>> st = tr.select(component="Z")
> >>> st.detrend(type="demean")
> >>> st_demean = st.copy()
> >>> st.detrend(type="linear")
> >>> st_detrend = st.copy()
> >>> st.filter('bandpass', freqmin=0.005, freqmax=1.6, corners=6,zerophase=True)
> >>> st_filt = st.copy()
> >>> t1 = st[0].stats.starttime
> >>> respfile = 'RESP.PO.YOSQ..HHZ'
> >>> seedresp = {'filename': respfile, 'date': t1, 'units':'VEL'}
> >>> st.simulate(paz_remove=None, remove_sensitivity=False,seedresp=seedresp)
> 
> error
>  WARNING: no response found for NET=CN,STA=YOSQ,LOCID=,CHAN=HHZ,DATE=2010,059
> Traceback (most recent call last):
>   File "<stdin>", line 1, in <module>
>   File "/home/akuponiyi/anaconda/lib/python2.7/site-packages/obspy-0.10.2-py2.7-linux-x86_64.egg/obspy/core/stream.py", line 1930, in simulate
>     simulate_sensitivity=simulate_sensitivity, **kwargs)
>   File "/home/akuponiyi/anaconda/lib/python2.7/site-packages/obspy-0.10.2-py2.7-linux-x86_64.egg/obspy/core/trace.py", line 231, in new_func
>     result = func(*args, **kwargs)
>   File "/home/akuponiyi/anaconda/lib/python2.7/site-packages/obspy-0.10.2-py2.7-linux-x86_64.egg/obspy/core/trace.py", line 1291, in simulate
>     simulate_sensitivity=simulate_sensitivity, **kwargs)
>   File "/home/akuponiyi/anaconda/lib/python2.7/site-packages/obspy-0.10.2-py2.7-linux-x86_64.egg/obspy/signal/invsim.py", line 538, in seisSim
>     channel=seedresp['channel'])
>   File "/home/akuponiyi/anaconda/lib/python2.7/site-packages/obspy-0.10.2-py2.7-linux-x86_64.egg/obspy/signal/invsim.py", line 286, in evalresp
>     raise ValueError(msg)
> ValueError: evalresp failed to calculate a response.
> 
> <RESP.PO.YOSQ..HHZ>


More information about the obspy-users mailing list