I've assembled a list of links to bring some order to the many astrometery codes I have been using. Most of these routines were then used to build a series of demostrations scripts (SimIm0, SimIm1, ...). Most of the scripts listed above are easy to use and well-documented (both online help via "-- help" and webpage documents). To demonstarte their use I am composing a set of scripts that allow some user interaction at runtime. The basic steps are to (1) build a simulated image, (2) apply a coordinate transformation to the image, (3) intsall WCS header in each of the two images that yield uniform astrometry.

1. Building a simulated image using ellradius.sh and make2dimg.sh (via model_images_1).
2. Play with coordinate transformation basics with the trs_ codes (Optional Demo Material)
3. Use clip_imshift_trs to transform the siumulated image.
4. Compose WCS header cards for a FITS image.
5. Install WCS header cards for a FITS image.
6. Use image_pa (and image_ps) to visualize sky compass.

Most of the scripts listed above are easy to use and well-documented (both online help via "-- help" and webpage documents). The goal here is to create a terse but clear sequence of commands that demontsrate the path to a WCS-calibrated image. The commands should closely follow the links above. The links (mostly to Big_Code_Table entries) give usage details.

The primary purose of these exercises is is t demsonstrate that I can transform images and properly compensate for these changes in the final WCS of the header, I manually made a region file of fk5 positions to aid in the demmo. If the region file (named "Test0.reg") is present in the directory where we run the SimIm* scripts, then the positions are overplotted in the original and the transformed images. Here is the region file I use:

% cat Test0.reg
# Region file format: DS9 version 4.1
# Filename: Modelwcs.fits
global color=green dashlist=8 3 width=1 font="helvetica 10 normal roman" select=1 highlite=1 dash=0 fixed=0 edit=1 move=1 delete=1 include=1 source=1
fk5
circle(12:00:04.763,+32:30:58.14,15.606") # color=black font="helvetica 12 normal roman" text={12:00:04.763  +32:30:58.14}
circle(11:59:57.291,+32:30:00.76,18.1465") # color=black font="helvetica 12 normal roman" text={11:59:57.291  +32:30:00.76}