Download + remove_response คลื่นไหวสะเทือน
Python
from obspy import UTCDateTime
from obspy.clients.fdsn import Client
############# load data ############
t1 = UTCDateTime("2025-03-28T06:20:00.000")
t2 = UTCDateTime("2025-03-28T06:40:00.000")
fdsn_client = Client('IRIS')
# Fetch waveforms and reponse from IRIS FDSN web service into ObsPy
# stream and inventory objects
st = fdsn_client.get_waveforms(
network='IU', station='CHTO', location='00', channel='BHZ',
starttime=t1, endtime=t2)
st.plot()
inv = fdsn_client.get_stations(
network='IU', station='CHTO', location='00', channel='BHZ',
starttime=t1, endtime=t2, level='response')
# define a filter band to prevent amplifying noise during the deconvolution
pre_filt = (0.005, 0.006, 30.0, 35.0)
st.remove_response(inventory=inv, output='DISP', pre_filt=pre_filt)
st.plot()
- 2
โฆษณา
- ดาวน์โหลดแอปพลิเคชัน
- © 2025 Blockdit
