procedure distcor(iname)

# By R. Elston
# Task to produce fit coords image to be used to 
# rectify images for astigmatic distortion. Use
# the database with the task transform to straighten
# night sky lines

string iname{prompt = 'Emission line  image name'}

begin

identify (iname,
section="m c", database="database", coordlist="irlines", nsum="10",
match=10., maxfeatures=50, zwidth=100., ftype="emission", fwidth=4.,
cradius=5., threshold=10., minsep=2., function="leg", order=2, sample="*",
niterate=0, low_reject=3., high_reject=3., grow=0., autowrite=no,
graphics="stdgraph", cursor="")

reidentify (iname,
iname, "yes", interactive="no", section="m c", newaps=yes, override=yes,
refit=no, trace=yes, step="5", nsum="10", shift="0.", nlost=0, cradius=5.,
threshold=10., addfeatures=no, coordlist="irlines", match=10.,
maxfeatures=50, minsep=2., database="database", logfiles="logfile",
plotfile="", verbose=yes, graphics="stdgraph", cursor="")

fitcoords (iname,
fitname="", interactive=yes, combine=no, database="database",
deletions="deletions.db", function="chebyshev", xorder=3, yorder=2,
logfiles="STDOUT,logfile", plotfile="plotfile", graphics="stdgraph",
cursor="")

imdel("test")

transform (iname,
"test", iname , database="database", interptype="poly5",
x1=INDEF, x2=INDEF, dx=INDEF, nx=INDEF, xlog=no, y1=INDEF, y2=INDEF, dy=INDEF,
ny=INDEF, ylog=no, flux=yes, logfiles="STDOUT,logfile")

display ("test", 1)



end