segy2rseis: Read in multiple segy files.

In RSEIS: Seismic Time Series Analysis Tools

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)

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

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.

Value

List of seismic traces.

Note

The Endian-ness of the input files is set by the system that created them. If the read1segy or read1sac does not make sense, try a different endian or BIGLONG setting.

Author(s)

Jonathan M. Lees<[e-mail]>

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

RSEIS documentation built on Sept. 13, 2024, 1:09 a.m.