# ALINEA IMAGENES DE FORMA RAPIDA
# By ANGEL R. 
# angelrls@ll.iac.es
# Adaptado de macro de ZVEZDA (Jose Caballero) para NIR en 1.5m CST
# zvezda@ll.iac.es
#
# Cargar paquete "tables" y subaquete "ttools"
#
# tables
# ttools
#
# fichero 'out' debe tener las imagenes de entrada
# en formato k1, k2, k3 ... 
#
# ejecutar con: cl < alinea.cl


int primera = 1     # primera imagen k1,k2,k3...
int posiciones = 2  # n'umero de imagenes totales

s2 = 'putita'      # cadena con el nombre del objeto y filtro

del ("p_ali")
for (i = primera; i < primera + posiciones; i=i+1)  {

print ("p"//i, >> "p_ali")

}



# presiona "a" en la * de referencia y despu'es pulsa "n"
# =)
del ("coor"); imexam ("@out", keeplog=yes, logfile="coor")


# Creaci'on del fichero de desplazamientos
# (c'odigo complejo)

del ("xey"); fields (files="coor",fields="1-2", > "xey")

head ("xey", n=1) | scan (x, y)  # first

del ("coo_ref"); print (x,y, >> "coo_ref")

del ("sx"); print (x, "-C1", >> "sx")
del ("sy"); print (y, "-C2", >> "sy")

tcalc (table="xey", outcol="C1", equals="@sx")
tcalc (table="xey", outcol="C2", equals="@sy")

del ("shifts")

head ("xey", n=1) | scan (x,y)  # second


 print (x, y, >> "shifts")


for (k = 2; k <= posiciones; k=k+1)  {


  tail ("xey", n=-k+1) | head n=1 | scan (x,y)
  print (x, y, >> "shifts")


}


# alineamiento

imalign (input="@out", reference="k"//primera,
    coords="coo_ref", output="@p_ali", shifts="shifts",
    boxsize=9, bigbox=11, interp_type="linear", niterat=6, trimimages=yes)


# combinaci'on final

imdel ("z"//s2)
imcombine ("@p_ali", "z"//s2, combine="median",
    reject="minmax", nlow=1, nhigh=1, zero="mode")


# imagen de salida

display ("z"//s2,1)


# borra las im'agenes de pasos intermedios

#imdel ("@p_ali")

flpr
flpr
flpr