convert2Rseis: Convert Seismic data to RSEIS

Description Usage Arguments Details Value Note Author(s) See Also Examples

View source: R/convert2Rseis.R

Description

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

Usage

1
2
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<jonathan.lees@unc.edu>

See Also

JGET.seis, JSAC.seis , Mine.seis

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
## 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. 2, 2021, 5:07 p.m.