Vlocate(Rquake)
Vlocate()所属R语言包：Rquake
Hypocenter Determination
震源深度测定
译者：生物统计家园网 机器人LoveR
描述----------Description----------
Hypocenter Determination with error checking and adjustments.
震源深度测定错误检查和调整。
用法----------Usage----------
Vlocate(Ldat,EQ,vel,
distwt = 10,
lambdareg =100,
REG = TRUE,
WTS = TRUE,
STOPPING = TRUE,
tolx = 0.1,
toly = 0.1,
tolz = 0.5,
RESMAX = c(.4,.5),
maxITER = c(7, 5, 7, 4),
PLOT=FALSE)
参数----------Arguments----------
参数：Ldat
list, must inlude: lat, lon ,err, sec, cor (see details)
列表，必须inlude：纬度，经度，错了，秒，肺心病（见详情）
参数：EQ
list, must inlude: lat,lon,z, t
列表，必须inlude：纬度，经度，Z，T
参数：vel
list, 1D velocity structure
列表中，一维速度结构
参数：distwt
distance weighting factor
距离权重系数
参数：lambdareg
regularization parameter for damping
正则化参数的阻尼
参数：REG
logical, TRUE=use regularization
逻辑，TRUE =使用正规化
参数：WTS
logical, TRUE==use weighting
逻辑，TRUE ==使用权
参数：STOPPING
logical, TRUE=use stopping criteria
逻辑，TRUE =使用停止准则
参数：tolx
numeric, tolerance in km in x direction
数字，宽容公里，在x方向
参数：toly
numeric, tolerance in km in y direction
数字，宽容公里，在y方向
参数：tolz
numeric, tolerance in km in z direction
数字，宽容公里，在z方向
参数：RESMAX
vector, residual max for P and S, default=c(4,5)
向量，剩余的最大P波和S，默认目录为C（4,5）
参数：maxITER
vector, Maximum number of iterations for each section of the location routine, default=c(7,5,7,4)
矢量，最大迭代次数为各部分的定位程序，默认= C（7,5,7,4）
参数：PLOT
logical, plot results during iterations
在迭代过程中的逻辑，图结果
Details
详细信息----------Details----------
This is a wrapper for XYlocate, only here the lat-lon of the stations is passed and the code does the projection internally.
这是一个包装为XYlocate，只有在这里站传递纬度，经度和代码内部的投影。
There are 3 main loops, each controled by differing input params: first event is located only in XY keeping the depth fixed (7 iterations). Then an initial free solution is estimated using robust elimination of residual based on RESMAX (5 iterations). Finally a set of 7 iterations is applied providing the final estimate, along with error bars, elliposids, etc.
有3个主要的循环，每个受控于不同的输入PARAMS：只有在第一个事件的深度保持固定的XY（7次迭代）。一个初步的估计免费的解决方案，使用可靠的消除残余的基础上RESMAX（5次迭代）。最后一组的7次迭代应用提供最终的估计，随着错误的条形，elliposids，等
In the event no good solution is derived, the regularization parameter is doubled and a loop with 4 iterations is applied, and the result returned.
在衍生事件没有很好的解决方案中，正则化参数的一倍，并施加一个循环中有4个迭代，并返回结果。
值----------Value----------
list:
列表：
参数：EQ
Hypocenter lcoation
震源lcoation
参数：ERR
Error Analysis
错误分析
参数：its
number of iteration
迭代次数
参数：Ksolutions
list of matrices, each with intermediate x,y,z,t locations
的矩阵，每一个与中间的x，Y，Z，吨位置列表
注意----------Note----------
The schedule may be adjusted by duplicating this function and changing the maxit parameters.
该时间表可通过复制此功能和改变麦克斯特参数的调整。
（作者）----------Author(s)----------
Jonathan M. Lees<jonathan.lees@unc.edu>
参考文献----------References----------
参见----------See Also----------
XYlocate, Klocate, DoRLocate
XYlocate，Klocate，DoRLocate
实例----------Examples----------
data(GH)
g1 = GH$pickfile
data(VELMOD1D)
vel= VELMOD1D
w1 = which(!is.na(g1$STAS$lat))
sec = g1$STAS$sec[w1]
N = length(sec)
Ldat = list(
name = g1$STAS$name[w1],
sec = g1$STAS$sec[w1],
phase = g1$STAS$phase[w1],
lat=g1$STAS$lat[w1],
lon = g1$STAS$lon[w1],
z = g1$STAS$z[w1],
err= g1$STAS$err[w1],
yr = rep(g1$LOC$yr , times=N),
jd = rep(g1$LOC$jd, times=N),
mo = rep(g1$LOC$mo, times=N),
dom = rep(g1$LOC$dom, times=N),
hr =rep( g1$LOC$hr, times=N),
mi = rep(g1$LOC$mi, times=N) )
wstart = which.min(Ldat$sec)
EQ = list(lat=Ldat$lat[wstart], lon=Ldat$lon[wstart], z=6, t=Ldat$sec[wstart] )
AQ = Vlocate(Ldat,EQ,vel,
distwt = 10,
lambdareg =100 ,
REG = TRUE,
WTS = TRUE,
STOPPING = TRUE,
tolx = 0.01,
toly = 0.01 ,
tolz = 0.05, maxITER = c(7,5,7,4) , RESMAX = c(0.1, 0.1), PLOT=FALSE)
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注：
注1：为了方便大家学习，本文档为生物统计家园网机器人LoveR翻译而成，仅供个人R语言学习参考使用，生物统计家园保留版权。
注2：由于是机器人自动翻译，难免有不准确之处，使用时仔细对照中、英文内容进行反复理解，可以帮助R语言的学习。
注3：如遇到不准确之处，请在本贴的后面进行回帖，我们会逐渐进行修订。
|
发表于 2012-9-28 20:30:42
