Initial commit

1 files changed, 83 insertions, 0 deletions

diff --git a/noao/astcat/src/agetcat/atwedit.x b/noao/astcat/src/agetcat/atwedit.x
new file mode 100644
index 00000000..7678eb35
--- /dev/null
+++ b/noao/astcat/src/agetcat/atwedit.x
@@ -0,0 +1,83 @@
+include <imhdr.h>
+include <pkg/cq.h>
+
+# Add a valid WCS to the image header if it does not already have one using
+# the image WCS status specified by the wcs keyword in the image survey
+# database. If the wcs keyord  is "fits", the image is assumed to have a
+# valid FITS WCS and no new wcs is computed, if it is "dss" the image is assumed
+# to have a valid DSS image header which will be transformed to a valid FITS
+# WCS if a FITS WCS is not already present, if it is "none" the image is
+# assumed to have no valid WCS and the code will attempt to insert one using
+# information in the image results structure. An error status is returned
+# only if there is no valid wcs code.
+
+procedure at_wedit (im, res, at, wcstype, update, verbose)
+
+pointer	im			#I the input image descriptor
+pointer	res			#I the image query results descriptor
+pointer	at			#I the astrometry package descriptor
+int	wcstype			#I the default wcs type
+bool	update			#I actually update the header ?
+bool	verbose			#I verbose mode ?
+
+int	cq_istati(), at_mkdss(), at_dbwcs(), at_parwcs()
+
+begin
+	# Update WCS from database
+	if (res != NULL) {
+
+	    switch (cq_istati (res, CQWCS)) {
+
+	    # Image surveys database indicates image already has a FITS WCS.
+	    case CQ_WFITS:
+	        ;
+
+	    # Image surveys database indicates image has a DSS WCS.
+	    case CQ_WDSS:
+	        if (at_mkdss (im, update, verbose) == ERR) {
+	            #if (update || verbose)
+	            if (verbose)
+		        call printf (
+			    "    Error converting DSS wcs to FITS wcs\n")
+	        }
+
+	    # Image surveys database indicates image has no WCS. If the proper
+	    # information is not in the image survey then default to awcspars.
+	    default:
+	        if (at_dbwcs (im, res, update, verbose) == ERR) {
+	            #if (update || verbose)
+	            if (verbose)
+		        call printf (
+		        "    Error creating FITS wcs using image survey db\n")
+	        }
+	    }
+
+	} else {
+
+	    switch (wcstype) {
+
+	    # User parameter indicates image already has a FITS WCS.
+	    case CQ_WFITS:
+		;
+
+	    # User parameter indicates image has a DSS WCS.
+	    case CQ_WDSS:
+	        if (at_mkdss (im, update, verbose) == ERR) {
+	            #if (update || verbose)
+	            if (verbose)
+		        call printf (
+			    "    Error converting DSS wcs to FITS wcs\n")
+	        }
+
+	    default:
+	        if (at == NULL)
+		    ;
+	        else if (at_parwcs (im, at, update, verbose) == ERR) {
+	            #if (update || verbose)
+	            if (verbose)
+		        call printf (
+		       "    Error creating FITS wcs using default parameters\n")
+	        }
+	    }
+	}
+end