# Macro de IRAF para reducir imagenes
# DEBEN estar los ficheros de texto "entrada" y "out"
# CAMBIA el nombre de la IMAGEN FINAL e imagen de FLATS para cada caso
#

# ejecutar con: cl < flats.cl

real m1
s1="biasr"    # BIAS
s2="flatr"   # FLAT

imarith (operand1="@entrada", op="-", operand2=s1, result="@out")

imcombine (input="@out", output=s2) 

imstat(s2)

print ("Introduce valor medio de pixel en flat")
scan(m1)

imarith (operand1=s2, op="/", operand2=m1, result="n"//s2)

imstat ("n"//s2)
display ("n"//s2,1)

imdel ("@out")
imdel (s2)

flpr
flpr