aux_splitcubechannels function

Create multiple single-component files from multi-component files

The function changes the meta data and file names of seismic data sets (organised in an eseis-compatible data structure). Data cubes in combination with splitter break-out boxes allow to record vertical component signals from three single-channel geophones, each connected by an individual cable. The logger interprets the data as E, N and Z components, regardless of the actual sensor and component connected to the respective cable. Hence, it is possible to record data from three independent (vertical) sensors with a single Data cube. However, the channels will be named according to the internal scheme (E, N, Z). To correct for that effect, the function will look up the right combination of station ID and component, and change meta data and name of each file in the input directory accordingly, before saving the updated file in the output directory.

aux_splitcubechannels(
  input,
  output,
  ID_in,
  component_in,
  ID_out,
  component_out,
  gipptools
)

Arguments

• input: Character value, path to directory where the Cube files to be processed are stored.
• output: Character value, path to directory where output data is written. It is not recommended to use the same directory for output and input, as files may be overwritten unintentionally.
• ID_in: Character vector with the IDs of the stations that will be changed. Note that for every channel of the input data, the corresponding ID must be repeated. See details for further information.
• component_in: Character vector with the components of the input stations that will be changed. The length must match the length of ID_in as well as ID_out and component_out.
• ID_out: Character vector with the IDs of the output files corresponding to the input stations and components.
• component_out: Character vector with the components of the output stations.
• gipptools: Character value, path to gipptools or cubetools directory. Only needed when mseed files are processed

Returns

A set of converted and organised seismic files written to disk.

Examples

## Not run:

## the example will convert the E, N and Z components of station RUEG1
## to the station IDs RUEGA, RUEGB, RUEGC, all with BHZ as component

aux_splitcubechannels(input =  paste0(system.file("extdata", 
                                      package="eseis"), "/"), 
                      output = tempdir(), 
                      ID_in = c("RUEG1", "RUEG1", "RUEG1"), 
                      component_in = c("BHE", "BHN", "BHZ"),
                      ID_out = c("RUEGA", "RUEGB", "RUEGC"),
                      component_out = c("BHZ", "BHZ", "BHZ"))
## End(Not run)

Author(s)

Michael Dietze