AGO Teleskop

Najprej aktiviram iraf paketa za obdelavo ccd slik:

ecl> imred
imred> ccdred

nato postavim vse parametre na privzete vrednosti. Toda pozor, če obdelave večkrat ponavljamo in imamo vse parametre že dobro nastavljene, ne izvedemo taska unlearn in tudi nastavitve lahko preskočimo. 

ccdred> unlearn ccdred
ccdred> unlearn ccdproc

Potem nastavim parametre za task ccdproc, če so moje slike posnete v binning 2 jmoram imeti v isti mapi tudi kalibracijske datoteke Zero-2.fts, Dark-2.fts, Flat-2B.fts, Flat-2V.fts, Flat-2R.fts, Flat-2W.fts:

ccdred> setinstrument ap16e directory="" site="" review=no
ccdred> ccdproc.fixpix=no
ccdred> ccdproc.overscan=no
ccdred> ccdproc.trim=no
ccdred> ccdproc.zero="Zero-2"
ccdred> ccdproc.dark="Dark-2"
ccdred> ccdproc.flat="Flat-2*"

Iraf si vedno zapomni nastavljene parametre, ozirma zadnje vnešene pri uspešnem izvajanju taska, če pa izvedemo unlearn se vse nastavitve izbrišejo in se postavijo na privzete vrednosti. Nastavitve hrani v datotekah v mapi ~/uparm

Zdal lahko preverim nastavitve:

ccdred> lpar ccdproc                                      
       images = ""              List of CCD images to correct
      (output = "")             List of output CCD images
     (ccdtype = "object")       CCD image type to correct
   (max_cache = 0)              Maximum image caching memory (in Mbytes)
      (noproc = no)             List processing steps only?\n
      (fixpix = no)             Fix bad CCD lines and columns?
    (overscan = no)             Apply overscan strip correction?
        (trim = no)             Trim the image?
     (zerocor = yes)            Apply zero level correction?
     (darkcor = yes)            Apply dark count correction?
     (flatcor = yes)            Apply flat field correction?
    (illumcor = no)             Apply illumination correction?
   (fringecor = no)             Apply fringe correction?
     (readcor = no)             Convert zero level image to readout correction?
     (scancor = no)             Convert flat field image to scan correction?\n
    (readaxis = "line")         Read out axis (column|line)
     (fixfile = "")             File describing the bad lines and columns
     (biassec = "")             Overscan strip image section
     (trimsec = "")             Trim data section
        (zero = "Zero-2")       Zero level calibration image
        (dark = "Dark-2")       Dark count calibration image
        (flat = "Flat-2*")      Flat field images
       (illum = "")             Illumination correction images
      (fringe = "")             Fringe correction images
  (minreplace = 1.)             Minimum flat field value
    (scantype = "shortscan")    Scan type (shortscan|longscan)
       (nscan = 1)              Number of short scan lines\n
 (interactive = no)             Fit overscan interactively?
    (function = "legendre")     Fitting function
       (order = 1)              Number of polynomial terms or spline pieces
      (sample = "*")            Sample points to fit
    (naverage = 1)              Number of sample points to combine
    (niterate = 1)              Number of rejection iterations
  (low_reject = 3.)             Low sigma rejection factor
 (high_reject = 3.)             High sigma rejection factor
        (grow = 0.)             Rejection growing radius
        (mode = "ql")           

Končno lahko naredim obdelavo za vse slike v tekoči mapi, ampak najprej samo preizkus, da vidimo, če vse deluje:

ccdred> ccdproc *.fts noproc=yes

in če je vse v redu, recimo najde vse kalibracijske datoteke, potem obdelamo vse slike hkrati:

ccdred> ccdproc *.fts

V glavo obdelane slike vpiše ključne besede:

ZEROCOR = 'Dec 11 21:19 Zero level correction image is Zero-2'
DARKCOR = 'Dec 11 21:19 Dark count correction image is Dark-2'
FLATCOR = 'Dec 11 21:19 Flat field image is Flat-2V.fts'
CCDPROC = 'Dec 11 21:19 CCD processing done'

Zapisi v glavi nam za vsako sliko dokujmentirajo kdaj in kaj je bilo narejeno pri obdelavi in tudi preprečujejo, da bi naprimer po pomoti še enkrat naredili kalibracije. Na računalniku dome.ago.fmf.uni-lj.si na katerem tudi zajemamo slike je že pripravljena skripta v pythonu, ki nam naredi vse navedeno. Skripto uporabimo za obdelavo ene ali več slik:

$ flatdark.py mojaslika.fts
$ flatdark.py m*.fts