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
|
include "../lib/daophotdef.h"
include "../lib/psfdef.h"
# DP_PLOTPSF -- Plot the psf using the default plot type.
procedure dp_plotpsf (dao, im, subrast, ncols, nlines, x1, y1, gd)
pointer dao # pointer to the daophot structure
pointer im # the input image descriptor
real subrast[ncols,nlines] # image subraster
int ncols, nlines # dimensions of the subraster
int x1, y1 # coordinates of the lower left corner
pointer gd # pointer to the graphics stream
real tx, ty
pointer sp, title
real dp_pstatr()
int dp_pstati()
begin
# Return if the graphics stream is undefined.
if (gd == NULL)
return
# Comvert the coordinates if necessary.
call dp_wout (dao, im, dp_pstatr(dao, CUR_PSFX), dp_pstatr(dao,
CUR_PSFY), tx, ty, 1)
# Construct the title.
call smark (sp)
call salloc (title, SZ_LINE, TY_CHAR)
call sprintf (Memc[title], SZ_LINE, "Star: %d X: %g Y: %g Mag: %g")
call pargi (dp_pstati (dao, CUR_PSFID))
call pargr (tx)
call pargr (ty)
call pargr (dp_pstatr (dao, CUR_PSFMAG))
# Initialize plot.
if (dp_pstati (dao, PLOTTYPE) == PSF_MESHPLOT)
call dp_surfpsf (dao, subrast, ncols, nlines, Memc[title], gd)
else if (dp_pstati (dao, PLOTTYPE) == PSF_CONTOURPLOT)
call dp_contpsf (dao, subrast, ncols, nlines, Memc[title], gd)
else if (dp_pstati (dao, PLOTTYPE) == PSF_RADIALPLOT)
call dp_radpsf (dao, subrast, ncols, nlines, x1, y1,
Memc[title], gd)
call gdeactivate (gd, 0)
call sfree (sp)
end
|
