viewCHAC: View Continuous Data

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

View source: R/viewCHAC.R

Description

Scroll through continuous data recorded in the field. Uses a database describing the locations and content of each file stored on disk.

Usage

1
2
viewCHAC(DBnov , gstas, gcomps,sched, stas, buts='GPIX', preFILT=list()
, replot=TRUE , kind=2, Iendian=1, BIGLONG=FALSE)

Arguments

DBnov

RSEIS Data Base (output of makeDB)

gstas

stations to extract

gcomps

components to extract

sched

schedule of start times for extraction

stas

station list

buts

buttons for swig

preFILT

Pre-Filter traces before plotting.

replot

logical, TRUE=rerun swig after done click

kind

kind of data, 0=nativeR, 1=segy, 2=sac

Iendian

endian

BIGLONG

big long or short long

Details

These are set up for the CHAC dataset.

Value

Graphics, and Side effects

Note

The preFILT argument is the standard way of assigning filters in RSEIS. For example, preFILT = list(ON=TRUE, fl=5 , fh=15, type="BP", proto="BU") will bandpass filter the traces between 5 and 15 Hz. If the logical ON flag of the filter is turned to FALSE, the filter is not applied.

Author(s)

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

See Also

makeDB, 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
## Not run: 

##########  set up data base:
fpath = "/home/lees/Site/CHAC/DATA"
fpat = "201111"
DBnov = makeDB(fpath, fpat, kind=2, Iendian=1, BIGLONG=FALSE)
###   get information:
IDB = infoDB(DBnov)

#####  select stations and components:
gstas = IDB$usta[-which(IDB$usta=="CHAC5")]
gcomps = IDB$ucomp[1:3]

####  extra buttons
buts = c("YPIX", "SPEC", "SGRAM", "WLET")
 fsta = "/home/lees/Site/CHAC/staLLZ.txt"
stas = scan(file=fsta,what=list(name="", lat=0, lon=0, z=0))
stas$z = stas$z/1000

###  set schedule
sched  =seq(from=325, to=335, by=1/24)

##########  open 2 windows
X11()
X11()

###  set main window to dev 2
dev.set(2)

###  set pre-filter (needs to be ON=TRUE to work

preFILT = list(ON=TRUE, fl=1/2 , fh=8, type="BP", proto="BU")

viewCHAC( DBnov , gstas, gcomps , sched, stas, buts =buts,
       preFILT = preFILT,kind = 2, Iendian = 1, BIGLONG = FALSE )



## End(Not run)

Rquake documentation built on Dec. 16, 2020, 5:06 p.m.