JSAC.seis {RSEIS}R Documentation

JSAC.seis

Description

Read SEGY/SAC format binary data

Usage

JSAC.seis(fnames, Iendian = 1 , HEADONLY=FALSE, BIGLONG=TRUE, PLOT = FALSE)
JSEGY.seis(fnames, Iendian = 1 , HEADONLY=FALSE, BIGLONG=TRUE, PLOT = FALSE)

Arguments

fnames vector of file names to be extracted and converted.
Iendian Endian-ness of the data: 1,2,3: "little", "big", "swap". Default = 1 (little)
HEADONLY logical, TRUE= header information only
BIGLONG logical, TRUE=long=8 bytes
PLOT logical, whether to plot the data after reading in

Details

Uses readBin to extract data in SAC format. user must know what kind of machine the data was created on for I/O purposes.

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

Note

SAC created on PC (windows) or LINUX machines typically will be in little endian format. SAC created on a SUN will be in big format. If you want to swap endian-ness , choose swap.

MAC uses different convention.

There are other issues regardin the size of long. These will be dealt with in future implementations.

Author(s)

Jonathan M. Lees<jonathan.lees@unc.edu>

See Also

Mine.seis, Package:Rsac

Examples


## Not run: 

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

S1  = JSAC.seis(Lname, Iendian = 1, PLOT = FALSE)

## End(Not run)




[Package RSEIS version 2.3-3 Index]