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
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
|
include <imhdr.h>
include <mii.h>
include <mach.h>
include "rpds.h"
# PDS_READ_HEADER -- Read a PDS header. EOT is detected by an EOF on the
# first read and EOF is returned to the calling routine. Errors are
# passed to the calling routine.
int procedure pds_read_header (pds_fd, im, parameters)
int pds_fd
pointer im
long parameters[LEN_PAR_ARRAY]
int nchars, sz_header, header[LEN_PDS_HEADER]
short temptext[LEN_PDS_TEXT]
char text[LEN_PDS_TEXT]
int read(), pds_roundup()
short pds_unpacks()
long pds_unpackl()
errchk read, miiupk
include "rpds.com"
begin
# Read the header record
sz_header = pds_roundup (LEN_PDS_HEADER, SZB_CHAR) / SZB_CHAR
nchars = read (pds_fd, header, sz_header)
if (nchars == EOF)
return (EOF)
else if (nchars != sz_header)
call error (1, "Error reading pds header.")
# Unpack ID string, convert to ASCII and copy to image header
call miiupk (header, temptext, LEN_PDS_TEXT, MII_SHORT, TY_SHORT)
call pds_apdp8s (temptext, temptext, LEN_PDS_TEXT)
call pds_apdp059 (temptext, text, LEN_PDS_TEXT)
text[LEN_PDS_TEXT + 1] = EOS
call strcpy (text, TITLE(im), SZ_TITLE)
# Unpack the remainder of the header
# If XCOORD or YCOORD are greater than TWO_TO_23 convert to
# -ve number using TWO_TO_24
P_DX(parameters) = DX(header)
P_DY(parameters) = DY(header)
P_NPTS_PER_SCAN(parameters) = NPTS_PER_SCAN(header)
P_NSCANS(parameters) = NSCANS(header)
P_SCANTYPE(parameters) = SCANTYPE(header)
P_SCANSPEED(parameters) = SCANSPEED(header)
P_SCANORIGIN(parameters) = SCANORIGIN(header)
P_CORNER(parameters) = CORNER(header)
P_NRECS_PER_SCAN(parameters) = NRECS_PER_SCAN(header)
P_XTRAVEL(parameters) = XTRAVEL(header)
P_YTRAVEL(parameters) = YTRAVEL(header)
P_NPTS_PER_REC(parameters) = NPTS_PER_REC(header)
P_XCOORD(parameters) = XCOORD(header)
if (P_XCOORD(parameters) >= TWO_TO_23)
P_XCOORD(parameters) = P_XCOORD(parameters) - TWO_TO_24
P_YCOORD(parameters) = YCOORD(header)
if (P_YCOORD(parameters) >= TWO_TO_23)
P_YCOORD(parameters) = P_YCOORD(parameters) - TWO_TO_24
# Write parameters to header
CT_VALID(im) = NO
NAXIS(im) = 2
NCOLS(im) = P_NPTS_PER_SCAN(parameters)
NLINES(im) = P_NSCANS(parameters)
# print the header
call pds_print_header (text, parameters)
return (OK)
end
# PDS_UNPACKS -- Procedure to unpack a short header value.
# The header value is stored in the integer array buffer, beginning
# at byte number offset. The value is unpacked into a temporary
# buffer, temp and converted to SPP format using the mii routines.
# Finally the PDS 10 or 12 bit bytes are converted to an SPP short value.
short procedure pds_unpacks (buffer, offset)
int buffer[ARB], offset
short value[1]
long temp[1]
errchk miiupk, miiupk
begin
call bytmov (buffer, offset, temp, 1, SZB_MIISHORT)
call miiupk (temp, value[1], 1, MII_SHORT, TY_SHORT)
call pds_apdp8s (value[1], value[1], 1)
return (value[1])
end
# PDS_UNPACKL -- Procedure to unpack a 24 bit long header value.
long procedure pds_unpackl (buffer, offset)
int buffer[ARB], offset
short temps[2]
long temp[1], value[1]
errchk miiupk, bytmov
begin
call bytmov (buffer, offset, temp, 1, SZB_MIILONG)
call miiupk (temp, temps, 2, MII_SHORT, TY_SHORT)
call pds_apdp8s (temps, temps, 2)
value[1] = temps[1] * 10000b + temps[2]
return (value[1])
end
# PDS_APDP8S -- Precedure to change a 12 or 10 bit PDP8 value to a short integer
# value.
procedure pds_apdp8s (a, b, npix)
short a[npix], b[npix]
int npix
int i
begin
for (i=1; i <= npix; i = i + 1)
b[i] = (a[i] / 400b) * 100b + mod (int (a[i]), 400b)
end
# PDS_APDP059 -- Procedure to convert PDP 059 code into ASCII
procedure pds_apdp059 (a, b, nchar)
char b[nchar]
int nchar
short a[nchar]
int i, j
char table[LEN_TABLE]
# Conversion table from PDS 059 code to ASCII
data table/ ' ', ' ', ' ', ' ', '$', ' ', ' ', '"', '(', ')',
'*', '+', ',', '-', '.', '/', '0', '1', '2', '3',
'4', '5', '6', '7', '8', '9', ' ', ' ', ' ', '=',
' ', ' ', ' ', 'A', 'B', 'C', 'D', 'E', 'F', 'G',
'H', 'I', 'J', 'K', 'L', 'M', 'N', 'O', 'P', 'Q',
'R', 'S', 'T', 'U', 'V', 'W', 'X', 'Y', 'Z', ' ',
' ', ' ', ' ', ' ', ' ', 'a', 'b', 'c', 'd', 'e',
'f', 'g', 'h', 'i', 'j', 'k', 'l', 'm', 'n', 'o',
'p', 'q', 'r', 's', 't', 'u', 'v', 'w', 'x', 'y',
'z', ' ', ' ', ' ', ' ', ' ', EOS/
begin
for (i = 1; i <= nchar; i = i + 1)
b[i] = a[i] - 159
for (i = 1; i <= nchar; i = i + 1) {
j = b[i]
if (j < 1 || j > LEN_TABLE)
j = 1
b[i] = table[j]
}
end
# PDS_PRINT_HEADER -- Procedure to print the header.
procedure pds_print_header (text, parameters)
char text[LEN_PDS_TEXT]
long parameters[LEN_PAR_ARRAY]
include "rpds.com"
begin
if (long_header == YES)
call pds_long_header (text, parameters)
if (short_header == YES && long_header == NO) {
call printf ("ID: %.30s ")
call pargstr (text)
call printf ("size =%d * %d \n")
call pargl (P_NPTS_PER_SCAN(parameters))
call pargl (P_NSCANS(parameters))
}
end
# PDS_LONG_HEADER -- Print a long header
procedure pds_long_header (text, parameters)
char text[LEN_PDS_TEXT]
long parameters[LEN_PAR_ARRAY]
begin
call printf ("ID:%s\n")
call pargstr (text)
call printf ("NPTS = %d ")
call pargl (P_NPTS_PER_SCAN(parameters))
call printf ("NSCANS = %d ")
call pargl (P_NSCANS(parameters))
call printf ("NRECS/SCAN = %d ")
call pargl (P_NRECS_PER_SCAN(parameters))
call printf ("PPERR = %d\n")
call pargl (P_NPTS_PER_REC(parameters))
call printf ("SCANTYPE = %s ")
if (P_SCANTYPE(parameters) == RASTER)
call pargstr ("RASTER")
else if (P_SCANTYPE(parameters) == EDGE)
call pargstr ("EDGE")
else
call pargstr ("FLIPPED")
call printf ("SCANSPEED = %d ")
call pargl (P_SCANSPEED(parameters))
call printf ("SCANORIGIN = %d ")
call pargl (P_SCANORIGIN(parameters))
call printf ("CORNER = %d\n")
call pargl (P_CORNER(parameters))
call printf ("DX = %d ")
call pargl (P_DX(parameters))
call printf ("XTRAVEL = %d ")
call pargl (P_XTRAVEL(parameters))
call printf ("XCOORD = %d\n")
call pargl (P_XCOORD(parameters))
call printf ("DY = %d ")
call pargl (P_DY(parameters))
call printf ("YTRAVEL = %d ")
call pargl (P_YTRAVEL(parameters))
call printf ("YCOORD = %d\n")
call pargl (P_YCOORD(parameters))
end
|
