getregionals(Rquake)
getregionals()所属R语言包:Rquake
Extract regional events
提取区域活动
译者:生物统计家园网 机器人LoveR
描述----------Description----------
Extract regional events from a hypocenter list (catalog)
提取区域活动从一个震源列表(目录)
用法----------Usage----------
getregionals(KAT, Mlat, Mlon, rad = 1000, t1 = 1, t2 = 2)
参数----------Arguments----------
参数:KAT
catalog list
目录列表
参数:Mlat
central latitude
中心纬度
参数:Mlon
central longitude
中央经度
参数:rad
radius (km)
半径(公里)
参数:t1
start time (julian days)
开始时间(儒略日)
参数:t2
end time (julian days)
结束时间(儒略日)
Details
详细信息----------Details----------
given an earthquake catalog from PDEs, for example, extract the events that are close to a network in a given time frame. The limited data set may be used to help predict arrival times for known hypocenter locations.
从偏微分方程的地震目录,例如,提取网络在给定的时间框架内的事件。有限的数据集,可以使用已知的震源位置,以帮助预测到达时间。
值----------Value----------
Catalog
目录
(作者)----------Author(s)----------
Jonathan M. Lees<jonathan.lees@unc.edu>
参见----------See Also----------
Mine.seis
Mine.seis
实例----------Examples----------
## Not run: [#不运行:]
fpde = '/Users/lees/Site/CHAC/pde_catalog_NOV_2011.txt'
ipde = prepPDE(fpde)
Mlat = median(stas$lat)
Mlon = median(stas$lon)
Lt1 = JtimL(upf$LOC)
t1 = Lt1-3/(24*6)
t2 = Lt1+3/(24*6)
localeqs = getregionals(KAT, Mlat, Mlon, rad=100 , t1=NULL, t2=NULL)
for(i in 1:length(localeqs))
{
j = localeqs[i]
at1 = KAT$jsec[j]
at2 = at1+20/(24*60)
GH = Mine.seis(at1, at2, DBnov , gstas, gcomps, kind=2, Iendian=1, BIGLONG=FALSE)
hord = which(GH$COMPS=="V")
gret = swig(GH, sel=hord, PADDLAB=buts)
#### save any picks to a file on disc for later use[###保存到一个文件中的任何选秀权在光盘上以备后用]
}
## End(Not run)[#(不执行)]
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注:
注1:为了方便大家学习,本文档为生物统计家园网机器人LoveR翻译而成,仅供个人R语言学习参考使用,生物统计家园保留版权。
注2:由于是机器人自动翻译,难免有不准确之处,使用时仔细对照中、英文内容进行反复理解,可以帮助R语言的学习。
注3:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
|