1 files changed, 153 insertions, 0 deletions
diff --git a/pkg/plot/crtpict/drawgraph.x b/pkg/plot/crtpict/drawgraph.x
new file mode 100644
index 00000000..5ee94045
--- /dev/null
+++ b/pkg/plot/crtpict/drawgraph.x
@@ -0,0 +1,153 @@
+# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
+
+include	<gset.h>
+include	<imhdr.h>
+include	<time.h>
+include "wdes.h"
+include	"crtpict.h"
+
+
+# CRT_DRAW_GRAPHICS -- Draw histogram plots and id information at the bottom of
+# the output print.
+
+procedure crt_draw_graphics (gp, im, cl, wdes)
+
+pointer	gp	# Pointer to graphics descriptor
+pointer	im	# Pointer to input image
+pointer	cl	# Pointer to cl parameter structure
+pointer	wdes	# Pointer to window descriptor
+
+pointer	w1, w0
+char	text[SZ_LINE],  system_id[SZ_LINE]
+int	ndev_rows, ndev_cols
+real	ndc_xs, ndc_xe, ndc_ys, ndc_ye 	 	# Graphics (NDC) viewport
+real	h_xs, h_xe, h_ys, h_ye			# Histogram viewport
+real	tx1_xc, tx1_ly, tx2_xs, tx2_ly, tx3_xs, tx3_ly, tx4_xs, tx4_ly
+real	px1, px2, py1, py2, pxcenter, pycenter, yres
+real	vx1, vx2, vy1, vy2, tx_start, xrf, yrf
+
+int	junk
+pointer	sp, buf
+int	envfind(), envputs()
+
+int	strlen(), ggeti()
+errchk	strlen, ggeti, crt_plot_hgrams, ggwind, ggview
+errchk	gtext
+
+begin
+	ndev_rows = ggeti (gp, "yr")
+	ndev_cols = ggeti (gp, "xr")
+	yres = (CRT_YE - CRT_YS) * ndev_rows
+
+	w0 = W_WC (wdes, 0)
+	w1 = W_WC (wdes, 1)
+
+	# The (NDC) device coordinates of the entire graphics viewport:
+	ndc_xs = CRT_XS 
+	ndc_xe = CRT_XE
+	ndc_ys = CRT_YS
+	ndc_ye = CRT_YS + real (yres * GRAPHICS_FRACTION(cl) / ndev_rows) 
+
+	# Working up from the bottom of the print, locations of various
+	# sections of the graphics are are calculated.  String TEXT1 will
+	# be centered in x at tx1_xc and have lower y coordinate tx1_ly:
+
+	tx1_xc = (ndc_xe + ndc_xs) / 2.0
+	tx1_ly = ndc_ys
+
+	# The three histograms occupy the space calculated next.  This
+	# space is broken into individual plots in a separate procedure.
+
+	h_xs = ndc_xs
+	h_xe = ndc_xe
+	h_ys = ndc_ys + ((ndc_ye - ndc_ys) * (TEXT1 + SPACE))
+	h_ye = h_ys + ((ndc_ye - ndc_ys) * HGRAMS)
+
+	# The left starting position of the text strings is calculated to
+	# line up with the leftmost histogram window:
+	tx_start = h_xs + ((h_xe - h_xs) / 6.0) - ((h_xe - h_xs) / 8.0)
+
+	# String TEXT2 has the following starting_x and lower_y coordinates:
+	tx2_xs = tx_start
+	tx2_ly = h_ye + ((ndc_ye - ndc_ys) * SPACE)
+
+	# String TEXT3 has these starting_x and lower_y coordinates:
+	tx3_xs = tx_start
+	tx3_ly = ndc_ys + ((TEXT1 + HGRAMS + TEXT2 + SPACE) * (ndc_ye - ndc_ys))
+
+	# String TEXT4 has these starting_x and lower_y coordinates:
+	tx4_xs = tx_start
+	tx4_ly = ndc_ys + ((TEXT1+HGRAMS+TEXT2+TEXT3+SPACE) * (ndc_ye - ndc_ys))
+
+	# Draw 3 plots describing transformation of image
+	call crt_plot_histograms (gp, cl, im, wdes, h_xs, h_xe, h_ys, h_ye)
+
+	# Set graphics WCS to WCS 0 for text plotting
+	call gseti (gp, G_WCS, 0)
+
+	# Text line 3 has the image filename and title string.
+	call sprintf (text, SZ_LINE, "%s: %s")
+	    call pargstr (W_IMSECT(wdes))
+	    call pargstr (IM_TITLE(im))
+	call gtext (gp, tx3_xs, tx3_ly, text, "s=0.5")
+
+	# Text line 2 contains image and transformation information; it
+	# is necessary to change to WCS_2 to retrieve the information:
+
+	call gseti (gp, G_WCS, 2)
+	call ggwind (gp, px1, px2, py1, py2)
+	call ggview (gp, vx1, vx2, vy1, vy2)
+	call gseti (gp, G_WCS, 0)
+
+	pxcenter = (px1 + px2) / 2.0
+	pycenter = (py1 + py2) / 2.0
+	xrf = ((vx2 * ndev_cols) - (vx1 * ndev_cols)) / (px2 - px1 + 1.0)
+	yrf = ((vy2 * ndev_rows) - (vy1 * ndev_rows)) / (py2 - py1 + 1.0)
+
+	call sprintf (text, SZ_LINE,
+	    "ncols=%d nrows=%d zmin=%g zmax=%g xc=%0.2f yc=%0.2f")  
+	        call pargi (IM_LEN(im,1))
+	        call pargi (IM_LEN(im,2))
+	        call pargr (IM_MIN(im))
+	        call pargr (IM_MAX(im))
+	        call pargr (pxcenter)
+	        call pargr (pycenter)
+	call sprintf (text[strlen(text)+1], SZ_LINE, " x_rep=%.2f y_rep=%.2f")
+	   call pargr (xrf)
+	   call pargr (yrf)
+	call gtext (gp, tx2_xs, tx2_ly, text, "s=0.35")
+
+	# Text line 1 gives the time and date the output was written
+	call sysid (system_id, SZ_LINE)
+	call gtext (gp, tx1_xc, tx1_ly, system_id, "h=c;s=0.45")
+
+	# Also output transformation information to STDOUT
+	call printf ("ncols=%d nrows=%d zmin=%g zmax=%g xc=%.2f yc=%.2f")
+	    call pargi (IM_LEN(im,1))
+	    call pargi (IM_LEN(im,2))
+	    call pargr (IM_MIN(im))
+	    call pargr (IM_MAX(im))
+	    call pargr (pxcenter)
+	    call pargr (pycenter)
+	call printf (" xrf=%.2f yrf=%.2f\n")
+	    call pargr (xrf)
+	    call pargr (yrf)
+
+	call printf ("%s \n")
+	    call pargstr (system_id)
+
+	# The following was added 17Dec85 at the request of the photo lab.
+	# It allows the negative to be identified easily by user name in
+	# addition to the sequence number written by the 11/23 program.
+
+	call smark (sp)
+	call salloc (buf, SZ_LINE, TY_CHAR)
+
+	if (envfind ("userid", Memc[buf], SZ_LINE) <= 0) {
+	    call getuid (Memc[buf], SZ_LINE)
+	    junk = envputs ("userid", Memc[buf])
+	}
+	call gtext (gp, CRT_XE, 0.001, Memc[buf], "h=r;v=b;s=1.2")
+
+	call sfree (sp)
+end