Initial commit

1 files changed, 250 insertions, 0 deletions

diff --git a/pkg/obsolete/generic/fixcol.x b/pkg/obsolete/generic/fixcol.x
new file mode 100644
index 00000000..ef64694a
--- /dev/null
+++ b/pkg/obsolete/generic/fixcol.x
@@ -0,0 +1,250 @@
+include	<imhdr.h>
+include	<imset.h>
+
+
+
+# FIXCOL -- Linearly interpolate columns across a region.
+
+procedure fixcols (image, x1, x2, y1, y2)
+
+pointer	image				# Image pointer
+int	x1, x2, y1, y2			# Region to be fixed
+
+int	i, j, nx, ny
+real	f1, f2, scale
+pointer	a, b, c
+pointer	imgs2s(), imps2s()
+
+begin
+	c = imps2s (image, x1, x2, y1, y2)
+
+	nx = x2 - x1 + 1
+	ny = y2 - y1 + 1
+	if (x1 == 1) {
+	    a = imgs2s (image, x2 + 1, x2 + 1, y1, y2)
+	    do i = 1, ny
+		call amovks (Mems[a + i - 1], Mems[c + (i - 1) * nx], nx)
+	} else if (x2 == IM_LEN (image, 1)) {
+	    a = imgs2s (image, x1 - 1, x1 - 1,  y1, y2)
+	    do i = 1, ny
+		call amovks (Mems[a + i - 1], Mems[c + (i - 1) * nx], nx)
+	} else {
+	    call imseti (image, IM_NBUFS, 2)
+	    a = imgs2s (image, x1 - 1, x1 - 1, y1, y2)
+	    b = imgs2s (image, x2 + 1, x2 + 1, y1, y2)
+	    do i = 1, ny {
+		f1 = Mems[a + i - 1]
+		f2 = Mems[b + i - 1]
+		scale = (f2 - f1) / (nx + 1)
+		do j = 1, nx
+		    Mems[c + (i - 1) * nx + j - 1] = j * scale + f1
+	    }
+	}
+end
+
+
+
+# FIXCOL -- Linearly interpolate columns across a region.
+
+procedure fixcoli (image, x1, x2, y1, y2)
+
+pointer	image				# Image pointer
+int	x1, x2, y1, y2			# Region to be fixed
+
+int	i, j, nx, ny
+real	f1, f2, scale
+pointer	a, b, c
+pointer	imgs2i(), imps2i()
+
+begin
+	c = imps2i (image, x1, x2, y1, y2)
+
+	nx = x2 - x1 + 1
+	ny = y2 - y1 + 1
+	if (x1 == 1) {
+	    a = imgs2i (image, x2 + 1, x2 + 1, y1, y2)
+	    do i = 1, ny
+		call amovki (Memi[a + i - 1], Memi[c + (i - 1) * nx], nx)
+	} else if (x2 == IM_LEN (image, 1)) {
+	    a = imgs2i (image, x1 - 1, x1 - 1,  y1, y2)
+	    do i = 1, ny
+		call amovki (Memi[a + i - 1], Memi[c + (i - 1) * nx], nx)
+	} else {
+	    call imseti (image, IM_NBUFS, 2)
+	    a = imgs2i (image, x1 - 1, x1 - 1, y1, y2)
+	    b = imgs2i (image, x2 + 1, x2 + 1, y1, y2)
+	    do i = 1, ny {
+		f1 = Memi[a + i - 1]
+		f2 = Memi[b + i - 1]
+		scale = (f2 - f1) / (nx + 1)
+		do j = 1, nx
+		    Memi[c + (i - 1) * nx + j - 1] = j * scale + f1
+	    }
+	}
+end
+
+
+
+# FIXCOL -- Linearly interpolate columns across a region.
+
+procedure fixcoll (image, x1, x2, y1, y2)
+
+pointer	image				# Image pointer
+int	x1, x2, y1, y2			# Region to be fixed
+
+int	i, j, nx, ny
+real	f1, f2, scale
+pointer	a, b, c
+pointer	imgs2l(), imps2l()
+
+begin
+	c = imps2l (image, x1, x2, y1, y2)
+
+	nx = x2 - x1 + 1
+	ny = y2 - y1 + 1
+	if (x1 == 1) {
+	    a = imgs2l (image, x2 + 1, x2 + 1, y1, y2)
+	    do i = 1, ny
+		call amovkl (Meml[a + i - 1], Meml[c + (i - 1) * nx], nx)
+	} else if (x2 == IM_LEN (image, 1)) {
+	    a = imgs2l (image, x1 - 1, x1 - 1,  y1, y2)
+	    do i = 1, ny
+		call amovkl (Meml[a + i - 1], Meml[c + (i - 1) * nx], nx)
+	} else {
+	    call imseti (image, IM_NBUFS, 2)
+	    a = imgs2l (image, x1 - 1, x1 - 1, y1, y2)
+	    b = imgs2l (image, x2 + 1, x2 + 1, y1, y2)
+	    do i = 1, ny {
+		f1 = Meml[a + i - 1]
+		f2 = Meml[b + i - 1]
+		scale = (f2 - f1) / (nx + 1)
+		do j = 1, nx
+		    Meml[c + (i - 1) * nx + j - 1] = j * scale + f1
+	    }
+	}
+end
+
+
+
+# FIXCOL -- Linearly interpolate columns across a region.
+
+procedure fixcolr (image, x1, x2, y1, y2)
+
+pointer	image				# Image pointer
+int	x1, x2, y1, y2			# Region to be fixed
+
+int	i, j, nx, ny
+real	f1, f2, scale
+pointer	a, b, c
+pointer	imgs2r(), imps2r()
+
+begin
+	c = imps2r (image, x1, x2, y1, y2)
+
+	nx = x2 - x1 + 1
+	ny = y2 - y1 + 1
+	if (x1 == 1) {
+	    a = imgs2r (image, x2 + 1, x2 + 1, y1, y2)
+	    do i = 1, ny
+		call amovkr (Memr[a + i - 1], Memr[c + (i - 1) * nx], nx)
+	} else if (x2 == IM_LEN (image, 1)) {
+	    a = imgs2r (image, x1 - 1, x1 - 1,  y1, y2)
+	    do i = 1, ny
+		call amovkr (Memr[a + i - 1], Memr[c + (i - 1) * nx], nx)
+	} else {
+	    call imseti (image, IM_NBUFS, 2)
+	    a = imgs2r (image, x1 - 1, x1 - 1, y1, y2)
+	    b = imgs2r (image, x2 + 1, x2 + 1, y1, y2)
+	    do i = 1, ny {
+		f1 = Memr[a + i - 1]
+		f2 = Memr[b + i - 1]
+		scale = (f2 - f1) / (nx + 1)
+		do j = 1, nx
+		    Memr[c + (i - 1) * nx + j - 1] = j * scale + f1
+	    }
+	}
+end
+
+
+
+# FIXCOL -- Linearly interpolate columns across a region.
+
+procedure fixcold (image, x1, x2, y1, y2)
+
+pointer	image				# Image pointer
+int	x1, x2, y1, y2			# Region to be fixed
+
+int	i, j, nx, ny
+real	f1, f2, scale
+pointer	a, b, c
+pointer	imgs2d(), imps2d()
+
+begin
+	c = imps2d (image, x1, x2, y1, y2)
+
+	nx = x2 - x1 + 1
+	ny = y2 - y1 + 1
+	if (x1 == 1) {
+	    a = imgs2d (image, x2 + 1, x2 + 1, y1, y2)
+	    do i = 1, ny
+		call amovkd (Memd[a + i - 1], Memd[c + (i - 1) * nx], nx)
+	} else if (x2 == IM_LEN (image, 1)) {
+	    a = imgs2d (image, x1 - 1, x1 - 1,  y1, y2)
+	    do i = 1, ny
+		call amovkd (Memd[a + i - 1], Memd[c + (i - 1) * nx], nx)
+	} else {
+	    call imseti (image, IM_NBUFS, 2)
+	    a = imgs2d (image, x1 - 1, x1 - 1, y1, y2)
+	    b = imgs2d (image, x2 + 1, x2 + 1, y1, y2)
+	    do i = 1, ny {
+		f1 = Memd[a + i - 1]
+		f2 = Memd[b + i - 1]
+		scale = (f2 - f1) / (nx + 1)
+		do j = 1, nx
+		    Memd[c + (i - 1) * nx + j - 1] = j * scale + f1
+	    }
+	}
+end
+
+
+
+# FIXCOL -- Linearly interpolate columns across a region.
+
+procedure fixcolx (image, x1, x2, y1, y2)
+
+pointer	image				# Image pointer
+int	x1, x2, y1, y2			# Region to be fixed
+
+int	i, j, nx, ny
+real	f1, f2, scale
+pointer	a, b, c
+pointer	imgs2x(), imps2x()
+
+begin
+	c = imps2x (image, x1, x2, y1, y2)
+
+	nx = x2 - x1 + 1
+	ny = y2 - y1 + 1
+	if (x1 == 1) {
+	    a = imgs2x (image, x2 + 1, x2 + 1, y1, y2)
+	    do i = 1, ny
+		call amovkx (Memx[a + i - 1], Memx[c + (i - 1) * nx], nx)
+	} else if (x2 == IM_LEN (image, 1)) {
+	    a = imgs2x (image, x1 - 1, x1 - 1,  y1, y2)
+	    do i = 1, ny
+		call amovkx (Memx[a + i - 1], Memx[c + (i - 1) * nx], nx)
+	} else {
+	    call imseti (image, IM_NBUFS, 2)
+	    a = imgs2x (image, x1 - 1, x1 - 1, y1, y2)
+	    b = imgs2x (image, x2 + 1, x2 + 1, y1, y2)
+	    do i = 1, ny {
+		f1 = Memx[a + i - 1]
+		f2 = Memx[b + i - 1]
+		scale = (f2 - f1) / (nx + 1)
+		do j = 1, nx
+		    Memx[c + (i - 1) * nx + j - 1] = j * scale + f1
+	    }
+	}
+end
+
+