RSEIS (version 4.2-4)

DISTxsec: Distance Cross section

Description

Plot time series vertically at specified distances. Produces a seismic cross section with correct spacing between traces.

Usage

DISTxsec(GH, dist, TIM.WIN = c(0, 3600), sel, trace.width = 10,
col = "black", text.col = "blue", text.font = 2, text.size = 0.8,
 add = FALSE, plot = TRUE)

Value

vector of x-y coordinates of the plot.

Arguments

GH

RSEIS seismic trace structure, output of prepSEIS used in swig

dist

distance for each station along x-axis

TIM.WIN

time window for cross section

sel

numeric, index of selected traces to plot.

trace.width

Width of each trace in plot. Should be in same units as x-axis

col

color for traces. If vector, each trace is plotted with assigned color.

text.col

color for text identifying each trace.

text.font

font for text identifying each trace.

text.size

size of text for identifying each trace.

add

logical, Whether to add traces, or just set up the figure

plot

logical, whether to plotthe traces.

Author

Jonathan M. Lees<[e-mail]>

Details

Distances should be a vector for each trace in the RSEIS list.

See Also

swig, prepSEIS

Examples

#### example using data in the RSEIS package
data(GH)
####  get the source location
lat.org  = GH$pickfile$LOC$lat
lon.org  = GH$pickfile$LOC$lon
####  get the station locations
g1 =GH$stafile
#### find the distance to each station
gd = rdistaz(lat.org, lon.org, g1$lat, g1$lon )

##### optional, filter the data
sel= which( GH$COMPS == 'V')
### filter traces
Fdef  <-  list(ON=TRUE, fl=1, fh=1, type="HP", proto="BU", RM=TRUE, zp=TRUE )
KF  <-  FILT.SEISN(GH, FILT=Fdef)

### match the stations in GH to the station distances
m1 = match(GH$STNS  , g1$name) 
dist.GH = gd$dist[m1]
TIM.WIN = range(GH$ex) 


#######  prepare plot, but do not add traces
A = DISTxsec(KF, dist.GH, TIM.WIN, sel, trace.width = 0.5 , add=FALSE,
plot=FALSE )
#####  add traces
B = DISTxsec(KF, dist.GH, TIM.WIN, sel, trace.width = 0.5 , add=TRUE,
plot=TRUE, col='black'  , text.col='red', text.size=1  )