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
|
# MOSPROC - Sky subtract, flat field and transpose images before mosaicing.
procedure mosproc (input, output, nxsub, nysub)
string input {prompt="Input images"}
string output {prompt="Output image"}
int nxsub {8, prompt="Number of subrasters in x"}
int nysub {8, prompt="Number of subrasters in y"}
bool skysubtract {yes, prompt="Sky subtract images before mosaicing"}
string sky {"", prompt="Sky image to subtract"}
string exclude {"", prompt="Input images excluded from sky frame"}
string expname {"EXPTIME", prompt="Image exposure time keywords"}
bool flatten {yes, prompt="Flatten images before mosaicing"}
string flat {"", prompt="Flat field image"}
bool transpose {no, prompt="Transpose images before mosaicing?"}
string trim_section {"[*,*]", prompt="Input image section to be extracted"}
string corner {"lr", prompt="Starting corner for the mosaic"}
string direction {"row", prompt="Starting direction for the mosaic"}
bool raster {no, prompt="Raster scan?"}
string median_section {"", prompt="Input subraster section for median ?"}
bool subtract {no, prompt="Substract median from each subraster?"}
real oval {-1.0, prompt="Mosaic border pixel values"}
bool delete {yes, prompt="Delete temporary images?"}
file logfile {"STDOUT", prompt="Log file name"}
struct *list1, *list2
begin
file tmpimg, tmptmp, tmpred, tmpexc
int nx, ny, i, nin, lo, hi
string skyframe, normframe, in, out, img, delim, junk
tmpimg = mktemp ("MOS")
tmptmp = mktemp ("MOS")
tmpred = mktemp ("MOS")
tmpexc = mktemp ("tmp$MOS")
# Get positional parameters
in = input
out = output
nx = nxsub
ny = nysub
# Expand input file name list removing the ".imh" extensions.
sections (in, option="fullname", > tmptmp)
list1 = tmptmp
for (nin = 0; fscan (list1, img) != EOF; nin += 1) {
i = strlen (img)
if (substr (img, i-3, i) == ".imh")
img = substr (img, 1, i-4)
print (img, >> tmpimg)
print (img // ".red", >> tmpred)
}
list1 = ""; delete (tmptmp, ver-, >& "dev$null")
# Expand the range of images to skip.
if (skysubtract && sky != "") {
skyframe = sky
imarith ("@"//tmpimg, "-", skyframe, "@"//tmpred, title="",
divzero=0., hparams="", pixtype="", calctype="", verbose+,
noact-, >> logfile)
} else if (skysubtract) {
print (exclude, ",") | translit ("", "^-,0-9", del+) |
translit ("", "-", "!", del-) | tokens (new-) |
translit ("", "\n,", " \n", del-, > tmpexc)
type (tmpexc, >> logfile)
list1 = tmpexc
while (fscan (list1, lo, delim, hi, junk) != EOF) {
if (nscan() == 0)
next
else if (nscan() == 1 && lo >= 1)
print (lo, >> tmptmp)
else if (nscan() == 3) {
lo = min (max (lo, 1), nin); hi = min (max (hi, 1), nin)
for (i = lo; i <= hi; i += 1)
print (i, >> tmptmp)
}
}
list1 = ""; delete (tmpexc, ver-, >& "dev$null")
if (access (tmptmp)) {
sort (tmptmp, col=0, ign+, num+, rev-) | unique (> tmpexc)
delete (tmptmp, ver-, >& "dev$null")
list1 = tmpimg; list2 = tmpexc; junk = fscan (list2, nin)
for (i = 1; fscan (list1, img) != EOF; i += 1) {
if (i == nin) {
junk = fscan (list2, nin)
next
}
print (img, >> tmptmp)
}
list1 = ""; list2 = ""; delete (tmpexc, ver-, >& "dev$null")
} else
tmptmp = tmpimg
skyframe = out // ".sky"
imcombine ("@"//tmptmp, skyframe, rejmask="", plfile="", sigma="",
logfile=logfile, combine="median", reject="none", project=no,
outtype="real", offsets="none", masktype="none", maskvalue=0.0,
blank=-1.0, scale="exposure", zero="none", weight="exposure",
statsec="", expname=expname, lthreshold=INDEF,
hthreshold=INDEF, nlow=1, nhigh=1, nkeep=1, mclip=yes,
lsigma=3.0, hsigma=3.0, rdnoise="0.0", gain="1.0", snoise="0.0",
sigscale=0.1, pclip=-0.5, grow=0)
print ("\n", >> logfile)
imarith ("@"//tmpimg, "-", skyframe, "@"//tmpred, title="",
divzero=0., hparams="", pixtype="", calctype="", verbose+,
noact-, >> logfile)
} else {
skyframe = ""
imcopy ("@"//tmpimg, "@"//tmpred, verbose-)
}
if (flatten) {
if (flat != "") {
print ("\n", >> logfile)
flatten ("@"//tmpred, flat, minflat=INDEF, pixtype="",
keeplog=yes, logfile=logfile)
} else if (skyframe != "") {
print ("\n", >> logfile)
normframe = out // ".norm"
imcopy (skyframe, normframe, verbose-)
bscale (normframe, normframe, bzero="0.0", bscale="mode",
section="", step=10, lower=INDEF, upper=INDEF,
verbose+, >>logfile)
print ("\n", >> logfile)
flatten ("@"//tmpred, normframe, minflat=INDEF, pixtype="",
keeplog=yes, logfile=logfile)
}
}
if (transpose) {
print ("\nTRANSPOSE: Transpose images", >> logfile)
time (, >> logfile)
imtrans ("@"//tmpred, "@"//tmpred)
time (, >> logfile)
print ("TRANSPOSE: done", >> logfile)
}
print ("\nIRMOSAIC: Mosaic images", >> logfile)
time (, >> logfile)
irmosaic ("@"//tmpred, out, "db"//out, nx, ny,
trim_section=trim_section, null_input="", corner=corner,
direction=direction, raster=raster, nxover=-1, nyover=-1,
nimcols=INDEF, nimrows=INDEF, oval=oval,
median_section=median_section, sub=subtract, opixtype="r",
verbose+, >> logfile)
time (, >> logfile)
print ("IRMOSAIC: done", >> logfile)
if (delete) {
if (access (tmpred))
imdelete ("@"//tmpred, ver-, >& "dev$null")
}
delete (tmpimg, ver-, >& "dev$null")
delete (tmptmp, ver-, >& "dev$null")
delete (tmpred, ver-, >& "dev$null")
end
|
