Skip to Content

Changes for page FDSN Guide

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

Summary

Details

Page properties
Content
... ... @@ -30,48 +30,21 @@
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 -== How to download event, station, instrument response ==
33 +== Sub-paragraph ==
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.
35 35  
36 -{{code language="python"}}
37 -import obspy
38 -from obspy.clients.fdsn import Client
37 +=== Sub-sub paragraph ===
39 39  
40 -# Use AusPass client for station, waveform, and earthquake information
41 -client = Client("AUSPASS")
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.
42 42  
43 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="*", level="response")
47 -
48 -# Inventory metadata is stored in a Inventory > Network > Station > Channel hierarchy
49 -
50 -print(inv) #inventory level
51 -
52 -print(inv[0]) # network level (the first network in the inventory)
53 -
54 -print(inv[0][0]) # station level (the first station of the first network in the inventory)
55 -
56 -print(inv[0][0][0]) # channel level (the first channel of the first station of the first network in the inventoy)
57 -
58 -# you can also select items directly
59 -
60 -print(inv.select(station='AUMTC',channel='HHZ')[0][0][0])
61 -
62 -# instrument response is attached to a channel object
63 -
64 -response = inv.select(station='AUMTC',channel='HHZ')[0][0][0].response
65 -{{/code}}
66 -
67 67  = Waveform Data =
68 68  
69 69  Waveform data (e.g. the actual seismic data) can be accessed directly (link) or via obspy's get_waveforms (link) tool. It can also be accessed via various tools such as seed-vault, pyweed, etc (add links).
70 70  
71 -== Downloading and Storing data ==
46 +== How to download waveform data ==
72 72  
73 -=== How to download waveform data ===
74 -
75 75  {{code language="python"}}
76 76  from obspy import UTCDateTime
77 77  from obspy.clients.fdsn import Client
... ... @@ -93,75 +93,9 @@
93 93  print("Downloaded and saved the MiniSEED file.")
94 94  {{/code}}
95 95  
96 -=== How to download a LOT of waveform data ===
69 +== How to remove instrument response ==
97 97  
98 98  {{code language="python"}}
99 -from obspy import UTCDateTime
100 -from obspy.clients.fdsn import Client
101 -import datetime
102 -
103 -# Initialize client and set parameters
104 -client = Client("IRIS")
105 -start_date = UTCDateTime("2023-07-01")
106 -end_date = UTCDateTime("2023-07-02")
107 -network = "S1"
108 -station = "AUMTS"
109 -channel = "BHZ"
110 -
111 -# Loop to download data one day at a time
112 -while start_date <= end_date:
113 - next_date = start_date + datetime.timedelta(days=1)
114 - try:
115 - st = client.get_waveforms(network, station, "*", channel, start_date, next_date)
116 - st.write(f"{start_date.date}.mseed", format="MSEED")
117 - except:
118 - print(f"Failed for {start_date} to {next_date}")
119 - start_date = next_date
120 -{{/code}}
121 -
122 -=== How to store and archive waveform data in SeisComP Data Structure (SDS) ===
123 -
124 -{{code language="python"}}
125 -import os
126 -from obspy import UTCDateTime, read
127 -from obspy.clients.fdsn import Client
128 -
129 -# Initialize the client
130 -client = Client("AUSPASS") # Replace with the correct client endpoint if different
131 -
132 -# Define event time and time window
133 -event_time = UTCDateTime("2021-09-21T23:15:53")
134 -starttime = event_time - 600 # 10 minutes before the event
135 -endtime = event_time + 1800 # 30 minutes after the event
136 -
137 -# Download the waveform data
138 -st = client.get_waveforms(network="S1", station="AUMTS", location="*", channel="*", starttime=starttime, endtime=endtime)
139 -
140 -# Create SDS structure: ROOT/YEAR/NET/STA/CHAN.TYPE/NET.STA.LOC.CHAN.YEAR.DAY
141 -sds_root = "." # Replace with your desired directory
142 -
143 -for tr in st:
144 - net = tr.stats.network
145 - sta = tr.stats.station
146 - loc = tr.stats.location
147 - chan = tr.stats.channel
148 - year = str(tr.stats.starttime.year)
149 - jday = str(tr.stats.starttime.julday).zfill(3)
150 -
151 - sds_path = os.path.join(sds_root, year, net, sta, f"{chan}.D", f"{net}.{sta}.{loc}.{chan}.{year}.{jday}")
152 -
153 - # Create directories if they don't exist
154 - os.makedirs(os.path.dirname(sds_path), exist_ok=True)
155 -
156 - # Save the trace as a MiniSEED file
157 - tr.write(sds_path, format="MSEED")
158 -{{/code}}
159 -
160 -== Common Data Operations ==
161 -
162 -=== How to remove instrument response ===
163 -
164 -{{code language="python"}}
165 165  from obspy import read
166 166  from obspy.core.util import AttribDict
167 167  
... ... @@ -183,7 +183,7 @@
183 183  st.write("Woodspoint_2021_corrected.mseed", format="MSEED")
184 184  {{/code}}
185 185  
186 -=== How to apply a bandpass filter ===
93 +== How to apply a bandpass filter ==
187 187  
188 188  {{code language="python"}}
189 189  from obspy import read
... ... @@ -203,7 +203,7 @@
203 203  st.write("Woodspoint_2021_filtered.mseed", format="MSEED")
204 204  {{/code}}
205 205  
206 -=== How to slice a waveform ===
113 +== How to slice a waveform ==
207 207  
208 208  {{code language="python"}}
209 209  from obspy import read, UTCDateTime, Stream # Importing Stream here
... ... @@ -225,7 +225,7 @@
225 225  sliced_st.write("Woodspoint_2021_filtered_sliced.mseed", format="MSEED")
226 226  {{/code}}
227 227  
228 -=== How to save a waveform ===
135 +== How to save a waveform ==
229 229  
230 230  {{code language="python"}}
231 231  # Save the sliced file as MiniSEED
... ... @@ -235,7 +235,7 @@
235 235  sliced_st.write("Woodspoint_2021_filtered_sliced.sac", format="SAC")
236 236  {{/code}}
237 237  
238 -=== How to convert miniSEED to SAC ===
145 +== How to convert miniseed to sac ==
239 239  
240 240  {{code language="python"}}
241 241  from obspy import read
... ... @@ -253,7 +253,6 @@
253 253  
254 254  = Earthquake Data =
255 255  
256 -Earthquake data can be accessed directly or via ObsPy's get_events code
257 257  
258 258  == How to download an Earthquake Catalog ==
259 259