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,12 +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 - 47 47 48 48 = Earthquake Data = 49 49 49 + 50 + 50 50 == How to download an Earthquake Catalog == 51 51 52 52 {{code language="python"}} ... ... @@ -74,4 +74,57 @@ 74 74 catalog.write("Woodspoint_earthquakes.xml", format="QUAKEML") 75 75 {{/code}} 76 76 78 +== How to plot (Global) Earthquakes == 79 + 80 +{{code language="python"}} 81 +from obspy import UTCDateTime 82 +from obspy.clients.fdsn import Client 83 + 84 +# Initialize FDSN client to connect to the IRIS data center 85 +client = Client("IRIS") 86 + 87 +# Set the time range for fetching earthquake data 88 +# Start time: January 1, 2023 89 +# End time: Current time 90 +starttime = UTCDateTime("2023-01-01") 91 +endtime = UTCDateTime() 92 + 93 +# Fetch earthquake events with a minimum magnitude of 7 94 +catalog = client.get_events(starttime=starttime, endtime=endtime, minmagnitude=7) 95 +#client.get_events(). This function returns a Catalog object that contains a list of Event objects. 96 +#Each Event object, in turn, has an Origins attribute that contains the depth information 97 + 98 +# Plot the fetched earthquake data using an orthographic projection 99 +catalog.plot(projection="ortho", title="Global Earthquakes with Magnitude >= 7 since 2023") 100 +#catalog.plot(), ObsPy automatically uses the depth information to color the events in the plot 101 +{{/code}} 102 + 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 + 77 77