NAME

batoccultmap - compute BAT exposure mask due to occultation

USAGE

batoccultmap infile outfile saofile

DESCRIPTION

batoccultmap computes the fractional exposure in the BAT field of view due to occultation by the earth, moon and/or sun. The BAT field of view is large enough that any of these bodies may occult cosmic sources. Since the bodies move with time, the task computes the exposure map as a function of time and merges the results. This task is used for sky images; if you want to data selection for the visibility of a single source, use batoccultgti.

The task can be run in two modes. If occultation="any", then any portion of the image which is blocked at any time is masked in the output. If occultation="fraction", then the fractional exposure is computed (a value between 0 and 1).

The user must supply the "SAO" file created by prefilter. This file contains the positions of the earth, moon and sun as a function of time. The sizes of the earth, moon and sun are configurable.

PARAMETERS

infile [filename]
The name (and optional extension) of the input image file name.

outfile [filename]
The name of the output map file. For occultation="fraction", the map will range continuously from 0.0 to 1.0, corresponding to the partial exposure. For occultation="any", the map will be either 0.0 (occulted at some time) or 1.0 (never occulted).

saofile [filename]
The name of the "SAO" file produced by prefilter. This file should contain the columns TIME, RA, DEC, ROLL, SAT_ALT, MOON_RA, MOON_DEC, SUN_RA, SUN_DEC, EARTH_RA, EARTH_DEC.

atmdepth = 5.0 [real]
Thickness of the earth "atmosphere" in kilometers. Points within this thickness are assumed to have zero exposure.

rearth = 6378.1 [real]
The effective radius of the earth in kilometers. The default value is the mean equatorial radius (a conservative number).

rmoon = 16.0 [real]
The effective radius of the moon in arcminutes. The default value is conservatively set at 0.5 arcminutes larger than the true value.

rsun = 16.5 [real]
The effective radius of the sun in arcminutes. The default value is conservatively set at 0.5 arcminutes larger than the true value.

method = "time" [string]
Time segmentation method, used for subdividing good time intervals. For method="time", good time intervals are divided into subintervals no larger than timesegtol seconds. For method="position", good time intervals are divided so that the spacecraft and earth/moon/sun move no more than timesegtol arcminutes.

timesegtol = 0.0 [real]
Time segmentation parameter (see "method" above). Generally smaller values will provide higher fidelity, but take longer to compute. A value of zero indicates the defaults of 3 arcmin / 10 seconds should be used.

occultation [string]
Either "any" or "fraction". See "outfile" above.

constraints = "EARTH" [string]
A comma-separated list of the constraints to check (EARTH, MOON and/or SUN). Note that automatic inclusion of the MOON and SUN is not recommended, because both are comparable to the size of a BAT point source and may confuse the source detection software.

gtifile = "INFILE" [string]
Name of file containing the good time interval extension. The file is searched for the first extension matching "*GTI*". A value of "INFILE" indicates to search the input file. A value of "NONE" indicates to use the TSTART/TSTOP image header keywords.

(multfiles = "NONE") [string]
Comma-separated list of image filenames to be corrected using the calculated correction image (or "@filename.lis" to supply a list of files in filename.lis). Images are corrected IN PLACE (i.e. no back-up is made) by MULTIPLYING the original image values by the correction.

(divfiles = "NONE") [string]
Comma-separated list of image filenames to be corrected using the calculated correction image (or "@filename.lis" to supply a list of files in filename.lis). Images are corrected IN PLACE (i.e. no back-up is made) by DIVIDING the original image values by the correction.

(algorithm = "CONTOUR") [string]
Name of algorithm to use, either IMAGE or CONTOUR. The IMAGE method computes visibility for each pixel; the CONTOUR method computes the contours of each object in pixel coordinates. By virtue of the algorithms, the two results are slightly different (CONTOUR produces slightly larger occultation regions).

(clobber = NO) [boolean]
If the output file already exists, then setting "clobber = yes" will cause it to be overwritten.

(chatter = 2) [integer, 0 - 5]
Controls the amount of informative text written to standard output. Setting chatter = 1 produces a basic summary of the task actions; chatter = 2 (default) additionally prints a summary of input parameters; chatter = 5 prints debugging information.

(history = YES) [boolean]
If history = YES, then a set of HISTORY keywords will be written to the header of the specified HDU in the output file to record the value of all the task parameters that were used to produce the output file.

EXAMPLES

1. Computes the fractional exposure map for the primary HDU of sky.img. Good time intervals are subdivided so that the spacecraft attitude/earth/moon/sun do not move more than 3 arcmin (the default).

     batoccultmap sky.img occult.img sw00111622000sao.fits.gz \
       method=position occultation=fraction

SEE ALSO

batoccultgti

LAST MODIFIED

Apr 2007