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
|
# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
include <imhdr.h>
include <mach.h>
# IC_EMASK -- Create exposure mask.
procedure ic_emask (pm, v, id, nimages, n, wts, w, npts)
pointer pm #I Pixel mask
long v[ARB] #I Output vector
pointer id[nimages] #I Image id pointers
int nimages #I Number of images
int n[npts] #I Number of good pixels
real wts[npts] #I Weights
pointer w[npts] #I Weight data pointers
int npts #I Number of output pixels per line
int i, j, k, impnli()
real exp
pointer buf
pointer exps # Exposure times
pointer ev # IMIO coordinate vector
real ezero # Integer to real zero
real escale # Integer to real scale
int einit # Initialization flag
common /emask/ exps, ev, ezero, escale, einit
begin
# Write scaling factors to the header.
if (einit == NO) {
if (ezero != 0. || escale != 1.) {
call imaddr (pm, "MASKZERO", ezero)
call imaddr (pm, "MASKSCAL", escale)
}
einit = YES
}
call amovl (v, Meml[ev], IM_MAXDIM)
i = impnli (pm, buf, Meml[ev])
call aclri (Memi[buf], npts)
if (w[1] == NULL) {
do i = 1, npts {
exp = 0.
do j = 1, n[i] {
k = Memi[id[j]+i-1]
if (wts[k] > 0.)
exp = exp + Memr[exps+k-1]
}
Memi[buf] = nint((exp-ezero)/escale)
buf = buf + 1
}
} else {
do i = 1, npts {
exp = 0.
do j = 1, n[i] {
k = Memi[id[j]+i-1]
if (Memr[w[id[j]+i-1]+i-1] > 0.)
exp = exp + Memr[exps+k-1]
}
Memi[buf] = nint((exp-ezero)/escale)
buf = buf + 1
}
}
end
# IC_EINIT -- Initialize exposure mask.
procedure ic_einit (in, nimages, key, default, maxval)
int in[nimages] #I Image pointers
int nimages #I Number of images
char key[ARB] #I Exposure time keyword
real default #I Default exposure time
int maxval #I Maximum mask value
int i
real exp, emin, emax, efrac, imgetr()
pointer exps # Exposure times
pointer ev # IMIO coordinate vector
real ezero # Integer to real zero
real escale # Integer to real scale
int einit # Initialization flag
common /emask/ exps, ev, ezero, escale, einit
begin
call malloc (ev, IM_MAXDIM, TY_LONG)
call malloc (exps, nimages, TY_REAL)
emax = 0.
emin = MAX_REAL
efrac = 0
do i = 1, nimages {
iferr (exp = imgetr (in[i], key))
exp = default
exp = max (0., exp)
emax = emax + exp
if (exp > 0.)
emin = min (exp, emin)
efrac = max (abs(exp-nint(exp)), efrac)
Memr[exps+i-1] = exp
}
# Set scaling.
ezero = 0.
escale = 1.
if (emin < 1.) {
escale = emin
emin = emin / escale
emax = emax / escale
} else if (emin == MAX_REAL)
emin = 0.
if (efrac > 0.001 && emax-emin < 1000.) {
escale = escale / 1000.
emin = emin * 1000.
emax = emax * 1000.
}
while (emax > maxval) {
escale = escale * 10.
emin = emin / 10.
emax = emax / 10.
}
einit = NO
end
|
