NAME

fgabor -- Perform a Gabor transform on a lightcurve and return the normalized power as a function of time and frequency.


USAGE

        fgabor infile fmin fmax nfreq window outfile timeunits ratecol extname maxgap 
          winfrac sigma normalize quiet clobber

DESCRIPTION

The Gabor transform takes a function of one variable, time, and maps it to a function of two variables, time and frequency:


f(t) -> G(t,k)

For each frequency k, the function of time only G(*,k) is the result of convolving the function f with a kernel which is the product of a fixed Gaussian window and the complex exponential of frequency k. So of G(*,k) is a filtered version of the original function (or signal) f.


PARAMETERS

infile [string - file name]
The name of the input FITS file containing the lightcurve to be analyzed. Any OGIP-compliant lightcurve should work, though you may need to set the ratecol and/or extname parameters appropriately (see below).

fmin [real - minimum frequency]
The minimum frequency to be considered in the transform. Should be greater than ndata/window. Measured in either Hertz or cycles per dataset, depending on the setting of the timeunits parameter (see below).

fmax [real - maximum frequency]
The maximum frequency to be considered in the transform. Should be less than than the Nyquist frequency. Measured in either Hertz or cycles per dataset, depending on the setting of the timeunits parameter (see below).

nfreq [integer - number of frequencies]
The number of frequencies to use between fmin and fmax. Must be at least two. Too high a value will result in a very large output file and will consume proportionally more system resources during the computation.

window = INDEF [real - window size]
The size (+/- 3 sigma) of the Gaussian window to use in the transform. Should be specified in either number of timesteps or in seconds, depending on the setting of the timeunits parameter (see below, also fmin and fmax). is used. If set to INDEF, the tool will choose a window size based on the size of the dataset and the winfrac parameter (see below).

outfile = gaborpwr.fits [string - filename]
The name for the output FITS file containing the normalized power array and a table of frequencies. Use the clobber parameter to allow the tool to overwrite a preexisting file.

(timeunits = yes) [boolean - use time units]
If "yes" then the fmin, fmax and window parameters will be interpreted as units of Hz and seconds, respectively. If "no" then they will be taken to be in cycles per dataset and number of rows instead.

(ratecol = RATE) [string - rate column name]
The name of the column in the input file containing the rate information.

(extname = RATE) [string - lightcurve extension name]
The name of the extension (EXTNAME keyword) which contains the lightcurve in the input FITS file.

(maxgap = 0.10) [real - maximum gapsize]
Defines the maximum number of consecutive rows which are missing or which contain NULL rates (as a fraction of the total size of the dataset). Up to this limit rows which contain NULL rate values or are missing will be replaced by random Gaussian noise prior to performing the Gabor transform.

(winfrac = 0.20) [real - default INDEF window]
Default window size (as a fraction of the total size of the dataset) used when the window parameter (see above) is set to INDEF.

(sigma = INDEF) [real - Standard deviation for random noise]
When set to INDEF (the default) the Gaussian used to fill NULL and/or missing values will be defined by the computed standard deviation of the existing rate values. This parameter allows users to instead set this standard deviation explicitly.

(normalize = yes) [boolean - Normalize output power?]
The default behavior is to normalize the output power to 1. If this parameter is set to "no" the normalization will not be done.

(quiet = no) [boolean]
If set to "yes" all standard output will be supressed.

(clobber = no) [boolean]
Overwrite existing output file having the same name?

(mode = ql) [string]
FTOOLS internal -- users should not need to change.


EXAMPLES


NOTES:

The output FITS file will include three data units. The first (primary array) will contain the result of the Gabor transform: a FITS image with dimensions of nfreq x nrows. The second data unit (first extension) is a binary table listing the frequencies used in the tranform, in units of either Hz or cycles per dataset, depending on the setting of the timeunits parameter. The third data unit (second extension) is a binary table which contains a copy of the data on which the transform was actually performed, ie with the mean value subtracted and with time gaps and NULL values filled with random, Gaussian-distributed noise.


BUGS

Please report problems to xtehelp@athena.gsfc.nasa.gov.


SEE ALSO

CATEGORY

Apr2001 ftools.xte