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
|
include <error.h>
include "refspectra.h"
# REFTABLE -- For each input image select reference spectrum list from a table.
# The table is read from the file and stored in a simple symbol table.
#
# The table consists of pairs of words. The first word is a list of spectra
# and the second word is the reference spectrum list to be used for each
# spectrum in the first list. Note that the first list is not an input
# list. As a convenience if a reference list is missing the preceding list
# is implied. Some examples follow.
#
# spec1 spec2,spec3,spec4
# spec5
# spec6,spec7 spect8,spec9
# spec10 spec11
# spec12 spec13
# spec14 spec15
procedure reftable (list, table, select)
pointer list # List of input spectra
char table[ARB] # Reference table
int select # Selection method
int i, fd, input, refs
pointer stp, sym
pointer sp, image, ref1, ref2
pointer stopen(), strefsbuf(), stenter(), stpstr(), stfind(), imtopen()
int imtgetim(), open(), fscan(), nscan()
errchk open
begin
# Read the table. Return an error if the file can't be opened.
# Read each table entry of spectrum list and reference list.
# Expand the input list to make a symbol table keyed on the
# spectrum with the reference list string as it's value.
# As a convenience if a reference list is missing the preceding
# list is implied.
fd = open (table, READ_ONLY, TEXT_FILE)
call smark (sp)
call salloc (image, SZ_FNAME, TY_CHAR)
call salloc (ref1, SZ_FNAME, TY_CHAR)
call salloc (ref2, SZ_FNAME, TY_CHAR)
stp = stopen ("table", 10, 10, 20*SZ_FNAME)
while (fscan (fd) != EOF) {
call gargwrd (Memc[image], SZ_FNAME)
call gargwrd (Memc[ref1], SZ_FNAME)
if (nscan() < 1)
next
if (nscan() < 2)
call strcpy (Memc[ref2], Memc[ref1], SZ_FNAME)
else
call strcpy (Memc[ref1], Memc[ref2], SZ_FNAME)
i = stpstr (stp, Memc[ref1], SZ_FNAME)
input = imtopen (Memc[image])
while (imtgetim (input, Memc[image], SZ_FNAME) != EOF) {
call refnoextn (Memc[image])
sym = stenter (stp, Memc[image], 1)
Memi[sym] = i
}
call imtclose (input)
}
call close (fd)
# For each input spectrum find the appropriate reference spectrum list.
# If no list is found print a message and continue. Switch on the
# selection method.
while (imtgetim (list, Memc[image], SZ_FNAME) != EOF) {
call refnoextn (Memc[image])
sym = stfind (stp, Memc[image])
if (sym == NULL) {
call refmsgs (NO_REFSPEC, Memc[image], "", "", "", 0, 0, "")
next
}
input = imtopen (Memc[image])
refs = imtopen (Memc[strefsbuf (stp, Memi[sym])])
switch (select) {
case MATCH:
call refmatch(input, refs)
case NEAREST:
call refnearest (input, refs)
case PRECEDING:
call refprecede (input, refs)
case FOLLOWING:
call reffollow (input, refs)
case INTERP:
call refinterp (input, refs)
case AVERAGE:
call refaverage (input, refs)
}
call imtclose (input)
call imtclose (refs)
}
call stclose (stp)
call sfree (sp)
end
|
