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
|
include <imhdr.h>
include <error.h>
include "ccdred.h"
include "ccdtypes.h"
define CACHEUNIT 1000000. # Units of max_cache parameter
# T_CCDPROC -- Process CCD images
#
# This is the main procedure for processing CCD images. The images are
# corrected for bad pixels, overscan levels, zero levels, dark counts,
# flat field response, illumination errors, and fringe response. They
# may also be trimmed. The input is a list of images to be processed.
# Each image must match any image type requested. The checking of
# whether to apply each correction, getting the required parameters, and
# logging the operations is left to separate procedures, one for each
# correction. The actual processing is done by a specialized procedure
# designed to be very efficient. These procedures may also process
# calibration images if necessary. There are two data type paths; one
# for short pixel types and one for all other pixel types (usually
# real).
procedure t_ccdproc ()
int list # List of CCD images to process
int outlist # LIst of output images
int ccdtype # CCD image type
int interactive # Fit overscan interactively?
int max_cache # Maximum image cache size
bool clgetb()
real clgetr()
int imtopenp(), imtgetim(), imtlen()
pointer sp, input, output, str, in, out, ccd
errchk set_input, set_output, ccddelete, cal_open
errchk set_fixpix, set_zero, set_dark, set_flat, set_illum, set_fringe
begin
call smark (sp)
call salloc (input, SZ_FNAME, TY_CHAR)
call salloc (output, SZ_FNAME, TY_CHAR)
call salloc (str, SZ_LINE, TY_CHAR)
# Get input and output lists and check they make sense.
list = imtopenp ("images")
outlist = imtopenp ("output")
if (imtlen (outlist) > 0 && imtlen (outlist) != imtlen (list))
call error (1, "Input and output lists do not match")
# Get instrument translation file. Open the translation
# file. Initialize the interactive flag and the calibration images.
call clgstr ("instrument", Memc[input], SZ_FNAME)
if (Memc[input] == EOS)
call error (1, "No 'instrument' translation file specified.")
call hdmopen (Memc[input])
call set_interactive ("", interactive)
call cal_open (list)
if (imtlen (list) < 3)
max_cache = 0.
else
max_cache = CACHEUNIT * clgetr ("max_cache")
call ccd_open (max_cache)
# Process each image.
while (imtgetim (list, Memc[input], SZ_FNAME) != EOF) {
if (clgetb ("noproc")) {
call printf ("%s:\n")
call pargstr (Memc[input])
}
call set_input (Memc[input], in, ccdtype)
if (in == NULL)
next
# Set output image.
if (imtlen (outlist) == 0)
call mktemp ("tmp", Memc[output], SZ_FNAME)
else if (imtgetim (outlist, Memc[output], SZ_FNAME) == EOF)
call error (1, "Premature end of output list")
call set_output (in, out, Memc[output])
# Set processing parameters applicable to all images.
call set_proc (in, out, ccd)
call set_sections (ccd)
call set_trim (ccd)
call set_fixpix (ccd)
call set_overscan (ccd)
# Set processing parameters for the standard CCD image types.
switch (ccdtype) {
case ZERO:
case DARK:
call set_zero (ccd)
case FLAT:
call set_zero (ccd)
call set_dark (ccd)
CORS(ccd, FINDMEAN) = YES
CORS(ccd, MINREP) = YES
case ILLUM:
call set_zero (ccd)
call set_dark (ccd)
call set_flat (ccd)
case OBJECT, COMP:
call set_zero (ccd)
call set_dark (ccd)
call set_flat (ccd)
iferr {
call set_illum (ccd)
call set_fringe (ccd)
} then
call erract (EA_WARN)
default:
call set_zero (ccd)
call set_dark (ccd)
call set_flat (ccd)
iferr {
call set_illum (ccd)
call set_fringe (ccd)
} then
call erract (EA_WARN)
CORS(ccd, FINDMEAN) = YES
}
# Do the processing if the COR flag is set.
if (COR(ccd) == YES) {
call doproc (ccd)
call set_header (ccd)
call imunmap (in)
call imunmap (out)
if (imtlen (outlist) == 0) {
# Replace the input image by the corrected image.
iferr (call ccddelete (Memc[input])) {
call imdelete (Memc[output])
call error (1,
"Can't delete or make backup of original image")
}
call imrename (Memc[output], Memc[input])
}
} else {
# Delete the output image.
call imunmap (in)
iferr (call imunmap (out))
;
iferr (call imdelete (Memc[output]))
;
}
call free_proc (ccd)
# Do special processing on certain image types.
if (imtlen (outlist) == 0) {
switch (ccdtype) {
case ZERO:
call readcor (Memc[input])
case FLAT:
call ccdmean (Memc[input])
}
} else {
switch (ccdtype) {
case ZERO:
call readcor (Memc[output])
case FLAT:
call ccdmean (Memc[output])
}
}
}
# Finish up.
call hdmclose ()
call imtclose (list)
call imtclose (outlist)
call cal_close ()
call ccd_close ()
call sfree (sp)
end
|
