RSEIS (version 4.2-4)

FAKEDATA: Fake Data for Examples.

Description

Create a list of artifical seismic traces to illustrate examples that require a database or long sequences.

Usage

FAKEDATA(amp, OLDdt = 0.01, newdt = 0.1, yr = 2000,
JD = 5, mi = 0, sec = 0, Ntraces = 48, seed = 200,
noise.est = c(1, 100), verbose = FALSE)

Value

List of data in a format similar to the output of GET.seis.

Arguments

amp

vector, some signal that will be repeated

OLDdt

Orignal sample rate

newdt

New sample rate, usually less than the original

yr

year

JD

starting Julian day

mi

starting minute

sec

starting second

Ntraces

number of traces

seed

random seed

noise.est

2-vector, starting and ending sample to estimate noise level of trace

verbose

logical, message feed back

Author

Jonathan M. Lees<[e-mail]>

Details

The input signal can be any time series, or even a made up signal. This is just to give the look of the result something like real data. The noise level is extracted from the man and std of the real data at the samples indicated by noise.est.

The sampling rate (dt, sec/sample ) is increased mainly for speed and plotting. This may be skipped for certain functions involving spectrum analysis.

The signal is distributed randomly in each hour along the total span of the requested period, i.e. each hour has one instance of the signal.

The date is arbitrary, of course.

See Also

GET.seis

Examples

#####  get a time series
data(KH)

amp = KH$JSTR[[1]]
OLDdt = KH$dt[1]
####  downsample to:
newdt = 0.1

JK = FAKEDATA(amp, OLDdt=OLDdt, newdt = 0.1, yr = 2000,
         JD = 4, mi = 12, sec = 0,  Ntraces = 3,
seed=200, noise.est=c(1, 100) , verbose=TRUE  )

 op <- par(no.readonly = TRUE)
par(mfrow=c(length(JK), 1) )
for(i in 1:length(JK) )
{
    DATTIM = paste(c(unlist(JK[[i]]$DATTIM), JK[[i]]$N), collapse=' ')
    
    plotGH( JK[[i]] )
    mtext(DATTIM, side=3, at=JK[[i]]$DATTIM$t2/2)
    
}
 par(op)