convert2Rseis: Convert Seismic data to RSEIS

Description

Convert Seismic in SAC or SEGY format to RSEIS native format.

Usage

convert2Rseis(FLS, NEWDIR = ".", kind = 1, Iendian = "little", BIGLONG =
FALSE, NEWsta = "", NEWcomp = "")

Arguments

FLS

array of File names

NEWDIR

Destination directory path

kind

an integer 1, 2, 3; 0=R(DAT) , 1 = segy, 2 = sac, 3 = AH.

Iendian

Endian-ness of the data: 1,2,3: "little", "big", "swap". Default = 1 (little)

BIGLONG

logical, TRUE=long=8 bytes

NEWsta

character vector, stations associated with the vector of files

NEWcomp

character vector, component name associated with the vector of files

Value

Side effects - creates new files on local system

Details

Converts the data to R format so it can be loaded with the load command. After this conversion, files should be loaded in subsequent calls by using kind=0.

Examples

Run this code

############   example of different scenarios
###  make list of the files
fls = list.files(path="/home/lees/SAC", pat="sac")

###  convert them to R

####   was created on a little endian machine and is being read on a little endian
convert2Rseis(fls, NEWDIR = "/home/lees/MY_RDATA", kind = 2, Iendian = "little", BIGLONG = FALSE)


####   was created on a little endian machine and is being read on a big endian
convert2Rseis(fls, NEWDIR = "/home/lees/MY_RDATA", kind = 2, Iendian = "swap", BIGLONG = FALSE)

####   was created on a little endian machine and is being read on a big
####                       endian, with  the 64 bit machine used LONG=8
convert2Rseis(fls, NEWDIR = "/home/lees/MY_RDATA", kind = 2, Iendian = "swap", BIGLONG = TRUE)


F90dir = "/mnt/flash/FUEGO/2009/F900"
kfl = list.files(path=F90dir, full.names=TRUE)

bname = basename(kfl)

###  now we want to convert the station names

###  extract the STA name from the file names:
STA = unlist( lapply(  strsplit(split="\.", bname), getmem , 3) )
###   extract comp name - in this case needs further extraction
COMP1 = unlist( lapply(  strsplit(split="\.", bname), getmem , 4) )
COMP = substr(COMP1, 3, 3)

convert2Rseis(kfl, NEWDIR="/home/lees/Site/Fuego/2009", kind = 2, Iendian="swap" , BIGLONG=FALSE, NEWsta=STA, NEWcomp=COMP )