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
103
|
# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
include <pmset.h>
include <plio.h>
# PM_PLLS -- Put a line segment input as a list list to a mask, applying the
# given ROP to combine the pixels with those of the output mask.
procedure pm_plls (pl, v, ll_raw, ll_depth, npix, rop)
pointer pl #I mask descriptor
long v[PL_MAXDIM] #I vector coords of line segment
short ll_raw[ARB] #I input line list
int ll_depth #I line list depth, bits
int npix #I number of pixels affected
int rop #I rasterop
pointer sp, ll_src, ll_dst, ll_stn, ll_out, px_src, im
int ll_len, step, xstep, temp, np, ip, i
int pl_l2pi(), pl_p2li()
pointer pl_access()
include "pmio.com"
begin
im = PM_REFIM(pl)
if (PM_MAPXY(pl) == NO) {
call pl_plls (pl, v, ll_raw, ll_depth, npix, rop)
return
}
call smark (sp)
call salloc (ll_src, LL_MAXLEN(pl), TY_SHORT)
# Determine physical coords of line segment.
call amovl (v, v3, PM_MAXDIM)
call imaplv (im, v3, v1, PM_MAXDIM)
v3[1] = v3[1] + npix - 1
call imaplv (im, v3, v2, PM_MAXDIM)
# Get line scaling parameters.
if (npix <= 1)
xstep = 1
else
xstep = (v2[1] - v1[1]) / (npix - 1)
step = xstep
if (xstep < 0) {
temp = v1[1]; v1[1] = v2[1]; v2[1] = temp
step = -step
}
np = (npix - 1) * step + 1
ll_stn = NULL
# Resample and flip the line list if necessary. Construct a stencil
# list if the step size is greater than 1.
if (xstep < 0 || step > 1) {
call salloc (px_src, np, TY_INT)
i = pl_l2pi (ll_raw, 1, Memi[px_src], npix)
call aclri (Memi[px_src+i], np - i)
# Flip data array.
if (xstep < 0)
call imaflp (Memi[px_src], npix, SZ_INT)
if (step > 1) {
# Resample data array.
ip = px_src + npix - 1
do i = np, 1, -step {
Memi[px_src+i-1] = Memi[ip]
ip = ip - 1
}
# Construct stencil.
call salloc (ll_stn, LL_MAXLEN(pl), TY_SHORT)
call aclri (Memi[px_src], np)
do i = 1, np, step
Memi[px_src+i-1] = 1
ll_len = pl_p2li (Memi[px_src], 1, Mems[ll_stn], np)
}
# Convert flipped and resampled data back to line list.
ll_len = pl_p2li (Memi[px_src], 1, Mems[ll_src], np)
} else {
ll_len = LL_LEN(ll_raw)
call amovs (ll_raw, Mems[ll_src], ll_len)
}
# Copy to or combine with destination.
if (ll_stn == NULL)
call pl_plls (pl, v1, Mems[ll_src], ll_depth, np, rop)
else {
call salloc (ll_out, LL_MAXLEN(pl), TY_SHORT)
ll_dst = pl_access (pl, v1)
call pl_linestencil (Mems[ll_src], 1, MV(ll_depth),
Mems[ll_dst], v1, PL_MAXVAL(pl), Mems[ll_stn], 1,
Mems[ll_out], np, rop)
call pl_update (pl, v1, Mems[ll_out])
}
call sfree (sp)
end
|
