AGO Teleskop

Če imajo slike narejen match, to je če imajo v glavo vpisane astrometrične WCS koordinate, jih lahko zelo enaostavno seštevamo, zlagamo ali sestavljamo v mozaik. Seštevanje slik gre takole, recimo, da imam 5 slik, ki se imenujejo Polje1_B_1.fts, Polje1_B_2.fts, ..., Polje1_B_5.fts in jih želim povprečiti v skupno sliko, ki se bo imenovala Polje1_B.fts. Slike so med lahko minimalno zamaknjene glede na referenčno sliko. Najprej si izberem referenčno sliko, običajno prvo ali srednjo v sekvenci slikanje. Moja referenčna bo kar Polje1_B_1.fts. Slike najprej registriram, to pomeni poravnam glede na referenčno sliko:

    # slike najprej poravnam po wcs
    images
    immatch
    unlearn wregister
    files Polje1_B_*.fts > imlist
    files Polje1_B_*//-wcs.fts > outlist
    wregister.function = "polynomial"
    wregister.xxorder = 4
    wregister.xyorder = 4
    wregister.xxterms = "half"
    wregister.yxorder = 4
    wregister.yyorder = 4
    wregister.yxterms = "half"
    wregister.reject = INDEF
    wregister.calctype = "real"
    wregister.geometry = "geometric"
    wregister.xsample = 1.
    wregister.ysample = 1.
    wregister.interpolant = "poly5"
    wregister.boundary = "nearest"
    wregister.constant = 0.0
    wregister.fluxconserve = yes
    wregister @imlist Polje1_B_1.fts @outlist

Dobim serijo poravnanih slik, ki se imenujejo Polje1_B_1-wcs.fts, Polje1_B_2-wcs.fts,...Imena so v tekstovni datoteki outlist, to je zdaj vhodna datoteka. Izberem tip pikslov izhodne slike real, za operacijo hočem povprečenje in utežim glede na čas osvetlitve:

    unlearn imcombine
    imcombine.outtype = "real"
    imcombine.lthreshold = 2
    imcombine.combine = "average"
    imcombine.expname = "EXPOSURE"
    imcombine.weight = "exposure"
    imcombine @outlist Polje1_B.fts 

in dobim:

 List of images to combine (@outlist):
 List of output images (Polje1_B.fts):
 Dec 29 19:58: IMCOMBINE
 combine = average, scale = none, zero = none, weight = none
 blank = 0.
       Images
        Polje1_B_1.fts
        Polje1_B_2.fts
        Polje1_B_3.fts
        Polje1_B_4.fts
        Polje1_B_5.fts
Output image = Polje1_B.fts, ncombine = 5
immatch>

Za metodo lahko izberemo različne operacije: seštevanje, povprečenje, mediano. Paziti moramo le, da ne prekoračimo obseg števil v katerem računamo. Privzeto je outtype=real, kar je v redu, so pa slike dvakrat večje kot pri tipu outtype=ushort (unsigned short). Kakorkoli že zložena slika podeduje WCS koordinate referenčne slike.

 ccdred> imstat *Polje1*B*fts fields=image,mean,stddev nclip=5
#               IMAGE      MEAN    STDDEV
       Polje1_B_1.fts     507.2     34.43
       Polje1_B_2.fts     499.3     32.64
       Polje1_B_3.fts      500.     32.44
       Polje1_B.fts       501.2     20.14

Povprečna vrednost neba je ostala enaka, standardna deviacija oziroma šum seštete slike pa se je zmanjšal. Če bi upoštavali samo fotonski šum in bi zanemarili šuma branja ccdja se šum zmanjša za \( \sqrt {n} \)  in je n število seštetih slik.