JGET.seis {RSEIS} | R Documentation |
This fuction calls c-routines to read in ‘segy’, ‘sac’.
JGET.seis(fnames, kind = 1, PLOT = FALSE) getseisinfo(fnames, kind = 1)
fnames |
list of file names. |
kind |
an integer 1, 2, 3; 1 = segy, 2 = sac, 3 = AH. |
PLOT |
logical, for if the traces should be plotted. |
JGET.seis extracts digital seismic data from binary files stored in the file system. The program uses C-code for I/O and passes data back to R. Currently SAC, SEGY formats are installed but it is easy to extend. AH format is available for LINUX systems, but there were problems compiling in WINDOWS and MACOS so this feature was removed. A filter for mseed format is currently being developed.
List containing the seismic data and header information. Each trace consists of a list with:
fn |
original file name |
sta |
station name |
comp |
compnent |
dt |
delta t in seconds |
DATTIM |
time list |
yr |
year |
jd |
julian day |
mo |
month |
dom |
day of month |
hr |
hour |
mi |
minute |
sec |
sec |
msec |
milliseconds |
dt |
delta t in seconds |
t1 |
time start of trace |
t2 |
time end of trace |
off |
off-set |
N |
number of points in trace |
units |
units |
amp |
vector of trace values |
getseisinfo gets only the header information from the trace. Information in the file names is ignored, so be sure to modify headers prior to using this method of extracting meta-data. (Or modify the meta data from the file names after reading in the data.)
For SEGY files, in LINUX-UNIX, use: rename, segymod (PASSCAL) to modify the headers
Jonathan M. Lees <jonathan.lees@unc.edu>
Mine.seis, Package:Rsac
## Not run: Lname = list.files(path='/data/wadati/bourbon/LaurelCanNC/R005.01' , pattern="08.005.01.41.23.9024", full.names=TRUE) S1 = JGET.seis(Lname, kind = 1, PLOT = FALSE) ## End(Not run)