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
|
include <mach.h>
include "../lib/apphotdef.h"
include "../lib/apphot.h"
include "../lib/noisedef.h"
include "../lib/photdef.h"
include "../lib/phot.h"
# APREMAG -- Procedure to recompute the magnitudes inside a set of apertures
# given that the sums and effective areas have already been computed.
int procedure apremag (ap, im, positive, skyval, skysig, nsky)
pointer ap # pointer to the apphot structure
pointer im # the input image descriptor
int positive # emission and absorption features
real skyval # sky value
real skysig # sigma of sky
int nsky # number of sky pixels
int nap
pointer nse, phot
real zmag
begin
# Initalize.
phot = AP_PPHOT(ap)
nse = AP_NOISE(ap)
call amovkr (INDEFR, Memr[AP_MAGS(phot)], AP_NAPERTS(phot))
call amovkr (INDEFR, Memr[AP_MAGERRS(phot)], AP_NAPERTS(phot))
if (IS_INDEFR(AP_PXCUR(phot)) || IS_INDEFR(AP_PYCUR(phot))) {
AP_OPXCUR(phot) = AP_PXCUR(phot)
AP_OPYCUR(phot) = AP_PYCUR(phot)
} else {
switch (AP_WCSOUT(ap)) {
case WCS_WORLD, WCS_PHYSICAL:
call ap_ltoo (ap, AP_PXCUR(phot), AP_PYCUR(phot),
AP_OPXCUR(phot), AP_OPYCUR(phot), 1)
case WCS_TV:
call ap_ltov (im, AP_PXCUR(phot), AP_PYCUR(phot),
AP_OPXCUR(phot), AP_OPYCUR(phot), 1)
default:
AP_OPXCUR(phot) = AP_PXCUR(phot)
AP_OPYCUR(phot) = AP_PYCUR(phot)
}
}
# Check for errors.
if (IS_INDEFR(AP_PXCUR(phot)) || IS_INDEFR(AP_PYCUR(phot)))
return (AP_APERT_NOAPERT)
if (IS_INDEFR(skyval))
return (AP_APERT_NOSKYMODE)
nap = min (AP_NMINAP(phot) - 1, AP_NMAXAP(phot))
# Compute the magnitudes and errors.
if (positive == YES)
call apcopmags (Memd[AP_SUMS(phot)], Memd[AP_AREA(phot)],
Memr[AP_MAGS(phot)], Memr[AP_MAGERRS(phot)], nap,
skyval, skysig, nsky, AP_ZMAG(phot), AP_NOISEFUNCTION(nse),
AP_EPADU(nse))
else
call apconmags (Memd[AP_SUMS(phot)], Memd[AP_AREA(phot)],
Memr[AP_MAGS(phot)], Memr[AP_MAGERRS(phot)], nap,
skyval, skysig, nsky, AP_ZMAG(phot), AP_NOISEFUNCTION(nse),
AP_EPADU(nse), AP_READNOISE(nse))
# Correct for itime.
zmag = 2.5 * log10 (AP_ITIME(ap))
call aaddkr (Memr[AP_MAGS(phot)], zmag, Memr[AP_MAGS(phot)], nap)
if (AP_NMAXAP(phot) < AP_NAPERTS(phot))
return (AP_APERT_OUTOFBOUNDS)
else if (AP_NMINAP(phot) <= AP_NMAXAP(phot))
return (AP_APERT_BADDATA)
else
return (AP_OK)
end
|
