R: Matrix Seismic Record

matsquiggle {RSEIS} | R Documentation

## Matrix Seismic Record

### Description

Plot a matrix of time series as a var-squiggle display (filled in half traces)

### Usage

    matsquiggle(XMAT, dt1, dist = NULL, thick = 1,
     FLIP = FALSE, filcol='blue', tracecol="black", add=FALSE, PLOT=TRUE)

### Arguments

`XMAT` | matrix of traces
`dt1` | sample interval, s
`dist` | distance for each trace in the matrix
`thick` | thickness for each trace to be plotted
`FLIP` | logical, FALSE (default) plot horizontal, TRUE=plot vertical
`filcol` | color for shading
`tracecol` | color for trace
`add` | add traces to existing plot
`PLOT` | whether to create a new plotting region

### Details

see varsquiggle for more details

### Value

side effects.

### Author(s)

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

### See Also

varsquiggle, varsquig

### Examples

    data(GH)
    m <- match( GH$STNS,    GH$stafile$name)
    LATS <- GH$stafile$lat[m]
    LONS <- GH$stafile$lon[m]
    dees <- rdistaz( GH$pickfile$LOC$lat, GH$pickfile$LOC$lon, LATS, LONS)

    sel <- which(GH$COMPS=="V")
    sel <- sel[order(dees$dist[sel])]

    ###  plot normal way:
    ### swig(GH, sel=sel, WIN=c(5,10), SHOWONLY=TRUE)

    ###  plot with varsquiggle
    ### varsquiggle(GH, sel=sel, WIN=c(5,10))

    ex <- seq(from=0, by=GH$dt[sel[1]], length=length(GH$JSTR[[sel[1]]]))
    wx <- ex>=5 & ex<=10
    XMAT <- matrix(ncol=length(sel), nrow=length(which(wx)))

    for(i in 1:length(sel))
    {
    XMAT[,i] <-  GH$JSTR[[sel[i]]][wx]

    }

    matsquiggle(XMAT, GH$dt[sel[1]] , dist = dees$dist[sel] , thick = 1,
    FLIP = FALSE)

    axis(1)
    axis(2)
    title(xlab="Time, s", ylab="Distance, km")

[Package RSEIS version 3.7-4 Index]