Initial commit

1 files changed, 159 insertions, 0 deletions

diff --git a/noao/nproto/ir/t_irmatch2d.x b/noao/nproto/ir/t_irmatch2d.x
new file mode 100644
index 00000000..38c6feb3
--- /dev/null
+++ b/noao/nproto/ir/t_irmatch2d.x
@@ -0,0 +1,159 @@
+include <imhdr.h>
+include <fset.h>
+include "iralign.h"
+
+# T_IRMATCHD2 -- Align the individual subraster elements in the input image.
+# In order to run this program the user should have created the output image
+# and the database file with the IRMOSAIC task. In addition the user should
+# supply a coordinate list consisting of pairs of coordinates of identical
+# objects or features in two adjacent subrasters.
+
+procedure t_irmatchd2 ()
+
+int	cl, interp, align, verbose, nimages, nmatch
+pointer	sp, inimage, outimage, database, coords, matchlist, trimlimits, ranges
+pointer	str, ir, im, outim, dt
+
+bool	clgetb()
+int	open(), clgwrd(), btoi(), ir_links(), ir_clinks(), ir_flinks()
+int	decode_ranges()
+pointer	immap(), dtmap()
+real	clgetr()
+
+begin
+	# Allocate space for the ir strucuture.
+	call ir_init (ir)
+
+	# Set the standard output to flush on a new line.
+	call fseti (STDOUT, F_FLUSHNL, YES)
+
+	# Allocate temporary working space.
+	call smark (sp)
+	call salloc (inimage, SZ_FNAME, TY_CHAR)
+	call salloc (outimage, SZ_FNAME, TY_CHAR)
+	call salloc (coords, SZ_FNAME, TY_CHAR)
+	call salloc (database, SZ_FNAME, TY_CHAR)
+	call salloc (matchlist, SZ_LINE, TY_CHAR)
+	call salloc (trimlimits, SZ_FNAME, TY_CHAR)
+	call salloc (ranges, 3 * MAX_NRANGES + 1, TY_INT)
+	call salloc (str, SZ_LINE, TY_CHAR)
+
+	# Get the input and output images and the coordinate list.
+	call clgstr ("input", Memc[inimage], SZ_FNAME)
+	call clgstr ("output", Memc[outimage], SZ_FNAME)
+	call clgstr ("database", Memc[database], SZ_FNAME)
+	call clgstr ("coords", Memc[coords], SZ_FNAME)
+	align = clgwrd ("alignment", Memc[str], SZ_LINE, ",coords,shifts,file,")
+	call clgstr ("match", Memc[matchlist], SZ_LINE)
+	call clgstr ("trimlimits", Memc[trimlimits], SZ_FNAME)
+
+	# Open the images and files.
+	im = immap (Memc[inimage], READ_ONLY, 0)
+	outim = immap (Memc[outimage], NEW_COPY, im)
+	dt = dtmap (Memc[database], READ_ONLY)
+
+	# Get the data base parameters.
+	call ir_dtrparams (dt, Memc[inimage], ir)
+
+	call ir_params (ir, im, outim)
+
+	interp = clgwrd ("interpolant", Memc[str], SZ_LINE,
+	    ",nearest,linear,poly3,poly5,spline3,")
+	verbose = btoi (clgetb ("verbose"))
+
+	# Decode the list of input images to be intensity matched.
+	nimages = IR_NXSUB(ir) * IR_NYSUB(ir)
+	if (Memc[matchlist] == EOS) {
+	    Memi[ranges] = NULL
+	} else if (Memc[matchlist] == '*') {
+	    Memi[ranges] = 1
+	    Memi[ranges+1] = nimages
+	    Memi[ranges+2] = 1
+	    Memi[ranges+3] = NULL
+	} else if (decode_ranges (Memc[matchlist], Memi[ranges], MAX_NRANGES,
+	    nmatch) == ERR) {
+		call error (0,
+		    "Cannot decode list of rasters to be intensity matched.")
+
+	}
+
+	# Allocate working space.
+	call ir_arrays (ir, nimages)
+
+	# Compute the shifts for each subraster.
+	switch (align) {
+	case IR_COORDS:
+	    cl = open (Memc[coords], READ_ONLY, TEXT_FILE)
+	    if (ir_links (cl, Memr[IR_XRSHIFTS(ir)], Memr[IR_YRSHIFTS(ir)],
+	        Memr[IR_XCSHIFTS(ir)], Memr[IR_YCSHIFTS(ir)],
+		Memi[IR_NRSHIFTS(ir)], Memi[IR_NCSHIFTS(ir)],
+		IR_NCOLS(ir), IR_NROWS(ir), IR_NXRSUB(ir), IR_NYRSUB(ir),
+		IR_NXSUB(ir), IR_NYSUB(ir), IR_NXOVERLAP(ir), IR_NYOVERLAP(ir),
+		IR_ORDER(ir)) > 0) {
+	        call ir_shifts (ir, im, outim, Memr[IR_XRSHIFTS(ir)],
+		    Memr[IR_YRSHIFTS(ir)], Memr[IR_XCSHIFTS(ir)],
+		    Memr[IR_YCSHIFTS(ir)], Memi[IR_IC1(ir)], Memi[IR_IC2(ir)],
+		    Memi[IR_IL1(ir)], Memi[IR_IL2(ir)], Memi[IR_OC1(ir)],
+		    Memi[IR_OC2(ir)], Memi[IR_OL1(ir)], Memi[IR_OL2(ir)],
+		    Memr[IR_DELTAX(ir)], Memr[IR_DELTAY(ir)])
+	        call ir_m2match (ir, im, Memi[ranges], Memi[IR_IC1(ir)],
+		    Memi[IR_IC2(ir)], Memi[IR_IL1(ir)], Memi[IR_IL2(ir)],
+		    Memr[IR_DELTAX(ir)], Memr[IR_DELTAY(ir)],
+		    Memr[IR_DELTAI(ir)])
+	    } else
+	        call error (0, "There are no legal shifts in the coords file.")
+	    call close (cl)
+
+	case IR_SHIFTS:
+	    if (ir_clinks (Memr[IR_XRSHIFTS(ir)], Memr[IR_YRSHIFTS(ir)],
+	        Memr[IR_XCSHIFTS(ir)], Memr[IR_YCSHIFTS(ir)], IR_NXRSUB(ir),
+		IR_NYRSUB(ir), IR_NXSUB(ir), IR_NYSUB(ir), clgetr ("xshift"),
+		clgetr ("yshift")) > 0) {
+	        call ir_shifts (ir, im, outim, Memr[IR_XRSHIFTS(ir)],
+		    Memr[IR_YRSHIFTS(ir)], Memr[IR_XCSHIFTS(ir)],
+		    Memr[IR_YCSHIFTS(ir)], Memi[IR_IC1(ir)], Memi[IR_IC2(ir)],
+		    Memi[IR_IL1(ir)], Memi[IR_IL2(ir)], Memi[IR_OC1(ir)],
+		    Memi[IR_OC2(ir)], Memi[IR_OL1(ir)], Memi[IR_OL2(ir)],
+		    Memr[IR_DELTAX(ir)], Memr[IR_DELTAY(ir)])
+	        call ir_m2match (ir, im, Memi[ranges], Memi[IR_IC1(ir)],
+		    Memi[IR_IC2(ir)], Memi[IR_IL1(ir)], Memi[IR_IL2(ir)],
+		    Memr[IR_DELTAX(ir)], Memr[IR_DELTAY(ir)],
+		    Memr[IR_DELTAI(ir)])
+	    } else
+	        call error (0, "There are no legal shifts in the coords file.")
+
+	case IR_FILE:
+	    cl = open (Memc[coords], READ_ONLY, TEXT_FILE)
+	    if (ir_flinks (cl, Memr[IR_DELTAX(ir)], Memr[IR_DELTAY(ir)],
+	        Memr[IR_DELTAI(ir)], nimages) >= nimages) {
+	        call ir_fshifts (ir, im, outim, Memr[IR_DELTAX(ir)],
+		    Memr[IR_DELTAY(ir)], Memi[IR_IC1(ir)], Memi[IR_IC2(ir)],
+		    Memi[IR_IL1(ir)], Memi[IR_IL2(ir)], Memi[IR_OC1(ir)],
+		    Memi[IR_OC2(ir)], Memi[IR_OL1(ir)], Memi[IR_OL2(ir)])
+	    } else
+	        call error (0, "There are fewer shifts than input subrasters.")
+	    call close (cl)
+
+	default:
+	    call error (0, "T_IRALIGN: Undefined alignment algorithm")
+	}
+
+
+	# Fill the output image with the unknown value.
+	call ir_imzero (outim, int (IM_LEN(outim,1)), int (IM_LEN(outim, 2)),
+	    IR_OVAL(ir))
+
+	# Shift and match all the subrasters.
+	call ir_subalign (ir, im, outim, Memc[trimlimits], Memi[IR_IC1(ir)],
+	    Memi[IR_IC2(ir)], Memi[IR_IL1(ir)], Memi[IR_IL2(ir)],
+	    Memi[IR_OC1(ir)], Memi[IR_OC2(ir)], Memi[IR_OL1(ir)],
+	    Memi[IR_OL2(ir)], Memr[IR_DELTAX(ir)], Memr[IR_DELTAY(ir)], 
+	    Memr[IR_DELTAI(ir)], YES, interp, verbose)
+
+	# Close up files.
+	call imunmap (im)
+	call imunmap (outim)
+	call dtunmap (dt)
+	call sfree (sp)
+	call ir_free (ir)
+end