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
|
include "../identify.h"
# AID_SHIFT -- Find a new shifted dispersion solution assuming (nearly) the
# same dispersion per pixel and the same dispersion direction. The shift is
# assumed to be less than or equal to the dispersion range of the input
# dispersion. The input is an ID pointer have the previous dispersion
# solution and features but with the new spectrum. If there are more than 10
# features then the list of user feature coordinates is used as the reference
# list. If there are not enough features or the previous search fails then
# the the coordinate list is used as the reference. The returned result is a
# new ID pointer if the algorithm succeeds or the original ID pointer if it
# fails along with an error status.
procedure aid_shift (id, crsearch, cdsearch)
pointer id #I ID object
double crsearch #I Search range
double cdsearch #I Search range
pointer aid
bool found, aid_autoid()
double crpix, crval, cdelt, id_fitpt()
begin
# Set approximate dispersion from input dispersion solution.
crpix = ID_NPTS(id) / 2 + 1
crval = id_fitpt (id, crpix)
cdelt = (FITDATA(id,ID_NPTS(id)) - FITDATA(id,1)) /
(ID_NPTS(id) - 1)
# Initialize AUTOID.
call aid_init (aid, "aidpars")
call aid_setd (aid, "crval", crval)
call aid_setd (aid, "cdelt", cdelt)
call aid_setd (aid, "crpix", crpix)
call aid_sets (aid, "cddir", "sign")
call aid_setd (aid, "crsearch", crsearch)
call aid_setd (aid, "cdsearch", cdsearch)
call aid_seti (aid, "nbest", 5)
found = false
if (ID_NFEATURES(id) > 10) {
# Try shift using features.
call aid_seti (aid, "ntarget", ID_NFEATURES(id))
call aid_seti (aid, "nreference", ID_NFEATURES(id))
call aid_setr (aid, "wrms", 0.5)
call aid_setr (aid, "wfmatch", 0.5)
call aid_setr (aid, "wftmatch", 0.)
call aid_sets (aid, "refspec", "FEATURES")
found = aid_autoid (id, aid)
}
if (!found) {
# Try shift using coordinate list.
call aid_seti (aid, "ntarget", max (ID_NFEATURES(id),20))
call aid_seti (aid, "nreference", max (ID_NFEATURES(id),40))
call aid_setr (aid, "wrms", 0.5)
call aid_setr (aid, "wfmatch", 0.25)
call aid_setr (aid, "wftmatch", 0.25)
call aid_sets (aid, "refspec", "COORDLIST")
found = aid_autoid (id, aid)
}
call aid_free (aid)
if (!found)
call error (1, "No solution not found")
end
|
