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
|
# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
include <imhdr.h>
# IC_RMASKS -- Set pixels for rejection mask.
procedure ic_rmasks (pm, v, id, nimages, n, npts)
pointer pm #I Pixel mask
long v[ARB] #I Output vector (input)
pointer id[nimages] #I Image id pointers
int nimages #I Number of images
int n[npts] #I Number of good pixels
int npts #I Number of output points per line
int i, j, k, ndim, impnls()
long v1[IM_MAXDIM]
pointer buf
begin
ndim = IM_NDIM(pm)
do k = 1, nimages {
call amovl (v, v1, ndim-1)
v1[ndim] = k
i = impnls (pm, buf, v1)
do j = 1, npts {
if (n[j] == nimages)
Mems[buf+j-1] = 0
else {
Mems[buf+j-1] = 1
do i = 1, n[j] {
if (Memi[id[i]+j-1] == k) {
Mems[buf+j-1] = 0
break
}
}
}
}
}
end
|
