Initial commit

1 files changed, 159 insertions, 0 deletions

diff --git a/pkg/images/imutil/src/generic/imanl.x b/pkg/images/imutil/src/generic/imanl.x
new file mode 100644
index 00000000..8ec958c4
--- /dev/null
+++ b/pkg/images/imutil/src/generic/imanl.x
@@ -0,0 +1,159 @@
+# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
+
+include	<imhdr.h>
+
+# IMA_NL -- For each line in the output image lines from the input images
+# are returned.  The input images are repeated as necessary.  EOF is returned
+# when the last line of the output image has been reached.  One dimensional
+# images are read only once and the data pointers are assumed to be unchanged
+# from previous calls.  The image line vectors must be initialized externally
+# and then left untouched.
+# 
+# This procedure is typically used when operations upon lines or pixels
+# make sense in mixed dimensioned images.  For example to add a one dimensional
+# image to all lines of a higher dimensional image or to subtract a
+# two dimensional image from all bands of three dimensional image.
+# The lengths of the common dimensions should generally be checked
+# for equality with xt_imleneq.
+
+
+int procedure ima_nls (im, data, v, nimages)
+
+pointer	im[nimages]		# IMIO pointers; the first one is the output
+pointer	data[nimages]		# Returned data pointers
+long	v[IM_MAXDIM, nimages]	# Line vectors
+int	nimages			# Number of images
+
+int	i
+
+int	impnls(), imgnls()
+
+begin
+	if (impnls (im[1], data[1], v[1,1]) == EOF)
+	    return (EOF)
+
+	for (i=2; i <= nimages; i=i+1) {
+	    if (imgnls (im[i], data[i], v[1,i]) == EOF) {
+	        if (IM_NDIM(im[i]) > 1) {
+	            call amovkl (long(1), v[1,i], IM_MAXDIM)
+                    if (imgnls (im[i], data[i], v[1,i]) == EOF)
+			call error (0, "Error reading image line")
+	        }
+	    }
+	}
+
+	return (OK)
+end
+
+int procedure ima_nli (im, data, v, nimages)
+
+pointer	im[nimages]		# IMIO pointers; the first one is the output
+pointer	data[nimages]		# Returned data pointers
+long	v[IM_MAXDIM, nimages]	# Line vectors
+int	nimages			# Number of images
+
+int	i
+
+int	impnli(), imgnli()
+
+begin
+	if (impnli (im[1], data[1], v[1,1]) == EOF)
+	    return (EOF)
+
+	for (i=2; i <= nimages; i=i+1) {
+	    if (imgnli (im[i], data[i], v[1,i]) == EOF) {
+	        if (IM_NDIM(im[i]) > 1) {
+	            call amovkl (long(1), v[1,i], IM_MAXDIM)
+                    if (imgnli (im[i], data[i], v[1,i]) == EOF)
+			call error (0, "Error reading image line")
+	        }
+	    }
+	}
+
+	return (OK)
+end
+
+int procedure ima_nll (im, data, v, nimages)
+
+pointer	im[nimages]		# IMIO pointers; the first one is the output
+pointer	data[nimages]		# Returned data pointers
+long	v[IM_MAXDIM, nimages]	# Line vectors
+int	nimages			# Number of images
+
+int	i
+
+int	impnll(), imgnll()
+
+begin
+	if (impnll (im[1], data[1], v[1,1]) == EOF)
+	    return (EOF)
+
+	for (i=2; i <= nimages; i=i+1) {
+	    if (imgnll (im[i], data[i], v[1,i]) == EOF) {
+	        if (IM_NDIM(im[i]) > 1) {
+	            call amovkl (long(1), v[1,i], IM_MAXDIM)
+                    if (imgnll (im[i], data[i], v[1,i]) == EOF)
+			call error (0, "Error reading image line")
+	        }
+	    }
+	}
+
+	return (OK)
+end
+
+int procedure ima_nlr (im, data, v, nimages)
+
+pointer	im[nimages]		# IMIO pointers; the first one is the output
+pointer	data[nimages]		# Returned data pointers
+long	v[IM_MAXDIM, nimages]	# Line vectors
+int	nimages			# Number of images
+
+int	i
+
+int	impnlr(), imgnlr()
+
+begin
+	if (impnlr (im[1], data[1], v[1,1]) == EOF)
+	    return (EOF)
+
+	for (i=2; i <= nimages; i=i+1) {
+	    if (imgnlr (im[i], data[i], v[1,i]) == EOF) {
+	        if (IM_NDIM(im[i]) > 1) {
+	            call amovkl (long(1), v[1,i], IM_MAXDIM)
+                    if (imgnlr (im[i], data[i], v[1,i]) == EOF)
+			call error (0, "Error reading image line")
+	        }
+	    }
+	}
+
+	return (OK)
+end
+
+int procedure ima_nld (im, data, v, nimages)
+
+pointer	im[nimages]		# IMIO pointers; the first one is the output
+pointer	data[nimages]		# Returned data pointers
+long	v[IM_MAXDIM, nimages]	# Line vectors
+int	nimages			# Number of images
+
+int	i
+
+int	impnld(), imgnld()
+
+begin
+	if (impnld (im[1], data[1], v[1,1]) == EOF)
+	    return (EOF)
+
+	for (i=2; i <= nimages; i=i+1) {
+	    if (imgnld (im[i], data[i], v[1,i]) == EOF) {
+	        if (IM_NDIM(im[i]) > 1) {
+	            call amovkl (long(1), v[1,i], IM_MAXDIM)
+                    if (imgnld (im[i], data[i], v[1,i]) == EOF)
+			call error (0, "Error reading image line")
+	        }
+	    }
+	}
+
+	return (OK)
+end
+