makeDB {RSEIS}R Documentation

Create a seismic Waveform Database

Description

Create a seismic Waveform Database

Usage

makeDB(path, pattern, kind = 1)

Arguments

path Path to directory where files and directories exist
pattern pattern for listing of files
kind kind of data: segy=1; sac=2

Details

The files are typically located in a directory structure created by programs like ref2segy, a PASSCAL program for downloading data in the field. Each file contains one seismogram, with a header. makeDB reads in all the headers and creates a list of meta-data for later use in RSEIS.

Value

list:

fn file name
yr year
jd julian day
hr hour
mi minute
sec second
dur duration, seconds
t1 time 1 in Epoch days
t2 time 2 in Epoch days
sta station name
comp component name
origyr origin year for epoch times

Note

Epoch times are used to accomodate problems where julian days cross year end boundaries, so that day 366 comes before day 1 of the next year.

Author(s)

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

See Also

setupDB, Mine.seis , getseis24, plotseis24, EPOCHday, PICK.GEN

Examples

## Not run: 

#############  set directory
path = '/home/lees/Site/Santiaguito/SG09'
pattern = "R0*"

###   get DB information
XDB  =  makeDB(path, pattern, kind =1)

##### select a station
usta = "CAL"
acomp = "V"

#####   extract 24 hours worht of data
JJ = getseis24(DB, 2009, 2, usta, acomp, kind = 1)

#####  plot 24 hours worth of data

pjj = plotseis24(JJ, dy=1/18, FIX=24, SCALE=0, 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))  )

##########   window a small portion on 24 hour display
w = winseis24(pjj)

###  open a new window
X11()
VNE = c("V", "N", "E")
###  IJK = c("I", "J", "K")

gsta=c("CAL", "KAM", "DOM",  "LAV")

ucomp=VNE

###  or:   ucomp=IJK

###   set epoch day and get the times
 eday = EPOCHday(w$yr, jd = w$jd, origyr = DB$origyr)

########  using epoch day, set times for beginning and end of window
        at1 = eday$jday + (w$hr[1])/24
        at2 = eday$jday + (w$hr[2])/24
########  extract data from the data base:
        GH = Mine.seis(at1, at2, DB, gsta, ucomp)

###  show data:
        PICK.GEN(GH)

## End(Not run)


[Package RSEIS version 2.1-4 Index]