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
|
#{ NORMFLAT -- Make a flat field by normalizing and replacing low values.
# image,f,a,,,,Calibration image
# flatfield,f,a,,,,Flat field image
# norm,r,h,INDEF,,,Normalization if not INDEF
# minflat,r,h,INDEF,,,Minimum data value to use in the flat field
# sample_section,s,h,"[]",,,Sample section for determining normalization
# keeplog,b,h,@generic.keeplog,,,Keep log of processing?
# logfile,f,h,@generic.logfile,,,Log file
# img,f,h
# flt,f,h
# tmp,f,h
# rlist,*s,h
# mean,r,h
# stat,i,h
{
# Get query parameters and set temporary parameters.
img = image
flt = flatfield
tmp = mktemp ("tmp$gec")
# Startup message.
if (keeplog) {
time (>> logfile)
print (" NORMFLAT: Create a flat field.\n", >> logfile)
print (" Calibration image: ", img, >> logfile)
print (" Flat field: ", flt, >> logfile)
if (minflat != INDEF)
print (" Minimum data value used in flat field = ", minflat,
>> logfile)
}
# Determine normalization.
if (norm == INDEF) {
# Determine the mean of the sample region.
imstatistics (img // sample_section, fields="mean",
lower=minflat, upper=INDEF, format=no, > tmp)
rlist = tmp
stat = fscan (rlist, mean)
rlist = ""
delete (tmp, verify=no)
} else
mean = norm
if (keeplog)
print (" Normalization = ", mean, >> logfile)
# Replace low values by the mean and normalize.
if (mean != 0.) {
if (minflat != INDEF) {
imcopy (img, flt, verbose=no)
imreplace (flt, mean, upper=minflat)
imarith (flt, "/", mean, flt, pixtype="real")
} else
imarith (img, "/", mean, flt, pixtype="real")
} else
print (" ERROR: Cannot normalize calibration image.")
# Set CCDMEAN to 1.
hedit (flt, "ccdmean", "1.", add=yes, update=yes, show=no, verify=no)
# Ending message.
if (keeplog) {
time (>> logfile)
print (" NORMFLAT: Done.", >> logfile)
}
}
|
