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
|
include "../lib/apphot.h"
include "../lib/display.h"
include "../lib/center.h"
include "../lib/fitsky.h"
define HELPFILE "apphot$phot/iphot.key"
# AP_RADSETUP -- Procedure to set up phot interactively using a radial profile
# plot of a bright star.
procedure ap_radsetup (ap, im, wx, wy, gd, out, stid)
pointer ap # pointer to apphot structure
pointer im # pointer to the IRAF image
real wx, wy # cursor coordinates
pointer gd # pointer to graphics stream
int out # output file descriptor
int stid # output file sequence number
int cier, sier, pier, key, wcs
pointer sp, cmd, str
real xcenter, ycenter, xc, yc, rmin, rmax, imin, imax
real u1, u2, v1, v2, x1, x2, y1, y2, rval
int apfitcenter(), apfitsky(), ap_wmag(), apstati(), clgcur(), ap_showplot()
real apstatr(), ap_cfwhmpsf(), ap_ccapert(), ap_cannulus(), ap_cdannulus()
real ap_csigma(), ap_crgrow(), ap_crclean(), ap_crclip()
real ap_cdatamin(), ap_cdatamax()
begin
# Check for open graphics stream
if (gd == NULL)
return
call greactivate (gd, 0)
call gclear (gd)
# Store the viewport and window coordinates.
call ggview (gd, u1, u2, v1, v2)
call ggwind (gd, x1, x2, y1, y2)
# Make the plot.
if (ap_showplot (ap, im, wx, wy, gd, xcenter, ycenter, rmin, rmax,
imin, imax) == ERR) {
call gdeactivate (gd, 0)
return
}
# Initialize.
call smark (sp)
call salloc (cmd, SZ_LINE, TY_CHAR)
call salloc (str, SZ_LINE, TY_CHAR)
call printf (
"Waiting for setup menu command (?=help, v=default setup, q=quit):\n")
while (clgcur ("gcommands", xc, yc, wcs, key, Memc[cmd],
SZ_LINE) != EOF) {
switch (key) {
case 'q':
break
case '?':
call gpagefile (gd, HELPFILE, "")
case 'f':
rval = ap_cfwhmpsf (ap, gd, out, stid, rmin, rmax, imin, imax)
case 's':
rval = ap_csigma (ap, gd, out, stid, rmin, rmax, imin, imax)
case 'l':
rval = ap_cdatamin (ap, gd, out, stid, rmin, rmax, imin, imax)
case 'u':
rval = ap_cdatamax (ap, gd, out, stid, rmin, rmax, imin, imax)
case 'c':
rval = ap_ccapert (ap, gd, out, stid, rmin, rmax, imin, imax)
case 'n':
rval = ap_crclean (ap, gd, out, stid, rmin, rmax, imin, imax)
case 'p':
rval = ap_crclip (ap, gd, out, stid, rmin, rmax, imin, imax)
case 'a':
rval = ap_cannulus (ap, gd, out, stid, rmin, rmax, imin, imax)
case 'd':
rval = ap_cdannulus (ap, gd, out, stid, apstatr (ap, ANNULUS),
rmin, rmax, imin, imax)
case 'g':
rval = ap_crgrow (ap, gd, out, stid, rmin, rmax, imin, imax)
case 'r':
call ap_caper (ap, gd, out, stid, Memc[str], rmin, rmax,
imin, imax)
case 'v':
rval = ap_cfwhmpsf (ap, gd, out, stid, rmin, rmax, imin, imax)
rval = ap_ccapert (ap, gd, out, stid, rmin, rmax, imin, imax)
rval = ap_csigma (ap, gd, out, stid, rmin, rmax, imin, imax)
#rval = ap_ccthresh (ap, gd, out, stid, rmin, rmax, imin, imax)
rval = ap_cannulus (ap, gd, out, stid, rmin, rmax, imin, imax)
rval = ap_cdannulus (ap, gd, out, stid, apstatr (ap, ANNULUS),
rmin, rmax, imin, imax)
call ap_caper (ap, gd, out, stid, Memc[str], rmin, rmax,
imin, imax)
default:
call printf ("Unknown or ambiguous keystroke command\007\n")
}
call printf (
"Waiting for setup menu command (?=help, v=default setup, q=quit):\n")
}
call printf (
"Interactive setup is complete. Type w to save parameters.\n")
# Restore the viewport and window coordinates.
call gsview (gd, u1, u2, v1, v2)
call gswind (gd, x1, x2, y1, y2)
call gdeactivate (gd, 0)
call sfree (sp)
# Compute and print the answer.
cier = apfitcenter (ap, im, xcenter, ycenter)
if (! IS_INDEFR (apstatr (ap, XCENTER)) &&
! IS_INDEFR (apstatr (ap, YCENTER)))
sier = apfitsky (ap, im, apstatr (ap, XCENTER), apstatr (ap,
YCENTER), NULL, gd)
if (! IS_INDEFR (apstatr (ap, SKY_MODE)))
pier = ap_wmag (ap, im, apstatr (ap, XCENTER), apstatr (ap,
YCENTER), apstati (ap, POSITIVE), apstatr (ap, SKY_MODE),
apstatr (ap, SKY_SIGMA), apstati (ap, NSKY))
call ap_pplot (ap, im, 0, gd, apstati (ap, RADPLOTS))
call ap_qpmag (ap, cier, sier, pier)
end
|
