[CCD-CR] Iris - skladani snimku a prace s FITS

Lukas Kral krall at troja....cvut.cz
Wed Mar 27 21:46:35 CET 2002


Ahoj,
asi nejsem ten nejvice znaly veci, ale aspon na neco zkusim reagovat...

>jaky je dynamicky rozsah vysledneho snimku?
16 bitu, stejne jako u vstupnich FITSu. Uklada se to zase ve formatu FITS.
Moznou saturaci vzniklou diky poskladani resi specialni "skladaci" funkce
COMPOSIT tak, ze pokud v prubehu automatickeho skladani dojde k saturaci,
vysledny slozeny snimek je vynasoben nejakym koeficientem (<1) urcenym tak,
aby ani nejjasnejsi pixely nepresahly 32767.
Algoritmy pouzivane pro posun a skladani snimku v manualu bohuzel detailne
popsane nejsou (jsou popsany nekde jinde??). Vim jen, ze posun snimku vuci
referencnimu je provaden se subpixelovou presnosti, poloha referencni hvezdy
na jednotl. snimcich je automaticky urcovana metodou PSF.

>jak je osetren vyskyt spatnych pixelu?
Je osetren, pokud pro slozeni pouzijete zminenou funkci COMPOSIT - ze
skladani jsou vyrazeny pixely, ktere se od prumeru lisi o vice nez zvoleny
nasobek standardni deviace. Tento proces muze byt navic iterativni s
nastavitelnym poctem iteraci.

Jak jsem vcelku experimentalne zjistil, Iris zachazi s FITSy trochu
zajimavym zpusobem. Jak asi vite, standardem pro ukladani binarnich dat v
16-bitovem FITS souboru je signed integer (16-bitove cislo, kde prvni bit je
znamenko), tedy rozsah -32768...+32767. Takto uklada FITSy i Iris. Vetsina
programu, ktera uklada FITSy, ovsem pouziva take nepovinny parametr BZERO s
hodnotou +32768. Skutecna hodnota pixelu se pak behem nacitani snimku ze
souboru vypocte jako
pix = bin_int * BSCALE + BZERO,
kde bin_int je ono binarni cislo zapsane ve FITSu (-32768...32767), BSCALE
je cislo normalne rovne jednicce (dalsi nepovinny parametr) a BZERO je tedy
konstanta, kterou je nutno pricist, abychom dostali skutecne hodnoty pixelu
(nezaporne).
Priklad: mame-li na snimku pixel s nulovou hodnotou (pix=0), pak je ve FITSu
normalne zapsan jako cislo -32768 (protoze -32768 * 1 + 32768 = 0).

Proc to vsecko rozpitvavam? Protoze prave v tomhle se Iris od jinych
programu lisi! Iris totiz parametr BZERO pri ukladani do FITS nepouziva a
tak je tento povazovan (podle definice FITS) za rovny nule. V Irisu jsou
tedy dovoleny hodnoty pixelu -32768...32767 a takto jsou take ukladany.
Nenechte se ale zmylit, tento zpusob ukladani vyhovuje specifikaci formatu
FITS, jen je u astronomickych snimku dost podivny - proc pracovat se
zapornymi pixely?

Iris tedy nema zadny problem pracovat se zapornymi hodnotami pixelu, zato
vam vsak nedovoli presahnout hodnotu +32767 (u jinych podobnych programu
jsme vetsinou zvykli na praci v kladnem rozsahu 0...65535). Kdyz do nej
zkusite nacist snimek FITS z jineho programu, ktery obsahuje pixely
"jasnejsi" nez +32767, Iris jim tuto hodnotu proste natvrdo priradi! Tzn.
vsemu jasnejsimu nez +32767 nastavi prave hodnotu +32767. Ostatni pixely
jsou nacteny spravne a je jim spravne prictena hodnota BZERO, pokud je
parametr pritomen (Iris tedy tento parametr zna a pouziva ho pri nacitani,
ale ne pri ukladani).

Tohle ujasneni cele situace me stalo nemalo casu a tak jsem vam to tu chtel
popsat, at nemusite experimentovat s bity tak jako ja :-( Zajimalo by me
ovsem, proc to autor programu Ch. Buil takhle divne vyresil. Taky nevim,
proc narozdil od ostatnich programu jsou snimky v Irisu zobrazovany vzhuru
nohama (pixel 0,0 je vlevo dole a ne vlevo nahore). Ale krome techhle
podivnosti se mi ten program opravdu moc libi a ma spoustu skvelych funkci!

Jasne nebe preje
Lukas






More information about the CCD mailing list