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
|
include "refspectra.h"
# REFMATCH -- Assign reference spectrum by match against reference list.
procedure refmatch (input, refs)
pointer input # List of input spectra
pointer refs # List of reference spectra
int ap
double sortval
pointer sp, image, refimage, gval
bool refgref(), refginput()
int imtgetim(), imtlen()
begin
if (imtlen (input) != imtlen (refs))
call error (0, "Input and reference list have different lengths")
call smark (sp)
call salloc (image, SZ_FNAME, TY_CHAR)
call salloc (refimage, SZ_FNAME, TY_CHAR)
# Assign reference spectra to each input spectrum.
# Skip spectra which are not of the appropriate aperture
# or have been assigned previously (unless overriding).
while ((imtgetim (input, Memc[image], SZ_FNAME) != EOF) &&
(imtgetim (refs, Memc[refimage], SZ_FNAME) != EOF)) {
call refnoextn (Memc[image])
call refnoextn (Memc[refimage])
if (!refginput (Memc[image], ap, sortval, gval))
next
if (!refgref (Memc[refimage], ap, sortval, gval))
next
call refspectra (Memc[image], Memc[refimage], 1.,
Memc[refimage], 0.)
}
call sfree (sp)
end
|
