Last data update: 2014.03.03

R: Get azimuthal particle motion
GAZIR Documentation

Get azimuthal particle motion

Description

Do particle motion analysis

Usage

GAZI(ADAT, dt = 1, ex = seq(0, 100),
comp = c(4, 5, 6), sta = "ZZZ", az = 0,
len = 50, shift = 10, prev = 1, fileid = "", picks = NA, labs = NA)

Arguments

ADAT

Matrix of 3 component seismic data

dt

delta T (s)

ex

x-axis vector

comp

component names

sta

station name

az

azimuth of station orientation

len

length of time series

shift

amount to shift per window

prev

length of buffer at beginning of trace

fileid

character string to put on plot

picks

arrival times for annotation

labs

labels for arrival times for annotation

Value

list(aex=aex[1:jall], rateig=rateig[1:jall], aaz=aaz[1:jall], ai=ai[1:jall], figaz=figaz, azpar=azpar, incpar=incpar )

Examples


data("GH")


temp  <-  cbind(GH$JSTR[[4]], GH$JSTR[[5]], GH$JSTR[[6]])

pmolabs <- c("Vertical", "North", "East")


 G   <-  GAZI(temp, dt =GH$dt[4]  , comp = pmolabs, sta = GH$STNS[4] ,
az = 0, len =75, shift = 10, prev = 1)

Results


R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(RSEIS)
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/RSEIS/GAZI.Rd_%03d_medium.png", width=480, height=480)
> ### Name: GAZI
> ### Title: Get azimuthal particle motion
> ### Aliases: GAZI
> ### Keywords: hplot
> 
> ### ** Examples
> 
> 
> data("GH")
> 
> 
> temp  <-  cbind(GH$JSTR[[4]], GH$JSTR[[5]], GH$JSTR[[6]])
> 
> pmolabs <- c("Vertical", "North", "East")
> 
> 
>  G   <-  GAZI(temp, dt =GH$dt[4]  , comp = pmolabs, sta = GH$STNS[4] ,
+ az = 0, len =75, shift = 10, prev = 1)
> 
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>