longfft: Long FFT Spectrogram

In RSEIS: Seismic Time Series Analysis Tools

Long FFT Spectrogram

Description

Creates hourly spectrograms, either alternating seismic and infrasound data or sequences of one component.

Usage

longfft(DB, DAYS = c(233, 234), HRS = 1:24, sta = "KR1", comp = c("V",
"I"), NPP = 6, CSCALE = FALSE, pal = rainbow(100), PS = FALSE, kind = 1,
Iendian = 1, BIGLONG = FALSE)

longreset(NPP, PS)

longpstart(NPP = 6, asta = "", acomp = "", theday = 1, hr = 0)

Arguments

DB	RSEIS Data base
DAYS	vector of Days to display
HRS	vector of hours to display
sta	stations to extract
comp	component to extract
NPP	Number of plot strips per page, default = 6
CSCALE	scaling
pal	palettes to use (given two will alternate these)
PS	logical, TRUE postscript output
kind	data type, an integer -1, 0, 1, 2 ; 0=R(DAT) , -1=RDS, 0=RDATA, 1 = segy, 2 = sac
Iendian	Endian-ness of binary data
BIGLONG	logical, TRUE=long is 8 bytes
asta	character, one station
acomp	character, one component
theday	one day
hr	one hour

Details

Extracts dats from the DB data base and plots strips of spectrograms for perusal.

longpstart, longreset are auxilliary codes used to set up the postscript files and initialize the plotting.

Value

Graphical Side effects

Note

Program is set for data being ready from external sources in binary (SAC, SEGY) format. If data is in R-format already, the code may not work.

Author(s)

Jonathan M. Lees<[e-mail]>

See Also

SPECT.drive

Examples

if(interactive()){
 #####  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 = 24,
seed=200, noise.est=c(1, 100) , verbose=TRUE  )
     
    
tdir = tempdir()
for(i in 1:length(JK) )
{
    sig = JK[[i]]
  d1 =  dateStamp(sig$DATTIM, sep='_')
   nam1 =  paste(d1,sig$sta, sig$comp, sep='_')
    nam2 = paste0(nam1, '.RDS')
    nam3 = paste(tdir, nam2, sep='/')
    saveRDS(file=nam3, sig) 
    }

LF  = list.files(path=tdir,pattern='.RDS',  full.names=TRUE)
DB = FmakeDB(LF, kind=-1)
IDB = infoDB(DB)

   p1 <- RPMG::Gcols(plow=5, phi=0,  N=100, pal="topo.colors", mingray=0.8)
     p2 <- RPMG::Gcols(plow=5, phi=0,  N=100, pal="rainbow", mingray=0.8)
    
longfft(DB, DAYS=5 , HRS=1:24 ,
     sta=IDB$usta, comp=IDB$ucomp , NPP=6 , CSCALE=FALSE, 
     pal = list(p1=p1, p2=p2), PS=FALSE ,  kind = -1,
     Iendian=1, BIGLONG=TRUE )


}

RSEIS documentation built on Sept. 13, 2024, 1:09 a.m.