我想从Tungurahua火山周围地区下载地震数据。在功能obspy.clients.fdsn的obspy文档中,描述了下载数据的步骤:
from obspy import UTCDateTime
t = UTCDateTime("2010-02-27T06:45:00.000")
st = client.get_waveforms("IU", "ANMO", "00", "LHZ", t, t + 60 * 60)
st.plot()
此处get_waveforms方法具有参数
使用函数client.get_stations我在我的案例中找到了感兴趣的网络和电台名称:
但是我不知道要传递给函数get_waveforms
所需的位置/通道字符串。
我尝试过的事情:
client.get_stations(network="YO", station="TBAG")
返回:
Inventory created at 2017-04-23T18:27:16.000000Z
Created by: IRIS WEB SERVICE: fdsnws-station | version: 1.1.25
http://service.iris.edu/fdsnws/station/1/query?station=TBAG&network...
Sending institution: IRIS-DMC (IRIS-DMC)
Contains:
Networks (1):
YO
Stations (1):
YO.TBAG (Banos, Tungurahua, Ecuador)
Channels (0):
所以似乎没有渠道存在。但是在尝试时
client.get_stations(network="IU", station="ANMO")
我得到了
Inventory created at 2017-04-23T18:40:22.000000Z
Created by: IRIS WEB SERVICE: fdsnws-station | version: 1.1.25
http://service.iris.edu/fdsnws/station/1/query?station=ANMO&network...
Sending institution: IRIS-DMC (IRIS-DMC)
Contains:
Networks (1):
IU
Stations (6):
IU.ANMO (Albuquerque, New Mexico, USA)
IU.ANMO (Albuquerque, New Mexico, USA)
IU.ANMO (Albuquerque, New Mexico, USA)
IU.ANMO (Albuquerque, New Mexico, USA)
IU.ANMO (Albuquerque, New Mexico, USA)
IU.ANMO (Albuquerque, New Mexico, USA)
Channels (0):
所以频道" LHZ"虽然很明显存在,但这里没有列出。
答案 0 :(得分:2)
您需要使用level
参数:
client.get_stations(network="YO", station="TBAG", level="channel")
然后你应该得到:
Inventory created at 2017-05-09T12:58:29.000000Z
Created by: IRIS WEB SERVICE: fdsnws-station | version: 1.1.25
http://service.iris.edu/fdsnws/station/1/query?format=xml&network=Y...
Sending institution: IRIS-DMC (IRIS-DMC)
Contains:
Networks (1):
YO
Stations (1):
YO.TBAG (Banos, Tungurahua, Ecuador)
Channels (5):
YO.TBAG..HDF, YO.TBAG..HHZ, YO.TBAG..HHN, YO.TBAG..HHE,
YO.TBAG..LOG
此参数在the documentation:
中描述level (str)
Specify the level of detail for the results (“network”, “station”, “channel”,
“response”), e.g. specify “response” to get full information including
instrument response for each channel.