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
104
105
|
# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
include <plio.h>
# PL_L2P -- Convert a line list to a pixel array. The number of pixels output
# (always npix) is returned as the function value.
int procedure pl_l2pl (ll_src, xs, px_dst, npix)
short ll_src[ARB] #I input line list
int xs #I starting index in ll_src
long px_dst[ARB] #O output pixel array
int npix #I number of pixels to convert
long pv
bool skipword
int opcode, data, ll_len, ll_first
int x1, x2, i1, i2, xe, np, ip, op, otop, i
define putpix_ 91
begin
# Support old format line lists.
if (LL_OLDFORMAT(ll_src)) {
ll_len = OLL_LEN(ll_src)
ll_first = OLL_FIRST
} else {
ll_len = LL_LEN(ll_src)
ll_first = LL_FIRST(ll_src)
}
# No pixels?
if (npix <= 0 || ll_len <= 0)
return (0)
xe = xs + npix - 1
skipword = false
op = 1
x1 = 1
pv = 1
do ip = ll_first, ll_len {
if (skipword) {
skipword = false
next
}
opcode = I_OPCODE(ll_src[ip])
data = I_DATA(ll_src[ip])
switch (opcode) {
case I_ZN, I_HN, I_PN:
# Determine inbounds region of segment.
x2 = x1 + data - 1
i1 = max (x1, xs)
i2 = min (x2, xe)
# Process segment if any region is inbounds.
np = i2 - i1 + 1
if (np > 0) {
otop = op + np - 1
if (opcode == I_HN) {
do i = op, otop
px_dst[i] = pv
} else {
do i = op, otop
px_dst[i] = 0
if (opcode == I_PN && i2 == x2)
px_dst[otop] = pv
}
op = otop + 1
}
# Advance the line index.
x1 = x2 + 1
case I_SH:
pv = (int(ll_src[ip+1]) * I_SHIFT) + data
skipword = true
case I_IH:
pv = pv + data
case I_DH:
pv = pv - data
case I_IS:
pv = pv + data
goto putpix_
case I_DS:
pv = pv - data
putpix_
if (x1 >= xs && x1 <= xe) {
px_dst[op] = pv
op = op + 1
}
x1 = x1 + 1
}
if (x1 > xe)
break
}
# Zero any remaining output range.
do i = op, npix
px_dst[i] = 0
return (npix)
end
|
