Initial commit

1 files changed, 74 insertions, 0 deletions

diff --git a/noao/onedspec/ecidentify/ecgdata.x b/noao/onedspec/ecidentify/ecgdata.x
new file mode 100644
index 00000000..1087d38c
--- /dev/null
+++ b/noao/onedspec/ecidentify/ecgdata.x
@@ -0,0 +1,74 @@
+include	<imhdr.h>
+include	<imio.h>
+include	<pkg/gtools.h>
+include	<smw.h>
+include	<units.h>
+include	"ecidentify.h"
+
+# EC_GDATA -- Get image data.
+
+procedure ec_gdata (ec)
+
+pointer	ec				# ID pointer
+
+int	i, j
+pointer	im, mw, sh, sp, str1, str2
+
+double	smw_c1trand()
+pointer	immap(), smw_openim(), smw_sctran()
+errchk	immap, smw_openim, shdr_open
+
+begin
+	# Map the image.
+	im = immap (Memc[EC_IMAGE(ec)], READ_ONLY, 0)
+
+	# Free previous data
+	do i = 1, EC_NLINES(ec)
+	    call shdr_close (SH(ec,i))
+	call mfree (EC_SHS(ec), TY_POINTER)
+	call mfree (EC_PIXDATA(ec), TY_DOUBLE)
+
+	# Set MWCS
+	mw = smw_openim (im)
+	EC_LP(ec) = smw_sctran (mw, "logical", "physical", 1)
+	EC_PL(ec) = smw_sctran (mw, "physical", "logical", 1)
+
+	# Allocate new vectors.
+	EC_NCOLS(ec) = IM_LEN(im, 1)
+	EC_NLINES(ec) = IM_LEN(im, 2)
+	call calloc (EC_SHS(ec), EC_NLINES(ec), TY_POINTER)
+	call malloc (EC_PIXDATA(ec), EC_NCOLS(ec)*EC_NLINES(ec), TY_DOUBLE)
+
+	# Set the coordinates.
+	sh = NULL
+	do j = 1, EC_NLINES(ec) {
+	    call shdr_open (im, mw, j, 1, INDEFI, SHDATA, sh)
+	    if (EC_UN(ec) != NULL)
+		iferr (call shdr_units (sh, UN_UNITS(EC_UN(ec))))
+		    ;
+	    if (j != EC_NLINES(ec))
+		call shdr_copy (sh, SH(ec,j), NO)
+	    else
+		SH(ec,j) = sh
+	    call ec_gline (ec, j)
+	    do i = 1, EC_NPTS(ec)
+	        PIXDATA(ec,i) = smw_c1trand (EC_LP(ec), double(i))
+	}
+	EC_LINE(ec) = 1
+	call ec_gline (ec, EC_LINE(ec))
+	EC_AP(ec) = APS(ec,EC_LINE(ec))
+	EC_ORDER(ec) = ORDERS(ec,EC_LINE(ec))
+
+	# Set graph title.
+	call smark (sp)
+	call salloc (str1, SZ_LINE, TY_CHAR)
+	call salloc (str2, SZ_LINE, TY_CHAR)
+
+	call sprintf (Memc[str1], SZ_LINE, "ecidentify %s: %s")
+	    call pargstr (Memc[EC_IMAGE(ec)])
+	    call pargstr (IM_TITLE(im))
+	call gt_sets (EC_GT(ec), GTTITLE, Memc[str1])
+
+	call imunmap (im)
+	call sfree (sp)
+end