NAME¶
funcen - find centroid (for binary tables)
SYNOPSIS¶
funcen [\-i] [\-n iter] [\-t tol] [\-v lev] <iname> <region>
OPTIONS¶
-i # use image filtering (default: event filtering)
-n iter # max number of iterations (default: 0)
-t tol # pixel tolerance distance (default: 1.0)
-v [0,1,2,3] # output verbosity level (default: 0)
DESCRIPTION¶
funcen iteratively calculates the centroid position within one or more
regions of a Funtools table (FITS binary table or raw event file). Starting
with an input table, an initial region specification, and an iteration count,
the program calculates the average x and y position within the region and then
uses this new position as the region center for the next iteration. Iteration
terminates when the maximum number of iterations is reached or when the input
tolerance distance is met for that region. A count of events in the final
region is then output, along with the pixel position value (and, where
available, WCS position).
The first argument to the program specifies the Funtools table file to process.
Since the file must be read repeatedly, a value of "stdin" is not
permitted when the number of iterations is non\-zero. Use Funtools Bracket
Notation to specify FITS extensions and filters.
The second required argument is the initial region descriptor. Multiple regions
are permitted. However, compound regions (accelerators, variable argument
regions and regions connected via boolean algebra) are not permitted. Points
and polygons also are illegal. These restrictions might be lifted in a future
version, if warranted.
The
\-n (iteration number) switch specifies the maximum number of
iterations to perform. The default is 0, which means that the program will
simply count and display the number of events in the initial region(s). Note
that when iterations is 0, the data can be input via stdin.
The
\-t (tolerance) switch specifies a floating point tolerance value. If
the distance between the current centroid position value and the last position
values is less than this value, iteration terminates. The default value is 1
pixel.
The
\-v (verbosity) switch specifies the verbosity level of the output.
The default is 0, which results in a single line of output for each input
region consisting of the following values:
counts x y [ra dec coordsys]
The last 3 WCS values are output if WCS information is available in the data
file header. Thus, for example:
[sh] funcen -n 0 snr.ev "cir 505 508 5"
915 505.00 508.00 345.284038 58.870920 j2000
[sh] funcen -n 3 snr.ev "cir 505 508 5"
1120 504.43 509.65 345.286480 58.874587 j2000
The first example simply counts the number of events in the initial region. The
second example iterates the centroid calculation three times to determine a
final "best" position.
Higher levels of verbosity obviously imply more verbose output. At level 1, the
output essentially contains the same information as level 0, but with keyword
formatting:
[sh] funcen -v 1 -n 3 snr.ev "cir 505 508 5"
event_file: snr.ev
initial_region: cir 505 508 5
tolerance: 1.0000
iterations: 1
events: 1120
x,y(physical): 504.43 509.65
ra,dec(j2000): 345.286480 58.874587
final_region1: cir 504.43 509.65 5
Level 2 outputs results from intermediate calculations as well.
Ordinarily, region filtering is performed using analytic (event) filtering, i.e.
that same style of filtering as is performed by
fundisp and
funtable. Use the
\-i switch to specify image filtering, i.e.
the same style filtering as is performed by
funcnts. Thus, you can
perform a quick calculation of counts in regions, using either the analytic or
image filtering method, by specifying the
\-n 0 and optional
\-i switches. These two method often give
different results because of how boundary events are processed:
[sh] funcen snr.ev "cir 505 508 5"
915 505.00 508.00 345.284038 58.870920 j2000
[sh] funcen -i snr.ev "cir 505 508 5"
798 505.00 508.00 345.284038 58.870920 j2000
See Region Boundaries for more information about how boundaries are calculated
using these two methods.
SEE ALSO¶
See
funtools(7) for a list of Funtools help pages