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
|
include <gset.h>
# List of colon commands
define CMDS "|show|function|xorder|yorder|corners|"
define SHOW 1 # Show parameters
define FUNCTION 2 # Set or show function type
define XORDER 3 # Set or show x order of function
define YORDER 4 # Set or show y order of function
define CORNERS 5 # Show corners
# IGS_COLON -- Processes colon commands.
procedure igs_colon (cmdstr, gp, sf)
char cmdstr[ARB] # Command string
pointer gp # GIO pointer
pointer sf # Surface pointer
char cmd[SZ_LINE]
int ncmd, ival
int nscan(), strdic()
real xgseval()
string funcs "|chebyshev|legendre|"
include "igsfit.com"
begin
# Use formated scan to parse the command string.
# The first word is the command and it may be minimum match
# abbreviated with the list of commands.
call sscan (cmdstr)
call gargwrd (cmd, SZ_LINE)
ncmd = strdic (cmd, cmd, SZ_LINE, CMDS)
switch (ncmd) {
case SHOW: # :show - Show the values of the fitting parameters.
call gdeactivate (gp, AW_CLEAR)
call printf ("function %s\n")
call pargstr (function)
call printf ("xorder %d\n")
call pargi (xorder)
call printf ("yorder %d\n")
call pargi (yorder)
call printf ("Fitted coordinates at the corners of the images:\n")
call printf (" (%d, %d) = %g (%d, %d) = %g\n")
call pargr (xmin)
call pargr (ymin)
call pargr (xgseval (sf, xmin, ymin))
call pargr (xmax)
call pargr (ymin)
call pargr (xgseval (sf, xmax, xmin))
call printf (" (%d, %d) = %g (%d, %d) = %g\n")
call pargr (xmin)
call pargr (ymax)
call pargr (xgseval (sf, xmin, ymax))
call pargr (xmax)
call pargr (ymax)
call pargr (xgseval (sf, xmax, ymax))
call printf ("rms %g\n")
call pargr (rms)
call greactivate (gp, AW_PAUSE)
case FUNCTION: # :function - List or set the fitting function.
call gargwrd (cmd, SZ_LINE)
if (nscan() == 1) {
call printf ("function = %s\n")
call pargstr (function)
} else {
if (strdic (cmd, cmd, SZ_LINE, funcs) > 0)
call strcpy (cmd, function, SZ_LINE)
else
call printf ("Unknown or ambiguous function\n")
}
case XORDER: # xorder: List or set the function order.
call gargi (ival)
if (nscan() == 1) {
call printf ("xorder %d\n")
call pargi (xorder)
} else if (ival < 2)
call printf ("xorder must be at least 2\n")
else
xorder = ival
case YORDER: # yorder: List or set the function order.
call gargi (ival)
if (nscan() == 1) {
call printf ("yorder %d\n")
call pargi (yorder)
} else if (ival < 2)
call printf ("yorder must be at least 2\n")
else
yorder = ival
case CORNERS: # corners: List coordinates at corners.
call printf ("(%d,%d)=%g (%d,%d)=%g (%d,%d)=%g (%d,%d)=%g\n")
call pargr (xmin)
call pargr (ymin)
call pargr (xgseval (sf, xmin, ymin))
call pargr (xmax)
call pargr (ymin)
call pargr (xgseval (sf, xmax, xmin))
call pargr (xmin)
call pargr (ymax)
call pargr (xgseval (sf, xmin, ymax))
call pargr (xmax)
call pargr (ymax)
call pargr (xgseval (sf, xmax, ymax))
default:
call printf ("Unrecognized or ambiguous command\007")
}
end
|
