找回密码
 注册
查看: 293|回复: 0

R语言 Rquake包 getGAP()函数中文帮助文档(中英文对照)

[复制链接]
发表于 2012-9-28 20:26:34 | 显示全部楼层 |阅读模式
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:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

手机版|小黑屋|生物统计家园 网站价格

GMT+8, 2024-11-26 18:37 , Processed in 0.036190 second(s), 15 queries .

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表