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
|
include <imhdr.h>
include "ms.h"
# MSGET -- Allocate memory and get data from the MULTISPEC database
# and associated image.
#
# MSGHDR -- Allocate memory and get MULTISPEC header information.
# MSGCOMMENTS -- Allocate memory and get MULTISPEC comments.
# MSGPARAM -- Allocate memory and get a line of MULTISPEC parameter data.
# MSGSAMPLE -- Allocate memory and get SAMPLE line array.
# MSGFIT -- Get parameter fit for a spectrum.
# MSGFITS -- Get parameter fit for all spectra.
# MSGGAUSS5 -- Get a line of GAUSS5 parameter data.
# MSGIMAGE -- Get a line of the image with possible averaging.
# MSGHDR -- Allocate memory and get MULTISPEC header information.
procedure msghdr (ms)
pointer ms # MULTISPEC data structure
int i
int dbread()
begin
if (MS_DATA(ms, HDR) == NULL)
call calloc (MS_DATA(ms, HDR), LEN_MS_HDR, TY_STRUCT)
i = dbread (MS_DB(ms), NAME(ms, HDR), HEADER(ms), 1)
end
# MSGCOMMENTS -- Allocate memory and get MULTISPEC comments.
procedure msgcomments (ms)
pointer ms # MULTISPEC data structure
int i
int dbread()
begin
if (MS_DATA(ms, COMMENTS) == NULL)
call calloc (MS_DATA(ms, COMMENTS), SZ_MS_COMMENTS, TY_CHAR)
i = dbread (MS_DB(ms), NAME(ms, COMMENTS), COMMENT(ms, 1), 1)
end
# MSGPARAM -- Allocate memory and get a line of MULTISPEC parameter data.
procedure msgparam (ms, parameter, line)
pointer ms # MULTISPEC data structure
int parameter # Parameter ID
int line # Sample line to be obtained
int i
char reference[SZ_MS_KEY]
bool is_param_id()
int dbread()
begin
# Check if the the requested parameter is valid.
if (!is_param_id (parameter))
call error (MS_ERROR, "Bad parameter identifier")
if (MS_DATA(ms, parameter) == NULL)
call calloc (MS_DATA(ms, parameter), MS_NSPECTRA(ms), TY_REAL)
# Make reference to the desired database record.
call sprintf (reference, SZ_MS_KEY, "%s[%d]")
call pargstr (NAME(ms, parameter))
call pargi (line)
i = dbread (MS_DB(ms), reference, PARAMETER(ms, parameter, 1), 1)
end
# MSGSAMPLE -- Allocate memory and get SAMPLE line array.
procedure msgsample (ms)
pointer ms # MULTISPEC data structure
int i
int dbread()
begin
if (MS_DATA(ms, SAMPLE) == NULL)
call malloc (MS_DATA(ms, SAMPLE), MS_NSAMPLES(ms), TY_INT)
i = dbread (MS_DB(ms), NAME(ms, SAMPLE), LINE(ms,1), 1)
end
# MSGFIT -- Get parameter fit for a spectrum.
procedure msgfit (ms, parameter, spectrum)
pointer ms # MULTISPEC data structure
int parameter # Parameter ID for desired fit
int spectrum # Spectrum
int i
char reference[SZ_MS_KEY]
pointer sp, fit
bool is_fit_id()
int dbread()
errchk cvrestore
begin
# Check if for valid parameter id.
if (!is_fit_id (parameter))
call error (MS_ERROR, "Bad fit identifier")
# Allocate memory for the curfit pointers.
if (MS_DATA(ms, parameter) == NULL)
call malloc (MS_DATA(ms, parameter), MS_NSPECTRA(ms), TY_INT)
# Allocate memory for the curfit coefficients.
call smark (sp)
call salloc (fit, 7 + MS_NSAMPLES(ms), TY_REAL)
# Reference appropriate data.
call sprintf (reference, SZ_MS_KEY, "%s[%d]")
call pargstr (NAME(ms, parameter))
call pargi (spectrum)
i = dbread (MS_DB(ms), reference, Memr[fit], 1)
iferr (call cvrestore (CV(ms, parameter, spectrum), Memr[fit]))
;
call sfree (sp)
end
# MSGFITS -- Get parameter fits.
procedure msgfits (ms, parameter)
pointer ms # MULTISPEC data structure
int parameter # Parameter ID for desired fit
int i
begin
do i = 1, MS_NSPECTRA(ms)
call msgfit (ms, parameter, i)
end
# MSGGAUSS5 -- Get a line of GAUSS5 parameter data.
procedure msggauss5 (ms, line)
pointer ms # MULTISPEC data structure
int line # Sample line to be obtained
begin
call msgparam (ms, I0, line)
call msgparam (ms, X0, line)
call msgparam (ms, S0, line)
call msgparam (ms, S1, line)
call msgparam (ms, S2, line)
end
# MSGIMAGE -- Get a line of the image with possible averaging.
procedure msgimage (im, line, naverage, data)
pointer im # Image descriptor
int line # Line to be gotten from the image
int naverage # Number of line to use in average
real data[ARB] # The output data array
int i, line_start, line_end
real nlines
pointer buf
pointer imgl2r()
begin
# If naverage is <= 1 copy the image line to the data array
# Else average the several lines.
if (naverage <= 1) {
call amovr (Memr[imgl2r (im, line)], data, IM_LEN(im,1))
} else {
# Determine starting and ending lines for the average.
line_start = max (1, line - naverage / 2)
line_end = min (IM_LEN(im, 2), line_start + naverage - 1)
# Clear data array for accumulating sum and then vector
# add the image lines.
call aclrr (data, IM_LEN(im, 1))
do i = line_start, line_end {
buf = imgl2r (im, i)
call aaddr (Memr[buf], data, data, IM_LEN(im, 1))
}
# Vector divide by the number of lines to form average.
nlines = line_end - line_start + 1
call adivkr (data, nlines, data, IM_LEN(im, 1))
}
end
|
