makeDB(RSEIS)
makeDB()所属R语言包：RSEIS
Create a seismic Waveform Database
创建一个地震波形数据库
译者：生物统计家园网 机器人LoveR
描述----------Description----------
Create a seismic Waveform Database
创建一个地震波形数据库
用法----------Usage----------
makeDB(path, pattern, kind = 1, Iendian=1, BIGLONG=FALSE)
参数----------Arguments----------
参数：path
Path to directory where files and directories exist
目录的路径，文件和目录是否存在
参数：pattern
pattern for listing of files
模式的文件列表
参数：kind
kind of data: R(DAT)=0, segy=1; sac=2
种数据：R（DAT）= 0，SEGY = 1;囊= 2
参数：Iendian
default=1Endian-ness of the data: 1,2,3: "little", "big", "swap". Default = 1 (little)
默认情况下，的= 1Endian性的数据：1,2,3：“小”，“大”，“交换”。默认值= 1（小）
参数：BIGLONG
logical, TRUE means long=8 bytes
长= 8个字节的逻辑，TRUE手段
Details
详细信息----------Details----------
The files are typically located in a directory structure created by programs like ref2segy, a PASSCAL program for downloading data in the field. Each file contains one seismogram, with a header. makeDB reads in all the headers and creates a list of meta-data for later use in RSEIS.
这些文件通常位于在目录结构中，像ref2segy，一个PASSCAL的程序字段中的数据下载的程序。每个文件都包含一个地震记录，一个标头。 makeDB中读取所有的标头和购买在RSEIS使用创建的元数据的列表。
Uses readBin to extract data in SAC format. user must know what kind of machine the data was created on for I/O purposes.
使用readBin在SAC格式中提取数据。用户必须知道什么样的机器上创建的数据I / O的目的。
If data was created on a little endian machine but is being read on big endian machine, need to call the endian "swap" for swapping.
如果数据创建一个little endian机器上，但正在读big endian机器上，需要调用的字节序的交换“交换”。
If data was created on a machine with LONG=4 bytes, be sure to call the program with BIGLONG=FALSE.
如果数据在计算机上创建LONG = 4个字节，BIGLONG = FALSE，一定要调用该程序。
值----------Value----------
list:
列表：
参数：fn
file name
文件名
参数：yr
year
年
参数：jd
julian day
儒略日
参数：hr
hour
小时
参数：mi
minute
分钟
参数：sec
second
第二
参数：dur
duration, seconds
时间，秒
参数：t1
time 1 in Epoch days
在大纪元天的时间1
参数： t2
time 2 in Epoch days
在大纪元天的时间2
参数：sta
station name
站名
参数：comp
component name
组件名称
参数：origyr
origin year for epoch times
大纪元的起源一年
注意----------Note----------
Epoch times are used to accomodate problems where julian days cross year end boundaries, so that day 366 comes before day 1 of the next year.
大纪元时报，以适应问题，朱利安天跨越年底边界的，所以明年1天前的那一天366。
（作者）----------Author(s)----------
Jonathan M. Lees<jonathan.lees@unc.edu>
参见----------See Also----------
setupDB, Mine.seis , getseis24, plotseis24, EPOCHday, swig
setupDB，Mine.seis，getseis24，plotseis24，EPOCHday，痛饮
实例----------Examples----------
## Not run: [＃不运行：]
############# set directory[＃＃＃＃＃＃＃＃＃＃＃＃设置目录]
path = '/home/lees/Site/Santiaguito/SG09'
pattern = "R0*"
### get DB information[＃＃DB信息]
XDB = makeDB(path, pattern, kind =1)
##### select a station[＃＃＃＃中选择一个电台]
usta = "CAL"
acomp = "V"
##### extract 24 hours worht of data[＃＃＃＃24小时worht的数据提取]
JJ = getseis24(DB, 2009, 2, usta, acomp, kind = 1)
##### plot 24 hours worth of data[＃＃＃＃图24个小时的数据]
pjj = plotseis24(JJ, dy=1/18, FIX=24, SCALE=0, FILT=list(ON=FALSE, fl=0.05 , fh=20.0, type="BP", proto="BU"), RCOLS=c(rgb(0.2, .2, 1), rgb(.2, .2, .2)) )
########## window a small portion on 24 hour display[＃＃＃＃＃＃＃＃＃窗口中一个小的24小时显示部上]
w = winseis24(pjj)
### open a new window[＃＃打开一个新窗口]
X11()
VNE = c("V", "N", "E")
### IJK = c("I", "J", "K")[＃＃IJK = C（“I”，“J”，“K”）]
gsta=c("CAL", "KAM", "DOM", "LAV")
ucomp=VNE
### or: ucomp=IJK[＃＃：ucomp = IJK]
### set epoch day and get the times[＃＃设置划时代的一天，时代]
eday = EPOCHday(w$yr, jd = w$jd, origyr = DB$origyr)
######## using epoch day, set times for beginning and end of window[＃＃＃＃＃＃＃划时代日，设置开始和结束的时间窗口]
at1 = eday$jday + (w$hr[1])/24
at2 = eday$jday + (w$hr[2])/24
######## extract data from the data base:[＃＃＃＃＃＃＃从数据库中提取数据：]
GH = Mine.seis(at1, at2, DB, gsta, ucomp)
### show data:[＃＃显示的数据：]
swig(GH)
## End(Not run)[＃（不执行）]
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注：
注1：为了方便大家学习，本文档为生物统计家园网机器人LoveR翻译而成，仅供个人R语言学习参考使用，生物统计家园保留版权。
注2：由于是机器人自动翻译，难免有不准确之处，使用时仔细对照中、英文内容进行反复理解，可以帮助R语言的学习。
注3：如遇到不准确之处，请在本贴的后面进行回帖，我们会逐渐进行修订。
|
发表于 2012-9-28 21:34:45
