Changes for page FDSN Guide

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		-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.
	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,10 +41,91 @@
41	41
42	42	= Waveform Data =
43	43
44		-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.
	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		-== Sub-paragraph ==
47	47
48		-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.
49		-
50	50	= Earthquake Data =
	48	+
	49	+
	50	+
	51	+== How to download an Earthquake Catalog ==
	52	+
	53	+{{code language="python"}}
	54	+from obspy.clients.fdsn import Client
	55	+from obspy import UTCDateTime
	56	+
	57	+# Initialize the AusPass FDSN client
	58	+client = Client("AUSPASS")
	59	+
	60	+# Define the time range for the earthquake catalog
	61	+start_time = UTCDateTime("2021-08-01")
	62	+end_time = UTCDateTime("2022-01-01") # End of year
	63	+
	64	+# Define the geographic region (latitude and longitude for Woodspoint, Victoria, Australia)
	65	+latitude = -37.47
	66	+longitude = 146.10
	67	+max_radius = 5 # in degrees
	68	+
	69	+# Download the earthquake catalog
	70	+catalog = client.get_events(starttime=start_time, endtime=end_time,
	71	+ minmagnitude=2, latitude=latitude, longitude=longitude,
	72	+ maxradius=max_radius)
	73	+
	74	+# Save the catalog to a file (e.g., QuakeML format)
	75	+catalog.write("Woodspoint_earthquakes.xml", format="QUAKEML")
	76	+{{/code}}
	77	+
	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	+
	131	+