Changes for page FDSN Guide
Last modified by robert on 2025/03/24 12:02
Summary
-
Page properties (1 modified, 0 added, 0 removed)
Details
- Page properties
-
- Content
-
... ... @@ -27,7 +27,7 @@ 27 27 28 28 = Station Metadata = 29 29 30 - Lorempsumdolorsitamet,consecteturadipiscing elit, sed do eiusmod temporincididuntutlaboreetdoloremagna aliqua. Utenimad minim veniam, quis nostrud exercitationullamcolaborisnisiut aliquipex eacommodo consequat.Duis aute irure dolor inreprehenderitin voluptatevelit essecillumdoloreeu fugiatnulla pariatur. Excepteur sintoccaecatcupidatatnonproident,suntin culpaquifficia deseruntmollit animid estlaborum.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 33 == Sub-paragraph == ... ... @@ -41,11 +41,13 @@ 41 41 42 42 = Waveform Data = 43 43 44 - Loremipsumdolor sitamet,consecteturadipiscing elit,sed do eiusmod temporincididuntut labore etdoloremagna aliqua. Ut enimadminimveniam,quis nostrudexercitation ullamcolaborisnisiutaliquip ex ea commodoconsequat.Duisauteirure dolorin reprehenderitin voluptatevelit essecillumdoloreeu fugiatnullapariatur.Excepteursintoccaecatcupidatatnonproident,suntinculpaqui officiaeseruntmollit anim id est laborum.44 +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). 45 45 46 46 47 47 = Earthquake Data = 48 48 49 + 50 + 49 49 == How to download an Earthquake Catalog == 50 50 51 51 {{code language="python"}} ... ... @@ -98,4 +98,32 @@ 98 98 #catalog.plot(), ObsPy automatically uses the depth information to color the events in the plot 99 99 {{/code}} 100 100 103 +== How to plot (Local) Earthquakes == 104 + 105 +{{code language="python"}} 106 +from obspy import UTCDateTime 107 +from obspy.clients.fdsn import Client 108 + 109 +# Initialize FDSN client 110 +client = Client("AUSPASS") 111 + 112 +# Define time range 113 +starttime = UTCDateTime("2023-01-01") 114 +endtime = UTCDateTime() 115 + 116 +# Latitude and longitude bounds for Australia 117 +minlatitude = -44.0 118 +maxlatitude = -10.0 119 +minlongitude = 113.0 120 +maxlongitude = 154.0 121 + 122 +# Fetch event data for Australia with a minimum magnitude 123 +catalog = client.get_events(starttime=starttime, endtime=endtime, minmagnitude=4, 124 + minlatitude=minlatitude, maxlatitude=maxlatitude, 125 + minlongitude=minlongitude, maxlongitude=maxlongitude) 126 + 127 +# Plot the earthquakes 128 +catalog.plot(projection="local", title="Australia Earthquakes", resolution="i") 129 +{{/code}} 130 + 101 101