getseis24: Get 24 Hours of Seismic Data

In RSEIS: Seismic Time Series Analysis Tools

Get 24 Hours of Seismic Data

Description

Get 24 Hours of Seismic Data

Usage

getseis24(DB, iyear = 2009, iday = 1, usta = "",
 acomp = "", kind = 1,  Iendian=1, BIGLONG=FALSE)

Arguments

DB	Data base of meta-data about the seismic trace files
iyear	Year for extraction
iday	Julian day for extraction
usta	station to show
acomp	component to show
kind	kind of data, default=1, 0="RDATA" , -1="RDS", 0="RDATA", 1 = "segy", 2 = "sac"
Iendian	Endian-ness of the data: 1,2,3: "little", "big", "swap". Default = 1 (little)
BIGLONG	logical, TRUE=long=8 bytes

Details

The DB file consists of a list of information on where to find the data and what times are covered. DB is

fn

full path to file

yr

year

jd

julian day

hr

hour

mi

minute

sec

second

dur

duration, seconds

origyr

origin time for epoch calculations

Value

yr	start year
jd	start julian day
t1	start t1 (with epoch)
t2	start t2 (with epoch day)
ed	epoch day
hr	start hour
mi	start minute
sec	start seconds
gamp	Amplitude of each trace
gdt	delta-t, sample interval, in seconds
gnam	station name
gfile	file information
sigs	List of time series
zna	List of NA values in each time series

Author(s)

Jonathan M. Lees<[e-mail]>

See Also

setupDB, plotseis24

Examples

data(KH)

amp = KH$JSTR[[1]]
OLDdt = KH$dt[1]
newdt = 0.1
yr = 2000
GIVE = FAKEDATA(amp, OLDdt=0.01, newdt = 0.1, yr = 2000,
         JD = 4, mi = 12, sec = 0,  Ntraces = 24*3,
seed=200, noise.est=c(1, 100) , verbose=TRUE  )
####   each trace in a separate file
tdir = tempdir()
for(i in 1:length(GIVE) )
{
    sig = GIVE[[i]]
  d1 =  dateStamp(sig$DATTIM,  sep='_')
   nam1 =  paste(d1,sig$sta, sig$comp, sep='_')
    nam2 = paste0(nam1, '.RDS')
    nam3 = paste(tdir, nam2, sep='/')
    saveRDS(file=nam3, sig) 
    }

########################  Now  read files and make the DataBase:
LF  = list.files(path=tdir, pattern='.RDS', full.names=TRUE)

DB = FmakeDB(LF, kind=-1)

IDB = infoDB(DB)


START =  list(yr =yr , jd= 5 , hr= 0 , mi= 0 ,sec= 0)

END = list(yr =yr , jd= 7 , hr= 0 , mi= 0 ,sec= 0)   


h = getseis24(DB, iyear = 2000, iday = 5, usta = IDB$usta,
                     acomp = IDB$ucomp, kind = -1,  Iendian=1, BIGLONG=FALSE)

  pjj <- plotseis24(h, dy=1/18, FIX=24, SCALE=1,
     FILT=list(ON=FALSE, fl=0.05 , fh=20.0, type="BP", proto="BU"),
     RCOLS=c(rgb(0.2, .2, 1), rgb(.2, .2, .2))  )

RSEIS documentation built on Sept. 13, 2024, 1:09 a.m.