Initial commit

1 files changed, 111 insertions, 0 deletions

diff --git a/noao/onedspec/identify/idgraph.x b/noao/onedspec/identify/idgraph.x
new file mode 100644
index 00000000..2c38efb4
--- /dev/null
+++ b/noao/onedspec/identify/idgraph.x
@@ -0,0 +1,111 @@
+include	<gset.h>
+include	<pkg/gtools.h>
+include	<smw.h>
+include	"identify.h"
+
+# ID_GRAPH -- Graph image vector in which features are to be identified.
+
+procedure id_graph (id, gtype)
+
+pointer	id		# ID pointer
+int	gtype		# Graph type
+
+begin
+	switch (gtype) {
+	case 1:
+	    call id_graph1 (id)
+	case 2:
+	    call id_graph2 (id)
+	default:
+	    call id_graph1 (id)
+	}
+end
+
+
+procedure id_graph1 (id)
+
+pointer	id				# ID pointer
+
+int	i, n
+real	xmin, xmax, ymin, ymax, dy, gt_getr()
+pointer	sh, x, y
+
+begin
+	sh = ID_SH(id)
+	call malloc (x, SN(sh), TY_REAL)
+	y = SY(sh)
+	n = SN(sh)
+
+	call achtdr (FITDATA(id,1), Memr[x], n)
+
+	call gclear (ID_GP(id))
+	xmin = min (Memr[x], Memr[x+n-1])
+	xmax = max (Memr[x], Memr[x+n-1])
+	ymin = gt_getr (ID_GT(id), GTXMIN)
+	ymax = gt_getr (ID_GT(id), GTXMAX)
+	if ((!IS_INDEF(ymin) && xmax<ymin) || (!IS_INDEF(ymax) && xmin>ymax)) {
+	    call gt_setr (ID_GT(id), GTXMIN, INDEF)
+	    call gt_setr (ID_GT(id), GTXMAX, INDEF)
+	}
+	call alimr (Memr[y], n, ymin, ymax)
+	dy = ymax - ymin
+	call gswind (ID_GP(id), xmin, xmax, ymin - .2 * dy, ymax + .2 * dy)
+	call gt_swind (ID_GP(id), ID_GT(id))
+	call gt_labax (ID_GP(id), ID_GT(id))
+	call gt_plot (ID_GP(id), ID_GT(id), Memr[x], Memr[y], n)
+
+	do i = 1, ID_NFEATURES(id)
+	    call id_mark (id, i)
+
+	call mfree (x, TY_REAL)
+end
+
+
+# ID_GRAPH2 -- Make review graph for current feature.
+
+procedure id_graph2 (id)
+
+pointer	id				# ID pointer
+
+int	i, j, k, n
+real	xmin, xmax, ymin, ymax, dy
+pointer	sh, x, y
+
+begin
+	sh = ID_SH(id)
+	call malloc (x, SN(sh), TY_REAL)
+	y = SY(sh)
+	n = SN(sh)
+
+	call achtdr (FITDATA(id,1), Memr[x], n)
+
+	xmin = real (FIT(id,ID_CURRENT(id))) - ID_ZWIDTH(id) / 2.
+	xmax = real (FIT(id,ID_CURRENT(id))) + ID_ZWIDTH(id) / 2.
+
+	i = 0
+	do k = 1, n {
+	    if ((Memr[x+k-1] < xmin) || (Memr[x+k-1] > xmax))
+		next
+	    if (i == 0)
+		i = k
+	    j = k
+	}
+	k = j - i + 1
+
+	call alimr (Memr[y+i-1], k, ymin, ymax)
+	dy = ymax - ymin
+
+	call gclear (ID_GP(id))
+	call gswind (ID_GP(id), xmin, xmax, ymin - .2 * dy, ymax + .2 * dy)
+#	if (ID_GT(id) != NULL) {
+#	    call gseti (ID_GP(id), G_XTRAN, GT_XTRAN(ID_GT(id)))
+#	    call gseti (ID_GP(id), G_YTRAN, GT_YTRAN(ID_GT(id)))
+#	}
+	call gt_labax (ID_GP(id), ID_GT(id))
+	call gt_plot (ID_GP(id), ID_GT(id), Memr[x], Memr[y], n)
+
+	do i = 1, ID_NFEATURES(id)
+	    call id_mark (id, i)
+
+	call mfree (x, TY_REAL)
+end