【Jupyter】基本操作

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格式写入工作区