INGV station instrument response error #170
Description
Dear INGV colleagues,
I am presently doing an analysis using many seismic stations across Europe, using the vertical component.
I am working in obspy, downloading inventory, then downloading data, then simple preprocessing.
All works well for INGV stations**, except for station CSOB.** The program fails when deconvolving with the instrument response.
Here is an example that doesn't work. I am working on this stream:
IV.CSOB..HHZ | 2022-02-01T22:43:07.300000Z - 2022-02-01T22:49:47.300000Z | 100.0 Hz, 40001 samples
There may be errors also on the other components and for other times.
The processing is simple : :
def preproc(st,fmin,fmax,inv) :
st_plot_filt=st.copy()
st_plot_filt=st_plot_filt.detrend("linear")
st_plot_filt.taper(max_percentage=0.04)
st_plot_filt.remove_response(output="VEL", water_level=10, inventory=inv)
st_plot_filt.filter("bandpass", freqmin=fmin, freqmax=fmax)
st_plot_filt.normalize()
return st_plot_filt
It fails at the line
----> 5 st_plot_filt_ingv=preproc(st_final_ingv,fmin,fmax, inventory_ingv)
and looking through the errors, comes with
ValueError: check_channel: Illegal RESP format
Would it be able to fix this issue? It makes it quite a lot more complex to download data from a lot of INGV stations into a single stream, then deconvolving.
I will of course make an exception in my present program, but I'd be grateful if you could drop me a note if/when the issue is fixed.
Just contact me if you have any questions
Best regards
Helle Pedersen
Summary
- IV.CSOB
- MN.BNI