A simple python example

This script downloads, prints and plots an hour of seismic data recorded at the AUAYR station (AUSIS array, network code: S) on the 29th March 2018. The obspy.read( ) function makes a standard FDSN data request to access AusPass data, see:

  • our FDSN tutorial for a summary of the FDSN webservices implemented on the AusPass data server, and,
  • the FDSN website for a complete description of the FDSN specification and webservices.

Requirements :

  • Having a working internet connection
  • Having at least Python 2.7 installed. All codes are tested in Python 3.
  • Having the ObsPy module installed. Use a package manager (such as conda or pip) to install all dependencies (numpy, SQLAlchemy, lxml, Cython, future, proj4, matplotlib and cartopy or basemap). Detailed instructions are given on https://www.obspy.org/.

Description of obspy.read( ):

ElementMeaning
http://auspass.edu.au:80/Server name (AusPass website serving data on port 80)
fdsnws/FSDN Web Services on that server
dataselect/1/Data retrieval services (distinct for station metadata and event metadata)
queryDo a download
?Marks the end of the website adress serving the data
&Separates the arguments to select AusPass data
net=SSeismic network "S" - AUstralian Seismometers In Schools - AUSIS
sta=AUAYRStation "AUAYR" located at Ayr, Queensland, Australia
cha=BH?All broadband channels available starting with "BH" ie. BHE, BHN and BHZ
start=2018-03-29T20:00:00.000Data start time using ISO standard date and time format
end=2018-03-29T21:00:00.000Data end time using ISO standard date and time format
In [1]:
import obspy
x = obspy.read('http://auspass.edu.au:80/fdsnws/dataselect/1/query?net=S&sta=AUAYR&cha=BH?&start=2018-03-29T20:00:00.000&end=2018-03-29T21:00:00.000')
print(x)
x.plot()
3 Trace(s) in Stream:
S.AUAYR..BHE | 2018-03-29T19:59:48.000000Z - 2018-03-29T21:00:14.500000Z | 10.0 Hz, 36266 samples
S.AUAYR..BHN | 2018-03-29T19:59:57.800000Z - 2018-03-29T21:00:12.700000Z | 10.0 Hz, 36150 samples
S.AUAYR..BHZ | 2018-03-29T19:59:41.000000Z - 2018-03-29T21:00:09.000000Z | 10.0 Hz, 36281 samples