R: Long FFT Spectrogram

longfft {RSEIS} | R Documentation

## Long FFT Spectrogram

### Description

Creates hourly spectrograms, either alternating seismic and infrasound data or sequences of one component.

### Usage

    longfft(DB, DAYS = c(233, 234), HRS = 1:24, sta = "KR1", comp = c("V",
    "I"), NPP = 6, CSCALE = FALSE, pal = rainbow(100), PS = FALSE, kind = 1,
    Iendian = 1, BIGLONG = FALSE)

    longreset(NPP, PS)

    longpstart(NPP = 6, asta = "", acomp = "", theday = 1, hr = 0)

### Arguments

`DB` | RSEIS Data base
`DAYS` | vector of Days to display
`HRS` | vector of hours to display
`sta` | stations to extract
`comp` | component to extract
`NPP` | Number of plot strips per page, default = 6
`CSCALE` | scaling
`pal` | palettes to use (given two will alternate these)
`PS` | logical, TRUE postscript output
`kind` | kind of seismic data, 1=SEGY, 2=SAC
`Iendian` | Endian-ness of binary data
`BIGLONG` | logical, TRUE=long is 8 bytes
`asta` | character, one station
`acomp` | character, one component
`theday` | one day
`hr` | one hour

### Details

Extracts dats from the DB data base and plots strips of spectrograms for perusal.

longpstart, longreset are auxilliary codes used to set up the postscript files and initialize the plotting.

### Value

Graphical Side effects

### Note

Program is set for data being ready from external sources in binary (SAC, SEGY) format. If data is in R-format already, the code may not work.

### Author(s)

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

### See Also

SPECT.drive

### Examples

    ## Not run:
    ###  output directory
    setwd("/home/lees/MARIO_TUNG/LONGFFT")

    ###  2 color palettes
    p1 <- Gcols(plow=5, phi=0,  N=100, pal="topo.colors", mingray=0.8)
    p2 <- Gcols(plow=5, phi=0,  N=100, pal="rainbow", mingray=0.8)

    ###   one station, vertical component
    longfft(DB, DAYS=c(148:156) , HRS=1:24 ,
    sta="BBIL", comp=c("BHZ") , NPP=6 , CSCALE=FALSE,
    pal = list(p1=p1, p2=p2), PS=TRUE ,  kind = 2,
    Iendian=1, BIGLONG=FALSE )

    ###   one station, infrasound
    longfft(DB, DAYS=c(148:156) , HRS=1:24 , sta="BBIL",
    comp=c("BDF") , NPP=6 , CSCALE=FALSE,
    pal = list(p1=p2, p2=p1), PS=TRUE ,  kind = 2,
    Iendian=1, BIGLONG=FALSE )

    ###  one station, both vertical and infrasound, alternating
    longfft(DB, DAYS=c(148:156) , HRS=1:24 , sta="BBIL",
    comp=c("BHZ", "BDF") , NPP=6 , CSCALE=FALSE,
    pal = list(p1=p2, p2=p1), PS=TRUE ,  kind = 2,
    Iendian=1, BIGLONG=FALSE )

    ## End(Not run)

[Package RSEIS version 3.7-4 Index]