Python et Obspy

Python et Obspy

Il est possible d’utiliser la librairie Obspy pour programmer la récupération des données dans des scripts  Python .

# Récupération de une heure de donnée d’une station dans le réseau FR, et affichage de la forme d’onde.

from obspy.fdsn import Client
from obspy import UTCDateTime
client = Client("RESIF")
t = UTCDateTime("2014-01-28T00:00:00.000")
stream = client.get_waveforms("FR", "OGDI", "00", "HHZ", t, t + 60 * 60)
stream.plot()

 

# Récupération de l'inventaire des channels du réseau FR

from obspy.fdsn import Client
client = Client("RESIF")
inventory = client.get_stations(network='FR',level='channel')
inventory.get_contents()

Il est également possible de requêter via les services web en utilisant d’autres librairies que Obspy.

# on exécute une requête de type wget, et on écrit les données sur le disque local

import urllib2
request = urllib2.Request('http://ws.resif.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() )

Techniquement , les accès aux données restreintes sont basée sur une authentification de type HTTP digest. Tout client qui implémente cette technologie sera capable de requêter des données protégées.

Rechercher