This simply named script performs most of the critical WCS calibration work:
% wcs_getA name.fits
The image "name.fits" is the original CCD image we have begun
to calibrate.
The wcs_getA and wcs_getB scripts are mostly the same with one important exception: the way they gather the initial set of paired X,Y and Ra,Dec sets. In the getA case, we presumably have only a rough RA,DEC image position in the name.fits header. Hence, in order to find matching sets of X,Y and Ra,Dec values for stars in the CCD image field, we use the wcs_prep script to collect star catalogs (from USNO, SDSS, and 2MASS) and a DSS image matching the position of our CCD image. In the getB case, we use the a ds9 gui inspection (with A.fits) to gather a new (larger) set of input star pair sets and derive a new calibrated image (B.fits).