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
|
include "../lib/apphotdef.h"
include "../lib/centerdef.h"
include "../lib/center.h"
# APCINIT - Procedure to initialize the centering structure.
procedure apcinit (ap, function, cbox, fwhmpsf, noise)
pointer ap # pointer to the apphot structure
int function # centering algorithm
real cbox # half width of centering box
real fwhmpsf # FWHM of the PSF
int noise # noise function
begin
# Allocate space.
call malloc (ap, LEN_APSTRUCT, TY_STRUCT)
# Set up the global apphot package parameters.
call ap_defsetup (ap, fwhmpsf)
# Setup the noise structure.
call ap_noisesetup (ap, noise)
# Initialize the centering parameters.
call ap_ctrsetup (ap, function, cbox)
# Display the options.
call ap_dispsetup (ap)
# Unused structures are set to null.
AP_PSKY(ap) = NULL
AP_PPSF(ap) = NULL
AP_PPHOT(ap) = NULL
AP_POLY(ap) = NULL
AP_RPROF(ap) = NULL
end
# AP_CTRSETUP -- Procedure to setup the centering array parameters.
procedure ap_ctrsetup (ap, function, cbox)
pointer ap # pointer to apphot structure
int function # centering function
real cbox # radius of centering aperture
pointer ctr
begin
# Allocate space for the centering structure.
call malloc (AP_PCENTER(ap), LEN_CENSTRUCT, TY_STRUCT)
ctr = AP_PCENTER(ap)
AP_CXCUR(ctr) = INDEFR
AP_CYCUR(ctr) = INDEFR
AP_CENTERFUNCTION(ctr) = function
switch (function) {
case AP_CENTROID1D:
call strcpy ("centroid", AP_CSTRING(ctr), SZ_FNAME)
case AP_GAUSS1D:
call strcpy ("gauss", AP_CSTRING(ctr), SZ_FNAME)
case AP_NONE:
call strcpy ("none", AP_CSTRING(ctr), SZ_FNAME)
case AP_OFILT1D:
call strcpy ("ofilter", AP_CSTRING(ctr), SZ_FNAME)
default:
AP_CENTERFUNCTION(ctr) = DEF_CENTERFUNCTION
call strcpy ("centroid", AP_CSTRING(ctr), SZ_FNAME)
}
AP_CAPERT(ctr) = cbox
AP_CTHRESHOLD(ctr) = DEF_CTHRESHOLD
AP_MINSNRATIO(ctr) = DEF_MINSNRATIO
AP_CMAXITER(ctr) = DEF_CMAXITER
AP_MAXSHIFT(ctr) = DEF_MAXSHIFT
AP_CLEAN(ctr) = DEF_CLEAN
AP_RCLEAN(ctr) = DEF_RCLEAN
AP_RCLIP(ctr) = DEF_RCLIP
AP_SIGMACLEAN(ctr) = DEF_CLEANSIGMA
AP_NCTRPIX(ctr) = 0
AP_CTRPIX(ctr) = NULL
AP_XCTRPIX(ctr) = NULL
AP_YCTRPIX(ctr) = NULL
AP_OXINIT(ctr) = INDEFR
AP_OYINIT(ctr) = INDEFR
AP_XCENTER(ctr) = INDEFR
AP_YCENTER(ctr) = INDEFR
AP_OXCENTER(ctr) = INDEFR
AP_OYCENTER(ctr) = INDEFR
AP_XERR(ctr) = INDEFR
AP_YERR(ctr) = INDEFR
end
|
