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
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
|
# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
include <imhdr.h>
include <pmset.h>
include <plio.h>
include "../mio.h"
# MIO_PLSEG -- Put a line segment to a masked image. A line segment is a
# region of the data image for which the corresponding region of the mask has
# the constant nonzero value MVAL. Line segments are returned for each line in
# the region VS to VE, returning the number of pixels in each line segment as
# the function value, or EOF when the region is exhausted. Once EOF is
# reached, repeated calls will continue to return EOF until the next call to
# MIO_SETRANGE. Repeated calls to MIO_SETRANGE may be used to access a series
# of distinct regions in the image. If a subregion of the image is being
# accessed with MIO_SETRANGE, the vector coordinates V returned below will
# be relative to the defined subregion (if this is not what is desired,
# the range should be set to the full image and a region mask used to mask
# off the subregion to be accessed).
int procedure mio_plseg$t (mp, ptr, mval, v, npix)
pointer mp #I MIO descriptor
pointer ptr #O pointer to a buffer containing the data
int mval #O mask value for the output line segment
long v[IM_MAXDIM] #U vector coordinates of first pixel
int npix #O number of pixels in output line segment
int x1, i
long ve[IM_MAXDIM]
pointer pm, im, rl, rp, bp
pointer impl2$t(), impl3$t(), impgs$t()
errchk impl2$t, impl3$t, impgs$t, pm_glri
bool pm_sectnotempty()
int plloop()
begin
pm = M_PM(mp)
rl = M_RLP(mp)
# Initialization performed for the first i/o on a new region.
if (M_ACTIVE(mp) == NO) {
call plsslv (pm, M_VS(mp,1), M_VN(mp,1), M_V(mp,1), M_VE(mp,1))
call pm_glri (pm,
M_V(mp,1), Memi[rl], M_DEPTH(mp), M_VN(mp,1), PIX_SRC)
M_RLI(mp) = RL_FIRST
M_ACTIVE(mp) = YES
}
# Get a new mask line?
while (M_RLI(mp) > RLI_LEN(rl))
if (plloop (M_V(mp,1), M_VS(mp,1), M_VE(mp,1),
M_NDIM(mp)) == LOOP_DONE) {
return (EOF)
} else {
call amovl (M_V(mp,1), ve, M_NDIM(mp))
ve[1] = M_VE(mp,1)
if (pm_sectnotempty (pm, M_V(mp,1), ve, M_NDIM(mp))) {
call pm_glri (pm,
M_V(mp,1), Memi[rl], M_DEPTH(mp), M_VN(mp,1), PIX_SRC)
M_RLI(mp) = RL_FIRST
}
}
# Get a new image line?
if (M_RLI(mp) == RL_FIRST) {
call amovl (M_V(mp,1), v, IM_MAXDIM)
im = M_IM(mp)
if (M_LINEIO(mp) == YES && M_NDIM(mp) == 2)
bp = impl2$t (im, v[2])
else if (M_LINEIO(mp) == YES && M_NDIM(mp) == 3)
bp = impl3$t (im, v[2], v[3])
else
bp = impgs$t (im, v, ve, M_NDIM(mp))
M_LBP(mp) = bp
} else
bp = M_LBP(mp)
# Return the next line segment.
rp = rl + (M_RLI(mp) - 1) * RL_LENELEM
M_RLI(mp) = M_RLI(mp) + 1
x1 = Memi[rp+RL_XOFF]
npix = Memi[rp+RL_NOFF]
mval = Memi[rp+RL_VOFF]
ptr = bp + x1 - M_VS(mp,1)
if (M_REGCOORDS(mp) == NO) {
v[1] = x1
do i = 2, M_NDIM(mp)
v[i] = M_V(mp,i)
} else {
v[1] = x1 - M_VS(mp,1) + 1
do i = 2, M_NDIM(mp)
v[i] = M_V(mp,i) - M_VS(mp,i) + 1
}
return (npix)
end
|
