May
1
Hello
I wanted to remove instrument response from my data in an offline manner.
I just simply used obspy
from obspy import read, read_inventory
st = read(‘AM.R5868.00.EHZ.D.2022.044’)
inv = read_inventory(“AM.R5868.00.xml”)
pre_filt = [1, 5, 40, 45]
tr.remove_response(inventory=inv, pre_filt=pre_filt, output=“VEL”,water_level=60, plot=True)
then I got a errormessage
ValueError: No matching response information found.
just wondering what could cause this? a format problem or?
Thank you very much for help.
May
Hello May, welcome to our community!
Possibly, these two links could be useful for the issue you are experiencing:
https://manual.raspberryshake.org/metadata.html?#example-obspy-code-for-removing-response
https://docs.obspy.org/packages/autogen/obspy.core.trace.Trace.remove_response.html
I think you are right, it is likely a format problem. We provide different response .xml files so one of these could be what you need. They are in:
-
sc3ml format: https://fdsnws.raspberryshakedata.com/fdsnws/station/1/query?network=AM&station=R5868&level=resp&format=sc3ml
-
fdsnxml format: https://fdsnws.raspberryshakedata.com/fdsnws/station/1/query?network=AM&station=R5868&level=resp
-
dataless format (available for some Shake models/versions): Metadata - Instrument Response Files — Instructions on Setting Up Your Raspberry Shake
It is possible that the one that you are using is not recognized by Obspy, and maybe another one from the list above will work.
May
3
Thank you very much. I solved the problem! Thanks! May
2 Likes
No problem at all, you’re welcome!
May
5
Just because the “endDate” in the xml file was later than my start time in the trace;)
2 Likes