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
|
procedure quadproc (images)
begin
string ims, in_list, cal_list, qp_list
int nims
struct buffer
# Freeze input image list
in_list = mktemp ("uparm$tmp")
ccdlist (images, ccdtype="", names+, > in_list)
ims = "@"//in_list
# Check that the input list contains some images of the specified type.
ccdssselect (ims, ccdtype=ccdtype, subset="") | count | scan (nims)
if (nims == 0) { # Nothing to do !
delete (in_list, ver-)
return
}
# Set initial values for the ccdproc parameters used for fitting the
# overscan. These parameters may be modified during the interactive
# fitting process, the new values being used for all subsequent fits.
qccdproc.function = function
qccdproc.order = order
qccdproc.sample = sample
qccdproc.naverage = naverage
qccdproc.niterate = niterate
qccdproc.low_reject = low_reject
qccdproc.high_reject = high_reject
qccdproc.grow = grow
qccdproc.interactive = interactive
if (overscan || trim || fixpix) {
# Only those images which must be treated specialy are processed
# with qproc:
# 1) Multiple readout
# 2) Not already trimmed
# First process [OT] any calibration images which WILL BE USED so
# that we can be sure they will be reduced explicitly by qproc
# rather than automaticaly within ccdproc.
qp_list = mktemp ("uparm$tmp")
if (noproc) {
cal_list = mktemp ("uparm$tmp")
qpcalimage (images=ims, only_param=yes, check=no, > cal_list)
qpselect ("@"//cal_list, ccdtype="", stop=no, > qp_list)
delete (cal_list, ver-)
count (qp_list) | scan (nims)
if (nims > 0) {
printf ("Calibration images which will be processed:\n")
qnoproc (qp_list, fixpix=fixpix, overscan=overscan,
trim=trim, fixfile=fixfile, xskip1=xskip1, xskip2=xskip2,
xtrim1=xtrim1, xtrim2=xtrim2, ytrim1=ytrim1, ytrim2=ytrim2)
printf ("\nImages from the input list:\n")
}
} else {
cal_list = mktemp ("uparm$tmp")
qpcalimage (images=ims, only_param=no, check=no, > cal_list)
qpselect ("@"//cal_list, ccdtype="", stop=no, > qp_list)
delete (cal_list, ver-)
count (qp_list) | scan (nims)
if (nims > 0) {
qproc (qp_list, fixpix=fixpix, overscan=overscan, trim=trim,
readaxis=readaxis, fixfile=fixfile, xskip1=xskip1,
xskip2=xskip2, xtrim1=xtrim1, xtrim2=xtrim2, ytrim1=ytrim1,
ytrim2=ytrim2)
}
}
delete (qp_list, ver-)
# Now process all the remaining images.
qpselect (ims, ccdtype=ccdtype, stop=no, >> qp_list)
if (noproc) {
qnoproc (qp_list, fixpix=fixpix, overscan=overscan,
trim=trim, fixfile=fixfile, xskip1=xskip1, xskip2=xskip2,
xtrim1=xtrim1, xtrim2=xtrim2, ytrim1=ytrim1, ytrim2=ytrim2)
} else {
qproc (qp_list, fixpix=fixpix, overscan=overscan, trim=trim,
readaxis=readaxis, fixfile=fixfile, xskip1=xskip1,
xskip2=xskip2, xtrim1=xtrim1, xtrim2=xtrim2, ytrim1=ytrim1,
ytrim2=ytrim2)
}
delete (qp_list, ver-)
}
# Do all other operations with the standard ccdproc.
if (noproc) {
# Handle those images that will be processed with qproc ....
qpselect (ims, ccdtype=ccdtype, stop=no, >> qp_list)
# We must also include the calibration images or ccdproc will
# complain about missing calibrations.
qpcalimage (images=ims, only_param=no, check=no, > cal_list)
qpselect ("@"//cal_list, ccdtype="", stop=yes, >> qp_list)
if (zerocor || darkcor || flatcor || illumcor || fringecor ||
readcor || scancor) {
qccdproc ("@"//qp_list, noproc=yes,
fixpix=no, overscan=no, trim=no, zerocor=zerocor,
darkcor=darkcor, flatcor=flatcor, illumcor=illumcor,
fringecor=fringecor, readcor=readcor, scancor=scancor,
ccdtype=ccdtype, max_cache=max_cache, readaxis=readaxis,
fixfile=fixfile, biassec="image", trimsec="image", zero=zero,
dark=dark, flat=flat, illum=illum, fringe=fringe,
minreplace=minreplace, scantype=scantype, nscan=nscan)
}
printf ("\n")
delete (qp_list, ver-)
# ..... and those that won't
qpselect (ims, ccdtype=ccdtype, stop=yes, >> qp_list)
qccdproc ("@"//qp_list, noproc=yes,
fixpix=fixpix, overscan=overscan, trim=trim, zerocor=zerocor,
darkcor=darkcor, flatcor=flatcor, illumcor=illumcor,
fringecor=fringecor, readcor=readcor, scancor=scancor,
ccdtype=ccdtype, max_cache=max_cache, readaxis=readaxis,
fixfile=fixfile, biassec="image", trimsec="image", zero=zero,
dark=dark, flat=flat, illum=illum, fringe=fringe,
minreplace=minreplace, scantype=scantype, nscan=nscan)
delete (qp_list, ver-)
} else {
# Validate fixfile
if (fixpix) {
match ("single_readout", fixfile) | scan (buffer)
if (stridx ("#", buffer) == 0) {
buffer = "fixfile " // fixfile //
" cannot be used with multi-readout images"
error (0, buffer)
}
}
qccdproc (ims, ccdtype=ccdtype, max_cache=max_cache, noproc=no,
fixpix=fixpix, overscan=overscan, trim=trim, zerocor=zerocor,
darkcor=darkcor, flatcor=flatcor, illumcor=illumcor,
fringecor=fringecor, readcor=readcor, scancor=scancor,
readaxis=readaxis, fixfile=fixfile, biassec="image",
trimsec="image", zero=zero, dark=dark, flat=flat, illum=illum,
fringe=fringe, minreplace=minreplace, scantype=scantype,
nscan=nscan, backup=backup, interactive=interactive,
verbose=verbose, logfile=logfile)
# Set task parameters used for overscan fitting to the ccdproc
# values which may have been adjusted interactively
function.p_value = qccdproc.function
order.p_value = qccdproc.order
sample.p_value = qccdproc.sample
naverage.p_value = qccdproc.naverage
niterate.p_value = qccdproc.niterate
low_reject.p_value = qccdproc.low_reject
high_reject.p_value = qccdproc.high_reject
grow.p_value = qccdproc.grow
}
delete (in_list, ver-)
end
|
