R: Veiw seismic data window

view.seis {RSEIS} | R Documentation

## Veiw seismic data window

### Description

Veiw seismic data (segy) window on an hourly basis.

### Usage

    view.seis(aday, ihour, inkhour, SAVEFILE, days,
    DB, usta, acomp,  STDLAB =c("QUIT",  "NEXT", "PREV", "HALF"),
    TZ=NULL)

### Arguments

`aday` | index of which day to use in vector days
`ihour` | hour to start
`inkhour` | increment in hours for viewing panel
`SAVEFILE` | file to save window picks in
`days` | vector of days to select from
`DB` | data base list of file names and start-times and durations
`usta` | stations to select
`acomp` | compnents to select
`STDLAB` | vector of buttons, DEFAULT = c("QUIT", "NEXT", "PREV", "HALF", "WPIX", "zoom out", "refresh", "restore", "SPEC", "SGRAM" ,"WLET", "FILT", "Pinfo", "WINFO")
`TZ` | Number of hours to add to GMT to get local time

### Details

The program view.seis assumes the data is stored in files accessable by the user and that the DB list has been scanned in and parsed.

### Value

Graphical side effects and save.wpix stores appended picks.

### Note

On LINUX systems I wrote these (non-R) programs to set up the data base for segy data:FLS.prl, segydatabase. To get these contact me directly. TZ is (-6) for Guatemala.

### Author(s)

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

### See Also

swig, save.wpix

### Examples

    ## Not run:

    ##############  method one:
    #######  read in data base from files:
    DB  <- scan(file="/home/lees/Site/Santiaguito/SG09/segyDB", list(fn="",
    yr=0, jd=0, hr=0, mi=0, sec=0, dur=0))

    ##  note: the file name includes the station SMI and the component J

    ##########  extract names of stations and components
    das1  <- substr(DB$fn, 58, 60)
    comp1 <- substr(DB$fn, 62,62))

    ###########  set up times (start and end) in julian days for selection
     DB$t1 <- DB$jd + DB$hr/24 + DB$mi/(24*60) + DB$sec/(24*3600)
    DB$t2 <- DB$t1 + DB$dur/(24*3600)

    #######   method 2:

    DB <- makeDB("/Users/lees/Site/Santiaguito/SG09", pattern = "R", kind = 1)
    udays <- unique(paste(DB$yr, DB$jd))
    sdays <-  as.numeric( unlist( strsplit(udays, split=" ") ) )

    ye <- sdays[seq(from=1, to=length(sdays), by=2)]
    d  <- sdays[seq(from=2, to=length(sdays), by=2)]
    o <- order(ye+d/366)

    days <- list(yr = ye[o], jd=d[o])

     das1  <- DB$sta
      comp1 <-  DB$comp

    #######   get unique stations from list:
    usta <- unique(das1)

    ###########   choose the desired components
    acomp <- c('I','V')

    ####   aday refers to one of the days listed in the days structure

    ###  days is a list of days (and associated years) that are in teh DB
    days <- list(jd=c(365,366, 1,2,3,4), yr=c(2008, 2008, rep(2009, times=4)))

    #### extract days from the DB list:
    udays <- unique(paste(DB$yr, DB$jd))
    sdays <-  as.numeric( unlist( strsplit(udays, split=" ") ) )
    ye <- sdays[seq(from=1, to=length(sdays), by=2)]
    d  <- sdays[seq(from=2, to=length(sdays), by=2)]
    o <- order(ye+d/366)

    days <- list(yr = ye[o], jd=d[o])

    ####

    aday <- 3
    SAVEFILE <- 'wpix.out'
    ihour <- 15
    inkhour <- .5

    ####   aday refers to one of the days listed in the days structure

    view.seis(aday, ihour, inkhour, SAVEFILE, days, DB, usta, acomp, TZ=(-6))

    ## End(Not run)

[Package RSEIS version 3.7-4 Index]