RSEIS (version 2.1-6)

JGET.seis: Reads various seismic file formats

Description

This fuction calls c-routines to read in segy, sac.

Usage

JGET.seis(fnames, kind = 1, PLOT = FALSE)
getseisinfo(fnames, kind = 1)

Arguments

fnames

list of file names.

kind

an integer 1, 2, 3; 1 = segy, 2 = sac, 3 = AH.

PLOT

logical, for if the traces should be plotted.

Value

• List containing the seismic data and header information. Each trace consists of a list with:
• fn
  original file name
• sta
  station name
• comp
  compnent
• dt
  delta t in seconds
• DATTIM
  time list
• yr
  year
• jd
  julian day
• mo
  month
• dom
  day of month
• hr
  hour
• mi
  minute
• sec
  sec
• msec
  milliseconds
• dt
  delta t in seconds
• t1
  time start of trace
• t2
  time end of trace
• off
  off-set
• N
  number of points in trace
• units
  units
• amp
  vector of trace values

Details

JGET.seis extracts digital seismic data from binary files stored in the file system. The program uses C-code for I/O and passes data back to R. Currently SAC, SEGY formats are installed but it is easy to extend. AH format is available for LINUX systems, but there were problems compiling in WINDOWS and MACOS so this feature was removed. A filter for mseed format is currently being developed.

See Also

Mine.seis, Package:Rsac

Examples

Lname = list.files(path='/data/wadati/bourbon/LaurelCanNC/R005.01' , pattern="08.005.01.41.23.9024", full.names=TRUE)

S1  = JGET.seis(Lname, kind = 1, PLOT = FALSE)