convert2Rseis: Convert Seismic data to RSEIS

In RSEIS: Seismic Time Series Analysis Tools

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

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.

Value

Side effects - creates new files on local system

Note

JGET.seis extracts digital seismic data from binary files stored in the file system. The program uses readBin for I/O and passes data back to R. Currently SAC, SEGY formats are installed but it is easy to extend. AH format is available for LINUX systems, but there were problems compiling in WINDOWS and MACOS so this feature was removed. A filter for mseed format is currently being developed.

Author(s)

Jonathan M. Lees<[e-mail]>

See Also

JGET.seis, JSAC.seis , Mine.seis

Examples

## Not run: 

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







## End(Not run)

RSEIS documentation built on Aug. 19, 2023, 5:07 p.m.