Initial commit

1 files changed, 67 insertions, 0 deletions

diff --git a/noao/onedspec/identify/autoid/aidshift.x b/noao/onedspec/identify/autoid/aidshift.x
new file mode 100644
index 00000000..1b910338
--- /dev/null
+++ b/noao/onedspec/identify/autoid/aidshift.x
@@ -0,0 +1,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