eqswath {GEOmap} | R Documentation |
Extract a set of eathquakes in swath along a cross sectional line
Description
Extract a set of eathquakes in swath along a cross sectional line
Usage
eqswath(x, y, z, L, width = 1, PROJ = NULL)
Arguments
x |
x-coordinates of earthquakes |
y |
y-coordinates of earthquakes |
z |
z-coordinates of earthquakes |
L |
list of x-y coordinates of cross section |
width |
width of swath (km) |
PROJ |
projection information |
Details
All units should be the same.
Value
r |
r-distance along cross section (x-coordinate) |
dh |
distance from cross seection |
depth |
depth in cross section (y-coordinate) |
flag |
index vector of which earthquakes fell in swath and depth range |
InvBox |
coordinates of swath for plotting on map |
Author(s)
Jonathan M. Lees<jonathan.lees@unc.edu>
See Also
XSECwin, XSECEQ
Examples
############# create data
x = runif(100, 1, 100)
y = runif(100, 1, 100)
z = runif(100, 1, 10)
plot(x,y, asp=1)
## L = locator()
L=list()
L$x=c( 5.42328560757,64.62879777806)
L$y=c(89.843266449785,-0.174423911329)
J = eqswath(x, y, z, L, width = 10, PROJ = NULL)
########## show box:
plot(x,y, asp=1)
lines(J$InvBox$x, J$InvBox$y)
############ show cross section with events plotted
plot(J$r, -J$depth)
[Package GEOmap version 2.5-11 Index]