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
|
include <gset.h>
include <pkg/gtools.h>
# AP_GTINIT -- Initialize the gtools package for the apphot routines.
pointer procedure ap_gtinit (image, wx, wy)
char image[ARB] # the image name
real wx, wy # center of sky subraster
pointer sp, gt, str
pointer gt_init()
begin
# Allocate working space.
gt = gt_init ()
call smark (sp)
call salloc (str, SZ_LINE, TY_CHAR)
# Set the plot title.
call sprintf (Memc[str], SZ_LINE, "Image: %s: %.2f %.2f\n")
call pargstr (image)
call pargr (wx)
call pargr (wy)
call gt_sets (gt, GTTITLE, Memc[str])
call sfree (sp)
return (gt)
end
# AP_GTFREE -- Free the gtools package.
procedure ap_gtfree (gt)
pointer gt # pointer to gtools structure
begin
call gt_free (gt)
end
# AP_PLOTRAD -- Plot the radial profile of a list of pixels.
procedure ap_plotrad (gd, gt, r, i, npts, polymark)
pointer gd # pointer to graphics stream
pointer gt # the GTOOLS pointer
real r[ARB] # the radii array
real i[ARB] # the intensity array
int npts # number of points
char polymark[ARB] # polyline type
begin
call gt_sets (gt, GTTYPE, "mark")
call gt_sets (gt, GTMARK, polymark)
call gt_plot (gd, gt, r, i, npts)
end
# AP_PLOTPTS -- Plot the radial profile of a list of pixels excluding points
# that are outside the plotting window altogether.
procedure ap_plotpts (gd, gt, r, i, npts, xmin, xmax, ymin, ymax, polymark)
pointer gd # pointer to graphics stream
pointer gt # the GTOOLS pointer
real r[ARB] # the radii array
real i[ARB] # the intensity array
int npts # number of points
real xmin, xmax # the x plot limits
real ymin, ymax # the y plot limits
char polymark[ARB] # polyline type
int j
begin
call gt_sets (gt, GTTYPE, "mark")
call gt_sets (gt, GTMARK, polymark)
do j = 1, npts {
if (r[j] < xmin || r[j] > xmax)
next
if (i[j] < ymin || i[j] > ymax)
next
call gt_plot (gd, gt, r[j], i[j], 1)
}
end
# AP_RSET -- Set up the parameters for the radial profile plot.
procedure ap_rset (gd, gt, xmin, xmax, ymin, ymax, xscale)
pointer gd # pointer to GRAPHICS stream
pointer gt # pointer to GTOOLS structure
real xmin, xmax # min and max of x vector
real ymin, ymax # min and max of y vector
real xscale # image scale
pointer sp, str, title
real vx1, vx2, vy1, vy2, aspect
real gstatr()
begin
call smark (sp)
call salloc (title, 3 * SZ_LINE, TY_CHAR)
call salloc (str, SZ_LINE, TY_CHAR)
# Reset the aspect ratio.
aspect = gstatr (gd, G_ASPECT)
call gsetr (gd, G_ASPECT, 0.75)
# Construct the title.
call sysid (Memc[title], 3 * SZ_LINE)
call strcat ("\n", Memc[title], 3 * SZ_LINE)
call gt_gets (gt, GTTITLE, Memc[str], SZ_LINE)
call strcat (Memc[str], Memc[title], 3 * SZ_LINE)
call strcat ("\n\n", Memc[title], 3 * SZ_LINE)
# Draw three axes.
call gseti (gd, G_XDRAWAXES, 2)
call gswind (gd, xmin / xscale, xmax / xscale, ymin, ymax)
call glabax (gd, Memc[title], "", "Intensity")
# Draw the bottom x axis
call gseti (gd, G_YDRAWAXES, 0)
call gseti (gd, G_XDRAWAXES, 1)
call ggview (gd, vx1, vx2, vy1, vy2)
call gsview (gd, vx1, vx2, vy1, vy2)
call gswind (gd, xmin, xmax, ymin, ymax)
call glabax (gd, "",
"Radial Distance (lower-pixels, upper-scale units)", "")
# Reset to standard gio parameters.
call gseti (gd, G_YDRAWAXES, 3)
call gseti (gd, G_XDRAWAXES, 3)
call gsetr (gd, G_ASPECT, aspect)
call sfree (sp)
end
|
