1 files changed, 29 insertions, 0 deletions
diff --git a/sys/vops/ak/agltl.x b/sys/vops/ak/agltl.x
new file mode 100644
index 00000000..3a416d37
--- /dev/null
+++ b/sys/vops/ak/agltl.x
@@ -0,0 +1,29 @@
+# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
+
+# AGLT -- Given a list of ranges, replace the value of each input pixel
+# which falls within a given range by applying the corresponding linear
+# transformation (b = a * kmul + kadd).  If KMUL is identically zero,
+# B is replaced by the constant KADD.
+
+procedure agltl (a, b, npix, low, high, kmul, kadd, nrange)
+
+long	a[ARB], b[ARB], pixval
+int	npix, i
+long	low[nrange], high[nrange]	# range limits
+double	kmul[nrange], kadd[nrange]	# linear transformation
+int	nrange, nr
+
+begin
+	do i = 1, npix {
+	    pixval = a[i]
+	    b[i] = pixval
+	    do nr = 1, nrange
+		if (pixval >= low[nr] && pixval <= high[nr]) {
+			if (kmul[nr] == 0.0D0)
+			    b[i] = kadd[nr]
+			else
+			    b[i] = (pixval * kmul[nr]) + kadd[nr]
+			break
+		}
+	}
+end