getGAP(Rquake)
getGAP()所属R语言包:Rquake
Get Seismic Gap
地震空区
译者:生物统计家园网 机器人LoveR
描述----------Description----------
Given an earthquake and a set of stations, return the maximum angle subtended between adjacent stations relative to the epicenter.
给定一个地震和一组工作站,返回相对震中相邻站之间的最大角。
用法----------Usage----------
getGAP(EQ, Ldat, PLOT = FALSE)
参数----------Arguments----------
参数:EQ
List, Earthequake location, elements (lat, lon) must be present
列表,Earthequake的位置,元素(纬度,经度)必须存在
参数:Ldat
List, station information, (lat, lon) must be present
列表,站内信息(纬度,经度)必须存在
参数:PLOT
logical, plot the stations and show the gap
逻辑,绘制站和显示的差距
Details
详细信息----------Details----------
Theangles are calculated in cartesian coordinates with the epicenter at the origin using a UTM projection.
Theangles计算在直角坐标系中的震中在使用UTM投影的起源。
值----------Value----------
numeric, gap in degrees
数字差距度
(作者)----------Author(s)----------
Jonathan M. Lees<jonathan.lees@unc.edu>
参见----------See Also----------
eqwrapup
eqwrapup
实例----------Examples----------
require(GEOmap)
set.seed(0)
N = 10
snames = paste(sep="", "A", as.character(1:N))
stas = list(name=snames, lat=runif(N, 35.9823, 36.1414), lon=runif(N, -118.0031, -117.6213))
NEQ = 3
WEQ = list(lat=runif(NEQ, 35.9823, 36.1414), lon=runif(NEQ, -118.0031, -117.6213))
MLAT = median(stas$lat)
MLON = median(stas$lon)
proj = setPROJ(type=2, LAT0=MLAT, LON0=MLON)
XYSTAS = GLOB.XY(stas$lat, stas$lon , proj)
eqxy = GLOB.XY(WEQ$lat, WEQ$lon, proj)
plot(range(c(XYSTAS$x, eqxy$x)), range(c(XYSTAS$y, eqxy$y)), type='n', asp=1, xlab="km", ylab="km" )
points(XYSTAS$x, XYSTAS$y, pch=6)
for(i in 1:NEQ)
{
EQ = list(lat=WEQ$lat[i], lon=WEQ$lon[i])
g = getGAP(EQ, stas, PLOT=FALSE)
points(eqxy$x[i], eqxy$y[i], pch=8, col='red')
text(eqxy$x[i], eqxy$y[i], labels=paste("gap=", format(g)), pos=3)
}
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注:
注1:为了方便大家学习,本文档为生物统计家园网机器人LoveR翻译而成,仅供个人R语言学习参考使用,生物统计家园保留版权。
注2:由于是机器人自动翻译,难免有不准确之处,使用时仔细对照中、英文内容进行反复理解,可以帮助R语言的学习。
注3:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
|