|
1. 读取示例地震数据
from obspy import read #导入read包
st = read() #read example
print(st) #print event ifo
st.plot(); #绘制地震数据
tr = st[0] #将st中第一个数据赋给tr
print(tr.stats) #输出tr这一道的头段信息
2. 从IRIS中下载指定发震时间的地震
from obspy.clients.fsdn import Client
client = Client('IRIS')
from obspy import UTCDateTime
time = UTCDateTime('2021-07-29T06:15:47') #指定心仪的地震事件，先从IRIS上选定
starttime = time + 600 #指定从哪个时间开始下载地震数据，s
endtime = starttime + 60*20 #数据结束相对时长，20min
net = 'IC' #指定台网
sta = 'ENH' #指定台站
loc = '00' #指定台站中的00台
chan = 'BH*' #指定地震道
st = client.get_waveforms(net,sta,loc,chan,starttime,endtime,attach_response = True) #最后一个为仪器响应
st.plot()
st_rem = st.copy()
st_rem.remove_response(output = 'VEL')
st[0].plot();
st_rem[0].plot(); #对比原始波形与去除仪器响应波形
st_rem = st.copy()
st_rem[0].remove_response(output = 'VEL',plot = True) #查看仪器响应，频谱与去除仪器响应的过程
filename = sta + '_' + chan + '.mseed'
st_rem.write(filename, formate ='MSEED')
#将内存中的地震数据以.mseed格式写入工作区
以上是08月02日练习的内容
|