Changes for page FDSN Guide

Last modified by robert on 2025/03/24 12:02

From version 1.2
edited by robert
on 2025/03/24 10:05
Change comment: There is no comment for this version
To version 1.3
edited by robert
on 2025/03/24 10:07
Change comment: There is no comment for this version

Summary

Details

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