aux_getIRISdata: Download seismic data from IRIS data base

In coffeemuggler/eseis: Environmental Seismology Toolbox

Download seismic data from IRIS data base

Description

This function accesses the IRIS internet data base of seismic signals and downloads seismic data based on the provided SNCL string and time information. The downloaded data is converted to the same structure as would be expected from read_sac or read_mseed.

Usage

aux_getIRISdata(
  start,
  duration,
  sncl,
  quality = "D",
  ID_iris = "IrisClient",
  eseis = TRUE
)

Arguments

start	POSIXct value, start time of the data to query.
duration	Numeric value, length of the data to query, in seconds.
sncl	Character vector, SNCL string used to identify station and component of interest. These strings should match the time criteria. Typically, the SNCL string can be taken from the output of the function aux_getirisstations.
quality	Character value, quality level of the data. One out of "D" (The state of quality control of the data is indeterminate), "R" (Raw Waveform Data with no Quality Control), "Q" (Quality Controlled Data, some processes have been applied to the data), "M" (Data center modified, time-series values have not been changed), "B". Default is "D".
ID_iris	Character value, IRIS ID. Default is "IrisClient".
eseis	Logical scalar, option to read data to an eseis object (recommended, see documentation of aux_initiateeseis), default is TRUE

Details

The function makes use of the package 'IRISSeismic'. It requires a working internet connection to perform the download.

Value

List with downloaded seismic data. For each element in sncl, a list element is created, which in turn contains a list with the typical seismic data organisation as, for example, created by read_sac.

Author(s)

Michael Dietze

Examples

## Not run: 

sncl <- aux_getIRISstation(start = as.POSIXct("2010-01-01 22:22:22", 
                                               tz = "UTC"), 
                            duration = 120, 
                            location = c(53, 13), 
                            radius = 0.7, 
                            component = "BHZ")

s <- aux_getIRISdata(start = as.POSIXct("2010-01-01 22:22:22", 
                                           tz = "UTC"), 
                        duration = 120,
                        sncl = sncl$sncl[1])
                        
plot_signal(data = s[[1]])

## End(Not run)
                     

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