Back

Ujemanje zvezd v IRAFu

Viewing page version #2
(Restore this version) 

Modified: 12 January 2016, 11:45 AM   User: Bojan Dintinjana  → Bojan

V irafu je procedura malo bolj komplicirana. Naredimo match takole:
Najprej si za našo sliko pripravimo vhodni zvezdni katalog. Za to lahko uporabimo kakšnega od večjih katalogov na VizieR.
Jaz uporabljam USNO_A2 katalog, ki ga imam na disku in program Wcstools za iskanje. Ukaz za generiranje je naslednji, recimo za sliko M27.fts, v ds9 odčitam koordinate centra slike 19:59:36 +22:43:11 in pripravim tabelo zvezd: 

 $ scat -c ua2 -n 170 -mx 16  -r 700x700 19:59:36 +22:43:11  J2000 > starlist.wcs
koordinate so približne koordinate centra slike. Omejitev nastavim za zvezde do 16 mag oziroma največ 170 zvezd. Iskalni radij je 700" (velikost slike na Vegi je 21'=1260"). Tabela se spravi v datoteko starlist.wcs, ki izgleda takole:
$ scat -c ua2 -n 170 -mx 16  -r 700x700 19:59:36 +22:43:11  J2000 
1125.14851851 20:00:06.101 +22:42:17.93   0.0   9.0     0  419.87
1125.14819594 19:59:20.774 +22:36:51.95   0.0  10.0     0  433.70
1125.14819606 19:59:20.793 +22:36:53.05  10.4  10.7   395  432.61
1125.14834854 19:59:42.848 +22:52:44.47  10.7   9.9   395  581.24
1125.14807327 19:59:04.105 +22:37:42.92  11.6   8.4   395  550.01
1125.14851717 20:00:05.909 +22:42:19.63  11.8  10.9   395  417.03
1125.14801286 19:58:56.336 +22:44:39.28  11.8  10.7   395  555.80
1125.14797044 19:58:50.839 +22:44:55.12  12.1  10.9   395  633.40
....
Prvi stolpec je katakoško ime zvezde, drugi in tretji stolpec sta rektascenzija in deklonacija zapisani v obliku ure:min:s. tretji in četrti stolpec pa sta magnituda B in R (pozor, te magnutude so fotografske in niso fotometrično dovolj natančne za kalibracijo fotometrije). 

Nato v irafu pripravimo datoteko, seznam zvezd na naši sliki: 

cl> imstat  M27.fts fields='stddev' nclip=5
cl> datapars.sigma=12
cl> daofind M27.fts
cl> psort  M27.fts.coo.1 field='MAG'
cl> !head -n 210  M27.fts.coo.1  > M27.fts.coo.2
Na koncu smo še odrezali datoteko, tako da vsebuje samo 170 svetlejših zvezd (40 vrstic je v datoteki glava/komentar na začetku). Zdaj naredimo match našega seznama in kataloga. Uporabili bomo metodo trikotnikov in tangentno projekcijo višjega reda, za urejanje parametrov uporabimo task epar ali pa tipkamo (kopiraj/prilepi) ukaze: 
cl> imcoords
cl> ccxymatch.input="M27.fts.coo.2"
cl> ccxymatch.reference = "starlist.wcs"
cl> ccxymatch.output = "_image_.mat.1"
cl> ccxymatch.tolerance = 7.
cl> ccxymatch.ptolerance = 10.
cl> ccxymatch.refpoints = ""
cl> ccxymatch.xin = 1025.
cl> ccxymatch.yin = 1025.
cl> ccxymatch.xmag = 0.636
cl> ccxymatch.ymag = 0.636
cl> ccxymatch.projection = "tnx"
cl> ccxymatch.lngref = 5.4759943
cl> ccxymatch.latref = 35.8256396
cl> ccxymatch.lngcolumn = 2
cl> ccxymatch.latcolumn = 3
cl> ccxymatch.xcolumn = 1
cl> ccxymatch.ycolumn = 2
cl> ccxymatch.lngunits = "hours"
cl> ccxymatch.latunits = "degrees"
cl> ccxymatch.separation = 20
cl> ccxymatch.pseparation = 20
cl> ccxymatch.matching = "triangles"
cl> ccxymatch.nmatch = 60
cl> ccxymatch.ratio = 10.
cl> ccxymatch.nreject = 10
cl> ccxymatch.verbose = "yes"
cl> ccxymatch M27.fts.coo.2  starlist.wcs _image_.mat.1 7. 10.
Za langref in latref vnesemo rektascenzijo v urah in deklinacijo centra slike. Če je vse v redu dobimo ~ 100 prepoznanih zvezd in naredimo še astrometrični fit, red transformacije nastavimo na 3 ali 4 : 
cl> ccmap.input = "_image_.mat.1"
cl> ccmap.database = "database.db"
cl> ccmap.solutions = "M27.fts"
cl> ccmap.images = "M27.fts"
cl> ccmap.results = ""
cl> ccmap.xcolumn = 3
cl> ccmap.ycolumn = 4
cl> ccmap.lngcolumn = 1
cl> ccmap.latcolumn = 2
cl> ccmap.refpoint = "user"
cl> ccmap.xref = "CRPIX1"
cl> ccmap.yref = "CRPIX2"
cl> ccmap.lngref = "CRVAL1"
cl> ccmap.latref = "CRVAL2"
cl> ccmap.refsystem = "J2000"
cl> ccmap.lngrefunits = "degrees"
cl> ccmap.latrefunits = "degrees"
cl> ccmap.projection = "tnx"
cl> ccmap.fitgeometry = "general"
cl> ccmap.function = "polynomial"
cl> ccmap.xxorder = 4
cl> ccmap.xyorder = 4
cl> ccmap.xxterms = "half"
cl> ccmap.yxorder = 4
cl> ccmap.yyorder = 4
cl> ccmap.yxterms = "half"
cl> ccmap.maxiter = 0
cl> ccmap.reject = 3.
cl> ccmap.update = "yes"
cl> ccmap.pixsystem = "logical"
cl> ccmap.verbose = "yes"
cl> ccmap.interactive = "yes"
cl> ccmap.graphics = "stdgraph"
cl> ccmap.cursor = ""
cl> ccmap  _image_.mat.1  database.db
Parameter update = "yes" pomeni, da se wcs rešitev vpiše v glavo fits slike.