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
|
include "../lib/apphot.h"
include "../lib/polyphot.h"
# AP_YWRITE -- Procedure to write a polygon to a file.
procedure ap_ywrite (py, im, cl, pl, x, y, nvertices, cid, pid, firstpoly,
newpoly)
pointer py # pointer to the apphot structure
pointer im # the input image descriptor
int cl # pointer to the coordinate list file
int pl # pointer to the polygon list file
real x[ARB] # x coordinates of the vertices
real y[ARB] # y coordinates of the vertices
int nvertices # number of vertices
int cid # coordinate list index
int pid # polygon list index
int firstpoly # first polygon measured
int newpoly # new polygon
real xtemp, ytemp
int i
real apstatr()
int apstati()
begin
# Make sure the output files are at EOF.
if (pl == NULL)
return
call seek (pl, EOF)
if (cl != NULL)
call seek (cl, EOF)
if (newpoly == YES) {
# Terminate the coord list that belongs with the first polygon.
if (firstpoly == NO && cl != NULL)
call fprintf (cl, ";\n")
# Write out the coordinates.
do i = 1, nvertices {
switch (apstati(py, WCSOUT)) {
case WCS_WORLD, WCS_PHYSICAL:
call ap_ltoo (py, x[i], y[i], xtemp, ytemp, 1)
case WCS_TV:
call ap_ltov (im, x[i], y[i], xtemp, ytemp, 1)
default:
xtemp = x[i]
ytemp = y[i]
}
call fprintf (pl, "%g %g\n")
call pargr (xtemp)
call pargr (ytemp)
}
if (nvertices > 0)
call fprintf (pl, ";\n")
pid = pid + 1
# Reset polygon parameters.
newpoly = NO
if (firstpoly == YES)
firstpoly = NO
}
# Write out the central coordinates of the polygon.
if (firstpoly == NO && cl != NULL) {
switch (apstati(py, WCSOUT)) {
case WCS_WORLD, WCS_PHYSICAL:
call ap_ltoo (py, apstatr(py,CWX), apstatr(py,CWY), xtemp,
ytemp, 1)
case WCS_TV:
call ap_ltov (im, apstatr(py,CWX), apstatr(py,CWY), xtemp,
ytemp, 1)
default:
xtemp = apstatr (py, CWX)
ytemp = apstatr (py, CWY)
}
call fprintf (cl, "%g %g\n")
call pargr (xtemp)
call pargr (ytemp)
cid = cid + 1
}
# Flush the output files.
if (pl != NULL)
call flush (pl)
if (cl != NULL)
call flush (cl)
end
|
