Projects image from one WCS to another
reproject_image infile matchfile outfile [resolution] [method]
[coord_sys] [lookupTab] [clobber] [verbose]
`reproject_image' maps an image in one WCS reference frame to
another. For each output pixel in the matchfile it maps an
n-point polygon
onto the input image via the WCS transforms. The area of each
pixel covered in the input image is then used to compute the
output pixel value.
`reproject_image' can be used to project an image w/ one tanget
point to another tanget point as is useful when merging multiple
observations or when matching data across missions.
Any Null/NaN valued pixels or pixels outside the data-subspace or
pixels outside the input image are assigned a value = 0.
reproject_image img.fits matchfile=new_tan.fits outfile=reproject.fits
The output image will have the same size and WCS as the 'matchfile'
file. Each pixel in 'matchfile' will be mapped to pixels in the
'infile.
reproject_image XMM.fits matchfile=Chandra.fits outfile=myout.fits
Projects an XMM file to the same WCS (tanget-point) as a Chandra image.
reproject_image @img.lis matchfile=new_tan.fits outfile=reproject.fits
Similar to first example but uses a stack as input. Each input
image is reprojected to the output coord system. The images are
then co-added. The 'method' normalization is applied before
the data are summed.
The reference image.
The 'infile' image is mapped to the same WCS as this image. The
'outfile' will have the same dimensions and WCS as this image.
The output file name
Contains the reprojected image.
Controls quality of projection; number of points per side of polygon
An n-sided polygon that outlines the output pixel is mapped to input
image via the WCS transforms. The number of points along each side
of the pixel is the resolution parameter. A value of '1' indicated
that just the corners of the pixel will be used. A value of '2'
indicates that the corner and the middle of the pixel-edge line
segments will be use. And so on. The more points on the polygon,
the better the polygon will approximate the possibly non-linear
transform between the images. However, the more points on the
polygon, the longer the run-time of the tool.
A value of '0' can also be used. This is a special quick mode that
simply maps the center output pixel to a single input pixel and uses
that value in the output image. This can be done very quickly;
however, when something other than a simple shift of the images is
need this can result in image artifacts (for example aliasing or
'dead' regions). Also the 'method' is not applicable in this
useage. The output image is essentially "interpolated", so if the
pixel scales are different the flux will not be preserved.
Controls output normalization
The output image can either represent a conservation of "sum" (
integral over an apperature on the input and output image would give
same value) or it can represent an "average" where the
output pixel value represents an average input value.
Typically users will use "sum" to reproject the COUNTS image and
will use "average" to reproject the EXPOSURE image when making
fluxed images.
Coordinate system to do the pixel mapping in
Currently only "world" coordinate system is supported.
Logical pixels in the output image are mapped to physical pixels
which are then mapped to world coordinate (RA,Dec). These are then
mapped back to physical pixels in the input image and then back to
image pixels in the input image. [TBD: implement "physical" and
"image" coord sys's]
The header merging table
Rules to merge the headers when more than one file supplied. If
set to NONE or a blank string then the header from the first file is used.
Remove output if it exists?
Used to specify whether or not to clobber existing file that
has the same name as the specified output file
The tool chatter level
Verbose can be from 0 to 5, generating different amounts
of debugging output.
|