Repatch (from linux) of OSX IRAF

1 files changed, 106 insertions, 0 deletions

diff --git a/noao/imred/quadred/src/ccdproc/ccdproc.x b/noao/imred/quadred/src/ccdproc/ccdproc.x
new file mode 100644
index 00000000..1b2a133c
--- /dev/null
+++ b/noao/imred/quadred/src/ccdproc/ccdproc.x
@@ -0,0 +1,106 @@
+include	<error.h>
+include	"ccdred.h"
+include	"ccdtypes.h"
+
+# CCDPROC -- Process a CCD image of a specified CCD image type.
+#
+# The input image is corrected for bad pixels, overscan levels, zero
+# levels, dark counts, flat field, illumination, and fringing.  It may also
+# be trimmed.  The checking of whether to apply each correction, getting the
+# required parameters, and logging the operations is left to separate
+# procedures, one for each correction.  The actual processing is done by
+# a specialized procedure designed to be very efficient.  These
+# procedures may also process calibration images if necessary.
+# The specified image type overrides the image type in the image header.
+# There are two data type paths; one for short data types and one for
+# all other data types (usually real).
+
+procedure ccdproc (input, ccdtype)
+
+char	input[ARB]		# CCD image to process
+int	ccdtype			# CCD type of image (independent of header).
+
+pointer	sp, output, str, in, out, ccd, immap()
+errchk	immap, set_output, ccddelete
+errchk	set_fixpix, set_zero, set_dark, set_flat, set_illum, set_fringe
+
+begin
+	call smark (sp)
+	call salloc (output, SZ_FNAME, TY_CHAR)
+	call salloc (str, SZ_LINE, TY_CHAR)
+
+	# Map the image, make a working output image and set the processing
+	# parameters.
+
+	in = immap (input, READ_ONLY, 0)
+	call mktemp ("tmp", Memc[output], SZ_FNAME)
+	call set_output (in, out, Memc[output])
+	call set_proc (in, out, ccd)
+	call set_sections (ccd)
+	call set_trim (ccd)
+	call set_fixpix (ccd)
+	call set_overscan (ccd)
+
+	# Set processing appropriate for the various image types.
+	switch (ccdtype) {
+	case ZERO:
+	case DARK:
+	    call set_zero (ccd)
+	case FLAT:
+	    call set_zero (ccd)
+	    call set_dark (ccd)
+	    CORS(ccd, FINDMEAN) = YES
+	    CORS(ccd, MINREP) = YES
+	case ILLUM:
+	    call set_zero (ccd)
+	    call set_dark (ccd)
+	    call set_flat (ccd)
+	case OBJECT, COMP:
+	    call set_zero (ccd)
+	    call set_dark (ccd)
+	    call set_flat (ccd)
+	    call set_illum (ccd)
+	    call set_fringe (ccd)
+	default:
+	    call set_zero (ccd)
+	    call set_dark (ccd)
+	    call set_flat (ccd)
+	    call set_illum (ccd)
+	    call set_fringe (ccd)
+	    CORS(ccd, FINDMEAN) = YES
+	}
+
+	# Do the processing if the COR flag is set.
+	if (COR(ccd) == YES) {
+	    call doproc (ccd)
+	    call set_header (ccd)
+
+	    # Replace the input by the output image.
+	    call imunmap (in)
+	    call imunmap (out)
+	    iferr (call ccddelete (input)) {
+		call imdelete (Memc[output])
+		call error (1,
+		    "Can't delete or make backup of original image")
+	    }
+	    call imrename (Memc[output], input)
+	} else {
+	    # Delete the temporary output image leaving the input unchanged.
+	    call imunmap (in)
+	    iferr (call imunmap (out))
+		;
+	    iferr (call imdelete (Memc[output]))
+		;
+	}
+	call free_proc (ccd)
+
+	# Do special processing for calibration images.
+	switch (ccdtype) {
+	case ZERO:
+	    call readcor (input)
+	case FLAT:
+	    call ccdmean (input)
+	}
+
+	call sfree (sp)
+end