JSAC.seis {RSEIS} | R Documentation |
Read SEGY/SAC format binary data
JSAC.seis(fnames, Iendian = 1 , HEADONLY=FALSE, BIGLONG=TRUE, PLOT = FALSE) JSEGY.seis(fnames, Iendian = 1 , HEADONLY=FALSE, BIGLONG=TRUE, PLOT = FALSE)
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 |
Uses readBin to extract data in SAC format. user must know what kind of machine the data was created on for I/O purposes.
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 |
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.
Jonathan M. Lees<jonathan.lees@unc.edu>
Mine.seis, Package:Rsac
## 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)