read_fdsn {eseis} | R Documentation |
Download and import seismic data from an FDSN service provider
Description
The function implements download and import of seismic data from FDSN data providers via the fdsnws-dataselect service (see https://www.fdsn.org/webservices/). It is basically a wrapper for the query approach.
Usage
read_fdsn(
start,
duration,
station,
network,
component = "BHZ",
url,
eseis = TRUE,
...
)
Arguments
start |
|
duration |
|
station |
|
network |
|
component |
|
url |
|
eseis |
|
... |
Additional query arguments sent to the FDSN data provider. See details for available argument names and conventions. |
Details
The FDSN (International Federation of Digital Seismograph Networks)
provides access to a large number of seismic stations worldwide. The data
are organised by network, station, component and further arguments. In
order to use the eseis function read_fdsn
, one must know at least
the former three criteria for the data of interest. A list of networks is
available here: https://www.fdsn.org/networks/. The function expects
the 2-digit network code, the 3- or 4-digit station code, a single seismic
component ID, and the URL to the data archive. Additional query arguments
can be added (and must be added to point at a single seismic trace to
download and import). A complete list of query arguments is available
here: https://www.fdsn.org/webservices/fdsnws-dataselect-1.1.pdf.
For each network listed there, one can find the URL that gives access to the data (if existing) under "Data Access". Note that the function only requires the first URL part, e.g., https://geofon.gfz-potsdam.de.
Value
An eseis
object or a list
with the time
($time
) and $signal
vectors as well as meta information.
Author(s)
Michael Dietze
Examples
## Not run:
## read and plot 10 min of data from Ecuador, specifying the component
s <- read_fdsn(start = "2020-05-16 22:42:00",
duration = 360,
station = "IMBA",
network = "EC",
component = "HHZ")
plot(s)
## read and plot 10 min of data from Germany, specifying the URL
s <- read_fdsn(start = "2017-03-21 04:38:00",
duration = 360,
station = "RGN",
network = "GE",
url = "http://geofon.gfz-potsdam.de")
plot(s)
## End(Not run)