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 "../lib/apphotdef.h"
include "../lib/fitskydef.h"
include "../lib/fitsky.h"
# APSINIT - Procedure to initialize the sky fitting structure and parameters.
procedure apsinit (ap, function, annulus, dannulus, fwhmpsf, noise)
pointer ap # pointer to the apphot structure
int function # sky fitting algorithm
real annulus # radius of sky annulus
real dannulus # width of sky annulus
real fwhmpsf # FWHM of the PSF
int noise # noise function
begin
# Set up the object parameters.
call malloc (ap, LEN_APSTRUCT, TY_STRUCT)
# Set up the global apphot package parameters.
call ap_defsetup (ap, fwhmpsf)
# Set up the noise model parameters.
call ap_noisesetup (ap, noise)
# Set up the display options.
call ap_dispsetup (ap)
# Initialize the sky fitting parameters.
call ap_skysetup (ap, function, annulus, dannulus)
# Unused structures are set to null.
AP_PCENTER(ap) = NULL
AP_PPHOT(ap) = NULL
AP_POLY(ap) = NULL
AP_PPSF(ap) = NULL
AP_RPROF(ap) = NULL
end
# AP_SKYSETUP -- Procedure to set up the sky fitting arrays and parameters.
procedure ap_skysetup (ap, function, annulus, dannulus)
pointer ap # pointer to apphot structure
int function # sky fitting function
real annulus # inner radius of sky annulus
real dannulus # outer radius of sky annulus
pointer sky
begin
call malloc (AP_PSKY(ap), LEN_SKYSTRUCT, TY_STRUCT)
sky = AP_PSKY(ap)
AP_SXCUR(sky) = INDEFR
AP_SYCUR(sky) = INDEFR
# Initialize the sky fitting parameters.
AP_SKYFUNCTION(sky) = function
switch (function) {
case AP_CONSTANT:
call strcpy ("constant", AP_SSTRING(sky), SZ_FNAME)
case AP_MODE:
call strcpy ("mode", AP_SSTRING(sky), SZ_FNAME)
case AP_CENTROID:
call strcpy ("centroid", AP_SSTRING(sky), SZ_FNAME)
case AP_SKYFILE:
call strcpy ("file", AP_SSTRING(sky), SZ_FNAME)
case AP_HISTPLOT:
call strcpy ("histplot", AP_SSTRING(sky), SZ_FNAME)
case AP_RADPLOT:
call strcpy ("radplot", AP_SSTRING(sky), SZ_FNAME)
case AP_MEDIAN:
call strcpy ("median", AP_SSTRING(sky), SZ_FNAME)
case AP_GAUSS:
call strcpy ("gauss", AP_SSTRING(sky), SZ_FNAME)
case AP_OFILT:
call strcpy ("ofilt", AP_SSTRING(sky), SZ_FNAME)
case AP_CROSSCOR:
call strcpy ("crosscor", AP_SSTRING(sky), SZ_FNAME)
case AP_MEAN:
call strcpy ("mean", AP_SSTRING(sky), SZ_FNAME)
default:
AP_SKYFUNCTION(sky) = DEF_SKYFUNCTION
call strcpy ("mode", AP_SSTRING(sky), SZ_FNAME)
}
AP_SKYBACKGROUND(sky) = DEF_SKYVALUE
AP_ANNULUS(sky) = annulus
AP_DANNULUS(sky) = dannulus
AP_K1(sky) = DEF_K1
AP_BINSIZE(sky) = DEF_BINSIZE
AP_SMOOTH(sky) = DEF_SMOOTH
AP_SLOCLIP(sky) = DEF_SLOCLIP
AP_SHICLIP(sky) = DEF_SHICLIP
AP_SMAXITER(sky) = DEF_SMAXITER
AP_RGROW(sky) = DEF_RGROW
AP_SNREJECT(sky) = DEF_SNREJECT
AP_SLOREJECT(sky) = DEF_SLOREJECT
AP_SHIREJECT(sky) = DEF_SHIREJECT
# Initialize the sky pixel buffers.
AP_LENSKYBUF(sky) = 0
AP_NSKYPIX(sky) = 0
AP_SKYPIX(sky) = NULL
AP_INDEX(sky) = NULL
AP_COORDS(sky) = NULL
AP_SWGT(sky) = NULL
# Initialize results parameters.
AP_SKY_MODE(sky) = INDEFR
AP_SKY_SIG(sky) = INDEFR
AP_SKY_SKEW(sky) = INDEFR
AP_NSKY(sky) = 0
AP_NSKY_REJECT(sky) = 0
end
|
