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
|
# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
include <plset.h>
include <plio.h>
# PL_REGIONROP -- Perform a rasterop operation on an irregular region in a
# 2-dimensional plane of a mask. The boundary of the region is defined by
# a user supplied function which, given the line number (Y) returns a range
# list defining the region for that line.
#
# rl_new = ufcn (ufd, y, rl_out, xs, npix)
#
# where
# rl_new = true if range list for line Y is different than for Y-1
# ufd = user region descriptor (parameters defining region)
# y = input y value
# rl_out = output range list
# xs = first pixel to be edited in dst mask
# npix = number of pixels in edit region
#
# If the dimensionality of the mask exceeds 2, the pl_setplane() procedure
# should be called first to define the plane of the mask to be modified.
procedure pl_regionrop (pl, ufcn, ufd, y1, y2, rop)
pointer pl #I mask descriptor
extern ufcn() #I user supplied region tracing procedure
pointer ufd #I user region descriptor
int y1 #I first mask line to be modified
int y2 #I last mask line to be modified
int rop #I rasterop defining operation
bool rl_new
long v[PL_MAXDIM]
int ll_len, xs, npix
pointer sp, ll_out, ll_reg, ll_dst, ol_dst, rl_out
pointer pl_access()
int pl_r2li()
bool ufcn()
errchk plvalid
begin
call plvalid (pl)
call smark (sp)
call salloc (ll_out, LL_MAXLEN(pl), TY_SHORT)
call salloc (ll_reg, LL_MAXLEN(pl), TY_SHORT)
call salloc (rl_out, RL_MAXLEN(pl), TY_INT)
call amovl (PL_PLANE(pl,1), v, PL_MAXDIM)
ol_dst = 1
for (v[2]=y1; v[2] <= y2; v[2]=v[2]+1) {
ll_dst = pl_access (pl, v)
rl_new = ufcn (ufd, v[2], Memi[rl_out], xs, npix)
if (rl_new)
ll_len = pl_r2li (Memi[rl_out], 1, Mems[ll_reg], npix)
if (ll_dst != ol_dst || rl_new) {
call pl_linestencil (Mems[ll_reg], 1, 1,
Mems[ll_dst], xs, PL_MAXVAL(pl),
Mems[ll_reg], 1,
Mems[ll_out], npix, rop)
ol_dst = ll_dst
}
# Update the affected line of the destination mask.
call pl_update (pl, v, Mems[ll_out])
}
# Compress the mask if excessive free space has accumulated.
if (PL_NEEDCOMPRESS(pl))
call pl_compress (pl)
call sfree (sp)
end
|
