RSEIS (version 4.2-4)

segy2rseis: Read in multiple segy files.

Description

Read in multiple segy files, and create a list of seismic traces.

Usage

segy2rseis(fnames, Iendian = 1, HEADONLY = FALSE, BIGLONG = FALSE, PLOT
= -1, RAW = FALSE)
sac2rseis(fnames, Iendian = 1, HEADONLY = FALSE,
BIGLONG = FALSE, PLOT = -1, RAW = FALSE)

Value

List of seismic traces.

Arguments

fnames

character vector of file names.

Iendian

Endian-ness of the files

HEADONLY

logical, TRUE=read only the header information. default=FALSE

BIGLONG

logical, indicating whether long is 8 or 4 bytes.

PLOT

logical, TRUE = plot traces

RAW

logical, TRUE=do not convert data to volts

Author

Jonathan M. Lees<[e-mail]>

Details

Segy format files are in integer format. The time series ususally represents counts recorded in a data acquisition system. The header includes meta-data and other identifying information.

See Also

read1sac, read1segy, sac2rseis, prepSEIS

Examples

#####  make some SAC files, then read them in
data(GH)
apath = tempdir()
##  setwd(apath)
## apath = 'TEMP'
J =  rseis2sac(GH, sel =1:5,  path = apath, BIGLONG =FALSE )
Iendian = .Platform$endian
####### next read them in
Lname  <-  list.files(path=J , pattern='SAC', full.names=TRUE)

H = sac2rseis(Lname , Iendian =Iendian , HEADONLY = FALSE,
BIGLONG = FALSE, PLOT = -1, RAW = FALSE)

####  should have 5 traces, look at elements of the first one:
names(H[[1]])

plotGH(H[[1]])