getseis24: Get 24 Hours of Seismic Data

getseis24R Documentation

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

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.