SYNOPSIS
mAddCube [-p imgdir] [-n(o-areas)] [-a mean|median|count] [-e(xact-size)] [-d level] [-s statusfile] images.tbl template.hdr out.fits
DESCRIPTION
Coadd the reprojected cubes in an input list to form an output mosaic with FITS header keywords specified in a header file. Creates two output files, one containing the coadded pixel values for each cube plane, and the other containing coadded pixel area values. The pixel area values can be used as a weighting function if the output cube values are themselves to be coadded with other projected cubes, and may also be used in validating the fidelity of the output pixel values.
OPTIONS
- -p imgdir
- Specifies path to directory containing reprojected images. If the -p switch is not included, mAddCube will look for the input images in the current working directory.
- -n
- Co-addition ignores weighting by pixel areas and performs coaddition based only on pixel postions. Will not output an area image for the output image.
- -a type
-
Select type of averaging to be done on accumulated pixel values (either mean or median). To generate a map showing counts of how many times each pixel was overlapped by the input images, use count.
- -e
- Enables exact size mode. The output image will match the header template exactly, instead of shrinking the output to fit the data.
- -d level
- Turns on debugging to the specified level (1-3).
- -s statusfile
-
mAddCube output and errors will be written to statusfile instead of stdout.
ARGUMENTS
- images.tbl
- ASCII table (generated by mImgtbl) containing metadata for all images to be coadded. For cubes, mImgtbl needs to be run with the -C option which picks up the third dimension parameters (and fourth if they exist).
- template.hdr
- FITS header template to be used in generation of output FITS cube.
- out.fits
-
Name of output FITS image.
RESULT
If successful, mAddCube creates a FITS datacube file out.fits that is a coadd of all the FITS cubes in the table images.tbl, according to the header template given. A corresponding out_area.fits is also created.
MESSAGES
- OK
- [struct stat = "OK", time=time]
- ERROR
- Invalid argument for -a flag
- ERROR
- Cannot open status file: statusfile
- ERROR
- Invalid image metadata file: filename
- ERROR
- Need columns: cntr,fname, crpix1, crpix2, cdelt1, cdelt2, naxis1, naxis2, crval1, crval2, naxis3 in image list
- ERROR
- Too many open files
- ERROR
- Input wcsinit() failed.
- ERROR
- Bad WCS in header template.
- ERROR
- FITS library error
- ERROR
- general error message
- ERROR
- Allocation failed for file info (realloc)
- ERROR
- Allocation failed for file info structs
- ERROR
- Allocation failed for area file info structs
- ERROR
- Allocation failed for input buffer
- ERROR
- Allocation failed for input area buffer
- ERROR
- Allocation failed for start/end info
- ERROR
- Allocation failed for data line pointers
- ERROR
- Allocation failed for data line
- ERROR
- Allocation failed for area line pointers
- ERROR
- Allocation failed for area line
- ERROR
- Allocation failed for data counts
- ERROR
- Allocation failed for output data line
- ERROR
- Allocation failed for output area line
- ERROR
- Allocation failed for data line (realloc)
- ERROR
- Allocation failed for area line (realloc)
- ERROR
- Allocation failed for median array
- ERROR
- Allocation failed for median array (realloc)
- ERROR
- Allocation failed for linked list structs
- ERROR
-
Allocation failed for linked list structs (additions)
EXAMPLES
The following example runs mAddCube on 4 FITS images, generating the output file mosaic.fits. Related files are images.tbl and template.hdr.
- $ mAddCube -p proj images.tbl template.hdr output/mosaic.fits
-
[struct stat="OK", time=1]
BUGS
The drizzle algorithm has been implemented but has not been tested in this release.If a header template contains carriage returns (i.e., created/modified on a Windows machine), the cfitsio library will be unable to read it properly, resulting in the error: [struct stat="ERROR", status=207, msg="illegal character in keyword"]
It is best for the background correction algorithms if the area described in the header template completely encloses all of the input images in their entirety. If parts of input images are "chopped off" by the header template, the background correction will be affected. We recommend you use an expanded header for the reprojection and background modeling steps, returning to the originally desired header size for the final coaddition. The default background matching assumes that there are no non-linear background variations in the individual images (and therefore in the overlap differences). If there is any uncertainty in this regard, it is safer to turn on the "level only" background matching (the "-l" flag in mBgModel.
COPYRIGHT
2001-2015 California Institute of Technology, Pasadena, CaliforniaIf your research uses Montage, please include the following acknowledgement: "This research made use of Montage. It is funded by the National Science Foundation under Grant Number ACI-1440620, and was previously funded by the National Aeronautics and Space Administration's Earth Science Technology Office, Computation Technologies Project, under Cooperative Agreement Number NCC5-626 between NASA and the California Institute of Technology."
The Montage distribution includes an adaptation of the MOPEX algorithm developed at the Spitzer Science Center.