Hi,
I would appreciate it if someone could help me out. I used the following lines to remove the instrument response, but I was wondering if my procedure is correct? I need to remove the instrument response of this record and check the correlation with my simulation up to 5Hz.
from obspy import read
from obspy.signal.invsim import corn_freq_2_paz
st = read(r"C:\Users\NAV\Desktop\Van018\Alfred_van2015-07-15\ALFS\ALFS_20150715_215948.seed")
st.detrend(type=“demean”)
st.detrend(type=“linear”)
st.detrend(type=“polynomial”, order=2) # Adjust the order as needed
paz_sts2 = {‘poles’: [-3.70080e-02-3.70080E-02j, -3.70080E-02+3.70080E-02j,
-1131+0.0j,
-5.02650E+02+0.0j, -1.00500E+03+0.0j],
‘zeros’: [0+0j, 0+0j],
‘gain’: 5.71508E8,
‘sensitivity’:7.99998E8}
pre_filt = (0.08, 0.1, 5.2, 6.0)
st.simulate(paz_remove=paz_sts2, pre_filt=pre_filt )
st.plot()
st.write(“ALFS_vel_removedInstrumentresponse”, format=“TSPAIR”, Units=‘m/s’)
I have used this line before, but it creates a significant discrepancy at low frequencies between my recordings and simulations.
paz_10hz = corn_freq_2_paz(10.0)
st.simulate(paz_remove=paz_sts2, paz_simulate=paz_10hz)
Thank you so much