TimeSeries FR EN
Ce service donne accès aux signaux temporels des stations sismologiques du réseau RESIF. Des traitements supplémentaires peuvent être réalisés sur les données.
Limitations
- La plage demandée ne peut excéder 31 jours.
- Seules les données non restreintes sont accessibles.
Fonctionnalités disponibles
- Filtres passe-bas, passe-haut ou passe-bande (optionnellement à phase nulle).
- Suppression de la moyenne ou de la tendance temporelle du signal.
- Mise à l’échelle du signal par un facteur constant.
- Déconvolution avec la réponse instrumentale (avec conversion d’unité et préfiltrage).
- Dérivée et intégration.
- Calcul de l’enveloppe du signal.
- Décimation (sous-échantillonnage).
- Sortie aux formats : miniSEED, SAC binary, ASCII, PNG ou JPEG plot.
Générateur d'URL
Documentation
Exemples
- ws.resif.fr/resifws/timeseries/1/query?net=RA&station=PYTO&cha=HN2&loc=02&demean&correct&start=2017-11-02T13:35:00&end=2017-11-02T13:40:00&format=ascii
- ws.resif.fr/resifws/timeseries/1/query?net=RA&station=PYTO&cha=HN2&loc=02&demean&correct&start=2017-11-02T13:35:00&end=2017-11-02T13:40:00&format=plot
Utilisation
/query? (channel-options) (date-range-options) (output-options) [filter-options] {plot-options} [nodata=404]
où :
channel-options :: (net=<network> sta=<station> loc=<location> cha=<channel>)
date-range-options :: (starttime=<date|durée>) (endtime=<date|durée>)
output-options :: (format=<ascii|mseed|sac|slist|tspair|plot>)
plot-options :: {showtitle=<TRUE|false>} {showscale=<TRUE|false>} {monochrome=<true|FALSE>}
plot-options :: {width=<400-2000>} {height=<200-2000>}
filter-options :: [demean=<TRUE|false>] [detrend=<true|FALSE>]
filter-options :: [taper=<0.0-0.5,HANNING|hamming|cosine>]
filter-options :: [decimate=<2-16>]
filter-options :: [correct=<true|FALSE>] {waterlevel=<nombre>} {freqlimits=<f1-f2-f3-f4>} {units=<AUTO|disp|vel|acc>}
filter-options :: [diff=<true|FALSE>] [int=<true|FALSE>]
filter-options :: [lpfilter=<fréquence>] [hpfilter=<fréquence>] [bpfilter=<fmin-fmax>] {zerophase=<true|FALSE>}
filter-options :: [envelope=<true|FALSE>]
filter-options :: [scale=<nombre>] [divscale=<nombre>]
(..) requis
[..] optionnel
{..} optionnel, mais toujours en complément d'une autre option
les valeurs par défaut sont en majusculesParamètres
Sélection de la source de données
| Paramètre | Alias | Description | Exemple |
|---|---|---|---|
| network | net | Code du réseau sismique. Accepte les jokers et les listes. | FR |
| station | sta | Code de la station. Accepte les jokers et les listes. | CIEL |
| location | loc | Code de localisation. Utilisez loc=-- pour les codes de localisations vides. Accepte les jokers et les listes. | 00 |
| channel | cha | Code du canal. Accepte les jokers et les listes. | HHZ |
Jokers et listes d’arguments
- Jokers : le point d’interrogation ? représente n’importe quel caractère unique, alors que l’astérisque * représente zéro caractère ou plus.
- Listes : plusieurs éléments peuvent être récupérés à l’aide d’une liste séparée par des virgules. Les jokers peuvent être inclus dans la liste.
Par exemple, pour le code des canaux : channel=EH?,BHZ
Sélection de la période temporelle
| Paramètre | Alias | Description | Exemple |
|---|---|---|---|
| starttime | start | Sélectionne les données à partir de l’heure spécifiée incluse. | 2010-01-10T00:00:00 |
| endtime | end | Sélectionne les données avant l’heure spécifiée incluse. | 2011-02-11T01:00:00 |
La définition de l’intervalle de temps avec starttime et endtime peut prendre différentes formes :
- une date, par exemple starttime=2015-08-12T01:00:00
- une durée en secondes, par exemple endtime=7200
- le mot-clé “currentutcday” qui signifie minuit de la date du jour (UTC), par exemple starttime=currentutcday
Formats des dates et des heures
YYYY-MM-DDThh:mm:ss[.ssssss] ex. 1997-01-31T12:04:32.123
YYYY-MM-DD ex. 1997-01-31 (une heure de 00:00:00 est supposée)
currentutcday (aujourd'hui à 00:00:00)
avec :
YYYY :: quatre chiffres de l'année
MM :: deux chiffres du mois (01=Janvier, etc.)
DD :: deux chiffres du jour du mois (01 à 31)
T :: séparateur date-heure
hh :: deux chiffres de l'heure (00 à 23)
mm :: deux chiffres des minutes (00 à 59)
ss :: deux chiffres des secondes (00 à 59)
ssssss :: un à six chiffres des microsecondes en base décimale (0 à 999999)Configuration du format de sortie
| Paramètre | Alias | Défaut | Description | Exemple |
|---|---|---|---|---|
| format | Format de sortie du fichier : ascii (équivalent à tspair), mseed, sac, slist ou plot (sortie graphique). | ascii | ||
| nodata | 204 | Code retour en cas d'absence de données | 404 |
Options de traitement des signaux temporels
Les paramètres suivants permettent de filtrer les signaux. L’ordre des paramètres compte puisque chaque opération sera effectuée dans l’ordre donné.
| Paramètre | Alias | Défaut | Unité | Description | Exemple | Documentation |
|---|---|---|---|---|---|---|
| demean | false | Retranche la valeur moyenne aux données. | true | ObsPy | ||
| detrend | false | Supprime la tendance temporelle du signal. | true | ObsPy | ||
| taper | Fenêtrage dans le domaine temporel. La taille de la fenêtre est donnée en fraction de la taille du signal (de 0 à 0.5). Elle peut être suivie du type de fenêtre : HANNING (défaut), HAMMING, COSINE. | 0.25 0.25, HANNING 0.35, COSINE |
ObsPy | |||
| decimate | deci | Facteur de décimation. Un filtre passe-bas est appliqué pour empêcher l’effet de repliement de spectre (antialiasing). | 2.0 | ObsPy | ||
| correct | false | Applique une correction instrumentale par déconvolution avec conversion en unités géophysiques. | true | ObsPy | ||
| waterlevel | 10 | dB | Niveau de l’eau utilisé pour la déconvolution. La valeur "none" correspond au filtrage inverse classique. |
1 none |
ObsPy | |
| freqlimits | Filtrage passe-bande appliqué avant la déconvolution dans le domaine fréquentiel1,2. Fenêtre rectangle d’amplitude 1 entre f2 et f3 et en cosinus jusqu’à 0 pour f1 < f < f2 et f3 < f < f4. | 0.01-0.04-0.5-0.6 | ObsPy | |||
| units | AUTO | Unité de sortie (AUTO, DISP, VEL, ACC). Actif uniquement avec l’option correct. | VEL | ObsPy | ||
| diff | false | Signal dérivé (approximation centrée). | true | ObsPy | ||
| int | false | Signal intégré (méthode des trapèzes). | true | ObsPy | ||
| lpfilter | lp | Hz | Filtre passe-bas (IIR d’ordre 4 / fenêtre de Butterworth)1. | 1.0 | ObsPy | |
| hpfilter | hp | Hz | Filtre passe-haut (IIR d’ordre 4 / fenêtre de Butterworth)1. | 2.0 | ObsPy | |
| bpfilter | bp | Hz | Filtre passe-bande (IIR d’ordre 4 / fenêtre de Butterworth)1,2. | 0.5-1.2 | ObsPy | |
| zerophase | false | Filtrage à phase nulle. Option à utiliser avec : lpfilter, hpfilter ou bpfilter. | true | ObsPy | ||
| envelope | false | Enveloppe du signal par transformée de Hilbert approchée. | true | ObsPy | ||
| scale | Multiplie les valeurs du signal par une constante. | 1.5 | ||||
| divscale | Divise les valeurs du signal par une constante. | 4.0 |
Notes :
- Les fréquences sont exprimées en Hertz.
- Les fréquences sont séparées par une virgule ou un tiret. (Exemple : 0.01,0.02)
Options graphiques
| Paramètre | Alias | Défaut | Description | Exemple |
|---|---|---|---|---|
| showtitle | true | Affiche ou non un titre montrant le canal et l’intervalle de dates. | false | |
| showscale | true | Affiche ou non une échelle sur le côté droit du graphique. | false | |
| monochrome | false | Image en niveaux de gris ou en couleur. | true | |
| width | 1200 | Largeur de l’image de sortie (pixels). De 400 à 2000. | 500 | |
| height | 400 | Hauteur de l’image de sortie (pixels). De 200 à 2000. | 400 |
Support
Le code source est distribué sous les termes de la licence GPLv3 et est disponible sur https://gricad-gitlab.univ-grenoble-alpes.fr/OSUG/RESIF/ws-timeseries
Pour toute question concernant l'accès aux données, la transmission d'un rapport de bug, les suggestions, les demandes d'aide, veuillez consulter notre service d'assistance ou envoyer un courriel à sismo-help@resif.fr