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
|
procedure qproc (image_list)
begin
struct buffer
string image, answr, imtype
int i, len, nampsx, nampsy
imtype = "." // envget ("imtype")
i = stridx (",", imtype)
if (i > 0)
imtype = substr (imtype, 1, i-1)
i = strlen (imtype)
cache ("quadsplit", "quadjoin", "qccdproc", "quadproc")
# 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)
}
}
# Initialise interactive query
if (qccdproc.interactive) {
answer.p_value = "yes"
answr = "yes"
} else {
answr = "NO"
}
fd = image_list
while (fscan (fd, image) != EOF) {
len = strlen (image)
if (substr(image, len-i+1, len) == imtype) {
image = substr (image, 1, len-i)
}
# Split out one image for each quadrant and set header sections
#quadsplit (image, output="",
#xskip1=xskip1, xskip2=xskip2, xtrim1=xtrim1, xtrim2=xtrim2,
#ytrim1=ytrim1, ytrim2=ytrim2, clobber=yes)
quadsplit (image, output="", clobber=yes)
# Find out of interactive fit is required for this image
if (answr == "yes" || answr == "no") {
printf ("Fit overscan vector for %s interactively\n", image) |
scan (buffer)
answer.p_prompt=buffer
answr = answer
}
# Overscan correct and trim
if (answr == "yes" || answr == "YES") {
qccdproc.interactive = yes
print ("YES") | qccdproc (image//".??"//imtype, fixpix=fixpix,
overscan=overscan, trim=trim, readaxis=readaxis,
fixfile=fixfile, biassec="image", trimsec="image",
ccdtype="", max_cache=0, noproc=no, zerocor=no, darkcor=no,
flatcor=no, illumcor=no, fringecor=no, readcor=no,
scancor=no, zero="", dark="", flat="", illum="", fringe="",
minreplace=1., scantype="shortscan", nscan=1, backup="",
logfile="", verbose=no, >> "dev$null")
# Set parameters of quadproc used for overscan fitting to match
# the ccdproc values which may have been adjusted interactively.
# We do this on every pass in case there is a later interupt
# of task execution.
quadproc.function.p_value = qccdproc.function
quadproc.order.p_value = qccdproc.order
quadproc.sample.p_value = qccdproc.sample
quadproc.naverage.p_value = qccdproc.naverage
quadproc.niterate.p_value = qccdproc.niterate
quadproc.low_reject.p_value = qccdproc.low_reject
quadproc.high_reject.p_value = qccdproc.high_reject
quadproc.grow.p_value = qccdproc.grow
# Force the parameter update
update ("quadproc")
} else {
qccdproc.interactive = no
qccdproc (image//".??"//imtype, fixpix=fixpix,
overscan=overscan, trim=trim, readaxis=readaxis,
fixfile=fixfile, biassec="image", trimsec="image",
ccdtype="", max_cache=0, noproc=no, zerocor=no, darkcor=no,
flatcor=no, illumcor=no, fringecor=no, readcor=no,
scancor=no, zero="", dark="", flat="", illum="", fringe="",
minreplace=1., scantype="shortscan", nscan=1, backup="",
logfile="", verbose=no)
}
# Combine processed quadrants into single image
quadjoin (image, output="", delete=yes)
}
# Reset interactive flag if we haven't recieved a definative NO
if (answr == "no") {
qccdproc.interactive = yes
}
end
|
