通过 fdsn 和 obspy 访问地震数据(使用网络、站点、位置、通道)
Accessing seismographic data via fdsn with obspy (using network, station, location, channel)
我想下载通古拉瓦火山周围地区的地震数据。在函数 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 方法有参数
- 网络="IU"
- 车站="ANMO"
- 位置=“00”
- 频道="LHZ"
使用函数 client.get_stations 我找到了我感兴趣的网络和站点名称:
- 网络="YO"
- 车站="TBAG"
但是我不知道需要传递给函数的位置/频道字符串 get_waveforms
。
我尝试过的:
client.get_stations(network="YO", station="TBAG")
哪个returns:
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",尽管它显然存在。
您需要使用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.
我想下载通古拉瓦火山周围地区的地震数据。在函数 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 方法有参数
- 网络="IU"
- 车站="ANMO"
- 位置=“00”
- 频道="LHZ"
使用函数 client.get_stations 我找到了我感兴趣的网络和站点名称:
- 网络="YO"
- 车站="TBAG"
但是我不知道需要传递给函数的位置/频道字符串 get_waveforms
。
我尝试过的:
client.get_stations(network="YO", station="TBAG")
哪个returns:
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",尽管它显然存在。
您需要使用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.