integratePower: Determine the total power within a given bandwidth

In astrochron: A Computational Tool for Astrochronology

Determine the total power within a given bandwidth

Description

Determine the total power within a given bandwidth, and also the ratio of this power to the total power in the spectrum (or up to a specified frequency). If bandwidth is not specified, generate interactive plots for bandwidth selection. For use with the function eha, integratePower can process spectrograms (time-frequency) or single spectra.

Usage

integratePower(spec,flow=NULL,fhigh=NULL,fmax=NULL,unity=F,f0=T,xmin=NULL,
               xmax=NULL,ymin=NULL,ymax=NULL,npts=NULL,pad=NULL,ydir=1,
               palette=6,ncolors=100,h=6,w=9,ln=F,genplot=T,verbose=T)

Arguments

spec	Spectral results to evaluate. If the data frame contains time-frequency results, it must have the following format: column 1=frequency; remaining columns (2 to n)=power; titles for columns 2 to n must be the location (depth or height). Note that this format is ouput by function eha. If the data frame contains one spectrum, it must have the following format: column 1=frequency, column 2=power.
flow	Low frequency cutoff for integration. If flow or fhigh are not specified, interactive plotting is activated.
fhigh	High frequency cutoff for integration. If flow or fhigh are not specified, interactive plotting is activated.
fmax	Integrate total power up to this frequency.
unity	Normalize spectra such that total variance (up to fmax) is unity. (T of F)
f0	Is f(0) included in the spectra? (T or F)
xmin	Minimum frequency for PLOTTING.
xmax	Maximum frequency for PLOTTING.
ymin	Minimum depth/height/time for PLOTTING. Only used if processing time-frequency results.
ymax	Maximum depth/height/time for PLOTTING. Only used if processing time-frequency results.
npts	The number of points in the processed time series window. This is needed for proper spectrum normalization.
pad	The total padded length of the processed time series window. This is needed for proper spectrum normalization.
ydir	Direction for y-axis in plots (depth or height). -1 = values increase downwards (slower plotting!), 1 = values increase upwards. Only used if processing time-frequency results.
palette	What color palette would you like to use? (1) rainbow, (2) grayscale, (3) blue, (4) red, (5) blue-white-red, (6) viridis
ncolors	Number of colors to use in plot. Only used if processing time-frequency results.
h	Height of plot in inches.
w	Width of plot in inches.
ln	Plot natural log of spectral results? (T or F)
genplot	Generate summary plots? (T or F)
verbose	Verbose output? (T or F)

Details

Depending on the normalization used, you may want to preprocess the power spectra prior to integration.

See Also

eha

Examples

# generate etp signal over past 10 Ma
ex=etp(tmax=10000)

# evolutive power
pwr=eha(ex,win=500,fmax=.1,pad=2000,output=2,pl=2)

# integrate power from main obliquity term
integratePower(pwr,flow=0.02,fhigh=0.029,npts=501,pad=2000)

astrochron documentation built on Sept. 30, 2024, 9:14 a.m.