eseis (version 0.8.0)

aux_getFDSNstation: 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)

Value

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

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.

Author

Michael Dietze

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".

See Also

aux_get_FDSNdata, aux_getIRISstation

Examples

if (FALSE) {

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)
}