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
|
# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
include <syserr.h>
include <imhdr.h>
include <imio.h>
include <fset.h>
include <mach.h>
include "stf.h"
define NBITS_CHAR (SZB_CHAR * NBITS_BYTE)
# STF_OPIX -- Open (or create) the pixel storage file. If the image header file
# is `image.hhh' the associated pixel storage file will be `image.hhd' in the
# same directory as the header file. STF_PFD is set if the pixel file is
# physically open. IM_PFD is not set until we have been called by IMIO, since
# we must be called to once set up all the descriptors, even if the pixel file
# was already opened to read the GPB.
#
# dlb, 18-may-88: added code to zero out gpb's in multi-group image for groups
# other than current; prevents strange numbers and when later programs try to
# read the gpb of an otherwise uninitialized group of the image.
# dlb, 29-dec-1988: added code to get default set of GPB parameters and
# correctly initialize STF-kernel descriptor.
procedure stf_opix (im, status)
pointer im # image descriptor
int status # return status
int compress, blklen
bool copy_of_stf_image
int pfd, sz_gpb, group, i
pointer stf, o_stf, o_im, ua, gpb
long sz_pixfile, pixoff, totpix, offset
int open()
errchk open, fseti, falloc, seek, syserrs, imioff, calloc
errchk write
include <szpixtype.inc>
begin
status = OK
if (IM_PFD(im) != NULL)
return
o_im = IM_OHDR(im)
stf = IM_KDES(im)
ua = IM_USERAREA(im)
pfd = STF_PFD(stf)
compress = YES
blklen = 1
pixoff = 1
switch (IM_ACMODE(im)) {
case READ_ONLY, READ_WRITE, WRITE_ONLY, APPEND:
if (pfd == NULL)
pfd = open (IM_PIXFILE(im), IM_ACMODE(im), BINARY_FILE)
case NEW_COPY, NEW_FILE:
# Initialize the IMIO and STF descriptors and allocate the pixel
# file.
if (STF_NEWIMAGE(stf) == YES) {
# Normalize IMIO header parameters for new image.
call imioff (im, pixoff, compress, blklen)
# Set up the required GPB parameters for the new image.
# Note - this call can change the STF pointer.
call stf_newimage (im)
stf = IM_KDES(im)
# Save the size of the old GPB user area header if we are
# making a new copy of an old STF format image.
copy_of_stf_image = false
if (IM_ACMODE(im) == NEW_COPY && o_im != NULL)
if (IM_KERNEL(o_im) == IM_KERNEL(im))
copy_of_stf_image = true
if (copy_of_stf_image) {
o_stf = IM_KDES(o_im)
STF_PCOUNT(stf) = STF_PCOUNT(o_stf)
STF_PSIZE(stf) = STF_PSIZE(o_stf)
}
# Since the stf_mergegpb code below has been deactivated,
# there is no need to do the complex and expensive spool/copy
# operation below. (dct 1/4/90)
# -------------------------------
# # We have to have space for the GPB data cards at the beginning
# # of the user area, so spool any existing user cards in a
# # buffer and truncate the user area at the end of the GPB.
#
# ua_fd = stropen (Memc[ua+sz_gpbhdr], ARB, READ_ONLY)
# spool = open ("opix_spool", READ_WRITE, SPOOL_FILE)
# call fcopyo (ua_fd, spool)
# call close (ua_fd)
# Memc[ua+sz_gpbhdr] = EOS
#
# # Merge any extra GPB parameters from the old image into the
# # GPB structure of the new image. The GPB data cards for
# # these parameters should already be in the user area.
# # Order the group parameters to match the ordering in the
# # old image. NOTE: since the STF now copies all relevant
# # GPB parameters from an old image into the new or
# # generates a default standard set (in stf_newimage),
# # the following is no longer necessary. Note that if we
# # eventually may add parameters to the GPB, these routines
# # will again be useful!
#
# #if (copy_of_stf_image) {
# # call stf_mergegpb (im, o_im)
# # call stf_ordergpb (o_stf, stf)
# #}
#
# # Now append the spooled user header cards to the new user
# # area following the GPB data cards, deleting any user cards
# # which redefine GPB cards in the process.
#
# call seek (spool, BOFL)
# ua_size = (IM_LENHDRMEM(im) - LEN_IMHDR) * SZ_STRUCT
# ua_fd = stropen (Memc[ua], ua_size, APPEND)
# call stf_copyfits (stf, spool, NULL, ua_fd)
# call close (ua_fd)
# call close (spool)
#
# # Compute the length of the new header
# IM_HDRLEN(im) = LEN_IMHDR +
# (strlen(Memc[ua]) + SZ_STRUCT-1) / SZ_STRUCT
# Open the new pixel storage file (preallocate space if
# enabled on local system). Save the physical pathname of
# the pixfile in the image header, in case "imdir$" changes.
sz_pixfile = STF_SZGROUP(stf) * STF_GCOUNT(stf)
call falloc (IM_PIXFILE(im), sz_pixfile)
# Zero out all remaining groups of the image
# Open pixel file if not already open
if (STF_PFD(stf) == NULL)
pfd = open (IM_PIXFILE(im), READ_WRITE, BINARY_FILE)
# Allocate a zeroed block of memory whose length is the same
# as that of the group parameter block
sz_gpb = STF_PSIZE(stf) / NBITS_BYTE / SZB_CHAR
call calloc (gpb, sz_gpb, TY_CHAR)
# Zero out every group except the current one.
do group = 1, STF_GCOUNT(stf) {
if (group != STF_GROUP(stf)) {
offset = (group * STF_SZGROUP(stf) + 1) - sz_gpb
call seek (pfd, offset)
call write (pfd, Memc[gpb], sz_gpb)
}
}
# Free the block of memory.
call mfree (gpb, TY_CHAR)
} else {
# If we are writing to a group of an existing multigroup image,
# verify that the important image parameters have not been
# changed.
if (STF_NAXIS(stf) != IM_NDIM(im))
call syserrs (SYS_IMGSZNEQ, IM_NAME(im))
do i = 1, IM_NDIM(im)
if (STF_LENAXIS(stf,i) != IM_LEN(im,i))
call syserrs (SYS_IMGSZNEQ, IM_NAME(im))
# Added 5/15/87--dlb to get correct size of each data portion
# of a group if image opened NEW_COPY and input file was a
# template of a different dimensionality used to get GPB.
# Compute the size of each group in the pixel file, in chars.
totpix = IM_LEN(im,1)
do i = 2, IM_NDIM(im)
totpix = totpix * IM_LEN(im,i)
STF_SZGROUP(stf) = totpix * pix_size[IM_PIXTYPE(im)] +
STF_PSIZE(stf) / (SZB_CHAR * NBITS_BYTE)
}
if (pfd == NULL)
pfd = open (IM_PIXFILE(im), READ_WRITE, BINARY_FILE)
# Tell IMIO where the pixels are.
pixoff = (STF_GROUP(stf) - 1) * STF_SZGROUP(stf) + 1
call imioff (im, pixoff, compress, blklen)
default:
call imerr (IM_NAME(im), SYS_IMACMODE)
}
STF_PFD(stf) = pfd
IM_PFD(im) = pfd
end
|
