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
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
|
include <mach.h>
include "../lib/apphotdef.h"
include "../lib/apphot.h"
include "../lib/noisedef.h"
include "../lib/centerdef.h"
include "../lib/center.h"
define CONVERT .424660900 # conversion factor from fwhmpsf to sigma
# APREFITCENTER -- Procedure to refit the centers assuming that the appropriate
# pixel buffer is in memory. See apfitcenter for further information.
int procedure aprefitcenter (ap, im, ier)
pointer ap # pointer to the apphot structure
pointer im # the input image descriptor
int ier # previous error code
int fier
pointer cen, nse
int ap_ctr1d(), ap_mctr1d(), ap_gctr1d(), ap_lgctr1d()
begin
cen = AP_PCENTER(ap)
nse = AP_NOISE(ap)
# Initialize
AP_XCENTER(cen) = AP_CXCUR(cen)
AP_YCENTER(cen) = AP_CYCUR(cen)
AP_XSHIFT(cen) = 0.0
AP_YSHIFT(cen) = 0.0
AP_OXSHIFT(cen) = 0.0
AP_OYSHIFT(cen) = 0.0
AP_XERR(cen) = INDEFR
AP_YERR(cen) = INDEFR
# Return if the center is undefined.
if (IS_INDEFR(AP_CXCUR(cen)) || IS_INDEFR(AP_CYCUR(cen))) {
AP_OXINIT(cen) = INDEFR
AP_OYINIT(cen) = INDEFR
AP_OXCENTER(cen) = INDEFR
AP_OYCENTER(cen) = INDEFR
return (AP_CTR_NOAREA)
}
# Convert the coordinates.
switch (AP_WCSOUT(ap)) {
case WCS_WORLD, WCS_PHYSICAL:
call ap_ltoo (ap, AP_CXCUR(cen), AP_CYCUR(cen), AP_OXINIT(cen),
AP_OYINIT(cen), 1)
call ap_ltoo (ap, AP_CXCUR(cen), AP_CYCUR(cen), AP_OXCENTER(cen),
AP_OYCENTER(cen), 1)
case WCS_TV:
call ap_ltov (im, AP_CXCUR(cen), AP_CYCUR(cen), AP_OXINIT(cen),
AP_OYINIT(cen), 1)
call ap_ltov (im, AP_CXCUR(cen), AP_CYCUR(cen), AP_OXCENTER(cen),
AP_OYCENTER(cen), 1)
default:
AP_OXINIT(cen) = AP_CXCUR(cen)
AP_OYINIT(cen) = AP_CYCUR(cen)
AP_OXCENTER(cen) = AP_CXCUR(cen)
AP_OYCENTER(cen) = AP_CYCUR(cen)
}
# Return input coordinates if no center fitting.
if (AP_CENTERFUNCTION(cen) == AP_NONE)
return (AP_OK)
# Choose the centering algorithm.
switch (AP_CENTERFUNCTION(cen)) {
case AP_CENTROID1D:
# Compute the x and y centroids.
if (AP_CTHRESHOLD(cen) <= 0.0) {
fier = ap_mctr1d (Memr[AP_CTRPIX(cen)], AP_CNX(cen),
AP_CNY(cen), AP_EPADU(nse), AP_XCENTER(cen),
AP_YCENTER(cen), AP_XERR(cen), AP_YERR(cen))
if (IS_INDEFR(AP_XERR(cen)))
AP_XCENTER(cen) = AP_CXC(cen)
if (IS_INDEFR(AP_YERR(cen)))
AP_YCENTER(cen) = AP_CYC(cen)
} else {
fier = ap_ctr1d (Memr[AP_CTRPIX(cen)], AP_CNX(cen),
AP_CNY(cen), AP_EPADU(nse), AP_XCENTER(cen),
AP_YCENTER(cen), AP_XERR(cen), AP_YERR(cen))
if (IS_INDEFR(AP_XERR(cen)))
AP_XCENTER(cen) = AP_CXC(cen)
if (IS_INDEFR(AP_YERR(cen)))
AP_YCENTER(cen) = AP_CYC(cen)
}
AP_XCENTER(cen) = AP_XCENTER(cen) + AP_CXCUR(cen) - AP_CXC(cen)
AP_YCENTER(cen) = AP_YCENTER(cen) + AP_CYCUR(cen) - AP_CYC(cen)
AP_XSHIFT(cen) = AP_XCENTER(cen) - AP_CXCUR(cen)
AP_YSHIFT(cen) = AP_YCENTER(cen) - AP_CYCUR(cen)
case AP_GAUSS1D:
# Compute the 1D Gaussian x and y centers.
fier = ap_gctr1d (Memr[AP_CTRPIX(cen)], AP_CNX(cen), AP_CNY(cen),
CONVERT * AP_FWHMPSF(ap) * AP_SCALE(ap), AP_CMAXITER(cen),
AP_XCENTER(cen), AP_YCENTER(cen), AP_XERR(cen), AP_YERR(cen))
AP_XCENTER(cen) = AP_XCENTER(cen) + AP_CXCUR(cen) - AP_CXC(cen)
AP_YCENTER(cen) = AP_YCENTER(cen) + AP_CYCUR(cen) - AP_CYC(cen)
AP_XSHIFT(cen) = AP_XCENTER(cen) - AP_CXCUR(cen)
AP_YSHIFT(cen) = AP_YCENTER(cen) - AP_CYCUR(cen)
case AP_OFILT1D:
# Compute the Goad 1D x and y centers.
fier = ap_lgctr1d (Memr[AP_CTRPIX(cen)], AP_CNX(cen), AP_CNY(cen),
AP_CXC(cen), AP_CYC(cen), CONVERT * AP_FWHMPSF(ap) *
AP_SCALE(ap), AP_CMAXITER(cen), AP_EPADU(nse),
AP_SKYSIGMA(nse), AP_XCENTER(cen), AP_YCENTER(cen),
AP_XERR(cen), AP_YERR(cen))
AP_XCENTER(cen) = AP_XCENTER(cen) + AP_CXCUR(cen) - AP_CXC(cen)
AP_YCENTER(cen) = AP_YCENTER(cen) + AP_CYCUR(cen) - AP_CYC(cen)
AP_XSHIFT(cen) = AP_XCENTER(cen) - AP_CXCUR(cen)
AP_YSHIFT(cen) = AP_YCENTER(cen) - AP_CYCUR(cen)
default:
# do nothing gracefully
}
# Convert the coordinates.
switch (AP_WCSOUT(ap)) {
case WCS_WORLD, WCS_PHYSICAL:
call ap_ltoo (ap, AP_XCENTER(cen), AP_YCENTER(cen),
AP_OXCENTER(cen), AP_OYCENTER(cen), 1)
call ap_ltoo (ap, AP_XCENTER(cen) - AP_XSHIFT(cen),
AP_YCENTER(cen) - AP_YSHIFT(cen), AP_OXINIT(cen),
AP_OYINIT(cen), 1)
AP_OXSHIFT(cen) = AP_OXCENTER(cen) - AP_OXINIT(cen)
AP_OYSHIFT(cen) = AP_OYCENTER(cen) - AP_OYINIT(cen)
case WCS_TV:
call ap_ltov (im, AP_XCENTER(cen), AP_YCENTER(cen),
AP_OXCENTER(cen), AP_OYCENTER(cen), 1)
call ap_ltov (im, AP_XCENTER(cen) - AP_XSHIFT(cen),
AP_YCENTER(cen) - AP_YSHIFT(cen), AP_OXINIT(cen),
AP_OYINIT(cen), 1)
AP_OXSHIFT(cen) = AP_OXCENTER(cen) - AP_OXINIT(cen)
AP_OYSHIFT(cen) = AP_OYCENTER(cen) - AP_OYINIT(cen)
default:
AP_OXCENTER(cen) = AP_XCENTER(cen)
AP_OYCENTER(cen) = AP_YCENTER(cen)
AP_OXINIT(cen) = AP_XCENTER(cen) - AP_XSHIFT(cen)
AP_OYINIT(cen) = AP_YCENTER(cen) - AP_YSHIFT(cen)
AP_OXSHIFT(cen) = AP_XSHIFT(cen)
AP_OYSHIFT(cen) = AP_YSHIFT(cen)
}
# Return appropriate error code.
if (fier != AP_OK) {
AP_XCENTER(cen) = AP_CXCUR(cen)
AP_YCENTER(cen) = AP_CYCUR(cen)
AP_XSHIFT(cen) = 0.0
AP_YSHIFT(cen) = 0.0
AP_XERR(cen) = INDEFR
AP_YERR(cen) = INDEFR
switch (AP_WCSOUT(ap)) {
case WCS_WORLD, WCS_PHYSICAL:
call ap_ltoo (ap, AP_CXCUR(cen), AP_CYCUR(cen), AP_OXINIT(cen),
AP_OYINIT(cen), 1)
call ap_ltoo (ap, AP_CXCUR(cen), AP_CYCUR(cen),
AP_OXCENTER(cen), AP_OYCENTER(cen), 1)
case WCS_TV:
call ap_ltov (im, AP_CXCUR(cen), AP_CYCUR(cen), AP_OXINIT(cen),
AP_OYINIT(cen), 1)
call ap_ltov (im, AP_CXCUR(cen), AP_CYCUR(cen),
AP_OXCENTER(cen), AP_OYCENTER(cen), 1)
default:
AP_OXCENTER(cen) = AP_CXCUR(cen)
AP_OYCENTER(cen) = AP_CYCUR(cen)
AP_OXINIT(cen) = AP_CXCUR(cen)
AP_OYINIT(cen) = AP_CYCUR(cen)
}
AP_OXSHIFT(cen) = 0.0
AP_OYSHIFT(cen) = 0.0
return (fier)
} else if (ier == AP_CTR_BADDATA) {
return (AP_CTR_BADDATA)
} else if (ier == AP_CTR_LOWSNRATIO) {
return (AP_CTR_LOWSNRATIO)
} else if (abs (AP_XSHIFT(cen)) > (AP_MAXSHIFT(cen) * AP_SCALE(ap))) {
return (AP_CTR_BADSHIFT)
} else if (abs (AP_YSHIFT(cen)) > (AP_MAXSHIFT(cen) * AP_SCALE(ap))) {
return (AP_CTR_BADSHIFT)
} else if (ier == AP_CTR_OUTOFBOUNDS) {
return (AP_CTR_OUTOFBOUNDS)
} else
return (AP_OK)
end
|
