view.seis(RSEIS)
view.seis()所属R语言包:RSEIS
Veiw seismic data window
Veiw地震数据窗口
译者:生物统计家园网 机器人LoveR
描述----------Description----------
Veiw seismic data (segy) window on an hourly basis.
Veiw以小时为单位的地震数据(SEGY)窗口。
用法----------Usage----------
view.seis(aday, ihour, inkhour, SAVEFILE, days, DB, usta, acomp, STDLAB =c("QUIT", "NEXT", "PREV", "HALF"), TZ=NULL)
参数----------Arguments----------
参数:aday
index of which day to use in vector days
指数每天使用的矢量天
参数:ihour
hour to start
小时开始
参数:inkhour
increment in hours for viewing panel
查看面板,以小时为单位递增
参数:SAVEFILE
file to save window picks in
文件保存窗口选秀权
参数:days
vector of days to select from
天向量的选择
参数:DB
data base list of file names and start-times and durations
数据库列表中的文件名和启动时间和持续时间
参数:usta
stations to select
站选择
参数:acomp
compnents to select
compnents选择
参数:STDLAB
vector of buttons, DEFAULT = c("QUIT", "NEXT", "PREV", "HALF", "WPIX", "zoom out", "refresh", "restore", "SPEC", "SGRAM" ,"WLET", "FILT", "Pinfo", "WINFO")
向量的按钮,DEFAULT = C(“退出”,“下一页”,“上一页”,“半壁江山”,“WPIX”,“缩小”,“刷新”,“恢复”,“ SPEC“,”SGRAM“中,”WLET“,”FILT“中,”PINFO“,”WINFO“)
参数:TZ
Number of hours to add to GMT to get local time
小时数添加到GMT以获取本地时间
Details
详细信息----------Details----------
The program view.seis assumes the data is stored in files accessable by the user and that the DB list has been scanned in and parsed.
的的程序view.seis假定的数据被存储在文件accessable由用户和数据库列表已经被扫描和解析。
值----------Value----------
Graphical side effects and save.wpix stores appended picks.
图形副作用,save.wpix的店附加选秀权。
注意----------Note----------
On LINUX systems I wrote these (non-R) programs to set up the data base for segy data:FLS.prl, segydatabase. To get these contact me directly. TZ is (-6) for Guatemala.
在Linux系统上,我写这些(非R)程序来设置数据的基础上,为segy数据:FLS.prl,segydatabase。为了获得这些直接与我联系。 TZ(-6)危地马拉。
(作者)----------Author(s)----------
Jonathan M. Lees<jonathan.lees@unc.edu>
参见----------See Also----------
swig, save.wpix
痛饮,save.wpix
实例----------Examples----------
## Not run: [#不运行:]
############## method one:[#############方法1:]
####### read in data base from files:[######从文件中读取数据的基础上:]
DB = scan(file="/home/lees/Site/Santiaguito/SG09/segyDB", list(fn="",
yr=0, jd=0, hr=0, mi=0, sec=0, dur=0))
## this file has first line:[#这个文件的第一行:]
#/home/lees/Site/Santiaguito/SG09/R001.02/09.001.11.57.09.SMI.J 2009 1 11 57 9.245000 3600.000000[2009年1 / home/lees/Site/Santiaguito/SG09/R001.02/09.001.11.57.09.SMI.J 11 57 9.245000 3600.000000]
#c...|....1....|....2....|....3....|....4....|....5....|....6....|....7.$[C ... | .... 1 .... | .... .... | .... 3 .... | .... 4 .... | .... 5 .... | .... 6 .... | ... 7。]
## note: the file name includes the station SMI and the component J[注意:文件名中包含站SMI和J成分]
########## extract names of stations and components[#########提取站和组件的名称]
das1 = substr(DB$fn, 58, 60)
comp1 = substr(DB$fn, 62,62))
########### set up times (start and end) in julian days for selection[##########建立时间(开始和结束)胡利安选择天]
DB$t1 = DB$jd + DB$hr/24 + DB$mi/(24*60) + DB$sec/(24*3600)
DB$t2 = DB$t1 + DB$dur/(24*3600)
####### method 2:[######方法2:]
DB = makeDB("/Users/lees/Site/Santiaguito/SG09", pattern = "R", kind = 1)
udays = unique(paste(DB$yr, DB$jd))
sdays = as.numeric( unlist( strsplit(udays, split=" ") ) )
ye = sdays[seq(from=1, to=length(sdays), by=2)]
d = sdays[seq(from=2, to=length(sdays), by=2)]
o = order(ye+d/366)
days = list(yr = ye[o], jd=d[o])
das1 = DB$sta
comp1 = DB$comp
####### get unique stations from list:[######独特的站从列表:]
usta = unique(das1)
########### choose the desired components [##########中选择所需的组件]
acomp = c('I','V')
#### aday refers to one of the days listed in the days structure[###ADAY指天结构中列出的日子之一]
### days is a list of days (and associated years) that are in teh DB[##天是天的列表(以及相关的),在格兰DB]
days = list(jd=c(365,366, 1,2,3,4), yr=c(2008, 2008, rep(2009, times=4)))
#### extract days from the DB list:[###从数据库中提取天列表:]
udays = unique(paste(DB$yr, DB$jd))
sdays = as.numeric( unlist( strsplit(udays, split=" ") ) )
ye = sdays[seq(from=1, to=length(sdays), by=2)]
d = sdays[seq(from=2, to=length(sdays), by=2)]
o = order(ye+d/366)
days = list(yr = ye[o], jd=d[o])
####[###]
aday = 3
SAVEFILE = 'wpix.out'
ihour = 15
inkhour = .5
#### aday refers to one of the days listed in the days structure[###ADAY指天结构中列出的日子之一]
view.seis(aday, ihour, inkhour, SAVEFILE, days, DB, usta, acomp, TZ=(-6))
## End(Not run)[#(不执行)]
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注:
注1:为了方便大家学习,本文档为生物统计家园网机器人LoveR翻译而成,仅供个人R语言学习参考使用,生物统计家园保留版权。
注2:由于是机器人自动翻译,难免有不准确之处,使用时仔细对照中、英文内容进行反复理解,可以帮助R语言的学习。
注3:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
|