Learn R Programming

RSEIS (version 3.3-3)

plotseis24: Plot 24 hours of seismic data

Description

Plot 24 hours of seismic data using output of getseis24.

Usage

plotseis24(JJ, dy = 1/18, FIX = 24, SCALE = 0, FILT = list(ON = FALSE,
fl = 0.05, fh = 20, type = "BP", proto = "BU"), RCOLS = c(rgb(0.2, 0.2,
1), rgb(0.2, 0.2, 0.2)), add=FALSE )

Arguments

JJ
output list of getseis24
dy
Delta-y in percentage of trace
FIX
Fix 24 hour plot. If FIX is less than 24, the plot will show only that number of hours.
SCALE
scale, 0=scale each trace, 1=scale window
FILT
filter data
RCOLS
colors
add
logical, if TRUE, add to existing plot (i.e. do not issue a plot command)

Value

  • list:
  • xx-axis
  • yy-axis
  • yryear
  • jdjulian day

Details

Plots full 24 hours of data. The list returned can be used by winseis24 to get picks and windows for zooming.

The FIX argument is currently not available.

See Also

getseis24, winseis24

Examples

Run this code
JJ <- getseis24(DB, 2008, 366, usta, acomp, kind = 1)

### show with no scaling
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(grey(0.8) , grey(0.8)), add=FALSE  )


pjj <- plotseis24(JJ, 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)), add=TRUE  )

Run the code above in your browser using DataLab