Fourier Spektroskopi Data Reduktion
Indlæsning og nulpunktskorrektion
Vi går ud fra, at I har Jeres billeder liggende som
FITS billeder, d.v.s filer med 'extension' .fit (eller .FIT).
De skal nu indlæses i IMSYS format med kommandoen
hvor image er filnavnet minus 'extension' og indholdet af
den kantede parentes kun skrives, hvis 'extension' er FIT.
En manual til IMSYS findes som en postscript fil,
som I kan printe ud med kommandoen
-
lpr -Phpnew ~srf/public_html/lab/imsys.ps
Når alle billeder er konverteret kan man med kommandoen
list checke, om de er med. list -l image giver
informationer om det enkelte billede.
Alle billederne skal nu nulpunkts-korrigeres med kommandoen
-
bias image image_b 2 10 12 xmax
hvor image_b er det nye billednavn og '2 10' er de søjler der
bruges i 'overscan' striben til beregning af nulpunktet og '12 xmax'
er de søjler der udtrækkes af billedet. Der er ingen
grund til at gemme 'overscan' arealet.
Fladtfeltet
Alle fladtfelt optagelser skal nu midles. Det kan man enten gøre
med et IDL program af egen konstruktion eller man kan bruge IMSYS.
IMSYS har to 'regne'-kommandoer til billeder: calc og
const. Fladtfelt billederne adderes med calc kommandoen,
middelværdien findes og middelbilldet normeres efter denne
opskrift
-
icp f1_b flat_lab3
-
calc flat_lab3 "+" f2_b
-
mean flat_lab3 200 250 200 250 (giver middelærdien = f)
-
const flat_lab3 "/" f/10000 (hvor f/10000 skal regnes ud og indsættes)
Her er f1_b og f2_b to nulpunktskorrigerede fladtfelts optagelser.
Man har nu et standard fladtfelt billede flat_lab3.
Næste trin er nu at kalibrere alle optagelserne af Mg lampens
interferens striber. Det gøres simpelt med
-
calc image_b "/" flat_lab3 10000
image_b skal nu bare bringes tilbage i FITS format med kommandoen
og dette billede er nu klar til at behandles (i IDL f.eks.)
Et simpelt program er kopieret til hvert data katalog med
navnet fringes.pro. Det kan bruges som udgangspunkt.