Initial commit

1 files changed, 145 insertions, 0 deletions

diff --git a/noao/onedspec/t_sflip.x b/noao/onedspec/t_sflip.x
new file mode 100644
index 00000000..b14e2ae0
--- /dev/null
+++ b/noao/onedspec/t_sflip.x
@@ -0,0 +1,145 @@
+include	<error.h>
+include	<imhdr.h>
+include	<smw.h>
+
+
+# SFLIP -- Flip data and/or coordinate system in spectra.
+
+procedure t_sflip ()
+
+pointer	inlist		# Input list
+pointer	outlist		# Output list
+bool	coord_flip	# Flip coordinates?
+bool	data_flip	# Flip data?
+
+bool	in_place
+int	i, j, k, n, axis
+pointer	sp, input, output, temp, a, b, c
+pointer	in, out, smw, mw, tmp, inbuf, outbuf
+
+bool	clgetb(), streq()
+int	imtgetim(), imtlen()
+pointer	imtopenp(), immap(), smw_openim(), imgl3r(), impl3r()
+errchk	immap, smw_openim
+
+begin
+	call smark (sp)
+	call salloc (input, SZ_FNAME, TY_CHAR)
+	call salloc (output, SZ_FNAME, TY_CHAR)
+	call salloc (temp, SZ_FNAME, TY_CHAR)
+	call salloc (a, 3*3, TY_DOUBLE)
+	call salloc (b, 3, TY_DOUBLE)
+	call salloc (c, 3, TY_DOUBLE)
+
+	# Get task parameters.
+	inlist = imtopenp ("input")
+	outlist = imtopenp ("output")
+	coord_flip = clgetb ("coord_flip")
+	data_flip = clgetb ("data_flip")
+
+	# Loop over all input images.
+	in = NULL
+	out = NULL
+	smw = NULL
+	while (imtgetim (inlist, Memc[input], SZ_FNAME) != EOF) {
+	    if (imtlen (outlist) > 0) {
+		if (imtgetim (outlist, Memc[output], SZ_FNAME) == EOF)
+		    break
+	    } else
+		call strcpy (Memc[input], Memc[output], SZ_FNAME)
+	    if (streq (Memc[input], Memc[output])) {
+		if (data_flip) {
+		    in_place = false
+		    call mktemp ("temp", Memc[temp], SZ_FNAME)
+		} else
+		    in_place = true
+	    } else {
+		in_place = false
+		call strcpy (Memc[output], Memc[temp], SZ_FNAME)
+	    }
+
+	    iferr {
+		# Map the images and WCS.
+		if (in_place) {
+		    tmp = immap (Memc[input], READ_WRITE, 0); in = tmp
+		    out = in
+		} else {
+		    tmp = immap (Memc[input], READ_ONLY, 0); in = tmp
+		    tmp = immap (Memc[temp], NEW_COPY, in); out = tmp
+		}
+		tmp = smw_openim (in); smw = tmp
+
+		# Flip coordinates.
+		if (coord_flip) {
+		    mw = SMW_MW(smw,0)
+		    n = SMW_PDIM(smw)
+		    axis = SMW_PAXIS(smw,1) - 1
+		    call mw_gltermd (mw, Memd[a], Memd[b], n)
+		    Memd[a+axis*(n+1)] = -Memd[a+axis*(n+1)]
+		    Memd[b+axis] = SMW_LLEN(smw,1) - Memd[b+axis] + 1
+		    call mw_sltermd (mw, Memd[a], Memd[b], n)
+		    call smw_saveim (smw, out)
+		}
+
+		# Flip data.
+		if (data_flip) {
+		    n = IM_LEN(in,1)
+		    do j = 1, IM_LEN(in,3) {
+			do i = 1, IM_LEN(in,2) {
+			    inbuf = imgl3r (in, i, j)
+			    switch (SMW_FORMAT(smw)) {
+			    case SMW_ND:
+				switch (SMW_LAXIS(smw,1)) {
+				case 1:
+				    outbuf = impl3r (out, i, j) + n - 1
+				    do k = 0, n-1
+					Memr[outbuf-k] = Memr[inbuf+k]
+				case 2:
+				    outbuf = impl3r (out, IM_LEN(in,2)-i+1, j)
+				    call amovr (Memr[inbuf], Memr[outbuf], n)
+				case 3:
+				    outbuf = impl3r (out, i, IM_LEN(in,3)-j+1)
+				    call amovr (Memr[inbuf], Memr[outbuf], n)
+				}
+			    case SMW_ES, SMW_MS:
+				outbuf = impl3r (out, i, j) + n - 1
+				do k = 0, n-1
+				    Memr[outbuf-k] = Memr[inbuf+k]
+			    }
+			}
+		    }
+		} else if (!in_place) {
+		    n = IM_LEN(in,1)
+		    do j = 1, IM_LEN(in,3) {
+			do i = 1, IM_LEN(in,2) {
+			    inbuf = imgl3r (in, i, j)
+			    outbuf = impl3r (out, i, j)
+			    call amovr (Memr[inbuf], Memr[outbuf], n)
+			}
+		    }
+		}
+	    } then {
+		if (!in_place && out != NULL) {
+		    call imunmap (out)
+		    call imdelete (Memc[temp])
+		}
+		call erract (EA_WARN)
+	    }
+
+	    if (smw != NULL)
+		call smw_close (smw)
+	    if (!in_place && out != NULL) {
+		call imunmap (out)
+		call imunmap (in)
+		if (streq (Memc[input], Memc[output])) {
+		    call imdelete (Memc[input])
+		    call imrename (Memc[temp], Memc[output])
+		}
+	    } else if (in != NULL)
+		call imunmap (in)
+	}
+
+	call imtclose (inlist)
+	call imtclose (outlist)
+	call sfree (sp)
+end