Initial commit

1 files changed, 55 insertions, 0 deletions

diff --git a/noao/digiphot/apphot/polyphot/apydraw.x b/noao/digiphot/apphot/polyphot/apydraw.x
new file mode 100644
index 00000000..aadcbed0
--- /dev/null
+++ b/noao/digiphot/apphot/polyphot/apydraw.x
@@ -0,0 +1,55 @@
+include "../lib/polyphot.h"
+
+# APYDRAW -- Procedure to draw polygons interactively on the  display.
+
+procedure ap_ydraw (py, im, cl, pl, ld, pd, id)
+
+pointer py			# pointer to apphot structure
+pointer	im			# theinput image descriptor
+int	cl			# coordinates file descriptor
+int	pl			# vertices list file descriptor
+int	ld			# coordinate list number
+int	pd			# polygon list number
+pointer	id			# pointer to image display stream
+
+int	req_num, prev_num, nvertices, delim
+pointer	sp, x, y
+int	ap_ynextobj()
+data	delim /';'/
+
+begin
+	# Allocate temporary space for arrays.
+	call smark (sp)
+	call salloc (x, MAX_NVERTICES + 1, TY_REAL)
+	call salloc (y, MAX_NVERTICES + 1, TY_REAL)
+
+	# Rewind the polygon and coordinate files.
+	if (pl != NULL)
+	    call seek (pl, BOF)
+	if (cl != NULL)
+	    call seek (cl, BOF)
+
+	# Initialize.
+	pd = 0
+	prev_num = 0
+	req_num = ld + 1
+	nvertices = ap_ynextobj (py, im, id, pl, cl, delim, Memr[x], Memr[y],
+	    MAX_NVERTICES, prev_num, req_num, ld, pd)
+
+	# Loop over the coordinate and polygon file.
+	while (nvertices != EOF) {
+	    
+	    # Draw the polygon if it has more than three vertices.
+	    if (id != NULL && nvertices >= 3)
+		call appymark (py, id, Memr[x], Memr[y], nvertices + 1,
+		    NO, NO, YES)
+
+	    # Setup for next polygon.
+	    prev_num = ld
+	    req_num = ld + 1
+	    nvertices = ap_ynextobj (py, im, id, pl, cl, delim, Memr[x],
+	        Memr[y], MAX_NVERTICES, prev_num, req_num, ld, pd)
+	}   
+
+	call sfree (sp)
+end