prepSEIS(RSEIS)
prepSEIS()所属R语言包:RSEIS
Prepare structure for RSEIS
准备为RSEIS结构
译者:生物统计家园网 机器人LoveR
描述----------Description----------
Takes list of traces and prepares new list for analysis in RSEIS
轨迹列表,并准备新的列表在RSEIS分析
用法----------Usage----------
prepSEIS(GG)
参数----------Arguments----------
参数:GG
Output list of Rsac function GET.seis
输出列表RSAC功能GET.seis
Details
详细信息----------Details----------
prepSEIS is offered to reformat the output of a list of seismic traces (or other time series) for inpout to program swig()
prepSEIS提供,重新格式化地震迹线(或其他的时间序列)的列表的输出inpout方案痛饮()
值----------Value----------
RSEIS list
RSEIS列表
(作者)----------Author(s)----------
Jonathan M. Lees<jonathan.lees.edu>
参见----------See Also----------
swig, JGET.seis, GET.seis(package="Rsac"), Package:Rsac
痛饮,JGET.seis,GET.seis(包=“RSAC”的),包装:RSAC
实例----------Examples----------
data(sunspots)
ES = prep1wig(wig=sunspots, dt=1/12, sta="STA", comp="CMP", units="UNITS" )
EH=prepSEIS(ES)
STDLAB = c("DONE", "zoom out", "refresh", "restore", "XTR", "SPEC", "SGRAM" ,"WLET")
xx = swig( EH, STDLAB = STDLAB)
###################################[##################################]
###################################[##################################]
#################### example with multiple signals[###################例如与多个信号]
dt = 0.01
t = seq(0, 60, by=0.01)
###### % sample rate = 100 Hz, 0.001 seconds 601 samples[#####%采样率= 100赫兹,0.001秒601样品]
### set up the fequencies and amplitudes for signals that have 2 frequencies[##成立fequencies和幅度的信号,有2个频率]
afreqs1 = c(8, 13,12, 5 )
amps1 = c(6, 2,3, 2 )
####[###]
afreqs2 = c(4,3,2, 3 )
amps2 = c(10,5, 9, 2 )
x = cbind( amps1[1]*sin(2*pi*afreqs1[1]*t) + amps2[1]* sin(2*pi*afreqs2[1]*t),
amps1[2]*sin(2*pi*afreqs1[2]*t) + amps2[2]* sin(2*pi*afreqs2[2]*t),
amps1[3]*sin(2*pi*afreqs1[3]*t) + amps2[3]* sin(2*pi*afreqs2[3]*t),
amps1[4]*sin(2*pi*afreqs1[4]*t) + amps2[4]* sin(2*pi*afreqs2[4]*t))
d = dim(x)
######## names of signals[#######的信号名称]
mysta=c("R1", "R2", "R3", "R4")
MYLIST = list()
############ set up the initial list of wiggles[###########设置的初始列表晃]
for(i in 1:d[2])
{
A = prep1wig(wig =x[,i], sta=mysta[i], dt=dt, comp="DO", units="amp")
A[[1]]$DATTIM$yr=2000
MYLIST = c(MYLIST, A)
}
### reorganize into RSEIS format:[##重组到RSEIS格式的:]
PH1=prepSEIS(MYLIST)
STDLAB = c("DONE", "zoom out", "refresh", "restore", "XTR", "SPEC", "SGRAM" ,"WLET")
swig(PH1, STDLAB = STDLAB)
## Not run: [#不运行:]
###### this example shows how to read data from SEGY files using Rsac[#####这个例子显示了如何读取数据从SEGY文件RSAC]
##### and entering RSEIS[####和进入RSEIS的]
KG4 = JGET.seis(c(fnames), kind = 1, PLOT = -1)
KH=prepSEIS(KG4)
STDLAB = c("DONE", "zoom out", "refresh", "restore", "XTR", "SPEC", "SGRAM" ,"WLET")
xx = swig( KH, sel=which(KH$COMPS == "V"), STDLAB = STDLAB)
## End(Not run)[#(不执行)]
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注:
注1:为了方便大家学习,本文档为生物统计家园网机器人LoveR翻译而成,仅供个人R语言学习参考使用,生物统计家园保留版权。
注2:由于是机器人自动翻译,难免有不准确之处,使用时仔细对照中、英文内容进行反复理解,可以帮助R语言的学习。
注3:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注:
注1:为了方便大家学习,本文档为生物统计家园网机器人LoveR翻译而成,仅供个人R语言学习参考使用,生物统计家园保留版权。
注2:由于是机器人自动翻译,难免有不准确之处,使用时仔细对照中、英文内容进行反复理解,可以帮助R语言的学习。
注3:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
|