The following image astrometry tasks have been installed in the new IRAF 2.11 imcoords package.
The main function of the new tasks is to compute plate solutions for images and store them in the image headers in a FITS compatible format. The new tasks can also be used for computing accurate pixel and world coordinates for objects in images and for determining geometric parameters for images such as the plate scale and orientation. These tasks do not comprise a fully featured astrometry package. However, they have been successfully integrated into the prototype astrometry package tfinder (see article in this Newsletter) and the NOAO Mosaic reduction package mscred (see article in this Newsletter).
The examples that follow demonstrate the basic functionality of the new tasks. To reproduce these examples on a local system copy dev$wpix to a local directory and edit in the missing EQUINOX keyword as shown below.
Get the Reference Object List
# if it exists
cl> ccfind wpix.coo wpix lngref=13:27:47 \
>>> latref=47:27:14 refsystem=B1950 \
>>> xmag=-0.77 ymag=0.77
Compute the Plate Solution
The ccmap task computes the plate solution using the matched celestial and pixel coordinate list produced by ccfind and the coordinates of the tangent point, then stores the computed solution in a database file for later use by the ccsetwcs and cctran tasks. The computed solution is defined in the coordinate system of the input coordinates, and J2000.0 is assumed by default. The task ccmap may be run interactively as shown below or non-interactively.
By default ccmap performs a standard 6-coefficient linear fit to the input coordinates. Higher order fits to the residuals from the linear fit are supported by ccmap and the plate solution evaluation task cctran . However, due to the lack of an FITS convention for representing higher order plate solutions, only the linear part of the computed plate solution is normally stored in the image header by the ccmap and ccsetwcs tasks. Users running IRAF V2.11.1 or later versions can force the entire plate solution to be recorded in the image header by setting the ccmap projection parameter to "tnx". All IRAF tasks which compute world coordinates, e.g., starfind , imexamine , wcsctran , etc., understand the "tnx" system. Users planning to export their IRAF images to another image processing system should avoid using the "tnx" projection if they are concerned about portable coordinate systems.
The ccsetwcs task creates the image world coordinate system from the computed plate solution as shown in the example below, or from a list of parameters supplied by the user.
Once the plate solution is computed and stored in the image header either the new image world coordinate system or the cctran task can be used to compute celestial coordinates for program objects. For 6-coefficient linear plate solutions the two methods are equivalent; for higher order plate solutions computed prior to the release of IRAF 2.11.1 the ccmap database record must be used to retrieve the full accuracy of the plate solution. In the following examples the starfind task is used to compute pixel and world coordinates for the objects that it detects in the image, and cctran uses the plate solution to convert the starfind pixel coordinates into accurate celestial coordinates.