source: trunk/external-alma/components/SpectralComponents/test/tSpectralFit.out@ 3023

Last change on this file since 3023 was 2980, checked in by Malte Marquarding, 10 years ago

Add a copy of casacore/components/SpectralComponents to external-alma directory to prepare for its removal from casacore-trunk. DOn;t activate in SConscript yet.

File size: 10.6 KB
RevLine 
[2980]1Test RecordInterface
2---------------------------------------------------
3toRecord
4fromRecord
5
6Test SpectralEstimate
7---------------------------------------------------
8Made spectrum with components:
9GAUSSIAN element:
10 Amplitude: 5, 0
11 Center: 45, 0
12 Sigma: 1.1, 0
13
14GAUSSIAN element:
15 Amplitude: 4, 0
16 Center: 50, 0
17 Sigma: 1.1, 0
18
19GAUSSIAN element:
20 Amplitude: 6, 0
21 Center: 55, 0
22 Sigma: 1.1, 0
23
24*** n 100
25*** y [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2.09099e-40, 3.34739e-35, 2.34499e-30, 7.18876e-26, 9.64375e-22, 5.66131e-18, 1.45434e-14, 1.63492e-11, 8.04274e-09, 1.73137e-06, 0.000163101, 0.00672359, 0.12129, 0.957476, 3.30757, 5.00013, 3.31295, 1.05451, 0.887271, 2.65278, 4.00036, 2.65413, 0.911529, 1.246, 3.97447, 6.00013, 3.96909, 1.14897, 0.145548, 0.00806831, 0.000195721, 2.07765e-06, 9.65129e-09, 1.9619e-11, 1.74521e-14, 6.79357e-18, 1.15725e-21, 8.62651e-26, 2.81398e-30, 4.01687e-35, 2.50918e-40, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
26Found (using rms, cutoff, minsig, q): (0.5, 1, 0.5, 2) 3 estimates
27Estimate 0: GAUSSIAN element:
28 Amplitude: 6.00013, 0
29 Center: 55.0715, 0
30 Sigma: 1.33398, 0
31
32Estimate 1: GAUSSIAN element:
33 Amplitude: 5.00013, 0
34 Center: 44.9152, 0
35 Sigma: 1.32309, 0
36
37Estimate 2: GAUSSIAN element:
38 Amplitude: 4.00036, 0
39 Center: 50.0043, 0
40 Sigma: 1.26335, 0
41
42*** res [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2.09099e-40, 3.34739e-35, 2.34499e-30, 7.18876e-26, 9.64375e-22, 5.66131e-18, 1.45434e-14, 1.63492e-11, 8.04274e-09, 1.73137e-06, 0.000163101, 0.00672359, 0.12129, 0.957476, 3.30757, 5.00013, 3.31295, 1.05451, 0.887271, 2.65278, 4.00036, 2.65413, 0.911529, 1.246, 3.97447, 6.00013, 3.96909, 1.14897, 0.145548, 0.00806831, 0.000195721, 2.07765e-06, 9.65129e-09, 1.9619e-11, 1.74521e-14, 6.79357e-18, 1.15725e-21, 8.62651e-26, 2.81398e-30, 4.01687e-35, 2.50918e-40, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
43Minimum, Maximum residuals: -0.961305, 0.0086807
44Window (using rms, cutoff, minsig, q): (0.5, 1, 0.5, 2) 0 estimates
45Found (using rms, cutoff, minsig, q): (0.5, 1, 0.5, 1) 3 estimates
46Estimate 0: GAUSSIAN element:
47 Amplitude: 6.00013, 0
48 Center: 54.9829, 0
49 Sigma: 1.17129, 0
50
51Estimate 1: GAUSSIAN element:
52 Amplitude: 5.00013, 0
53 Center: 45.0202, 0
54 Sigma: 1.17355, 0
55
56Estimate 2: GAUSSIAN element:
57 Amplitude: 4.00036, 0
58 Center: 50.0182, 0
59 Sigma: 1.38534, 0
60
61Found (using rms, cutoff, minsig, q): (0, 0, 0, 1) 3 estimates
62Estimate 0: GAUSSIAN element:
63 Amplitude: 6.00013, 0
64 Center: 54.9829, 0
65 Sigma: 1.17129, 0
66
67Estimate 1: GAUSSIAN element:
68 Amplitude: 5.00013, 0
69 Center: 45.0202, 0
70 Sigma: 1.17355, 0
71
72Estimate 2: GAUSSIAN element:
73 Amplitude: 4.00036, 0
74 Center: 50.0182, 0
75 Sigma: 1.38534, 0
76
77Found (using rms, cutoff, minsig, q): (0, 0, 0, 2) 3 estimates
78Estimate 0: GAUSSIAN element:
79 Amplitude: 6.00013, 0
80 Center: 55.0715, 0
81 Sigma: 1.33398, 0
82
83Estimate 1: GAUSSIAN element:
84 Amplitude: 5.00013, 0
85 Center: 44.9152, 0
86 Sigma: 1.32309, 0
87
88Estimate 2: GAUSSIAN element:
89 Amplitude: 4.00036, 0
90 Center: 50.0043, 0
91 Sigma: 1.26335, 0
92
93Test SpectralFit
94---------------------------------------------------
95Spectral elements:
96GAUSSIAN element:
97 Amplitude: 1, 0
98 Center: 1402.5, 0
99 Sigma: 0.0195312, 0
100
101GAUSSIAN element:
102 Amplitude: 2.5, 0
103 Center: 1405, 0
104 Sigma: 0.0390625, 0
105
106GAUSSIAN element:
107 Amplitude: 0.5, 0
108 Center: 1407.5, 0
109 Sigma: 0.078125, 0
110
111GAUSSIAN element:
112 Amplitude: 5, 0
113 Center: 1403, 0
114 Sigma: 0.0292969, 0
115
116---------------------------------------------------
117Data for frequencies 1402.93 - 1403.08
1181402.93: 0.219133
1191402.94: -0.480645
1201402.95: -0.281389
1211402.96: 0.0423511
1221402.97: 0.0837382
1231402.98: 1.05514
1241402.99: 3.14053
1251403: 4.93017
1261403.01: 4.53475
1271403.02: 1.54477
1281403.03: 0.247699
1291403.04: -0.50962
1301403.05: 0.0446798
1311403.06: -0.0518962
1321403.07: -0.750935
1331403.08: 0.499586
134---------------------------------------------------
135Specify fitter on the 4 gaussians:
136GAUSSIAN element:
137 Amplitude: 1, 0
138 Center: 1402.5, 0
139 Sigma: 0.0195312, 0
140
141GAUSSIAN element:
142 Amplitude: 2.5, 0
143 Center: 1405, 0
144 Sigma: 0.0390625, 0
145
146GAUSSIAN element:
147 Amplitude: 0.5, 0
148 Center: 1407.5, 0
149 Sigma: 0.078125, 0
150
151GAUSSIAN element:
152 Amplitude: 5, 0
153 Center: 1403, 0
154 Sigma: 0.0292969, 0
155
156---------------------------------------------------
157Execute the fitter, can do (1:Y): 1
158The results:
159GAUSSIAN element:
160 Amplitude: 0.917797, 0.310674
161 Center: 1402.5, 0.00347673
162 Sigma: 0.00882989, 0.00343021
163Parameters: [0.917797, 1402.5, 0.00882989]
164Errors: [0.310674, 0.00347673, 0.00343021]
165
166GAUSSIAN element:
167 Amplitude: 2.76117, 0.238196
168 Center: 1405, 0.00151172
169 Sigma: 0.0151761, 0.00151172
170Parameters: [2.76117, 1405, 0.0151761]
171Errors: [0.238196, 0.00151172, 0.00151172]
172
173GAUSSIAN element:
174 Amplitude: 0.612907, 0.18968
175 Center: 1407.48, 0.0085527
176 Sigma: 0.0239303, 0.0085527
177Parameters: [0.612907, 1407.48, 0.0239303]
178Errors: [0.18968, 0.0085527, 0.0085527]
179
180GAUSSIAN element:
181 Amplitude: 5.24752, 0.272235
182 Center: 1403, 0.000695977
183 Sigma: 0.0116179, 0.000695945
184Parameters: [5.24752, 1403, 0.0116179]
185Errors: [0.272235, 0.000695977, 0.000695945]
186
187>>>
188Number of iterations: 29
189<<<
190---------------------------------------------------
191GAUSSIAN element:
192 Amplitude: 1, 0
193 Center: 1402.5, 0
194 Sigma: 0.0195312, 0
195
196GAUSSIAN element:
197 Amplitude: 2.5, 0
198 Center: 1405, 0
199 Sigma: 0.0390625, 0
200
201GAUSSIAN element:
202 Amplitude: 0.5, 0
203 Center: 1407.5, 0
204 Sigma: 0.078125, 0
205
206GAUSSIAN element:
207 Amplitude: 5, 0 (fixed)
208 Center: 1403, 0
209 Sigma: 0.0292969, 0
210
211---------------------------------------------------
212Execute the fitter again with a fixed ampl, can do (1:Y): 1
213The results:
214GAUSSIAN element:
215 Amplitude: 0.917772, 0.310662
216 Center: 1402.5, 0.00347618
217 Sigma: 0.00883038, 0.00342961
218Parameters: [0.917772, 1402.5, 0.00883038]
219Errors: [0.310662, 0.00347618, 0.00342961]
220
221GAUSSIAN element:
222 Amplitude: 2.76117, 0.23818
223 Center: 1405, 0.00151162
224 Sigma: 0.0151761, 0.00151162
225Parameters: [2.76117, 1405, 0.0151761]
226Errors: [0.23818, 0.00151162, 0.00151162]
227
228GAUSSIAN element:
229 Amplitude: 0.61293, 0.18968
230 Center: 1407.48, 0.00855088
231 Sigma: 0.0239286, 0.00855088
232Parameters: [0.61293, 1407.48, 0.0239286]
233Errors: [0.18968, 0.00855088, 0.00855088]
234
235GAUSSIAN element:
236 Amplitude: 5, 0 (fixed)
237 Center: 1403, 0.000741387
238 Sigma: 0.0119707, 0.00060533
239Parameters: [5, 1403, 0.0119707]
240Errors: [0, 0.000741387, 0.00060533]
241
242>>>
243Number of iterations: 26
244<<<
245---------------------------------------------------
246Different start values:
247GAUSSIAN element:
248 Amplitude: 0.9, 0
249 Center: 1402.5, 0
250 Sigma: 0.0185547, 0
251
252GAUSSIAN element:
253 Amplitude: 2.25, 0
254 Center: 1405, 0
255 Sigma: 0.0371094, 0
256
257GAUSSIAN element:
258 Amplitude: 0.45, 0
259 Center: 1407.5, 0
260 Sigma: 0.0742188, 0
261
262GAUSSIAN element:
263 Amplitude: 4.5, 0
264 Center: 1403, 0
265 Sigma: 0.027832, 0
266
267Execute the fitter, can do (1:Y): 1
268The results:
269GAUSSIAN element:
270 Amplitude: 0.917785, 0.310679
271 Center: 1402.5, 0.00347653
272 Sigma: 0.00883013, 0.00342999
273
274GAUSSIAN element:
275 Amplitude: 2.76117, 0.238196
276 Center: 1405, 0.00151172
277 Sigma: 0.0151761, 0.00151172
278
279GAUSSIAN element:
280 Amplitude: 0.61293, 0.189693
281 Center: 1407.48, 0.00855143
282 Sigma: 0.0239285, 0.00855143
283
284GAUSSIAN element:
285 Amplitude: 5.24752, 0.272235
286 Center: 1403, 0.000695977
287 Sigma: 0.0116179, 0.000695945
288
289>>>
290Number of iterations: 32
291<<<
292---------------------------------------------------
293Differences:
294Min difference: -1.09994, max: 1.01271, average: 0.00978336, sigma: 0.0100325
295---------------------------------------------------
296---------------------------------------------------
297Estimates:
298Found 20 estimates
299Estimate 0: GAUSSIAN element:
300 Amplitude: 4.93017, 0
301 Center: 307.195, 0
302 Sigma: 3.05174, 0
303
304Estimate 1: GAUSSIAN element:
305 Amplitude: 2.79525, 0
306 Center: 511.819, 0
307 Sigma: 3.29148, 0
308
309Estimate 2: GAUSSIAN element:
310 Amplitude: 0.906691, 0
311 Center: 766.653, 0
312 Sigma: 2.31423, 0
313
314Estimate 3: GAUSSIAN element:
315 Amplitude: 0.90563, 0
316 Center: 501.372, 0
317 Sigma: 3.60214, 0
318
319Estimate 4: GAUSSIAN element:
320 Amplitude: 0.861827, 0
321 Center: 839.949, 0
322 Sigma: 3.43708, 0
323
324Estimate 5: GAUSSIAN element:
325 Amplitude: 0.846625, 0
326 Center: 112.878, 0
327 Sigma: 2.28534, 0
328
329Estimate 6: GAUSSIAN element:
330 Amplitude: 0.751603, 0
331 Center: 871.157, 0
332 Sigma: 2.78939, 0
333
334Estimate 7: GAUSSIAN element:
335 Amplitude: 0.699624, 0
336 Center: 36.4464, 0
337 Sigma: 3.65743, 0
338
339Estimate 8: GAUSSIAN element:
340 Amplitude: 0.682235, 0
341 Center: 395.917, 0
342 Sigma: 3.5893, 0
343
344Estimate 9: GAUSSIAN element:
345 Amplitude: 0.663379, 0
346 Center: 450.12, 0
347 Sigma: 0.548263, 0
348
349Estimate 10: GAUSSIAN element:
350 Amplitude: 0.640783, 0
351 Center: 827.254, 0
352 Sigma: 3.35359, 0
353
354Estimate 11: GAUSSIAN element:
355 Amplitude: 0.616448, 0
356 Center: 245.555, 0
357 Sigma: 3.28022, 0
358
359Estimate 12: GAUSSIAN element:
360 Amplitude: 0.578568, 0
361 Center: 120.043, 0
362 Sigma: 1.45528, 0
363
364Estimate 13: GAUSSIAN element:
365 Amplitude: 0.569486, 0
366 Center: 663.971, 0
367 Sigma: 1.81295, 0
368
369Estimate 14: GAUSSIAN element:
370 Amplitude: 0.56377, 0
371 Center: 611.848, 0
372 Sigma: 2.99822, 0
373
374Estimate 15: GAUSSIAN element:
375 Amplitude: 0.561451, 0
376 Center: 937.91, 0
377 Sigma: 2.32834, 0
378
379Estimate 16: GAUSSIAN element:
380 Amplitude: 0.54666, 0
381 Center: 336.559, 0
382 Sigma: 1.8761, 0
383
384Estimate 17: GAUSSIAN element:
385 Amplitude: 0.524891, 0
386 Center: 443.644, 0
387 Sigma: 0.809077, 0
388
389Estimate 18: GAUSSIAN element:
390 Amplitude: 0.523377, 0
391 Center: 645.266, 0
392 Sigma: 3.01063, 0
393
394Estimate 19: GAUSSIAN element:
395 Amplitude: 0.515985, 0
396 Center: 446.61, 0
397 Sigma: 1.18317, 0
398
399---------------------------------------------------
400Limited number of estimates:
401Found 4 estimates
402Estimate 0: GAUSSIAN element:
403 Amplitude: 4.93017, 0
404 Center: 307.195, 0
405 Sigma: 3.05174, 0
406
407Estimate 1: GAUSSIAN element:
408 Amplitude: 2.79525, 0
409 Center: 511.819, 0
410 Sigma: 3.29148, 0
411
412Estimate 2: GAUSSIAN element:
413 Amplitude: 0.906691, 0
414 Center: 766.653, 0
415 Sigma: 2.31423, 0
416
417Estimate 3: GAUSSIAN element:
418 Amplitude: 0.90563, 0
419 Center: 501.372, 0
420 Sigma: 3.60214, 0
421
422---------------------------------------------------
Note: See TracBrowser for help on using the repository browser.