grotseis(RSEIS)
grotseis()所属R语言包：RSEIS
Get seismic rotation matrix
地震旋转矩阵
译者：生物统计家园网 机器人LoveR
描述----------Description----------
Set up a rotation matrix for a seismic trace. Rotation matrix is 3D, although this rotation only creates a rotation for conversion to radial-transverse orientation.
设置一个旋转矩阵的地震轨迹。旋转矩阵是三维的，虽然只旋转创建一个旋转转换为径向，横向方向。
用法----------Usage----------
grotseis(ang, flip = FALSE)
参数----------Arguments----------
参数：ang
Angle to rotate horizontal components, degrees
要旋转水平分量度
参数：flip
Logical, TRUE=flip the vertical axis, default=FALSE
逻辑，TRUE =翻转，垂直轴，默认= FALSE
Details
详细信息----------Details----------
Returns a 3 by 3 matrix used for rotationg a 3-component seismic record, usually stored as an N by 3 matrix. Only the N-E components are rotated, North component is assumed to be in position 2 and the east compnent is assumed to be in position 3. The vertical component can be flipped.
返回一个3×3的矩阵，用于回转时的3分量地震记录，通常被存储为一个N 3矩阵。只有网元组件转动，北分量被假定为在位置2和东部compnent被假定为在位置3。可翻转的垂直分量。
值----------Value----------
3 by 3 rotation matrix.
3×3旋转矩阵。
注意----------Note----------
Positive radial is away from the source (direction of wave propagation). Positive transverse is to the right when facing the direction of wave propagation.
正径向远离源（波传播方向）。正横时所面临的波传播方向的右侧。
（作者）----------Author(s)----------
Jonathan M. Lees<jonathan.lees@unc.edu>
参见----------See Also----------
rdistaz
rdistaz
实例----------Examples----------
STAXY=list()
STAXY$'x'=c(-2.9162198461534,-2.49599248511068,-2.85909405321704,-1.96135073099434,
-6.50413342506259,2.64026676599765,
-3.95701139503518,-2.84082134537436,-0.0457817300378462,-2.74214190991955)
STAXY$'y'=c(-7.83435541676815,-4.46180337254565,-6.46036190991833,-5.01212763828746,
-2.56091416028758,
5.31173503708142,2.10545324503380,-0.87490923667824,-0.172422188354707,-1.52055218789877)
STAXY$'lat'=c(14.685621984127,14.7159182222222,14.6979647030651,14.710975070028,
14.7329873333333,14.8037143111518
,14.7749104943935,14.7481391460905,14.7544511215933,14.7423394025875)
STAXY$'lon'=c(268.420918730159,268.424817925926,268.421447725096,268.429783940243,268.387586722222,
268.472531954619,268.41123843527,268.421611351166,268.447574716981,268.422528671994)
STAXY$'z'=c(0.92522857142857,1.48225333333333,1.14740517241379,1.4423781512605,1.51148,
2.53268681318681,2.70014678899083,2.04094444444444,2.90827547169811,2.31817123287671)
STAXY$'cen'=c(14.756,-91.552)
STAXY$name=c('OBS','CAR','MAR','CAS','MTB','STA','STE','MOT','SUM','DOM')
sguitoXY=list()
sguitoXY$'x'=c(-1.78551922571555)
sguitoXY$'y'=c(-1.80850340813817)
sguitoXY$'lat'=c(14.7397535236)
sguitoXY$'lon'=c(268.4314147874)
sguitoXY$'z'=c(2.501)
DAZ = rdistaz( sguitoXY$lat, sguitoXY$lon , STAXY$lat, STAXY$lon)
STAXY$az = DAZ$baz
plot(STAXY$x, STAXY$y, asp=1, xaxs="r" , yaxs="r" )
text(STAXY$x, STAXY$y,STAXY$name, pos=3)
points(0,0, pch=3)
points(sguitoXY$x,sguitoXY$y , pch=8)
segments(sguitoXY$x, sguitoXY$y, STAXY$x, STAXY$y, col="green", lty=2)
vecs = rbind(c(0,0,1), c(0,1,0))
for( i in 1:length(STAXY$x))
{
rbaz = grotseis(STAXY$az[i], flip=FALSE)
bvec = vecs
############## red is north, blue east[＃＃＃＃＃＃＃＃＃＃＃＃＃红，蓝东是北]
######## red is radial positive away from source, blue is transverse[径向正远离源＃＃＃＃＃＃＃红色，蓝色的是横向]
########## blue is positive to the right, if facing away from source[＃＃＃＃＃＃＃＃＃蓝色右侧为正，如果面向远离源]
## arrows(STAXY$x[i],STAXY$y[i], STAXY$x[i]+vecs[,2], STAXY$y[i]+vecs[,3], col=c("red", "blue"), lty=2)[＃箭头（STAXY $ x [i]的，STAXY美元的Y [i]，STAXY $ x [我] +血管内皮单元[2]，STAXY $ Y [我] +血管内皮单元[3]，列= C（“红” ，“蓝色”），LTY = 2）]
arrows(STAXY$x[i],STAXY$y[i], STAXY$x[i]+bvec[,2], STAXY$y[i]+bvec[,3],
col=c("red", "blue"), length=.08)
}
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注：
注1：为了方便大家学习，本文档为生物统计家园网机器人LoveR翻译而成，仅供个人R语言学习参考使用，生物统计家园保留版权。
注2：由于是机器人自动翻译，难免有不准确之处，使用时仔细对照中、英文内容进行反复理解，可以帮助R语言的学习。
注3：如遇到不准确之处，请在本贴的后面进行回帖，我们会逐渐进行修订。
|
发表于 2012-9-28 21:31:40
