From fa080de7afc95aa1c19a6e6fc0e0708ced2eadc4 Mon Sep 17 00:00:00 2001
From: Joseph Hunkeler <jhunkeler@gmail.com>
Date: Wed, 8 Jul 2015 20:46:52 -0400
Subject: Initial commit

---
 sys/plio/plloadim.x | 99 +++++++++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 99 insertions(+)
 create mode 100644 sys/plio/plloadim.x

(limited to 'sys/plio/plloadim.x')

diff --git a/sys/plio/plloadim.x b/sys/plio/plloadim.x
new file mode 100644
index 00000000..6fc45adc
--- /dev/null
+++ b/sys/plio/plloadim.x
@@ -0,0 +1,99 @@
+# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
+
+include	<imhdr.h>
+include	<imset.h>
+include	<plset.h>
+include	<plio.h>
+
+# PL_LOADIM -- Load a mask stored as a conventional image, i.e., convert an
+# image to a mask.
+
+procedure pl_loadim (pl, imname, title, maxch)
+
+pointer	pl			#I mask descriptor
+char	imname[ARB]		#I image name or section
+char	title[ARB]		#O mask title
+int	maxch			#I max chars out
+
+bool	sampling
+pointer	im, px, im_pl, bp
+long	vs_l[PL_MAXDIM], vs_p[PL_MAXDIM]
+long	ve_l[PL_MAXDIM], ve_p[PL_MAXDIM]
+int	npix, naxes, maxdim, maxval, depth, sz_buf, i
+long	v_in[PL_MAXDIM], v_out[PL_MAXDIM], vn[PL_MAXDIM]
+
+pointer	immap()
+int	imgnli(), imstati(), im_pmsvhdr()
+errchk	immap, imgnli
+
+begin
+	# Open the input image.
+	im = immap (imname, READ_ONLY, 0)
+
+	# Encode and output the image header.
+	bp = NULL;  i = im_pmsvhdr (im, bp, sz_buf)
+	call strcpy (Memc[bp], title, maxch)
+	call mfree (bp, TY_CHAR)
+
+	# Determine the mask depth in bits.
+	maxval = IM_MAX(im)
+	if (maxval <= 0)
+	    maxval = I_PVMAX
+	do i = 0, ARB
+	    if (2**i > min (I_PVMAX, maxval)) {
+		depth = i
+		break
+	    }
+
+	# Initialize the mask to the size of the image section.
+	npix = IM_LEN(im,1)
+	naxes = IM_NDIM(im)
+	maxdim = min (IM_MAXDIM, PL_MAXDIM)
+	call amovl (IM_LEN(im,1), vn, maxdim)
+	call pl_ssize (pl, naxes, vn, depth)
+
+	# If the image is already a mask internally, check whether any
+	# subsampling, axis flipping, or axis mapping is in effect.
+	# If so we can't use PLIO to copy the mask section.
+
+	im_pl = imstati (im, IM_PLDES)
+	sampling = false
+
+	if (im_pl != NULL) {
+	    call amovkl (long(1), vs_l, maxdim)
+	    call amovl (IM_LEN(im,1), ve_l, maxdim)
+	    call imaplv (im, vs_l, vs_p, maxdim)
+	    call imaplv (im, ve_l, ve_p, maxdim)
+
+	    do i = 1, maxdim {
+		vn[i] = ve_l[i] - vs_l[i] + 1
+		if (vn[i] != ve_p[i] - vs_p[i] + 1) {
+		    sampling = true
+		    break
+		}
+	    }
+	}
+
+	# If the source image is already a mask internally and no axis
+	# geometry is in effect in the image section (if any), then we can
+	# use a PLIO rasterop to efficiently copy the mask subsection.
+
+	if (im_pl != NULL && !sampling) {
+	    # Copy a mask subsection (or entire image if no section).
+	    call pl_rop (im_pl, vs_p, pl, vs_l, vn, PIX_SRC)
+
+	} else {
+	    # Copy image pixels.  Initialize the vector loop indices.
+	    call amovkl (long(1), v_in, maxdim)
+	    call amovkl (long(1), v_out, maxdim)
+
+	    # Copy the image.
+	    while (imgnli (im, px, v_in) != EOF) {
+		call pl_plpi (pl, v_out, Memi[px], 0, npix, PIX_SRC)
+		call amovl (v_in, v_out, maxdim)
+	    }
+	}
+
+	call pl_compress (pl)
+	call imunmap (im)
+end
-- 
cgit