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
|
include <mach.h>
include <pkg/gtools.h>
include "identify.h"
# T_IDENTIFY -- Identify features
procedure t_identify ()
int list, clscan(), clgeti(), clgwrd(), nscan(), imtopenp(), imtgetim()
real clgetr()
pointer sp, str, id, gt_init(), un_open()
begin
call smark (sp)
call salloc (str, SZ_LINE, TY_CHAR)
# Allocate the basic data structure.
call id_init (id)
# Get task parameters.
list = imtopenp ("images")
if (clscan ("nsum") != EOF) {
call gargi (ID_NSUM(id,1))
call gargi (ID_NSUM(id,2))
if (nscan() == 0)
call error (1, "Error in 'nsum' parameter")
if (nscan() == 1)
ID_NSUM(id,2) = ID_NSUM(id,1)
ID_NSUM(id,1) = max (1, ID_NSUM(id,1))
ID_NSUM(id,2) = max (1, ID_NSUM(id,2))
}
ID_MAXFEATURES(id) = clgeti ("maxfeatures")
ID_MINSEP(id) = clgetr ("minsep")
ID_MATCH(id) = clgetr ("match")
ID_ZWIDTH(id) = clgetr ("zwidth")
ID_FTYPE(id) = clgwrd ("ftype", Memc[str], SZ_LINE, FTYPES)
ID_FWIDTH(id) = clgetr ("fwidth")
ID_CRADIUS(id) = clgetr ("cradius")
ID_THRESHOLD(id) = clgetr ("threshold")
call clgstr ("database", ID_DATABASE(id), ID_LENSTRING)
call clgstr ("coordlist", ID_COORDLIST(id), ID_LENSTRING)
call clgstr ("units", Memc[str], SZ_LINE)
call xt_stripwhite (Memc[str])
if (Memc[str] != EOS)
ID_UN(id) = un_open (Memc[str])
ID_LABELS(id) = 1
# Initialize features data structure.
ID_GT(id) = gt_init()
call gt_sets (ID_GT(id), GTTYPE, "line")
ID_CV(id) = NULL
ID_CURRENT(id) = 0
ID_SHIFT(id) = 0.
# Initialize ICFIT
call ic_open (ID_IC(id))
call clgstr ("function", Memc[str], SZ_LINE)
call ic_pstr (ID_IC(id), "function", Memc[str])
call ic_puti (ID_IC(id), "order", clgeti ("order"))
call clgstr ("sample", Memc[str], SZ_LINE)
call ic_pstr (ID_IC(id), "sample", Memc[str])
call ic_puti (ID_IC(id), "naverage", 1)
call ic_puti (ID_IC(id), "niterate", clgeti ("niterate"))
call ic_putr (ID_IC(id), "low", clgetr ("low_reject"))
call ic_putr (ID_IC(id), "high", clgetr ("high_reject"))
call ic_putr (ID_IC(id), "grow", clgetr ("grow"))
call ic_pstr (ID_IC(id), "xlabel", "Feature positions")
call ic_pstr (ID_IC(id), "xunits", "pixels")
call ic_pstr (ID_IC(id), "ylabel", "")
call ic_pkey (ID_IC(id), 1, 'y', 'x')
call ic_pkey (ID_IC(id), 2, 'y', 'v')
call ic_pkey (ID_IC(id), 3, 'y', 'r')
call ic_pkey (ID_IC(id), 4, 'y', 'd')
call ic_pkey (ID_IC(id), 5, 'y', 'n')
call ic_puti (ID_IC(id), "key", 3)
# Get the line list.
call id_mapll (id)
# Expand the image template and identify features in each image.
while (imtgetim (list, ID_IMAGE(id), ID_LENSTRING) != EOF)
call id_identify (id)
# Finish up.
call smw_daxis (NULL, NULL, 0, 0, 0)
call id_free (id)
call imtclose (list)
call sfree (sp)
end
|
