Plot Hodogram and show seismic particle motion
PMOT.drive(temp, dt, pmolabs = c("Vertical", "North", "East"), STAMP = "", baz = 0)
Graphical Side Effect.
matrix of 3-component seismic signal
sample interval (delta-T, seconds)
labels for traces
Character string Identification stamp
Back Azimuth, degrees
Jonathan M. Lees<jonathan.lees.edu>
Input matrix should V, N, E. Baz is not implemented yet.
data("GH")
sel <- which(GH$STNS == "CE1")
YMAT <- cbind(GH$JSTR[[sel[1]]][1168:1500],
GH$JSTR[[sel[2]]][1168:1500],
GH$JSTR[[sel[3]]][1168:1500])
dt <- GH$dt[ sel[1] ]
ftime <- Zdate(GH$info, sel[1], 1)
if (FALSE) {
PMOT.drive(YMAT, dt, pmolabs = c("Vertical", "North", "East"),
STAMP =ftime )
}
Run the code above in your browser using DataLab