Wavesdownloader User's Page
Nera Project « Back to table of contents
NERA: Network of European Research Infrastructures for Earthquake Risk Assessment and Mitigation
NERA (2010-2014) is an EC infrastructure project that integrates key research infrastructures in Europe for monitoring earthquakes and assessing their hazard and risk. NERA integrates and facilitates the use of these infrastructures and access to data for research, provides services and access to earthquake data and parameters, and hazard and risk products and tools. NERA activities are coordinated with other relevant EC-projects and European initiatives, and contribute to the ESFRI EPOS infrastructure and the OECD GEM program. Grants providing access to highly specialised research facilities are also available within NERA.
JRA2/WP12 - JRA2: Tools for real time seismology, acquisition and mining « Back to table of contents
Off-line tool-box for data download
JRA2 aim is to interface existing research analysis software with the data to the end of improving seismic real-time monitoring. This allows the testing of the software in a semi-productive mode without having to insert it within existing specialised seismic monitoring software packages such as EarthWorm or SeisComP.
WavesDownloader « Back to table of contents
Wavesdownloader is a python script package. The current version performs waveform and metadata request to IRIS and EIDA datacenters using a unique command line request. Addittionaly wavesdownloader select data without gaps, remove trend and mean; perform filtering, instrument removal, plotting and compute PGMs parameters suitable for ShakeMap.
Download here wavesdownloader and user examples.
Package requirements « Back to table of contents
- Obspy
- numpy
- scipy
- rdseed-5.0 or higher
- Eida registration: Fill form for automatic registration. This step is required in order to access to EIDA services
Note for Mac OS users: numpy and scipy are already included into dmg installing file
Usage « Back to table of contents
>>> wavesdownloader --help or wavesdownloader -h
>>> Usage: wavesdownloader.py [-h] [--beg BEG] [--end END] [--len LEN] [--usr USR]
[--pas PAS] [--sta STA] [--net NET] [--loc LOC]
[--cha CHA] [--rot ROT] [--res RES] [--mode MODE]
[--center CENTER] [--radius RADIUS]
[--supCor SUPCOR] [--infCor INFCOR]
[--format FORMAT] [--outdir OUTDIR]
[--server SERVER] [--fsfile FSFILE]
[--rmgaps RMGAPS] [--mingap MINGAP]
[--maxgap MAXGAP] [--reject REJECT] [--cfreq CFREQ]
[--demean DEMEAN] [--bandpass BANDPASS]
[--highpass HIGHPASS] [--lowpass LOWPASS]
[--wfiltr WFILTR] [--deci DECI] [--deco DECO]
[--flim FLIM] [--slta SLTA] [--wcf WCF] [--pgm PGM]
[--pgmfile PGMFILE] [--shake SHAKE] [--sa SA]
[--pltmode PLTMODE] [--pltchan PLTCHAN]
[--pltNERT PLTNERT] [--pltazi PLTAZI]
[--summary SUMMARY]
Options « Back to table of contents
| -h, --help | show this help message and exit |
| --beg BEG | Begin Time. !! No defaults!! Format YYYY-MM-DDThh:mm:ss (UTC Time) (e.g.: 2011-07-25T12:30:00) |
| --end END | End Time. !! No defaults!! F\ormat YYYY-MM-DDThh:mm:ss (UTC Time) (e.g.: 2011-07-25T12:35:00) |
| --len LEN | Length of signal in seconds. !! No defaults!! Mandatory option if --end not specified |
| --usr USR | User name for eida !! MANDATORY OPTION, No defaults!! |
| --pas PAS | Passwd for eida \!! MANDATORY OPTION, No defaults!! |
| --sta STA | Station list. Default=* |
| --net NET | Network list. Default=* |
| --loc LOC | Station information location ID. Default=* |
| --cha CHA | Channel list. Default=BHZ |
| --rot ROT | Rotate horizontal components from North-East -> Radial-Trasversal. Only with --mode "center". Default [Y]/N. |
| --res RES | Instrument response extraction: 0=none,1=RESP,2=PAZ(default),3=RESP/PAZ |
| --mode MODE | Area selection mode: [circular|rectangular]. Default=circular |
| --center CENTER | Lat Lon inner position for circular request. Default="41.9 12.5" |
| --radius RADIUS | Radius in DEGREE from innerPos to outerPos for --mode circular. Default="0 10" |
| --supCor SUPCOR | Max latitude and longitude for --mode rectangular. Default="60 60" |
| --infCor INFCOR | Min latitude and longitude for --mode rectangular. Default="10 10" |
| --format FORMAT | file format extraction storage [SAC,SACXY,GSE1,GSE2,SH_ASC,WAV]. Default=SAC |
| --outdir OUTDIR | directory for data extraction. Default=data |
| --server SERVER | servers [EIDA,IRIS,WEBDC,LOCAL]. LOCAL look for fseed files stored on your local machine. Only fseed files are allowed. --server "LOCAL needs --fsfile to be specified. LOCAL and external server are allowed within the same request. Default="EIDA IRIS WEBDC" |
| --fsfile FSFILE | fseed file name inclusive of path if path different than \. Default=None |
| --rmgaps RMGAPS | Remove traces with gaps. N/[Y] |
| --mingap MINGAP | Minimum gap allowed in seconds. Default=0 |
| --maxgap MAXGAP | Maximum gap allowed in seconds. Default=0 |
| --reject REJECT | Minimum length in percent for trace rejection. Default=100 |
| --cfreq CFREQ | Get dominant period [N]/Y. |
| --demean DEMEAN | Remove mean and trend. N/[Y] |
| --bandpass BANDPASS | Bandpass filter "corners fimn fmax". E.g.: 2 0.01 0.1. Default=None |
| --highpass HIGHPASS | Highpass filter "corners fimn". E.g.: 2 0.01. Default=None |
| --lowpass LOWPASS | Lowpass filter "corners fmax". E.g.: 2 0.1. Default=None |
| --wfiltr WFILTR | Write new filtered files into --wfiltr path. [N]/Y |
| --deci DECI | Decimation factor for sampling rate. Only integer decimation factor allowed. Default=None |
| --deco DECO | Deconvolution from instrument response [N]/Y. Requires --res=[2|3]. |
| --flim FLIM | Corner frequency for deconvolution filtering. Default="0.002 0.005 0.5 1" |
| --slta SLTA | Make Short-term/long-term average and trigs picks "STA LTA ON OFF". STA: Short-term average LTA: Long-term average ON: level trigger on OFF: level trigger off. E.g.: "0.5 5 7 1.5". Default="None" |
| --wcf WCF | Write characteristic function from slta into sac file format [N]/Y. Only binary sac file format allowed. No need to use --format option |
| --pgm PGM | Peaks Ground Motion parameters. max_displacement [m], max_velocity[m/s], max_acceleration[m/s^2]. Output into summary3.log file. Automatically enable --sa option (Spectral acceleration response). Default [N]/Y |
| --pgmfile PGMFILE | Name file for shakemap. User file-name MUST end with "_dat.xml" in order to be accepted by ShakeMap. Default="shakeList_dat.xml" |
| --shake SHAKE | Write PGMs for ShakeMap. [N]/Y | --sa SA | Modify defaults Spectral Acceleration Response parameters (only if –pgm="Y"): damping and cornerfrequencies in Hz. Example: --sa "0.1 1 10": damping factor (0.1) and corner frequencies in Hz. Spectral Acceleration Response in [m/s^2]. Defaults: sa="0.05 3.33 1.00 3.33". |
| --pltmode PLTMODE | Plot traces. 0=No plot,1=y_axe regular, 2=y_axe distance from epicenter. Plots saved in pdf format. Default=0. Plot names: plotWavesZ.pdf for vertical, plotWavesNS.pdf and plotWavesEW.pdf with --pltNERT "NE"; plotWavesT.pdf and plotWavesR.pdf with --pltNERT "RT" | --pltchan PLTCHAN | Channel to plot. Default=None |
| --pltNERT PLTNERT | Horizontal components [NE|RT]. Default="None" !! ONLY with --rot=Y |
| --pltazi PLTAZI | Plot traces within azimuth. Default="0 360" |
| --summary SUMMARY | Print data request summary on screen. [N]/Y. |
Output files « Back to table of contents
summary1.log, summary2.log and summary3.log files are automatically saved into the --outdir directory.- summary1.log
- summary2.log
- summary3.log
Includes a summary of all option used to run wavesdownloader
List downloaded stations and metadata. Download an explained file format example
List downloaded stations, metadata, PGMs, central frequency and picker values. Download an explained file format example