aux_getFDSNstation: Query FDSN data base for stations

In coffeemuggler/eseis: Environmental Seismology Toolbox

Query FDSN data base for stations

Description

This function queries as series of data bases for seismic stations that match a set of criteria for seismic data. The criteria include signal time stamp and location, and component. The returned data can be used to download data using the function aux_FDSNdata.

Usage

aux_getFDSNstation(centre, radius, start, access, url)

Arguments

centre	Numeric vector of length two, center coordinates of the location to search data for (c(latitude, longitude)). Units must be decimal degrees.
radius	Numeric value, radius within which to search for seismic stations. Unit must be decimal degrees.
start	POSIXct value, start time of the data to query. If omitted, stations are queried for the full time available.
access	Logical value, access type of the data. If omitted, all data sets are returned, if set TRUE, only data with access flag "open" are returned.
url	Character vector, optional other FDSN base web addresses to search for stations. See details for default addresses and their format.

Details

The function requires a working internet connection to perform the query. It uses the following FDSN data bases by default:

• orfeus "http://www.orfeus-eu.org"

• geofon "http://geofon.gfz-potsdam.de/"

• bgr "http://eida.bgr.de"

• sss "http://eida.ethz.ch"

Other FDSN data base addresses can be provided in the same way as the addresses in the above list. They need to be provided as character vector. For a list of addresses see "http://www.fdsn.org/webservices/datacenters/" and "http://docs.obspy.org/packages/obspy.clients.fdsn.html#module-obspy.clients.fdsn".

Value

Data frame with query results. The data frame contains information for all seismic stations fulfilling the defined criteria.

Author(s)

Michael Dietze

See Also

aux_get_FDSNdata, aux_getIRISstation

Examples

## Not run: 

x <- aux_getFDSNstation(start = as.POSIXct(x = "2010-01-01 22:22:22", 
                                           tz = "UTC"), 
                        centre = c(45, 10), 
                        radius = 1)
                           
## optionally plot station locations on a map (requires RgoogleMaps)
center <- c(mean(x$station_latitude), 
            mean(x$station_longitude))

zoom <- min(RgoogleMaps::MaxZoom(range(x$station_latitude), 
                                 range(x$station_longitude)))
                                 
Map <- RgoogleMaps::GetMap(center = center,
                           zoom = zoom, 
                           maptype = "terrain")
                           
RgoogleMaps::PlotOnStaticMap(MyMap = Map, 
                             lat = x$station_latitude, 
                             lon = x$station_longitude, 
                             pch = 15, 
                             col = 4)

## End(Not run)
                     

coffeemuggler/eseis documentation built on Aug. 19, 2023, 9:57 p.m.