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
|
include <tbset.h>
include "../lib/daophotdef.h"
include "../lib/apseldef.h"
# DP_RESTARS -- Read in the IDS of the stars to be excluded, find these stars
# in the photometry list, and set their magnitudes to INDEF.
int procedure dp_restars (dao, im, ext, text_file)
pointer dao # pointer to the daophot structure
pointer im # the input image descriptor
int ext # the exclude list file descriptor
bool text_file # text or table file ?
real tx, ty, rjunk
pointer apsel, sp, fields, indices, key
int i, nrow, idno, nexcl, starno
int tbpsta(), dp_apsel(), dp_exfind()
begin
# Get some pointers.
apsel = DP_APSEL(dao)
# Get some working space.
call smark (sp)
call salloc (fields, SZ_LINE, TY_CHAR)
call salloc (indices, 1, TY_INT)
# Initialize the read.
if (text_file) {
call pt_kyinit (key)
Memi[indices] = DP_PAPID
call dp_gappsf (Memi[indices], Memc[fields], 1)
} else {
call dp_tptinit (ext, Memi[indices])
nrow = tbpsta (ext, TBL_NROWS)
}
i = 1
nexcl = 0
repeat {
# Read the next star.
if (text_file) {
if (dp_apsel (key, ext, Memc[fields], Memi[indices], idno,
rjunk, rjunk, rjunk, rjunk) == EOF)
break
} else {
if (i > nrow)
break
call dp_tptread (ext, Memi[indices], idno, rjunk, rjunk, rjunk,
i)
}
# Subtract star from the photometry list.
if (idno > 0) {
starno = dp_exfind (Memi[DP_APID(apsel)],
Memr[DP_APXCEN(apsel)], Memr[DP_APYCEN(apsel)],
Memr[DP_APMAG(apsel)], DP_APNUM(apsel), idno)
if (starno > 0) {
if (DP_VERBOSE(dao) == YES) {
call dp_wout (dao, im, Memr[DP_APXCEN(apsel)+starno-1],
Memr[DP_APYCEN(apsel)+starno-1], tx, ty, 1)
call printf (
"EXCLUDING - Star:%5d X =%8.2f Y =%8.2f Mag =%8.2f\n")
call pargi (Memi[DP_APID(apsel)+starno-1])
call pargr (tx)
call pargr (ty)
call pargr (Memr[DP_APMAG(apsel)+starno-1])
}
nexcl = nexcl + 1
} else if (DP_VERBOSE(dao) == YES) {
call printf ("EXCLUDING - Star:%5d not found\n")
call pargi (idno)
}
}
i = i + 1
}
if (text_file)
call pt_kyfree (key)
call sfree (sp)
return (nexcl)
end
# DP_EXFIND -- Find the star to be exclude in the photometry list.
int procedure dp_exfind (ids, xcen, ycen, mags, nstars, idex)
int ids[ARB] # array of stellar ids
real xcen[ARB] # array of x coordinates
real ycen[ARB] # array of y coordinates
real mags[ARB] # array of magnitudes
int nstars # number of stars in photometry list
int idex # id of star to be excluded
int i, found
begin
found = 0
do i = 1, nstars {
if (ids[i] != idex)
next
found = i
break
}
if (found > 0)
mags[i] = INDEFR
return (found)
end
|
