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
|
include <mach.h>
include <imhdr.h>
include <error.h>
include "cyber.h"
# READ_IPPS_ROWS -- reads the ipps raster image row by row from the tape and
# writes the output image. At the completion of READ_IMAGE, the tape is
# positioned to the header record of the next ipps raster.
procedure read_ipps_rows (rd, out_fname, data_type, dt)
int rd, data_type
char out_fname[SZ_FNAME]
pointer dt
pointer sp, im, spp_buf
int npru_skip, i, stat
long clktime()
int read_dumpf()
pointer immap(), impl2r()
errchk cy_skip_pru, immap, ipps_to_iraf, read_dumpf()
begin
# Allocate buffer for ipps raster pixels
call smark (sp)
call salloc (spp_buf, WRDS_PER_ROW(dt) * NINT_CYBER_WRD, TY_INT)
# Map new iraf image and set up image header
im = immap (out_fname, NEW_IMAGE, 0)
IM_LEN(im, 1) = NCOLS(dt)
IM_LEN(im, 2) = NROWS(dt)
call strcpy (IPPS_ID(dt), IM_TITLE(im), SZ_IMTITLE)
IM_MIN(im) = DATA_MIN(dt)
IM_MAX(im) = DATA_MAX(dt)
# Set optimum image pixel type
if (data_type == NOT_SET) {
switch (BITS_PIXEL(dt)) {
case 12:
IM_PIXTYPE(im) = TY_SHORT
case 20:
IM_PIXTYPE(im) = TY_REAL
case 30:
IM_PIXTYPE(im) = TY_REAL
case 60:
IM_PIXTYPE(im) = TY_REAL
default:
call error (3, "IPPS BITS_PIXEL is incorrect")
}
} else
IM_PIXTYPE(im) = data_type
IM_LIMTIME(im) = clktime (long(0))
# Loop over rows to read, reorder and convert pixels.
for (i=1; i <= NROWS(dt); i=i+1) {
iferr {
stat = read_dumpf (rd, Memi[spp_buf], WRDS_PER_ROW(dt))
} then {
call imunmap (im)
call sfree (sp)
call erract (EA_WARN)
return
}
if (stat == EOF) {
call imunmap (im)
call sfree (sp)
call eprintf ("Premature EOF in image at row # %d\n")
call pargi (i)
return
}
call ipps_to_iraf (Memi[spp_buf], Memr[impl2r(im,i)], NCOLS(dt),
BITS_PIXEL(dt))
}
# Skip from present position to end of rcopy raster
npru_skip = PRU_EOR(dt) - PRU_ROW_ONE(dt) - (NPRU_ROW(dt)*NROWS(dt))
call cy_skip_pru (rd, npru_skip)
call imunmap (im)
call sfree (sp)
end
# IPPS_TO_IRAF -- performs the conversion from Cyber pixels to IRAF pixels.
# Each row of the ipps image is required to occupy an integral number of Cyber
# PRU's, so the input buffer contains pixels plus filler. The number of
# pixels unpacked from this buffer will always fill an integral number of
# Cyber words; a maximum of 4 extraneous pixels will be unpacked.
procedure ipps_to_iraf (in_buf, iraf_real, npix, nbits_pixel)
int in_buf[ARB]
real iraf_real[npix]
pointer iraf_int, sp
int nbits_pixel, npix, offset, npix_unpk, npix_cyber_wrd
errchk unpk_12, unpk_20
begin
# Allocate (maximum) space needed on the stack.
call smark (sp)
call salloc (iraf_int, npix + 4, TY_INT)
offset = 1
switch (nbits_pixel) {
case 12:
npix_cyber_wrd = 5
npix_unpk = ((npix + 4) / npix_cyber_wrd) * npix_cyber_wrd
call unpk_12 (in_buf, offset, Memi[iraf_int], npix_unpk)
call achtir (Memi[iraf_int], iraf_real, npix)
case 20:
npix_cyber_wrd = 3
npix_unpk = ((npix + 2) / npix_cyber_wrd) * npix_cyber_wrd
call unpk_20 (in_buf, offset, Memi[iraf_int], npix_unpk)
call achtir (Memi[iraf_int], iraf_real, npix)
case 30:
call unpk_30 (in_buf, offset, iraf_real, npix)
case 60:
call unpk_60r (in_buf, offset, iraf_real, npix)
default:
call error (5, "Illegal IPPS #B/P")
}
call sfree (sp)
end
# CY_SKIP_IMAGE -- skips over an IPPS raster once the header has been
# read. When SKIP_IMAGE returns, the tape is positioned to the first
# record of the next tape image.
procedure cy_skip_image (rd, dt)
int rd
pointer dt
int npru_skip
errchk cy_skip_pru
begin
# Calculate number of PRU's in image
npru_skip = PRU_EOR(dt) - PRU_ROW_ONE(dt)
call cy_skip_pru (rd, npru_skip)
end
# CY_SKIP_PRU -- positions the tape by skipping the requested number of PRU's.
procedure cy_skip_pru (rd, npru_skip)
int rd, npru_skip
pointer sp, dummy
int read_dumpf()
begin
call smark (sp)
call salloc (dummy, NINT_CYBER_WRD * LEN_PRU * npru_skip, TY_INT)
if (read_dumpf (rd, Memi[dummy], npru_skip * LEN_PRU) == EOF)
call eprintf ("Unexpected EOF when skipping image\n")
call sfree (sp)
end
|
