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
|
include <error.h>
define TYPES "|string|real|integer|"
define SVAL 1 # String value
define RVAL 2 # Real value
define IVAL 3 # Integer value
# T_CCDHEDIT -- Add, delete, or change CCD image header parameters.
# This task differs from HEDIT in that it uses the CCD instrument translation
# file.
procedure t_ccdhedit ()
int list # List of CCD images
pointer param # Parameter name
int type # Parameter type
pointer sval # Parameter value
pointer instrument # Instrument file
int ip, ival, imtopenp(), imtgetim(), clgwrd(), ctoi(), ctor()
real rval
bool streq()
pointer sp, im, immap()
errchk hdmpstr, hdmputr, hdmputi
begin
call smark (sp)
call salloc (param, SZ_LINE, TY_CHAR)
call salloc (sval, SZ_LINE, TY_CHAR)
call salloc (instrument, SZ_FNAME, TY_CHAR)
# Get the task parameters.
list = imtopenp ("images")
call clgstr ("parameter", Memc[param], SZ_LINE)
type = clgwrd ("type", Memc[sval], SZ_LINE, TYPES)
call clgstr ("value", Memc[sval], SZ_LINE)
call clgstr ("instrument", Memc[instrument], SZ_FNAME)
call xt_stripwhite (Memc[sval])
# Open the instrument translation file.
call hdmopen (Memc[instrument])
# If the parameter is IMAGETYP then change the parameter value from
# the package form to the image form using the inverse mapping in the
# translation file.
if (streq (Memc[param], "imagetyp"))
call hdmparm (Memc[sval], Memc[sval], SZ_LINE)
# Edit each image in the input list.
while (imtgetim (list, Memc[instrument], SZ_FNAME) != EOF) {
iferr (im = immap (Memc[instrument], READ_WRITE, 0)) {
call erract (EA_WARN)
next
}
# If the parameter value is null then delete the entry.
if (Memc[sval] == EOS) {
iferr (call hdmdelf (im, Memc[param]))
call erract (EA_WARN)
# Otherwise add the parameter of the specified type.
} else {
switch (type) {
case SVAL:
call hdmpstr (im, Memc[param], Memc[sval])
case RVAL:
ip = 1
if (ctor (Memc[sval], ip, rval) == 0)
call error (0, "Parameter value is not a number")
call hdmputr (im, Memc[param], rval)
case IVAL:
ip = 1
if (ctoi (Memc[sval], ip, ival) == 0)
call error (0, "Parameter value is not a number")
call hdmputi (im, Memc[param], ival)
}
}
call imunmap (im)
}
# Finish up.
call hdmclose ()
call imtclose (list)
call sfree (sp)
end
|
