Skip to Content

Changes for page FDSN Guide

Last modified by robert on 2025/03/24 12:02
From version 2.3
edited by robert
on 2025/03/24 10:21
Change comment: There is no comment for this version
To version 2.4
edited by robert
on 2025/03/24 10:30
Change comment: There is no comment for this version

Summary

Details

Page properties
Content
... ... @@ -30,10 +30,43 @@
30 30  Information such as site locations, sensor and data logger types, response information, etc are in the station metadata. This can be accessed directly(link) or via the obspy get_stations (link) tool.
31 31  
32 32  
33 -== Sub-paragraph ==
33 +== How to download event, station, instrument response ==
34 34  
35 -Lorem ipsum dolor sit amet, consectetur adipiscing elit, sed do eiusmod tempor incididunt ut labore et dolore magna aliqua. Ut enim ad minim veniam, quis nostrud exercitation ullamco laboris nisi ut aliquip ex ea commodo consequat. Duis aute irure dolor in reprehenderit in voluptate velit esse cillum dolore eu fugiat nulla pariatur. Excepteur sint occaecat cupidatat non proident, sunt in culpa qui officia deserunt mollit anim id est laborum.
36 36  
36 +{{code language="python"}}
37 +import obspy
38 +from obspy.clients.fdsn import Client
39 +
40 +# Use AusPass client for station, waveform, and earthquake information
41 +client = Client("AUSPASS")
42 +
43 +
44 +# Download station information for AUMTC station in S1 network at the response level
45 +inv = client.get_stations(network="S1", station="AUMTC", location="*",
46 + channel="*", starttime=event_time - 60,
47 + endtime=event_time + 1000, level="response")
48 +print(inv)
49 +
50 +# Inventory metadata is stored in a Inventory > Network > Station > Channel hierarchy
51 +
52 +print(inv) #inventory level
53 +
54 +print(inv[0]) # network level (the first network in the inventory)
55 +
56 +print(inv[0][0]) # station level (the first station of the first network in the inventory)
57 +
58 +print(inv[0][0][0]) # channel level (the first channel of the first station of the first network in the inventoy)
59 +
60 +# you can also select items directly
61 +
62 +print(inv.select(station='AUMTC',channel='HHZ')[0][0][0])
63 +
64 +# instrument response is attached to a channel object
65 +
66 +response = inv.select(station='AUMTC',channel='HHZ')[0][0][0].response
67 +{{/code}}
68 +
69 +
37 37  === Sub-sub paragraph ===
38 38  
39 39  Lorem ipsum dolor sit amet, consectetur adipiscing elit, sed do eiusmod tempor incididunt ut labore et dolore magna aliqua. Ut enim ad minim veniam, quis nostrud exercitation ullamco laboris nisi ut aliquip ex ea commodo consequat. Duis aute irure dolor in reprehenderit in voluptate velit esse cillum dolore eu fugiat nulla pariatur. Excepteur sint occaecat cupidatat non proident, sunt in culpa qui officia deserunt mollit anim id est laborum.
... ... @@ -241,4 +241,3 @@
241 241  catalog.plot(projection="local", title="Australia Earthquakes", resolution="i")
242 242  {{/code}}
243 243  
244 -