1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
|
include <imhdr.h>
define VPTR Memi[$1+$2-1] # Array of axis vector pointers
# IMJOIN -- Join the set of input images into an output image along the
# specified axis, any dimension up to 7 (NOT necessarily IM_MAXDIM!).
$for (silrdx)
procedure imjoin$t (inptr, nimages, out, joindim, outtype)
pointer inptr[nimages] # Input IMIO pointers
int nimages # Number of input images
pointer out # Output IMIO pointer
int joindim # Dimension along which to join images
int outtype # Output datatype
pointer in, inbuf, outbuf, sp, vin
int stat, image, cum_len, line, band, dim4, dim5, dim6, dim7
long vout[IM_MAXDIM]
pointer imgnl$t()
pointer impnl$t()
begin
call smark (sp)
call salloc (vin, nimages, TY_INT)
call amovkl (long(1), vout, IM_MAXDIM)
do image = 1, nimages {
call salloc (VPTR(vin,image), IM_MAXDIM, TY_LONG)
call amovkl (long(1), Meml[VPTR(vin,image)], IM_MAXDIM)
}
# Join input images along specified dimension. Joins along columns
# and lines require processing in special order, all others in the
# same order. In the first two cases we process all input images
# in inner loops, so we have to keep all those imdescriptors open.
switch (joindim) {
case 1: # join columns
do band = 1, IM_LEN(out,3)
do line = 1, IM_LEN(out,2) {
stat = impnl$t (out, outbuf, vout)
cum_len = 0
do image = 1, nimages {
in = inptr[image]
stat = imgnl$t (in, inbuf, Meml[VPTR(vin,image)])
call amov$t (Mem$t[inbuf], Mem$t[outbuf+cum_len],
IM_LEN(in,1))
cum_len = cum_len + IM_LEN(in,1)
}
}
case 2: # join lines
do band = 1, IM_LEN(out,3)
do image = 1, nimages {
in = inptr[image]
do line = 1, IM_LEN(in,2) {
stat = impnl$t (out, outbuf, vout)
stat = imgnl$t (in, inbuf, Meml[VPTR(vin,image)])
call amov$t (Mem$t[inbuf], Mem$t[outbuf], IM_LEN(in,1))
}
}
case 3,4,5,6,7: # join bands or higher
do image = 1, nimages {
in = inptr[image]
do dim7 = 1, IM_LEN(in,7)
do dim6 = 1, IM_LEN(in,6)
do dim5 = 1, IM_LEN(in,5)
do dim4 = 1, IM_LEN(in,4)
do band = 1, IM_LEN(in,3)
do line = 1, IM_LEN(in,2) {
stat = impnl$t (out, outbuf, vout)
stat = imgnl$t (in, inbuf,
Meml[VPTR(vin,image)])
call amov$t (Mem$t[inbuf],
Mem$t[outbuf], IM_LEN(in,1))
}
# Unmap last image to free resources.
call imunmap (in)
}
}
call sfree (sp)
end
$endfor
|
