AGO Teleskop

Arhivske posnetke lahko iščemo v raznih podatkovnih zbirkah. Tu bom prikazal kako lahko uporabimo program Aladin. V glavnem meniju izberem File->Open, nato izberem zavihek Aladin Images ter navedem še ime zvezde (tarče) naprimer Cyg 61, imenski strežnik mi ime takoj prevede v nebesne koordinate, za velikost iskanja pa je smiseno recimo 20'.

Izberem sliko, ki jo želim, nato še format FITS in kliknem SUBMIT. Izbiram slike ki imajo vsaj približno enako ločljivost na piksel in seveda različne datume opazovanja. Izbral sem sliko iz POSSI  (posneto okoli leta 1951), POSSII dve sliki, posneto 1989 in 1993 in še sliko iz 2MASS iz leta 1998. Nastavim format prikaza štiri slike, centriram prvo sliko in kliknem match zdaj se mi vse slike poravnajo na enake wcs (nebesne) koordinate in enako velikost:

Še preverim datume, v nekaterih podatkovnih bazah, če so slike sestavljene, težko dobiš pravi datum opazovanja. Najbolje je zato odpreti originalni FITS header:

in iščem ključne besede, kot DATE-OBS, MJD-OBS ali EPOCH, ki povejo kadaj je bila slika posneta. Izberem MJD-OBS:

Zdaj se lotim meritve koordinat obeh zvezd. Najlaže to naredim kar v Aladinu, najprej si nastavim za vsako sliko ustrezno skaliranje, tako da sta centra zvezde kar se da razločno vidni:

in slike tudi močno povečam; grem s kurzorjem-križcem na center zvezde, kar se da natančno, tu se zanesem na vidno presojo:

Desni klik na miški in izberem "Copy recticle location in clipboard" in potem koordinate prilepim v tekstovno datoteko, spodnja zvezda je komponenta B, zgornja A, meritve za vsako zvezdo v svojo datoteko:

# ime     opazovanje MJD     REKTAS.   DEKLINAC.
Cyg_61_B  50984.3575189815   316.72822 +38.74082
Cyg_61_B  48446.38625        316.71682 +38.73520
Cyg_61_B  47775.1979166667   316.71508 +38.73349
Cyg_61_B  33835.4345833333   316.65885 +38.69999

Cyg_g1_A  50984.3575189815   316.72235 +38.74805
Cyg_61_A  48446.38625        316.71158 +38.74182
Cyg_61_A  47775.1979166667   316.70920 +38.74038
Cyg_61_B  33835.4345833333   316.65294 +38.70574

za kontrolo narišem v gnuplotu, najprej rektascenzijo:

set grid 0
f(x)=a*x+b
fit f(x) 'cyg61B.txt' using 2:3  via a,b
plot 'cyg61B.txt' u 2:3 w t "Cyg 61 B" ls 1 , f(x) ls 2

 Final set of parameters            Asymptotic Standard Error
=======================            ==========================
a               = 4.02285e-06      +/- 4.727e-08    (1.175%)
b               = 316.523          +/- 0.002163     (0.0006832%)

in za deklinacijo:

set xlabel "Cas opazovanja MJD /dan"
set ylabel "Deklinacija /stopinje"
g(x)=c*x+d
fit g(x) 'cyg61B.txt' using 2:4  via c,d
plot [*:*][38.6:38.8] 'cyg61B.txt' u 2:4 t "Cyg 61 B" ls 1 , g(x) ls 2

Final set of parameters            Asymptotic Standard Error
=======================            ==========================
c               = 2.39275e-06      +/- 1.793e-08    (0.7492%)
d               = 38.6191          +/- 0.0008202    (0.002124%)

in graf:

Končna vrednost lastno gibanje:

 v rektascenziji  = a * 365.25 * 3600 * cos(38.72) * 1000 = 4125.3 mas/yr (1 +- 0.012)

v deklinaciji       = c * 365.25 * 3600 * 1000 =  3146.2 mas/yr (1 +- 0.075)

V Simbadu: