BETA LYRAE

Jedna od retkih oblasti amaterske astronomije koja u jednoj priličnoj meri može dati podatke dovoljno kvalitetne za izvlačenje naučnih zaključaka je fotometrija.
Po definiciji astrometrija se odnosi na merenje položaja nebeskih tela a fotometrija na merenje sjaja, i obe metode se u praksi lepo dopunjuju. Osim toga, zadnjih nekoliko decenija su unele pravu malu revoluciju pre svega zbog pojave CCD kamera, pa je to itekako uticalo na činjenicu da mnogi ljudi iz svog dvorišta mogu da snime tranzite vansolarnih planeta. Pre toga se još od 1861.g. za merenje sjaja koristio fotometar, neka vrsta detektora-brojača fotona; a još pre toga se za merenje sjaja koristila metoda vizuelne komparacije.

Današnje CCD kamere se mogu  posmatrati kao skup nekoliko miliona fotometara od kojih svaki za sebe meri intenzitet svetla jednog segmenta slike. To onda u praksi znači veliku prednost, jer sa starim fotometrom vi ste merili samo jednu zvezdu a sa kamerom možete recimo imati u kadru neko zvezdano jato i samim tim merenje sjaja stotina zvezda u isto vreme postaje moguće.

Metoda koja je nama najinteresantnija, a ujedno spada i u najprostije je diferencijalna fotometrija. Radi se o postupku merenja koji podrazumeva da kroz isto vidno polje optičkog sistema snimamo tj merimo ciljani objekt a i kontrolni objekt koji služi za upoređivanje sjaja. Ovim se izbacuje gomila promenljivih iz jednačine a i softver koji se pritom koristi veoma lako izlazi na kraj sa merenjem razlike. Tako se dobije krivulja razlike sjaja između promenljive i kontrolne zvezde, bez navođenja konkretne magnitude oba ova objekta. Ako znamo tačnu magnitudu kontrolne zvezde onda je relativno lako odrediti magnitudu promeljive zvezde, ali rekoh, ovde to uopšte nije toliko primarno koliko sama promena u jedinici vremena. Ako je u pitanju promenljiva, tranzit vansolarne planete, tranzit zvezde pratioca, rotacija asteroida ili bilo šta drugo, nama je po definiciji bitnija sama ta promena i njena magnituda od aposlutne magnitude.

Pošto za uplovljavanje u fotometrijske vode nije potrebna samo dobra volja već još ponešto, zacrtao sam da odradim jednostavnu vežbu - diferencijalnu fotometriju nečeg osnovnog. Tako valjda neću mnogo da zabrljam a možda i dobijem neke rezultate.
Optički sistem: EOS 40D i objektiv EF 50 1.8 zavrnut na 2.8; pravilo je da se objektivi u fotometriji koriste širom otvoreni ali je ovaj objektiv širom otvoren ne baš briljantan primer optičko-tehničkih dostignuća japanske industrije. Ova preporuka da se objektiv koristi širom otvoren je iz razloga što blende raznih objektiva imaju različite oblike, ovaj je recimo petougao, ali i zato da se ne bi bespotrebno rasipali resursi.
Jako bitno je pre početka snimanja podesiti vreme u aparatu po UT precizno u sekund. Ovo iz razloga jer softver uzima vreme iz aparata za merenje.
Snimci su sa tripoda i ekspozicija je bila po deset sekundi. Sve je blago defokusirano i pošto sam takve snimke dobio greškom - sinula mi je ideja da bih mogao da ih iskoristim za fotometriju. Jer, aperturna fotometrija se i radi sa blago defokusiranom slikom.

Evo isečka u originalnoj rezoluciji; Lirin trougao se lepo raspoznaje. Fokusa, dakle, ni u tragovima i to je ono što je odlično.


Zapazite da sem Vege nijedna zvezda nije pregorela. Double-double ni u nijednom kanalu nije, a Zeta Lyrae (donja zvezda u trouglu) je presaturirana samo u G i B kanalu. Ko ne veruje - Photoshop i pipetu u šake, podesite na point sample i gledajte gore desno u info boksu. Sve vrednosti ispod 255 su prihvatljive za merenje. U tome je i poenta - sve snimke za aperturnu fotometriju treba malo defokusirati kako bi ovi svetli kružići mogli da posluže za precizno merenje; u slučaju da imamo potpuno fokusiran snimak - šta možemo da merimo na pregoreloj beloj tačkici od dva-tri piksela? Ovako je razmažemo na deset piksela i njen sjaj se automatski distribuira na većoj površini i nema pregorevanja.

Što se tiče softvera koji sam koristio, svi znaju da ja volim da se mučim sa Irisom. Program vrlo moćan, napisan da može kompletno da radi i iz komandne linije, dozlaboga kompleksan i jednom rečju savršen. Vidi se da ga je napisao čovek veoma upućen istovremeno i u astronomiju i u informatiku - Cristian Buil. Što je najlepše program je potpuno džabe, besplatan, koštaju vas samo vaši živci dok naučite.

Na terenu treba snimiti proizvoljan broj light frejmova, ja sam sa goreopisanim parametrima snimio ukupno 23. Zatim sa tim istim podešavanjima po pitanju ISO i dužine izlaganja, ali poklopcem na objektivu treba snimiti dark frejmove. Koliko - to zavisi od vas; ja sam snimio 16 a preporuka je što više (mada više od tridesetak već gubi poentu, neće manji broj uneti šum toliko koliko ljudi ponekad misle). Na kraju a najbolje je kući, treba snimiti flat frejmove na ISO100. Oni će služiti za rešavanje prašine i vinjetiranja optičkog sistema, to je veoma bitno kod fotometrije - mnogo bitnije nego kod klasične astrofotografije gde u obzir mogu da dođu i sintetički fletovi (jer se po rubovima ionako ne nalazi ništa bitno). U fotometriji na samom rubu mogu biti zvezde koje se mere a takav slučaj ćemo upravo imati - bez ispravne kalibracije sve pada u vodu.
Flet se najbolje uzima na ravnomernom plavom nebu, ili ekranu monitora. Ja sam imao oblačno nebo i morao sam da u PS-u generišem ravnomerno sivu, objektiv maksimalno defokusiram i približim se monitoru. Možda nebitno pošto je u RAW-u, ali podesite i WB na daylight, nikako na auto. Neki kažu da nema nikakve veze a neki kažu da ima, ne košta ništa da izbegnemo dilemu pošto će auto opcija na Canonu (a možda i Nikonu) garantovano da promaši. Fletova treba da ima isto kao i darkova, ja sam snimio 16.

Čega sve ima možete videti na jednom fletu gde je kontrast pritegnut do maksimuma. Neko će imati prašinu, neko pečurke, a neko (kao ja) čist senzor i vinjetiranje.


Ovo deluje jako lepo, a to je iz razloga što je pedesetica objektiv za fulfrejm, odnosno moj aparat vidi samo optički najbolji, centralni deo slike koju objektiv pravi. Na FF aparatu zasigurno ne bismo bili tako sigurni u performanse ovog objektiva, bar kad je vinjetiranje u pitanju.

Za sledeći postupak najbolji opis je dat u veoma korisnom Vodiču kroz digitalnu astrofotografiju koji je napisala grupa hrvatskih autora, odakle sam ukratko prepričao (uz minimalne izmene radi uprošćavanja) jedan deo koji je pisao Vid Nikolić, a koji vodi poreklo sa stranice Saratogaskies. Ukratko, radi se iz komandne linije da bi bilo brže.

Odabere se opcija Digital Photo > Decode RAW files i zatim se mišem prevuku light frejmovi. Da im se neko ime, ja recimo koristim LFA (Light Frames A) a zatim se klikne na CFA. Zatim se sve to obriše opcijom Erase list i prebace se dark frejmovi; da se neko ime da asocira na to, recimo DFA (Dark Frames A) i opet CFA. Na kraju tako isto prebacite flat frejmove i imenujete (FLA). Kad je sve to gotovo onda se ide na opciju Done.
Prvo što treba uraditi je da dobijemo master dark, tj kombinaciju svih dark fajlova. Ako bismo koristili samo jedan dark za oduzimanje onda bi njegov šum izazvao neujednačenosti na finalnom steku. Dakle, komanda je smedian DFA 16, gde 16 označava broj dark frejmova u mom konkretnom slučaju. Nakon toga se snimi - save dark.

Sa finalnim dark frejmom u memoriji otkucate komandu stat koja daje očitavanja parametara koji su vezani za hot-piksele. Zatim se kuca find_hot cosmetic N gde je N broj koji varira u zavisnosti od tipa DSLR-a, starosti senzora i ambijentalne temperature. Meni se na 40D taj broj koji treba ukucati kreće oko 15000-16000 i slobodno eksperimentišite jer je cilj da Iris prijavi da je "hot pixel number" između 10 i 50. Naravno, leti na +30C se i očekuje da broj hot-piksela bude veći, na 20D je išao i preko 200.

Sad ide opet komandna linija:

>load DFA1
>fill 0
>save offset
>fill 1
>save flat

Nakon ovoga ide Digital photo > Make a flat-field a tu se upisuje:

Generic name: FLA
Offset image: offset
Normalisation value: 20000
Number: 16

Komandna linija: >save flat


Sledeće je iz menija Digital photo > Preprocesing.

Input generic name: LFA
Offset: offset
Dark: dark
Flat-field:  flat
Cosmetic file:  cosmetic
Output generic name: LFB
Number: 23 (broj LFB tj light-frejmova)
Optimize - isključeno.

Sad dolazi opet malo teorije. Svi DSLR senzori (sa izuzetkom nekih egzotičnih Sigma Foveon i Fuji senzora) imaju Bayer-ovu matricu zasnovanu na pikselima tri osnovna tipa (R, G i B). Za osnovnu interpretaciju nečega što je u fotometriji snimljeno moramo između ostalog da znamo i na kojoj je talasnoj dužini snimljeno (odnosno koji je filter korišćen pri snimanju). Za nas je najbitnije da od nekih par stotina filtera koji su se koristili (i koriste se još uvek) jedan praktično perfektno odgovara: Johnson V-band sa pikom na 551nm je najbliži G kanalu na DSLR-u.
Dakle, potrebno je samo izvaditi G kanal i raditi sa njim kao sa crnobelim snimkom. Razlog zašto je to tako je činjenica da mnoge zvezde imaju različite boje (plava zvezda je primetno veća u B nego u R kanalu) tako da se na ovaj način uzima neki etalon, koji se otprilike nalazi na sredini opsega vidljive svetlosti. Zato se V filter ponekad i naziva vizuelni - mada je činjenica da bilo koji kanal možemo koristiti za fotometriju.
Postoje dve mogućnosti za dobijanje zelenog kanala. Jedna je da u Photoshopu napravite akciju kojom bi se iz svakog pojedinačnog RAW-a selektovao samo G kanal i kopirao kao layer u novi fajl, a zatim se rezultat snima kao TIF. Jeste da bi bilo mukotrpno jedan po jedan RAW učitavati i primenjivati akciju, postoji opcija još veće automatizacije ali je ja zasad ne znam a ne troši mi se vreme na traganje za tim, ionako je ostala obrada sasvim dovoljno zahtevna. Tu akciju zasad prepuštam onima koji bi da eksperimentišu, možda neki drugi put bude stajao opis. Jedino ostaje dilema kako u Irisu otvoriti tolike TIF fajlove, ali i tu postoji odgovor - iz menija File se odabere Select files, odabere TIF i ubace svi potrebni fajlovi.
Druga solucija je da to umesto nas uradi Iris. Pošto postoje dva zelena piksela uz po jedan plavi i crveni, samim tim imamo i dva zelena kanala.

Komanda koja "razbija" CFA fajl na kanale - jedan crveni, jedan plavi i dva zelena; je komanda split. Zašto dva... pojma nemam, zato što je takva zapravo arhitektura DSLR senzora, imamo duplo više zelenih nego plavih i crvenih piksela. Dakle treba mukotrpno učitavati jedan po jedan LFB fajl, razbiti na četiri kanala a zatim spojiti dva zelena u jedan i snimiti takav fajl.
Dat je primer za prvi LFB, cifra 1 se dalje menja u 2 itd do kraja liste od 23 ili koliko već imate snimaka. Prečicu zasad nisam otkrio; ovu sekvencu ćete upisati i izvršiti 23 puta.

>load LFB1 * cifra 1 se menja
>split_cfa s1 s2 s3 s4
>load s1
>load s4
>add_norm s 2 * cifra 2 se NE menja
>save LFB1 * cifra 1 se menja

Registracija je neophodna da bismo dalje radili. Registrovani fajlovi imaju posmatranu zvezdu uvek na istom mestu, i to je neophodno za stekiranje (mada mi ne idemo u tom pravcu jer nam i ne treba - treba nam samo sekvenca snimaka za fotometriju):

>setspline 1
>coregister2 LFB LFC  23

I to bi sve bilo kad su u pitanju pripreme za ono glavno što sledi. Da bismo regularno nastavili sa merenjem moramo da izvršimo jednu proveru. Nađite zvezdu na snimku koju merite, nacrtajte levim tasterom miša mali kvadratić oko nje, a zatim desnim tasterom miša odaberite shape (neke verzije imaju opciju profile) iz malog padajućeg menija koji se pritom pojavljuje. Dobićete grafikon na kome je prestavljena distribucija sjaja zvezde idući od centra ka periferiji. Na apscisi je predstavljena udaljenost u desetim delovima piksela, a na ordinati je intenzitet signala zabeležen na senzoru.


Na gornjem grafiku vidite profil zvezde onako kako treba da izgleda, a na donjem je profil Vege gde se jasno vidi da je zvezda debelo presaturirana, i to u poluprečniku od 3 piksela. Mada je moguće i tu raditi fotometriju ipak preciznost će do jedne mere biti žrtvovana a ovde se sve i vrti oko preciznosti. Drugim rečima, ako se ne potrudite da otklonite svaku nepreciznost koju je moguće otkloniti - nemojte se ni baviti fotometrijom. Jer, kome su potrebni neprecizni rezultati? Kakve je zaključke moguće dobiti na osnovu nepreciznih merenja? U estetici je preciznost merenja često od manjeg značaja u odnosu na sposobnost improvizacije od strane astronoma, ali ovde je tačnost u postupku praktično imperativ.

Takođe iz grafika se vidi da su zvezde tako defokusirane da se zona gornje zvezde prostire tri piksela van centra, to je faktički poluprečnik; a da je donja zvezda poluprečnika 4.5 piksela.

Nakon što smo se uverili da ni promenljiva ni referentna zvezda nisu izgorele u centru, sledi fotometrija iz sekvence snimaka koji su obeleženi sa LFC. Iz menija se odabere Analyze > Select object. Sad sledi obeležavanje prvo promenljive zvezde, a zatim referentne (možete odabrati i više referentnih zvezdi, ja sam odabrao jednu zbog jednostavnosti ali je preporuka da ih bude bar dve-tri, zbog neuporedivo veće tačnosti). Nakon što se završi obeležavanje isključi se Select object. Sledeće je iz menija Analysis > Automatic photometry. U prozoru se redom upisuje:

Input generic name: LFC
Number: 23
Output data file: phot
Magnitude output: čekirano
No matching: čekirano
Aperture photometry: čekirano

Interpretacija rezultata može da bude na prvi pogled malo konfuzna (kao da ovo dosad nije bilo dovoljno konfuzno), i rezultati su smešteni pre svega u okviru fajla phot a takođe se vide i u Output prozoru. Leva kolona je UT vreme (jeste podesili UT na početku? Niste? A lepo sam vam rekao), druga kolona je izmerena magnituda prve zvezde, treća kolona je magnituda druge zvezde itd.
Par dodatnih objašnjenja. Prvo, magnituda koja je ovde izmerena nije ujedno i ona magnituda na koju svi mislimo (eh, da je tako prosto...), ovo je faktički magnituda kojom Iris izražava razliku između centra zvezde i onog selektovanog prstena neposredno oko zvezde. Drugim rečima ovo je razlika između zvezde i neba (koje takođe ima neku svoju magnitudu). Sve što nas ovom prilikom interesuje je diferencijalna magnituda odnosno razlika između promenljive i referentne zvezde i ona je vrlo tačno izmerena. Zato referentna zvezda ne sme da bude promenljiva, jer bi automatski unela grešku u merenje.
Druga stvar je da možemo lepo da ispratimo tok merenja pomoću fajla VERIF.DAT. Iz menija se ide na View > Plot > File > Load > verif. Ovde se prikazuje varijacija u magnitudi promenljive zvezde svakog fajla posebno - u poređenju sa prosekom magnitude referentne zvezde. Ovaj fajl služi za orijentacioni pregled profila promenljive zvezde, nikako za naučne komparacije.


Da sad ovde imamo eklipsnu promenljivu - profil bi bio sasvim drugačiji. Ovde vidimo random grešku u suštini.

Sve je ovo bilo neka vrsta uvoda; vreme je da se pređe na merenje. Zvezda koja je označena kao Beta Lyrae spada u promenljive eklipsnog tipa sa periodom od nešto manje od 13 dana. Varijacije su 3.2-4.4mag. Interesantno je što je ova zvezda kontinuirano zaklonjena pratiocem, ovo je binarni sistem polukontaktnog tipa (jedna zvezda prevazilazi Rošeov limit). Manja zvezda je nekad bila veća i evoluirala je u giganta koji je gubio masu u korist druge zvezde. Druga zvezda, danas masivnija, je okružena akrecionim diskom iz perioda usisavanja mase sa svog pratioca. Rezultat postojanja ovog gigantskog diska je činjenica da mi danas ne možemo da vidimo drugu zvezdu, upravo jer je skrivena unutar diska. A sjaj akrecionog diska (koji i dalje vredno radi) je otprilike jedna petina sjaja celokupnog binarnog sistema. Ovo je sve moguće iz jednog prostog razloga: dve zvezde orbitiraju na relativno maloj međusobnoj udaljenosti.

Ono što je izuzetno bitno u ovoj priči je da se na 46 sekundi (skoro jedan najveći Jupiterov prečnik) nalazi zvezda magnitude 7.2, a na 86" se nalazi zvezda magnitude 10. Zašto je to bitno, zato što prilikom defokusiranja na objektivu 50mm a i merenja spoljni prsten Irisa bez problema obuhvata i ove zvezde i to sve utiče na tačnost merenja. Kad bi se merenje obavilo na teleskopu preklapanja ne bi bilo.

Za kontrolnu zvezdu sam odabrao nešto tamniju μ Lyrae, magnitude 5.12. Podatke koji su u Irisu bili u Output prozoru sam ručno upisivao u Excel, tačnije drugu i treću kolonu, a onda sam ih oduzimao i to je bila još jedna kolona - tačno ono što mi treba, razlika između posmatranih zvezda izražena u magnitudama.
Ovo je grafikon iz Excela koji pokazuje Irisove izmerene nagnitude obe zvezde, primetite da to nema veze sa pravim magnitudama tih zvezda. Ovde su vrednosti kao prvo negativne, a kao drugo one predstavljaju razliku između zvezdine okoline i same zvezde.


Primetno je da postoje neke oscilacije u sjaju obe zvezde, npr na frejmu 11 obe zvezde su podjednako pretrpele pad sjaja. Moguće je da se radi o promeni atmosferskih uslova u smislu cirusa ili nečeg sličnog - ali ne mogu da isključim ni da je greška samog aparata iz nekog razloga bila 0.1 magnitude. U svakom slučaju nebitno je, jer mi ovde (upravo iz tih fantomskih razloga) merimo samo razliku u sjaju ove dve zvezde.
Oduzeo sam magnitude ove dve zvezde i ta razlika je na sledećem grafikonu:


Pošto je ovolika razlika malo isuviše, moram da priznam da je Beta Lyrae bila na visini od 16 stepeni nad horizontom u momentu merenja. Neko iskusniji u fotometriji bi rekao da sve ovo jednostavno treba obrisati, meriocu opaliti šamar i meriti u momentu kad je Beta bliža zenitu - tačno, ali ovo je bila samo vežba. Iz tog razloga i oscilacije: turbulencija i scintilacija itekako utiču na merenje, što se i vidi iz izmerene razlike. Oscilacije su išle do 0.13 magnitude što se praktično isključivo može pripisati visini od 16 stepeni.
Još malo diskusije: sudeći po donjoj tabeli, ulovio sam momenat kad je Beta Lyrae u minimumu. Međutim, zajedno sa njom su i pratioci koje sam gore naveo, sasvim je sigurno da je iris u sjaj Betae uračunao i njihov doprinos. To bi onda moglo da pokaže da je Beta apsurdno ispod svog minimuma - a sa druge strane u kalkulaciju nije uvedena apsorpcija koja na ovako malim visinama sasvim sigurno igra veliku ulogu... Tako da ovo merenje onda očigledno ne spada u red najtačnijih.

Na kraju imamo grafikon samo sjaja Betae Lyrae izračunat na osnovu poznate magnitude kontrolne zvezde, što je i bio cilj merenja. Jedina korisna sugestija bi bila da se sve ovo obavi tako što bi se merio sjaj jednom dnevno i tako u periodu od 13 dana. Onda bismo dobili pravu krivulju sjaja - ovde su date oscilacije sjaja u periodu od nepuna 4 minuta.


I na kraju da vidimo kako bi izgledao fokusiran snimak tog istog regiona neba sa istim parametrima: isto 23 snimka sa tripoda, isto f2.8...
Ovde je uzet isečak iz originalne rezolucije, odnosno možete smatrati da je fotka uvećana 2.5 puta.


U gornjem levom uglu se vidi referentna zvezda (μ Lyrae) a u donjem desnom Beta Lyrae. Čak ako bolje pogledate u donjem levom uglu tog kvadrata se vidi tirkizna tačkica - M57 - ali bez bilo kakvih detalja. Pa mnogo bi i bilo ako bih sa 50mm tu očekivao detalje.

Коментари