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
|
include <mach.h>
include "refspectra.h"
# REFPRECEDE -- Assign preceding reference spectrum based on sort key.
# If there is no preceding spectrum assign the nearest following spectrum.
procedure refprecede (input, refs)
pointer input # List of input spectra
pointer refs # List of reference spectra
bool ignoreaps # Ignore aperture numbers?
int i, i1, i2, nrefs, ap
double sortval, d, d1, d2
pointer sp, image, gval, refimages, refaps, refvals, refgvals
bool clgetb(), streq(), refginput(), refgref()
int imtgetim(), imtlen()
begin
call smark (sp)
call salloc (image, SZ_FNAME, TY_CHAR)
# Task parameters
ignoreaps = clgetb ("ignoreaps")
# Tabulate reference spectra. This expands the reference list,
# checks the spectrum is a reference spectrum of the appropriate
# aperture.
call salloc (refimages, imtlen (refs), TY_INT)
call salloc (refaps, imtlen (refs), TY_INT)
call salloc (refvals, imtlen (refs), TY_DOUBLE)
call salloc (refgvals, imtlen (refs), TY_INT)
nrefs = 0
while (imtgetim (refs, Memc[image], SZ_FNAME) != EOF) {
call refnoextn (Memc[image])
if (!refgref (Memc[image], ap, sortval, gval))
next
for (i=0; i<nrefs; i=i+1)
if (streq (Memc[image], Memc[Memi[refimages+i]]))
break
if (i == nrefs) {
call salloc (Memi[refimages+nrefs], SZ_FNAME, TY_CHAR)
call salloc (Memi[refgvals+nrefs], SZ_FNAME, TY_CHAR)
call strcpy (Memc[image], Memc[Memi[refimages+i]], SZ_FNAME)
Memi[refaps+i] = ap
Memd[refvals+i] = sortval
call strcpy (Memc[gval], Memc[Memi[refgvals+i]], SZ_FNAME)
nrefs = i + 1
}
}
if (nrefs < 1)
call error (0, "No reference images specified")
# Assign preceding reference spectra to each input spectrum.
# Skip input spectra which are not of the appropriate aperture
# or have been assigned previously (unless overriding).
while (imtgetim (input, Memc[image], SZ_FNAME) != EOF) {
call refnoextn (Memc[image])
if (!refginput (Memc[image], ap, sortval, gval))
next
i1 = 0
i2 = 0
d1 = MAX_REAL
d2 = -MAX_REAL
do i = 1, nrefs {
if (!streq (Memc[gval], Memc[Memi[refgvals+i-1]]))
next
if (!ignoreaps && ap != Memi[refaps+i-1])
next
d = sortval - Memd[refvals+i-1]
if ((d >= 0.) && (d < d1)) {
i1 = i
d1 = d
}
if ((d < 0.) && (d < d2)) {
i2 = i
d2 = d
}
}
if (i1 > 0) # Nearest preceding spectrum
call refspectra (Memc[image], Memc[Memi[refimages+i1-1]], 1.,
Memc[Memi[refimages+i1-1]], 0.)
else if (i2 > 0) # Nearest following spectrum
call refspectra (Memc[image], Memc[Memi[refimages+i2-1]], 1.,
Memc[Memi[refimages+i2-1]], 0.)
else { # No reference spectrum found
call refprint (STDERR, NO_REFSPEC, Memc[image], "", "", "",
ap, 0, "")
do i = 1, nrefs {
if (!streq (Memc[gval], Memc[Memi[refgvals+i-1]])) {
call refprint (STDERR, REF_GROUP, Memc[image],
Memc[Memi[refimages+i-1]], Memc[gval],
Memc[Memi[refgvals+i-1]], ap, Memi[refaps+i-1], "")
next
}
if (!ignoreaps && ap != Memi[refaps+i-1])
call refprint (STDERR, REF_AP, Memc[image],
Memc[Memi[refimages+i-1]], Memc[gval],
Memc[Memi[refgvals+i-1]], ap, Memi[refaps+i-1], "")
next
}
}
}
call sfree (sp)
end
|
