Začnem z datotekami v mapi graf, irafu naredim: ecl> cd /home/bojan/data/2014/V2014-01-07/graf ecl> ll total 32872 drwxrwxr-x 2 bojan bojan 4096 Mar 26 13:20 . drwxrwxr-x 14 bojan bojan 20480 Mar 26 13:19 .. -rw-rw-r-- 1 bojan bojan 8406720 Jan 20 2014 sLO_M67B.fts -rw-rw-r-- 1 bojan bojan 8406720 Jan 20 2014 sLO_M67V.fts -rw-rw-r-- 1 bojan bojan 8406720 Jan 20 2014 sSH_M67B.fts -rw-rw-r-- 1 bojan bojan 8406720 Jan 20 2014 sSH_M67V.fts Najprej določim parametre slike sigma in fwhm: ecl> imstat *.fts # IMAGE NPIX MEAN STDDEV MODE sLO_M67B.fts 3630142 352.5 10.29 352.1 sLO_M67V.fts 3592145 517.9 11.57 514.1 sSH_M67B.fts 3623131 71.71 4.859 71.96 sSH_M67V.fts 3552577 97.95 5.146 98.29 ecl> display sLO_M67V.fts ecl> imexamine z miško ozačim nakaj lepih zvezd in stisnem , Log file imex.log open # COL LINE RMAG FLUX SKY N RMOM ELLIP PA PEAK GFWHM 535.45 753.29 13.45 41714.0 364.00 315 5.86 0.066 -42.2 1420.19 4.80 1237.41 723.82 11.84 184376.0 351.00 316 5.62 0.092 -52.6 6807.05 4.62 498.68 1536.05 11.77 195639.0 359.00 312 6.31 0.291 -73.0 5291.24 5.37 1979.14 1228.85 12.00 158957.0 337.00 312 5.84 0.195 79.8 4416.85 5.63 1443.69 1111.00 12.43 106797.0 348.00 310 5.58 0.161 -69.7 3621.31 4.93 826.84 962.25 11.85 181676.0 363.00 315 5.74 0.096 -48.8 6275.39 4.78 676.20 719.79 12.79 76740.0 364.00 314 5.78 0.066 -36.0 2645.07 4.78 Fwhm je potem okoli 4.7 Grem v paket apphot in si nastavim update in verbose na yes ecl> noao.digiphot.apphot apphot> epar apphot PACKAGE = digiphot TASK = apphot (version= May00) (wcsin = logical) The input coordinates wcs (wcsout = logical) The output coordinates wcs (cache = no) Cache image in memory ? (verify = yes) Verify critical parameters ? (update = yes) Update critial parameters ? (verbose= yes) Print verbose output ? (graphic= stdgraph) Default graphics device (display= stdimage) Default display device (mode = ql) ($nargs = 0) Zdaj nastavim sigma in fwhm v datapars apphot> datapars (scale = 1.) Image scale in units per pixel (fwhmpsf = 4.7) FWHM of the PSF in scale units (emission = yes) Features are positive ? (sigma = 11.) Standard deviation of background in counts (datamin = 0.) Minimum good data value (datamax = 57000.) Maximum good data value (noise = "poisson") Noise model (ccdread = "RDNOISE") CCD readout noise image header keyword (gain = "GAIN") CCD gain image header keyword (readnoise = 0.) CCD readout noise in electrons (epadu = 1.) Gain in electrons per count (exposure = "EXPOSURE") Exposure time image header keyword (airmass = "AIRMASS") Airmass image header keyword (filter = "FILTER") Filter image header keyword (obstime = "JD") Time of observation image header keyword (itime = 1.) Exposure time (xairmass = INDEF) Airmass (ifilter = "FILTER") Filter (otime = "INDEF") Time of observation (mode = "ql") apphot> unlearn daofind apphot> daofind sLO_M67B.fts apphot> unlearn phot apphot> epar phot PACKAGE = apphot TASK = phot image = The input image(s) skyfile = The input sky file(s) (coords = default) The input coordinate files(s) (default: image.coo. (output = default) The output photometry file(s) (default: image.mag. (plotfil= ) The output plots metacode file (datapar= ) Data dependent parameters (centerp= ) Centering parameters (fitskyp= ) Sky fitting parameters (photpar= ) Photometry parameters (interac= no) Interactive mode ? (radplot= no) Plot the radial profiles in interactive mode ? (icomman= ) Image cursor: [x y wcs] key [cmd] (gcomman= ) Graphics cursor: [x y wcs] key [cmd] (wcsin = )_.wcsin) The input coordinate system (logical,tv,physical,w (wcsout = )_.wcsout) The output coordinate system (logical,tv,physical) (cache = )_.cache) Cache the input image pixels in memory ? (verify = )_.verify) Verify critical parameters in non-interactive mode (update = )_.update) Update critical parameters in non-interactive mode (verbose= )_.verbose) Print messages in non-interactive mode ? (graphic= )_.graphics) Graphics device apphot> phot sLO_M67B.fts Za radij zaslonke vzamem malo manj kot je fwhm torej 4.6 Zdaj nastavim še parametre sigma in fwhm za sliko V in ponovim postopek: apphot> daofind sLO_M67V.fts apphot> phot sLO_M67V.fts Dobljene .mag datoteke najprej sortiram po magnitudi, pride prav kasneje: apphot> psort sLO_M67B.fts.mag.1 field=mag apphot> psort sLO_M67V.fts.mag.1 field=mag in naredim izvleček v tekstovno datoteko: apphot> pdump sLO_M67V.fts.mag.1 fields="XC,YC,MAG,MER,IFILTER" expr="MAG!=INDEF & MERR!=INDEF" > LO_M67V.txt apphot> pdump sLO_M67B.fts.mag.1 fields="XC,YC,MAG,MER,IFILTER" expr="MAG!=INDEF & MERR!=INDEF" > LO_M67B.txt Grem ven iz irafa in si instaliram program match iz http://spiff.rit.edu/match/ Verzijo 10 že prevedeno za ubuntu 14.4 si lahko naložite tudi iz naše strani: http://dome.ago.fmf.uni-lj.si/podatki/razno/ Zdaj lahko uparim tekstovne datoteke: $ match LO_M67B.txt 0 1 2 LO_M67V.txt 0 1 2 TRANS: a=-33.283209127 b=0.999599801 c=-0.000380664 d=46.051714449 e=-0.000353179 f=0.999370929 sig=5.8306e-02 Nr=505 Nm=529 sx=1.7727e-01 sy=2.0889e-01 kar mi izpiše so koeficienti transformacijske matrika, velja x' = A + B*x + C*y y' = D + E*x + F*y Program mi tudi generira datoteke matched.mtA matched.mtB kjer so zvezde iz orve in druge slike ki se ujemajo in datoteki matched.unA in matched.unB kjer so zvezde ki se ne ujemajo (te me ne zanimajo) Odprem datoteki v tekstovnem urejevalniku (ki ima možnost kopirati po stolpcih): $ kate matched.mtA matched.mtB in označim "blok selection mode" da izbere stolpce, zadnji stolpec magntude iz matched.mtA prenesem v datoteko matched.mtB, zbrišem še prvi stolpec, ki je zaporedna številka, ker jo ne rabim in spravim v datoteko LO_M67VB.txt Datoteka (začetek) zdaj tako: # x y V B 1300.7380000 1464.6680000 14.14 15.15 840.2750000 880.6450000 14.21 15.37 1260.5790000 1245.0460000 14.29 15.26 902.0340000 1426.3010000 14.35 16.08 1569.8330000 59.9290000 14.39 15.77 988.5450000 1680.3650000 14.46 15.60 1486.8310000 1485.0380000 14.48 16.20 590.6520000 884.6670000 14.50 15.55 1118.2950000 770.3840000 14.51 15.40 .... Grem nazaj v iraf in naredim še pretvorbo koordinat x,y v nebesne (world) koordinate rektascenzija, deklinacija apphot> imcoords imcoords> epar wcsctran PACKAGE = imcoords TASK = wcsctran input = LO_M67VB.txt The input coordinate files output = LO_M67VB.wcs The output coordinate files image = sLO_M67V.fts The input images inwcs = logical The input coordinate system outwcs = world The output coordinate system (columns= 1 2 3 4 5 6 7) List of input file columns (units = h n) List of input coordinate units (formats= ) List of output coordinate formats (min_sig= 3) Minimum precision of output coordinates (verbose= yes) Write comments to the output file ? (mode = ql) Začetek datoteke je sedaj, koordinati sta rektascenzija in deklinacija v decimalnih stopinjah # Image: sLO_M67V.fts Wcsin: logical Wcsout: world # rek dek V B 132.89311342 11.852570599 14.14 15.15 132.8100906 11.74987854 14.21 15.37 132.88588056 11.813980279 14.29 15.26 132.8211742 11.845753505 14.35 16.08 132.94170592 11.6058267 14.39 15.77 132.8367605 11.890405852 14.46 15.6 Vse se preverim, xy koordinate iz LO_M67VB.txt narišem na sliko: imcoords> display sLO_M67V.fts imcoords> tvmark 1 LO_M67VB.txt mark=circle color=204 radii=7 label=no Datoteko LO_M67VB.wcs, vzamem ven recimo zvezdo 132.8100906 11.74987854 in grem poiskat v simbad, če najdem. Če Simbad nima točno tam zvezde je nekaj narobe. Še narišem za kontrolo v gnupoltu: gnuplot> plot "LO_M67VB.wcs" u ($4-$3):($3) graf je v redu je pa precej premaknjen in narobe obrnjen. Kalibnracija magnitude glede na standardne magnitude B,V Uporabili bomo Vezier in iskali naše zvezde po koordinatah, za referenčni katalog izberemo: Montgomery K.A., Marschall L.A., Janes K.A. (1993) Astron. J. 106, 181 CCD photometry of the old open cluster M67 1993AJ....106..181M To je katalog 54 iz http://www.univie.ac.at/webda/ Dostop do magnitud pa je pri webdi je malo bolj neroden ker moramo zvezde iskat v dveh tabelah. Zato raje uporabimo http://vizier.u-strasbg.fr/viz-bin/VizieR-3 za katalog vpišemo "Montgomery K.A., Marschall L.A., Janes K.A." ali pa kar "J/AJ/106/181" in izberemo table3. Zvezde lahko iščemo posamič, ali pa izberemo zavihek "List Of Tagret" tu pa označimo "single table" in "add your input as first column" ter pripravimo datoteko, iz LO_M67VB.wcs vzamemo po prvih 30 vrstic samo koordinate. Iz zadtkov in sestavimo tabelo v takšni obliki: # ----------------izmerjeno------------ ---katalog---- -----izracun---- # rek dek v b V B-V v-V b-B 132.89311342 11.852570599 14.14 15.15 10.940 0.220 132.8100906 11.74987854 14.21 15.37 11.078 0.430 132.88588056 11.813980279 14.29 15.26 11.063 0.190 132.8211742 11.845753505 14.35 16.08 11.200 1.080 132.8367605 11.890405852 14.46 15.60 11.251 0.415 132.92669177 11.856178819 14.48 16.20 11.330 1.070 .... Iz tabele izračunamo, recimo v preglednici, aditivni konstanti za V in B V B 'POVPREČJE 3,179279 3,899375 'STD.DEV. 0,038882 0,029192 Če kakšna zvezda močno izstopa jo izbrišemo iz tabele. Lako tudi preverimo ker so nekatere zvezde znane kot spremenjive, teh seveda ne vzamemo za kalibracijo. Zgaj v gnuplotu lahko narišem HR diagram: gnuplot> plot [0:2][18:10] "LO_M67VB.wcs" u ($4-3.89-$3+3.179):($3-3.179)