Python and obspy examples – EPOS-FRANCE Seismology

You may use Obspy library to natively handle seismic data into your Python scripts.

# this will download 1 hour of data from a FR network station, and display a wafeform plotfrom obspy.fdsn import Client
from obspy import UTCDateTime
client = Client("EPOSFR")
t = UTCDateTime("2026-01-28T00:00:00.000")
stream = client.get_waveforms("FR", "OGDI", "00", "HHZ", t, t + 60 * 60)
stream.plot()

# this will show an inventory for FR network from obspy.fdsn import Client
client = Client("EPOSFR")
inventory = client.get_stations(network='FR',level='channel')
inventory.get_contents()

You may also request web services with one of the numerous Python grabber library [5] Unlike Obspy, these libraries are seismic-data agnostic (like wget is):

# this will make a wget-like request and write resulting file on local disk import urllib2

request = urllib2.Request("https://seisdata.epos-france.fr/fdsnws/dataselect/1/query?network=FR&station=OGDI&channel=HHZ&starttime=2014-01-10T00:00:00&endtime=2014-01-11T00:00:00")

response = urllib2.urlopen(request)
newfile = open('./mydata.mseed', 'w')
newfile.write ( response.read() )

Technically, restricted access works with HTTP digest authentication. Any HTTP client implementing this technology should be able to request restricted-access data.

Search