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
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
|
include <gset.h>
include <error.h>
include <smw.h>
include "identify.h"
# List of colon commands.
define CMDS "|show|features|image|nsum|database|read|write|add|coordlist|match\
|maxfeatures|minsep|zwidth|labels|fwidth|ftype|cradius|threshold|"
define SHOW 1 # Show parameters
define FEATURES 2 # Show list of features
define IMAGE 3 # Set new image
define NSUM 4 # Set the number of lines or columns to sum
define DATABASE 5 # Set new database
define READ 6 # Read database entry
define WRITE 7 # Write database entry
define ADD 8 # Add features from database
define COORDLIST 9 # Set new coordinate list
define MATCH 10 # Set coordinate list matching distance
define MAXFEATURES 11 # Set maximum number of features for auto find
define MINSEP 12 # Set minimum separation distance
define ZWIDTH 13 # Set zoom window width
define LABEL 14 # Set label type
define WIDTH 15 # Set centering width
define TYPE 16 # Set centering type
define RADIUS 17 # Set centering radius
define THRESHOLD 18 # Set the centering threshold
# ID_COLON -- Respond to colon command.
procedure id_colon (id, cmdstr, newimage, prfeature)
pointer id # ID pointer
char cmdstr[ARB] # Colon command
char newimage[ARB] # New image name
int prfeature # Print current feature on status line
char cmd[SZ_LINE]
int i, ncmd, ival[2]
real rval[2]
pointer im
int nscan(), strdic()
pointer immap()
errchk immap, id_dbread, id_dbwrite, id_log
begin
# Scan the command string and get the first word.
call sscan (cmdstr)
call gargwrd (cmd, SZ_LINE)
ncmd = strdic (cmd, cmd, SZ_LINE, CMDS)
switch (ncmd) {
case SHOW: # :show - show values of parameters
call gargwrd (cmd, SZ_LINE)
if (nscan() == 1) {
call gdeactivate (ID_GP(id), AW_CLEAR)
call id_show (id, "STDOUT")
call greactivate (ID_GP(id), AW_PAUSE)
} else {
iferr (call id_show (id, cmd)) {
call erract (EA_WARN)
prfeature = NO
}
}
case FEATURES: # :features - list features
call gargwrd (cmd, SZ_LINE)
if (nscan() == 1) {
call gdeactivate (ID_GP(id), AW_CLEAR)
call id_log (id, "STDOUT", NULL)
call greactivate (ID_GP(id), AW_PAUSE)
} else {
iferr (call id_log (id, cmd, NULL)) {
call erract (EA_WARN)
prfeature = NO
}
}
case IMAGE: # :image - set image to identify
call gargwrd (cmd, SZ_LINE)
if (nscan() == 1) {
call printf ("image %s\n")
call pargstr (Memc[ID_IMAGE(id)])
prfeature = NO
} else {
call strcpy (cmd, newimage, SZ_FNAME)
iferr {
im = immap (newimage, READ_ONLY, 0)
call imunmap (im)
} then {
newimage[1] = EOS
call erract (EA_WARN)
prfeature = NO
}
}
case NSUM: # :nsum - set number of lines or columns to sum in image
call gargi (ival[1])
if (nscan() == 1) {
call printf ("nsum %d %d\n")
call pargi (ID_NSUM(id,1))
call pargi (ID_NSUM(id,2))
prfeature = NO
} else {
ID_NSUM(id,1) = ival[1]
call gargi (ival[2])
if (nscan() == 3)
ID_NSUM(id,2) = ival[2]
call smw_daxis (NULL, NULL, SMW_PAXIS(MW(ID_SH(id)),1),
ID_NSUM(id,1), ID_NSUM(id,2))
}
case DATABASE: # :database - set database
call gargwrd (cmd, SZ_LINE)
if (nscan() == 1) {
call printf ("database %s\n")
call pargstr (Memc[ID_DATABASE(id)])
prfeature = NO
} else {
call strcpy (cmd, Memc[ID_DATABASE(id)], SZ_FNAME)
ID_NEWDBENTRY(id) = YES
}
case READ: # :read - read database entry
prfeature = NO
iferr {
call gargwrd (cmd, SZ_LINE)
if (nscan() == 1)
call id_dbread (id, Memc[ID_IMAGE(id)], ID_AP(id,1),
NO, YES)
else {
call gargi (ival[1])
if (nscan() < 3)
ival[1] = ID_AP(id,1)
call gargi (ival[2])
if (nscan() < 4)
ival[2] = ID_AP(id,2)
call id_dbread (id, cmd, ival, NO, YES)
}
} then
call erract (EA_WARN)
case WRITE: # :write - write database entry
prfeature = NO
iferr {
ival[1] = ID_AP(id,1)
ival[2] = ID_AP(id,2)
call gargwrd (cmd, SZ_LINE)
if (nscan() == 1)
call id_dbwrite (id, Memc[ID_IMAGE(id)], ival, YES)
else {
call gargi (ival[1])
if (nscan() < 3)
ival[1] = ID_AP(id,1)
call gargi (ival[2])
if (nscan() < 4)
ival[2] = ID_AP(id,2)
call id_dbwrite (id, cmd, ival, YES)
}
} then
call erract (EA_WARN)
case ADD: # :add - add features from database entry
prfeature = NO
iferr {
call gargwrd (cmd, SZ_LINE)
if (nscan() == 1)
call id_dbread (id, Memc[ID_IMAGE(id)], ID_AP(id,1),
YES, YES)
else {
call gargi (ival[1])
if (nscan() < 3)
ival[1] = ID_AP(id,1)
call gargi (ival[2])
if (nscan() < 4)
ival[2] = ID_AP(id,2)
call id_dbread (id, cmd, ival, YES, YES)
}
} then
call erract (EA_WARN)
case COORDLIST: # :coordlist - set coordinate list
call gargwrd (cmd, SZ_LINE)
if (nscan() == 1) {
call printf ("coordlist %s\n")
call pargstr (Memc[ID_COORDLIST(id)])
prfeature = NO
} else {
call strcpy (cmd, Memc[ID_COORDLIST(id)], SZ_FNAME)
call id_unmapll (id)
call id_mapll (id)
}
case MATCH: # :match - set matching distance for coordinate list
call gargr (rval[1])
if (nscan() == 1) {
call printf ("match %g\n")
call pargr (ID_MATCH(id))
prfeature = NO
} else
ID_MATCH(id) = rval[1]
case MAXFEATURES: # :maxfeatures - set max num features for auto find
call gargi (ival[1])
if (nscan() == 1) {
call printf ("maxfeatures %d\n")
call pargi (ID_MAXFEATURES(id))
prfeature = NO
} else
ID_MAXFEATURES(id) = ival[1]
case MINSEP: # :minsep - set minimum feature separation allowed
call gargr (rval[1])
if (nscan() == 1) {
call printf ("minsep %g\n")
call pargr (ID_MINSEP(id))
prfeature = NO
} else
ID_MINSEP(id) = rval[1]
case ZWIDTH: # :zwidth - set zoom window width
call gargr (rval[1])
if (nscan() == 1) {
call printf ("zwidth %g\n")
call pargr (ID_ZWIDTH(id))
prfeature = NO
} else {
ID_ZWIDTH(id) = rval[1]
if (ID_GTYPE(id) == 2)
ID_NEWGRAPH(id) = YES
}
case LABEL: # :labels - set label type
call gargwrd (cmd, SZ_LINE)
if (nscan() == 1) {
switch (ID_LABELS(id)) {
case 2:
call printf ("labels index\n")
case 3:
call printf ("labels pixel\n")
case 4:
call printf ("labels coords\n")
case 5:
call printf ("labels user\n")
case 6:
call printf ("labels both\n")
default:
call printf ("labels none\n")
}
prfeature = NO
} else {
ID_LABELS(id) = strdic (cmd, cmd, SZ_LINE, LABELS)
do i = 1, ID_NFEATURES(id)
call id_mark (id, i)
}
case WIDTH: # :fwidth - set centering width
call gargr (rval[1])
if (nscan() == 1) {
call printf ("fwidth %g\n")
call pargr (ID_FWIDTH(id))
prfeature = NO
} else
ID_FWIDTH(id) = rval[1]
case TYPE: # :ftype - set centering type
call gargwrd (cmd, SZ_LINE)
if (nscan() == 1) {
switch (ID_FTYPE(id)) {
case EMISSION:
call printf ("ftype emission\n")
case ABSORPTION:
call printf ("ftype absorption\n")
case GEMISSION:
call printf ("ftype gemission\n")
case GABSORPTION:
call printf ("ftype gabsorption\n")
}
prfeature = NO
} else
ID_FTYPE(id) = strdic (cmd, cmd, SZ_LINE, FTYPES)
case RADIUS: # :cradius - set centering radius
call gargr (rval[1])
if (nscan() == 1) {
call printf ("cradius %g\n")
call pargr (ID_CRADIUS(id))
prfeature = NO
} else
ID_CRADIUS(id) = rval[1]
case THRESHOLD: # :threshold - set centering threshold
call gargr (rval[1])
if (nscan() == 1) {
call printf ("threshold %g\n")
call pargr (ID_THRESHOLD(id))
prfeature = NO
} else
ID_THRESHOLD(id) = rval[1]
default:
call printf ("Unrecognized or ambiguous command\007")
prfeature = NO
}
end
|
